Source author record

Victor P. Debattista

Victor P. Debattista 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

67works
7topics
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

67 published item(s)

preprint2026arXiv

Composite Bulges -- V. Detecting signatures of gas inflows in IFU data: The MUSE view of ionised gas kinematics in nearby galaxies

Using VLT/MUSE data, we study the ionised-gas kinematics in a mass- and volume-limited ($M_* \geq 10^{10} M_\odot$, $D \leq 20$\,Mpc) sample of 21 nearby galaxies to identify signatures of extended shocks within their inner kiloparsec, which appear as coherent velocity jumps in kinematic maps. By removing angular momentum, shocks in gas cause inflows, which can trigger nuclear star formation and fuel AGN activity. To identify the signatures of extended shocks, we examine residual velocity fields after subtracting a modelled rotating disc, and we study velocity difference between various gas tracers. Combining our kinematic analysis with BPT ionisation diagnostic maps and dust morphology, we find that 11 of 21 galaxies ($\sim$52%) show extended shock signatures with velocity jumps consistent with models of bar-driven shocks. This is likely a lower limit, as three additional galaxies ($\sim$15%) exhibit shocks along bars, potentially reaching the nucleus but obscured by AGN outflows. We trace shock signatures inwards close to the resolution limit, which suggests that shocks may be the prevailing mechanism of inflow in the central kpc of galaxies. The only two unbarred galaxies in our sample are also the only systems with unperturbed kinematics and no shocks, strongly linking the perturbed gas dynamics in centres of galaxies to the presence of bars. All galaxies with inner bars show LINER- or Seyfert-like nuclear emission, whereas galaxies without inner bars exhibit all emission types, indicating that regardless of gas supply, inner bars suppress star formation in galactic nuclei.

preprint2022arXiv

Age dissection of the vertical breathing motions in Gaia DR2: evidence for spiral driving

Gaia DR2 has revealed breathing motions in the Milky Way, with stars on both sides of the Galactic mid-plane moving coherently towards or away from it. The generating mechanism of these breathing motions is thought to be spiral density waves. Here we test this hypothesis. Using a self-consistent, high-resolution simulation with star formation, and which hosts prominent spirals, we first study the signatures of breathing motions excited by spirals. In the model, the breathing motions induced by the spiral structure have an increasing amplitude with distance from the mid-plane, pointing to an internal cause for them. We then show that, at fixed height, the breathing motion amplitude decreases with age. Next, we investigate the signature of the breathing motions in the Gaia DR2 dataset. We demonstrate that, at the location with a consistently large breathing motion, the corresponding amplitude increases monotonically with distance from the mid-plane, in agreement with the model. Furthermore, we show that at the same location, the breathing motion amplitude decreases with age, again similar to what we find in the model. This strengthens the case that the observed breathing motions are driven by spiral density waves.

preprint2022arXiv

Bending waves excited by irregular gas inflow along warps

Gaia has revealed clear evidence of bending waves in the vertical kinematics of stars in the Solar Neighbourhood. We study bending waves in two simulations, one warped, with the warp due to misaligned gas inflow, and the other unwarped. We find slow, retrograde bending waves in both models, with the ones in the warped model having larger amplitudes. We also find fast, prograde bending waves. Prograde bending waves in the unwarped model are very weak, in agreement with the expectation that these waves should decay on short, ~ crossing, timescales, due to strong winding. However, prograde bending waves are much stronger for the duration of the warped model, pointing to irregular gas inflow along the warp as a continuous source of excitation. We demonstrate that large amplitude bending waves that propagate through the Solar Neighbourhood give rise to a correlation between the mean vertical velocity and the angular momentum, with a slope consistent with that found by Gaia. The bending waves affect populations of all ages, but the sharpest features are found in the young populations, hinting that short wavelength waves are not supported by the older, kinematically hotter, populations. Our results demonstrate the importance of misaligned gas accretion as a recurrent source of vertical perturbations of disc galaxies, including in the Milky Way.

preprint2022arXiv

Excitation of vertical breathing motion in disc galaxies by tidally-induced spirals in fly-by interactions

