Source author record

E. Athanassoula

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

116works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

116 published item(s)

preprint2026arXiv

Gasflows in Barred Galaxies with Big Orbital Loops-A Comparative Study of Two Hydrocodes

We study the flow of gas in a barred-galaxy model, in which a considerable part of the underlying stable periodic orbits have loops where, close to the ends of the bar, several orbital families coexist and chaos dominates. Such conditions are typically encountered in a zone between the 4:1 resonance and corotation. The purpose of our study is to understand the gaseous flow in the aforementioned environment and trace the morphology of the shocks that form. We use two conceptually different hydrodynamic schemes for our calculations, namely, the mesh-free Lagrangian SPH method and the adaptive mesh refinement code RAMSES. This allows us to compare responses by means of the two algorithms. We find that the big loops of the orbits, mainly belonging to the x1 stable periodic orbits, do not help the shock loci to approach corotation. They deviate away from the regions occupied by the loops, bypass them and form extensions at an angle with the straight-line shocks. Roughly at the distance from the center at which we start to observe the big loops, we find characteristic "tails" of dense gas streaming towards the straight-line shocks. The two codes give complementary information for understanding the hydrodynamics of the models.

preprint2026arXiv

Morphologies arising from the gas flow in the innermost kiloparsec of barred galaxy models

Context. We study a series of response models to investigate the formation of specific morphological features in the central 1 kpc region of the gas component in barred spiral galaxies. Aims. We aim to understand how structures, such as nuclear rings and spirals, form by varying the parameters of a general gravitational potential and gas properties. Our goal is to determine how much the shape of these structures is driven by the orbital dynamics of the models compared to the influence of the hydrodynamics of the gas. In particular, we examine the effects of the bar strength, bar shape, pattern speed, and central density, as well as their mutual interdependence. Methods. We modeled the gas flow using hydrodynamical simulations run with the Eulerian RAMSES code. The underlying gravitational potential was a two-dimensional Ferrers bar and the gas was considered to be isothermal. Alongside analyzing the gas response to the imposed gravitational potentials, we carried out orbital studies for all models. This involved assessing the shapes and stability of periodic orbits and analyzing the distribution of regular versus chaotic regions within the systems. Results. The parameters of the gravitational potential alone are insufficient to accurately predict the gas dynamics in a system. The morphology of the gaseous response varies substantially with changes in sound speed, emphasizing the fundamental role of hydrodynamic processes in determining the structure of the gas within the central region. We identify the factors that affect the morphology of nuclear rings and trailing and leading nuclear spirals. The best alignment between our models and structures observed in local barred galaxies is achieved by assuming a sound speed of $c_s=20\,\rm{km\,s^{-1}}$.

preprint2022arXiv

Stellar masses, sizes, and radial profiles for 465 nearby early-type galaxies: an extension to the Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G)

The Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G) is a detailed study of over 2300 nearby galaxies in the near-infrared (NIR), which has been critical to our understanding of the detailed structures of nearby galaxies. Because the sample galaxies were selected only using radio-derived velocities, however, the survey favored late-type disk galaxies over lenticulars and ellipticals. A follow-up Spitzer survey was conducted to rectify this bias, adding 465 early-type galaxies (ETGs) to the original sample, to be analyzed in a manner consistent with the initial survey. We present the data release of this ETG extension, up to the third data processing pipeline (P3): surface photometry. We produce curves of growth and radial surface brightness profiles (with and without inclination corrections) using reduced and masked Spitzer IRAC 3.6$μ$m and 4.5$μ$m images produced through Pipelines 1 and 2, respectively. From these profiles, we derive the following integrated quantities: total magnitudes, stellar masses, concentration parameters, and galaxy size metrics. We showcase NIR scaling relations for ETGs among these quantities. We examine general trends across the whole S$^{4}$G and ETG extension among our derived parameters, highlighting differences between ETGs and late-type galaxies (LTGs). ETGs are, on average, more massive and more concentrated than LTGs, and also show subtle distinctions among ETG morphological sub-types. We also derive the following scaling relations and compare with previous results in visible light: mass--size (both half-light and isophotal), mass--concentration, mass--surface brightness (central, effective, and within 1 kpc), and mass--color. We find good agreement with previous works, though some relations (e.g., mass--central surface brightness) will require more careful multi-component decompositions to be fully understood.

preprint2020arXiv

Basis function expansions for galactic dynamics: Spherical versus cylindrical coordinates

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)

preprint2020arXiv

Disc galaxies formed from major mergers in Illustris

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.

preprint2020arXiv