It is now clear that the stars in the Solar neighbourhood display large-scale coherent vertical breathing motions. At the same time, Milky Way-like galaxies experience tidal interactions with satellites/companions during their evolution. While these tidal interactions can excite vertical oscillations, it is still not clear whether vertical breathing motions are excited \textit{directly} by the tidal encounters or are driven by the tidally-induced spirals. We test whether excitation of breathing motions are directly linked to tidal interactions by constructing a set of $N$-body models (with mass ratio 5:1) of unbound, single fly-by interactions with varying orbital configurations. We first reproduce the well-known result that such fly-by interactions can excite strong transient spirals (lasting for $\sim 2.9-4.2$ Gyr) in the outer disc of the host galaxy. The generation and strength of the spirals are shown to vary with the orbital parameters (the angle of interaction, and the orbital spin vector). Furthermore, we demonstrate that our fly-by models exhibit coherent breathing motions whose amplitude increases with height. The amplitudes of breathing motions show characteristic modulation along the azimuthal direction, with compressing breathing motions coinciding with the peaks of the spirals and expanding breathing motions falling in the inter-arm regions -- a signature of a spiral-driven breathing motion. These breathing motions in our models end when the strong tidally-induced spiral arms fade away. Thus, it is the tidally-induced spirals which drive the large-scale breathing motions in our fly-by models, and the dynamical role of the tidal interaction in this context is indirect.

preprint2022arXiv

The secular growth of bars revealed by flat (peak + shoulders) density profiles

The major-axis density profiles of bars are known to be either exponential or 'flat'. We develop an automated non-parametric algorithm to detect flat profiles and apply it to a suite of simulations (with and without gas). We demonstrate that flat profiles are a manifestation of a bar's secular growth, producing a 'shoulder' region (an overdensity above an exponential) in its outskirts. Shoulders are not present when bars form, but develop as the bar grows. If the bar does not grow, shoulders do not form. Shoulders are often accompanied by box/peanut bulges, but develop separately from them and are independent tracers of a bar's growth. They can be observed at a wide range of viewing orientations with only their slope varying significantly with inclination. We present evidence that shoulders are produced by looped x1 orbits. Since the growth rate of the bar moderately correlates with the growth rate of the shoulder strength, these orbits are probably recently trapped. Shoulders therefore are evidence of bar growth. The properties of the shoulders do not, however, establish the age of a bar, because secondary buckling or strong spirals may destroy shoulders, and also because shoulders do not form if the bar does not grow much. In particular, our results show that an exponential profile is not necessarily an indication of a young bar.

preprint2021arXiv

Co-formation of the thin and thick discs revealed by APOGEE-DR16 and Gaia-DR2

Since thin disc stars are younger than thick disc stars on average, the thin disc is predicted by some models to start forming after the thick disc had formed, around 10 Gyr ago. Accordingly, no significant old thin disc population should exist. Using 6-D coordinates from Gaia-DR2 and age estimates from Sanders & Das (2018), we select $\sim 24000$ old stars (${τ> 10}$ Gyr, with uncertainties $\lesssim 15\%$) within 2 kpc from the Sun (full sample). A cross-match with APOGEE-DR16 ($\sim 1000$ stars) reveals comparable fractions of old chemically defined thin/thick disc stars. We show that the full sample pericenter radius ($r_\mathrm{per}$) distribution has three peaks, one associated with the stellar halo and the other two having contributions from the thin/thick discs. Using a high-resolution $N$-body+Smooth Particle Hydrodynamics simulation, we demonstrate that one peak, at $r_\mathrm{per}\approx 7.1$ kpc, is produced by stars from both discs which were born in the inner Galaxy and migrated to the Solar Neighbourhood. In the Solar Neighbourhood, $\sim 1/2$ ($\sim 1/3$) of the old thin (thick) disc stars are classified as migrators. Our results suggest that thin/thick discs are affected differently by radial migration inasmuch as they have different eccentricity distributions, regardless of vertical scale heights. We interpret the existence of a significant old thin disc population as evidence for an early co-formation of thin/thick discs, arguing that clump instabilities in the early disc offer a compelling explanation for the observed trends.

preprint2021arXiv

Composite Bulges -- II. Classical Bulges and Nuclear Discs in Barred Galaxies: The Contrasting Cases of NGC 4608 and NGC 4643