Suppressed or enhanced central star formation rates in late-type barred galaxies

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.

preprint2020arXiv

The orbital content of bars: The origin of 'non-x1-tree', bar-supporting orbits

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.

preprint2020arXiv

Tidally induced warps of spiral galaxies in IllustrisTNG

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.

preprint2019arXiv

Fundamental Physics with the Square Kilometre Array

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.

preprint2019arXiv

The sequence of spiral arm classes: Observational signatures of persistent spiral density waves in grand-design galaxies

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$μ$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.

preprint2016arXiv

A close look at secular evolution: Boxy/peanut bulges reduce gas inflow to the central kiloparsec

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.

preprint2016arXiv

Evidence of Bar-induced Secular Evolution in the Inner Regions of Stellar Discs in Galaxies: What Shapes Disc Galaxies?

We present evidence of bar-induced secular evolution in galactic discs using 3.6 ${μ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 $Θ$-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.

preprint2016arXiv

Forming disk galaxies in wet major mergers. I. Three fiducial examples

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.

preprint2016arXiv

Forming first-ranked early-type galaxies through hierarchical dissipationless merging

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.

preprint2016arXiv

Halpha Kinematics of S4G Spiral Galaxies - III. Inner rotation curves

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 (μ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 μ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.

preprint2016arXiv

Morpho-kinematics of z~1 galaxies probe the hierarchical scenario

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 ΛCDM models, supporting a merger-dominated hierarchical scenario as being the main driver of galaxy formation at least during the last 8 billion years.

preprint2016arXiv

Morphology and environment of galaxies with disc breaks in the S4G and NIRS0S

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.

preprint2016arXiv

On the nature of the barlens component in barred galaxies: what do boxy/peanut bulges look like when viewed face-on?

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.

preprint2016arXiv

Orbital classification in an N-body bar

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.

preprint2016arXiv

The orbital structure of a tidally induced bar

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.

preprint2016arXiv

Tidally induced bars of galaxies in clusters

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.

preprint2016arXiv

Tidally Induced Offset Disks in Magellanic Spiral Galaxies

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.

preprint2015arXiv

2D kinematic signatures of boxy/peanut bulges

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.

preprint2015arXiv

A Classical Morphological Analysis of Galaxies in the Spitzer Survey of Stellar Structure in Galaxies (S4G)

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.

preprint2015arXiv

APOGEE Kinematics I: Overview of the Kinematics of the Galactic Bulge as Mapped by APOGEE

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.

preprint2015arXiv

Applying Schwarzschild's orbit superposition method to barred or non-barred disc galaxies

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.

preprint2015arXiv

ARRAKIS: Atlas of Resonance Rings As Known In the S4G

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.

preprint2015arXiv

Boxy/peanut/X bulges, barlenses and the thick part of galactic bars: What are they and how did they form?

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.

preprint2015arXiv

Environmental regulation of cloud and star formation in galactic bars

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.

preprint2015arXiv

Evolution of the Milky Way with radial motions of stars and gas I. The solar neighborhood and the thin and thick disk

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.

preprint2015arXiv

Globular Cluster Populations: Results Including S$^4$G Late-Type Galaxies

Using 3.6 and 4.5$μ$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.

preprint2015arXiv

H-alpha kinematics of S4G spiral galaxies-II. Data description and non-circular motions

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ées 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.

preprint2015arXiv

Spectral analysis of stellar orbits in a tidally induced bar

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.

preprint2015arXiv

Stellar Populations of Barred Quiescent Galaxies

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.

preprint2015arXiv

The Effects of Boxy/Peanut Bulges on Galaxy Models

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.

preprint2015arXiv

The Odd Offset between the Galactic Disk and Its Bar in NGC 3906

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.

preprint2015arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S$^4$G): Precise Stellar Mass Distributions from Automated Dust Correction at 3.6 microns

The mid-infrared is an optimal window to trace stellar mass in nearby galaxies and the 3.6$μ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$μ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$μ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$μ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$μ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.

preprint2015arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S4G): Stellar Masses, Sizes and Radial Profiles for 2352 Nearby Galaxies

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

preprint2014arXiv

A ~ 12 kpc HI extension and other HI asymmetries in the isolated galaxy CIG 340 (IC 2487)

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.

preprint2014arXiv

Adventures of a tidally induced bar

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.

preprint2014arXiv

Evolution of the Milky Way with radial motions of stars and gas II. The evolution of abundance profiles from H to Ni

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.

preprint2014arXiv