We present detailed morphological, photometric, and stellar-kinematic analyses of the central regions of two massive, early-type barred galaxies with nearly identical large-scale morphologies. Both have large, strong bars with prominent inner photometric excesses that we associate with boxy/peanut-shaped (B/P) bulges; the latter constitute ~ 30% of the galaxy light. Inside its B/P bulge, NGC 4608 has a compact, almost circular structure (half-light radius R_e approx. 310 pc, Sérsic n = 2.2) we identify as a classical bulge, amounting to 12.1% of the total light, along with a nuclear star cluster (R_e ~ 4 pc). NGC 4643, in contrast, has a nuclear disc with an unusual broken-exponential surface-brightness profile (13.2% of the light), and a very small spheroidal component (R_e approx. 35 pc, n = 1.6; 0.5% of the light). IFU stellar kinematics support this picture, with NGC 4608's classical bulge slowly rotating and dominated by high velocity dispersion, while NGC 4643's nuclear disc shows a drop to lower dispersion, rapid rotation, V-h3 anticorrelation, and elevated h4. Both galaxies show at least some evidence for V-h3 correlation in the bar (outside the respective classical bulge and nuclear disc), in agreement with model predictions. Standard 2-component (bulge/disc) decompositions yield B/T ~ 0.5-0.7 (and bulge n > 2) for both galaxies. This overestimates the true "spheroid" components by factors of four (NGC 4608) and over 100 (NGC 4643), illustrating the perils of naive bulge-disc decompositions applied to massive barred galaxies.

preprint2021arXiv

The COMBS Survey -- II. Distinguishing the Metal-Poor Bulge from the Halo Interlopers

The metal-poor stars in the bulge are important relics of the Milky Way's formation history, as simulations predict that they are some of the oldest stars in the Galaxy. In order to determine if they are truly ancient stars, we must understand the origins of this population. Currently, it is unclear if the metal-poor stars in the bulge ([Fe/H] < -1 dex) are merely halo interlopers, a unique accreted population, part of the boxy/peanut-shaped (B/P) bulge or a classical bulge population. In this work, we use spectra from the VLT/FLAMES spectrograph to obtain metallicity estimates using the Ca-II triplet (CaT) of 473 bulge stars (187 of which have [Fe/H]<-1 dex), targeted using SkyMapper photometry. We also use Gaia DR2 parallaxes and proper motions to infer the Galactic positions and velocities along with orbital properties for 523 bulge stars. We employ a probabilistic orbit analysis and find that about half of our sample has a > 50\% probability of being bound to the bulge, and half are halo interlopers. We also see that the occurrence rate of halo interlopers increases steadily with decreasing metallicity across the full range of our sample (-3 < [Fe/H] < 0.5). Our examination of the kinematics of the confined compared to the unbound stars indicates the metal-poor bulge ([Fe/H] < -1 dex) comprises at least two populations; those confined to the boxy/peanut bulge and halo stars passing through the inner galaxy. We conclude that an orbital analysis approach, as we have employed, is important to uncover and understand the composite nature of the metal-poor stars in the inner region.

preprint2021arXiv

Tracing birth properties of stars with abundance clustering

To understand the formation and evolution of the Milky Way disk, we must connect its current properties to its past. We explore hydrodynamical cosmological simulations to investigate how the chemical abundances of stars might be linked to their origins. Using hierarchical clustering of abundance measurements in two Milky Way-like simulations with distributed and steady star formation histories, we find that abundance clusters of stars comprise different groups in birth place ($R_\text{birth}$) and time (age). Simulating observational abundance errors (0.05 dex), we find that to trace discrete groups of ($R_\text{birth}$, age) requires a large vector of abundances. Using 15-element abundances (Fe, O, Mg, S, Si, C, P, Mn, Ne, Al, N, V, Ba, Cr, Co), up to $\approx$ 10 clusters can be defined with $\approx$ 25% overlap in ($R_\text{birth}$, age). We build a simple model to show that it is possible to infer a star's age and $R_\text{birth}$ from abundances with precisions of $\pm$0.06 Gyr and $\pm$1.17 kpc respectively. We find that abundance clustering is ineffective for a third simulation, where low-$α$ stars form distributed in the disc and early high-$α$ stars form more rapidly in clumps that sink towards the galactic center as their constituent stars evolve to enrich the interstellar medium. However, this formation path leads to large age-dispersions across the [$α$/Fe]-[Fe/H] plane, which is inconsistent with the Milky Way's observed properties. We conclude that abundance clustering is a promising approach toward charting the history of our Galaxy.

preprint2020arXiv

Geometric properties of galactic discs with clumpy episodes

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

preprint2020arXiv

Kinematic decomposition of IllustrisTNG disk galaxies: morphology and relation with morphological structures

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

preprint2020arXiv

Relations among structural parameters in barred galaxies with a direct measurement of bar pattern speed

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

preprint2020arXiv

The Splash without a merger

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

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

A New Channel of Bulge Formation via The Destruction of Short Bars

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

preprint2016arXiv

Caught in the Act: Direct Detection of Galactic Bars in the Buckling Phase

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

preprint2016arXiv

Imprints of radial migration on the Milky Way's metallicity distribution functions

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

preprint2016arXiv

Kinematic properties of double-barred galaxies: simulations vs. integral-field observations

Using high resolution $N$-body simulations, we recently reported that a dynamically cool inner disk embedded in a hotter outer disk can naturally generate a steady double-barred (S2B) structure. Here we study the kinematics of these S2B simulations, and compare them to integral-field observations from ATLAS$^{3D}$ and SAURON. We show that S2B galaxies exhibit several distinct kinematic features, namely: (1) significantly distorted isovelocity contours at the transition region between the two bars, (2) peaks in $σ_\mathrm{LOS}$ along the minor axis of inner bars, which we term "$σ$-humps", that are often accompanied by ring/spiral-like features of increased $σ_\mathrm{LOS}$, (3) $h_3-\bar{v}$ anti-correlations in the region of the inner bar for certain orientations, and (4) rings of positive $h_4$ when viewed at low inclinations. The most impressive of these features are the $σ$-humps, these evolve with the inner bar, oscillating in strength just as the inner bar does as it rotates relative to the outer bar. We show that, in cylindrical coordinates, the inner bar has similar streaming motions and velocity dispersion properties as normal large-scale bars, except for $σ_z$, which exhibits peaks on the minor axis, i.e., humps. These $σ_z$ humps are responsible for producing the $σ$-humps. For three well-resolved early-type S2Bs (NGC 2859, NGC 2950, and NGC 3941) and a potential S2B candidate (NGC 3384), the S2B model qualitatively matches the integral-field data well, including the "$σ$-hollows" previously identified. We also discuss the kinematic effect of a nuclear disk in S2Bs.

preprint2016arXiv

The origin of type-I profiles in cluster lenticulars: An interplay between ram pressure stripping and tidally-induced spiral migration

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

preprint2015arXiv

A Kiloparsec-Scale Nuclear Stellar Disk in the Milky Way as a Possible Explanation of the High Velocity Peaks in the Galactic Bulge

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

preprint2015arXiv

A unified framework for the orbital structure of bars and triaxial ellipsoids

We examine a large random sample of orbits in self-consistent simulations of N-body bars. Orbits in the bars are classified both visually and with a new automated orbit classification method based on frequency analysis. The well known prograde x1 orbit family originates from the same parent orbit as the box orbits in stationary and rotating triaxial ellipsoids. However only a small fraction of bar orbits ~4% have predominately prograde motion like their periodic parent orbit. Most bar orbits arising from the x1 orbit have little net angular momentum in the bar frame making them equivalent to box orbits in rotating triaxial potentials. A small fraction of bar orbits (~7%) are long axis tubes that behave exactly like those in triaxial ellipsoids:they are tipped about the intermediate-axis due to the Coriolis force, with the sense of tipping determined by the sign of their angular momentum about the long axis. No orbits parented by prograde periodic x2 orbits are found in the pure bar model, but a tiny population (~2%) of short axis tube orbits parented by retrograde x4 orbits are found. When a central point mass representing a supermassive black hole (SMBH) is grown adiabatically at the center of the bar, those orbits that lie in the immediate vicinity of the SMBH are transformed into precessing Keplerian orbits (PKOs) which belong the same major families (short axis tubes, long axis tubes and boxes) occupying the bar at larger radii. During the growth of a SMBH the inflow of mass and outward transport of angular momentum transforms some x1 and long axis tube orbits into prograde short axis tubes. This study has important implications for future attempts to constrain the masses of SMBHs in barred galaxies using orbit based methods like the Schwarzschild orbit superposition scheme and for understanding the observed features in barred galaxies.

preprint2015arXiv

Forming Double-barred Galaxies From Dynamically Cool Inner Disks