Galaxy Zoo: Are Bars Responsible for the Feeding of Active Galactic Nuclei at 0.2 < z < 1.0?

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.

preprint2014arXiv

Globular Cluster Populations: First Results from S$^4$G Early-Type Galaxies

Using 3.6$μ$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.

preprint2014arXiv

IRAC Near-Infrared Features in the Outer Parts of S4G Galaxies

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.

preprint2014arXiv

Milky Way mass galaxies with X-shaped bulges are not rare in the local Universe

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.

preprint2014arXiv

On the Correlation Between Metallicity and the X-Shaped Morphology of the Milky Way Bulge

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.

preprint2014arXiv

Reconstructing the stellar mass distributions of galaxies using S4G IRAC 3.6 and 4.5 micron images: II. The conversion from light to mass

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.

preprint2014arXiv

The Baryonic Tully-Fisher Relationship for S$^4$G Galaxies and the "Condensed" Baryon Fraction of Galaxies

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$ ($Δ$ log $M_{baryon}/Δ$ 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.

preprint2014arXiv

The Mass Profile and Shape of Bars in the Spitzer Survey of Stellar Structure in Galaxies (S4G): Search for an Age Indicator for Bars

We have measured the radial light profiles and global shapes of bars using two-dimensional 3.6 $μ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.

preprint2014arXiv

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

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

preprint2013arXiv

ARGOS IV: The Kinematics of the Milky Way Bulge

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.

preprint2013arXiv

Bar slowdown and the distribution of dark matter in barred galaxies

`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).

preprint2013arXiv

Embedded Star Formation in S4G Galaxy Dust Lanes

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.

preprint2013arXiv

Galaxy Zoo: Observing Secular Evolution Through Bars

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$).

preprint2013arXiv

Hα Kinematics of S4G spiral galaxies-I. NGC 864

We present a study of the kinematics of the isolated spiral galaxy NGC 864, using Hα Fabry-Perot data obtained with the GHαFaS instrument at the William Herschel Telescope in La Palma, complemented with images at 3.6 μm, in the R band and in Hα 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.

preprint2013arXiv

Morphological Parameters of Spitzer Survey of Stellar Structure in Galaxies (S4G)

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.

preprint2013arXiv

On the Origin of Lopsidedness in Galaxies as Determined from the Spitzer Survey of Stellar Structure in Galaxies (S4G)

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.

preprint2013arXiv

Radial migration in a bar-dominated disk galaxy I: Impact on chemical evolution

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

preprint2013arXiv

Statistics of the structure components in S0s: implications for bar induced secular evolution

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.

preprint2013arXiv

The effect of softening on dynamical simulations of galaxies

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.

preprint2013arXiv

The impact of bars on disk breaks as probed by S4G imaging

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.

preprint2013arXiv

Unveiling the Structure of Barred Galaxies at 3.6 $μ{\rm m}$ with the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G): I. Disk Breaks

We have performed two-dimensional multicomponent decomposition of 144 local barred spiral galaxies using 3.6 $μ{\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.

preprint2013arXiv

X-ray nuclear activity in S4G barred galaxies: No link between bar strength and co-occurrent supermassive black hole fueling

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.

preprint2012arXiv

A unified picture of breaks and truncations in spiral galaxies from SDSS and S^{4}G imaging

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$μ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.

preprint2012arXiv

ARGOS II: The Galactic Bulge Survey

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.

preprint2012arXiv

Bar formation and evolution in disc galaxies with gas and a triaxial halo: Morphology, bar strength and halo properties

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.

preprint2012arXiv

Bars and secular evolution in disk galaxies: Theoretical input

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.

preprint2012arXiv

Breaks in thin and thick discs of edge-on galaxies imaged in the Spitzer Survey of Stellar Structure in Galaxies (S4G)

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.

preprint2012arXiv

Chaotic mixing and the secular evolution of triaxial cuspy galaxy models built with Schwarzschild's method

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.

preprint2012arXiv

Early type galaxies with tidal debris and their scaling relations in the Spitzer Survey of Stellar Structure in Galaxies (S4G)

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

preprint2012arXiv

Galaxy disks do not need to survive in the L-CDM paradigm: the galaxy merger rate out to z~1.5 from morpho-kinematic data

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.

preprint2012arXiv

HI asymmetry in the isolated galaxy CIG 85 (UGC 1547)

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.

preprint2012arXiv

Hot Disks And Delayed Bar Formation

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.

preprint2012arXiv

Loops formed by tidal tails as fossil records of a major merger

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.

preprint2012arXiv

Manifold-driven spirals in N-body barred galaxy simulations

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.

preprint2012arXiv

Simulations of Barred Galaxies in Triaxial Dark Matter Haloes: The Effects of Gas

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.

preprint2012arXiv

The origin of the split red clump in the Galactic bulge of the Milky Way

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

preprint2012arXiv

The S4G perspective on circumstellar dust extinction of AGB stars in M100

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.

preprint2012arXiv

Towards understanding the dynamics of the bar/bulge region in our Galaxy

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.

preprint2011arXiv

Extensions and applications of the iterative method

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.

preprint2011arXiv

HI asymmetries in the isolated galaxy CIG 292

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.

preprint2011arXiv

Modelling the inner disc of the Milky Way with manifolds. I - A first step

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.

preprint2011arXiv

Reconstructing the stellar mass distributions of galaxies using S4G IRAC 3.6 and 4.5 micron images: I. Correcting for contamination by PAH, hot dust, and intermediate age stars

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.

preprint2011arXiv

Regular and chaotic orbits in barred galaxies - I. Applying the SALI/GALI method to explore their distribution in several models

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.

preprint2011arXiv

Star cluster evolution in barred disc galaxies. I. Planar periodic orbits

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.

preprint2011arXiv

Star formation in the extended gaseous disk of the isolated galaxy CIG 96

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.

preprint2011arXiv

The Formation of Large Galactic Disks: Revival or Survival?

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.

preprint2011arXiv

The unusual vertical mass distribution of NGC 4013 seen through the Spitzer Survey of Stellar Structure in Galaxies (S4G)

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.

preprint2011arXiv

Thick disks of edge-on galaxies seen through the Spitzer Survey of Stellar Structure in Galaxies (S4G): Lair of missing baryons?

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

preprint2010arXiv

Dark Matter Direct Detection Signals inferred from a Cosmological N-body Simulation with Baryons

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.

preprint2010arXiv

Dynamical models of the elliptical galaxy NGC 4494

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.

preprint2010arXiv

Loss of halo triaxiality due to bar formation

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]

preprint2010arXiv

Mid-Infrared Galaxy Morphology from the Spitzer Survey of Stellar Structure in Galaxies (S^4G): The Imprint of the de Vaucouleurs Revised Hubble-Sandage Classification System at 3.6 microns

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.

preprint2010arXiv

Rings and spirals in barred galaxies. III. Further comparisons and links to observations

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.

preprint2010arXiv

Star Clusters in M31: Old Clusters with Bar Kinematics

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.

preprint2010arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S^4G)

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.

preprint2010arXiv

The thick disk in the galaxy NGC 4244 from S4G imaging

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.

preprint2009arXiv

Rings and spirals in barred galaxies. I Building blocks

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.

preprint2009arXiv

Rings and spirals in barred galaxies. II. Ring and spiral morphology

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.

preprint2009arXiv

The Hubble sequence: just a vestige of merger events?

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.

preprint2008arXiv

Boxy/peanut bulges : formation, evolution and properties

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.

preprint2008arXiv

Dynamical study of 2D and 3D barred galaxy models

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é $(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.

preprint2008arXiv

Invariant manifolds as building blocks for the formation of spiral arms and rings in barred galaxies

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.

preprint2008arXiv

The role of the unstable equilibrium points in the transfer of matter in galactic potentials

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.

preprint2007arXiv

Evolution of the Bar Fraction in COSMOS: Quantifying the Assembly of the Hubble Sequence

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.

preprint2007arXiv

Globular clusters and dwarf galaxies in Fornax - I. Kinematics in the cluster core from multi-object spectroscopy

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.

preprint2007arXiv

Studying the global dynamics of conservative dynamical systems using the SALI chaos detection method

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.

preprint2006arXiv

The effect of giant molecular clouds on star clusters

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.

preprint2005arXiv

Halo properties and secular evolution in barred galaxies

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.

preprint2005arXiv

On the nature of bulges in general and of box/peanut bulges in particular. Input from $N$-body simulations

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.

preprint2004arXiv

Bar-Driven Evolution and 2D Spectroscopy of Bulges

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]

preprint2003arXiv

Angular momentum redistribution and the evolution and morphology of bars

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.

preprint2002arXiv

Bar-Halo Interaction and Bar Growth

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.

preprint1997arXiv

Hickson 96: a physical compact group

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.

preprint1996arXiv

Evolution of Compact Groups of Galaxies I. Merging Rates

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.

preprint1996arXiv

Formation of rings in galactic discs by infalling small companions

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.