About one third of early-type barred galaxies host small-scale secondary bars. The formation and evolution of such double-barred galaxies remain far from being well understood. In order to understand the formation of such systems, we explore a large parameter space of isolated pure-disk simulations. We show that a dynamically cool inner disk embedded in a hotter outer disk can naturally generate a steady secondary bar while the outer disk forms a large-scale primary bar. The independent bar instabilities of inner and outer disks result in long-lived double-barred structures whose dynamical properties are comparable with observations. This formation scenario indicates that the secondary bar might form from the general bar instability, the same as the primary bar. Under some circumstances, the interaction of the bars and the disk leads to the two bars aligning or single, nuclear, bars only. Simulations that are cool enough of the center to experience clump instabilities may also generate steady double-barred galaxies. In this case, the secondary bars are "fast", i.e., the bar length is close to the co-rotation radius. This is the first time that double-barred galaxies containing a fast secondary bar are reported. Previous orbit-based studies had suggested that fast secondary bars are not dynamically possible.

preprint2015arXiv

Internal Alignments of Red Versus Blue Discs in Dark Matter Halos

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

preprint2015arXiv

Measuring the mass of the central black hole in the bulgeless galaxy NGC 4395 from gas dynamical modeling

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

preprint2015arXiv

Nuclear Star Clusters and Bulges

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

preprint2015arXiv

Predicting the stellar and non-equilibrium dust emission spectra of high-resolution simulated galaxies with DART-Ray

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

preprint2015arXiv

Reinforcing the link between the double red clump and the X-shaped bulge of the Milky Way

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

preprint2015arXiv

Stellar orbital properties as diagnostics of the origin of the stellar halo

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

preprint2015arXiv

The effect of stellar feedback on a Milky Way-like galaxy and its gaseous halo

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

preprint2014arXiv

N-body simulation insights into the X-shaped bulge of the Milky Way: Kinematics and distance to the Galactic Centre

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

preprint2014arXiv

Re-interpretation of "Bar slowdown and the distribution of dark matter in barred galaxies" by Athanassoula

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

preprint2014arXiv

Secular- and merger-built bulges in barred galaxies

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

preprint2014arXiv

Systematic Problems With Using Dark Matter Simulations to Model Stellar Halos

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

preprint2014arXiv

The formation of stellar nuclear discs in bar-induced gas inflows

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

preprint2014arXiv

The Vertical Structure and Kinematics of Grand Design Spirals

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

preprint2013arXiv

3D kinematics through the X-shaped Milky Way bulge

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

preprint2013arXiv

Disk Assembly and the M_BH-sigma Relation of Supermassive Black Holes

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

preprint2013arXiv

Halo orbits in cosmological disk galaxies: tracers of formation history

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

preprint2013arXiv

On the Offset of Barred Galaxies From the Black Hole $M_{BH}$-$σ$ Relationship

We use collisionless $N$-body simulations to determine how the growth of a supermassive black hole (SMBH) influences the nuclear kinematics in both barred and unbarred galaxies. In the presence of a bar, the increase in the velocity dispersion $σ$ (within the effective radius) due to the growth of an SMBH is on average $\lesssim 10%$, whereas the increase is only $\lesssim 4%$ in an unbarred galaxy. In a barred galaxy, the increase results from a combination of three separate factors (a) orientation and inclination effects; (b) angular momentum transport by the bar that results in an increase in the central mass density; (c) an increase in the vertical and radial velocity anisotropy of stars in the vicinity of the SMBH. In contrast the growth of the SMBH in an unbarred galaxy causes the velocity distribution in the inner part of the nucleus to become less radially anisotropic. The increase in $σ$ following the growth of the SMBH is insensitive to a variation of a factor of 10 in the final mass of the SMBH, showing that it is the growth process rather than the actual SMBH mass that alters bar evolution in a way that increases $σ$. We argue that using an axisymmetric stellar dynamical modeling code to measure SMBH masses in barred galaxies could result in a slight overestimate of the derived $M_{BH}$, especially if a constant M/L ratio is assumed. We conclude that the growth of a black hole in the presence of a bar could result in an increase in $σ$ which is roughly of $4-8%$ larger than the increase that occurs in an axisymmetric system. While the increase in $σ$ due to SMBH growth in a barred galaxy might partially account for the claimed offset of barred galaxies and pseudo bulges from the $M_{BH}-σ$ relation obtained for elliptical galaxies and classical bulges in unbarred galaxies, it is inadequate to account for all of the offset.

preprint2013arXiv

Peanuts at an Angle: Detecting and Measuring the Three-Dimensional Structure of Bars in Moderately Inclined Galaxies

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

preprint2013arXiv

Searching for nuclear stellar discs in simulations of star cluster mergers

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

preprint2013arXiv

Stellar Ages through the Corners of the Boxy Bulge

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

preprint2013arXiv

The effect of bars on the M-sigma relation: offset, scatter and residuals correlations

We analyse a set of collisionless disc galaxy simulations to study the consequences of bar formation and evolution on the M-sigma relation of supermassive black holes. The redistribution of angular momentum driven by bars leads to a mass increase within the central region, raising the velocity dispersion of the bulge, sigma, on average by ~12% and as much as ~20%. If a disc galaxy with a SMBH satisfying the M-sigma relation forms a bar, and the SMBH does not grow in the process, then the increase in sigma moves the galaxy off the M-sigma relation. We explore various effects that can affect this result including contamination from the disc and anisotropy. The displacement from the M-sigma relation for individual model barred galaxies correlates with both M(B)/M(B+D) and beta_phi(B+D) measured within the effective radius of the bulge. Overall, this process leads to an M-sigma for barred galaxies offset from that of unbarred galaxies, as well as an increase in its scatter. We assemble samples of unbarred and barred galaxies with classical bulges and find tentative hints of an offset between the two consistent with the predicted. Including all barred galaxies, rather than just those with a classical bulge, leads to a significantly larger offset.

preprint2013arXiv

The Effects of Radial Migration on the Vertical Structure of Galactic Discs

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

preprint2013arXiv

What's Up in the Milky Way? The Orientation of the Disc Relative to the Triaxial Halo

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

preprint2012arXiv

Outer-Disk Populations in NGC 7793: Evidence for Stellar Radial Migration

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

preprint2012arXiv

Spatially Resolved Spectroscopic Star Formation Histories of Nearby Disks: Hints of Stellar Migration

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

preprint2011arXiv

Constraining the formation of inner bars. Photometry, kinematics and stellar populations in NGC 357

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

preprint2011arXiv

Constraining the role of star cluster mergers in nuclear cluster formation: Simulations confront integral-field data

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

preprint2011arXiv

Probing the shape and history of the Milky Way halo with orbital spectral analysis

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

preprint2010arXiv

IFU Spectroscopy of the Stellar Disk Truncation Region of NGC 6155

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

preprint2010arXiv

Misaligned Angular Momentum in Hydrodynamic Cosmological Simulations: Warps, Outer Discs, and Thick Discs

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

preprint2010arXiv

New tools for probing the phase space structure of dark matter halos

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

preprint2010arXiv

The ACS Nearby Galaxy Survey Treasury V. Radial Star Formation History of NGC 300

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

preprint2010arXiv

The Genesis of the Milky Way's Thick Disk via Stellar Migration

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

preprint2010arXiv

The NGC 404 Nucleus: Star Cluster and Possible Intermediate Mass Black Hole

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

preprint2009arXiv

LSST Science Book, Version 2.0

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

preprint2009arXiv

Orbital Eccentricity as a probe of Thick Disk Formation Scenarios

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

preprint2009arXiv

Stochasticity in N-body Simulations of Disc Galaxies

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

preprint2008arXiv

Beyond Inside-Out Growth: Formation and Evolution of Disk Outskirts

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

preprint2008arXiv

The Rotating Nuclear Star Cluster in NGC 4244

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

preprint2007arXiv

The bar pattern speed of NGC 4431

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

preprint2006arXiv

The Binary Nucleus in VCC 128: A Candidate Supermassive Black Hole in a Dwarf Elliptical Galaxy

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

preprint2005arXiv

Anomalously Weak Dynamical Friction in Halos

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

preprint2005arXiv

Fate of Supermassive Black Holes and Evolution of the M(BH)-sigma Relation in Merging Galaxies: The Effect of Gaseous Dissipation

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

preprint1997arXiv

Dynamical Friction and the Distribution of Dark Matter in Barred Galaxies

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

preprint1996arXiv

Secular Evolution in Barred Galaxies

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