We model for the first time the complete orbital evolution of a pair of Supermassive Black Holes (SMBHs) in a 1:10 galaxy merger of two disk dominated gas-rich galaxies, from the stage prior to the formation of the binary up to the onset of gravitational wave emission when the binary separation has shrunk to 1 milli parsec. The high-resolution smoothed particle hydrodynamics (SPH) simulations used for the first phase of the evolution include star formation, accretion onto the SMBHs as well as feedback from supernovae explosions and radiative heating from the SMBHs themselves. Using the direct N-body code \phi-GPU we evolve the system further without including the effect of gas, which has been mostly consumed by star formation in the meantime. We start at the time when the separation between two SMBHs is ~ 700 pc and the two black holes are still embedded in their galaxy cusps. We use 3 million particles to study the formation and evolution of the SMBH binary till it becomes hard. After a hard binary is formed, we reduce (reselect) the particles to 1.15 million and follow the subsequent shrinking of the SMBH binary due to 3-body encounters with the stars. We find approximately constant hardening rates and that the SMBH binary rapidly develops a high eccentricity. Similar hardening rates and eccentricity values are reported in earlier studies of SMBH binary evolution in the merging of dissipation-less spherical galaxy models. The estimated coalescence time is ~ 2.9 Gyr, significantly smaller than a Hubble time. We discuss why this timescale should be regarded as an upper limit. Since 1:10 mergers are among the most common interaction events for galaxies at all cosmic epochs, we argue that several SMBH binaries should be detected with currently planned space-borne gravitational wave interferometers, whose sensitivity will be especially high for SMBHs in the mass range considered here.

Galaxy centers are residing places for Super Massive Black Holes (SMBHs). Galaxy mergers bring SMBHs close together to form gravitationally bound binary systems which, if able to coalesce in less than a Hubble time, would be one of the most promising sources of gravitational waves for the Laser Interferometer Space Antenna (LISA). In spherical galaxy models, SMBH binaries stall at a separation of approximately one parsec, leading to the "final parsec problem" (FPP). On the other hand, it has been shown that merger-induced triaxiality of the remnant in equal-mass mergers is capable of supporting a constant supply of stars on so-called centrophilic orbits that interact with the binary and thus avoid the FPP. In this paper, using a set of direct N-body simulations of mergers of initially spherically symmetric galaxies with different mass ratios, we show that the merger-induced triaxiality is able to drive unequal-mass SMBH binaries to coalescence. The binary hardening rates are high and depend only weakly on the mass ratios of SMBHs for a wide range of mass ratios q. The hardening rates are significantly higher for galaxies having steep cusps in comparison with those having shallow cups at centers. The evolution of the binary SMBH leads to relatively shallower inner slopes at the centers of the merger remnants. The stellar mass displaced by the SMBH binary on its way to coalescence is ~ 1-5 times the combined mass of binary SMBHs. The coalescence times for SMBH binary with mass ~ million solar masses are less than 1 Gyr and for those at the upper end of SMBH masses (~ billion solar masses) are 1-2 Gyr for less eccentric binaries whereas less than 1 Gyr for highly eccentric binaries. SMBH binaries are thus expected to be promising sources of gravitational waves at low and high redshifts.

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

We investigate a purely stellar dynamical solution to the Final Parsec Problem. Galactic nuclei resulting from major mergers are not spherical, but show some degree of triaxiality. With $N$-body simulations, we show that massive black hole binaries (MBHB) hosted by them will continuously interact with stars on centrophilic orbits and will thus inspiral---in much less than a Hubble time---down to separations at which gravitational wave (GW) emission is strong enough to drive them to coalescence. Such coalescences will be important sources of GWs for future space-borne detectors such as the {\it Laser Interferometer Space Antenna} (LISA). Based on our results, we expect that LISA will see between $\sim 10$ to $\sim {\rm few} \times 10^2$ such events every year, depending on the particular MBH seed model as obtained in recent studies of merger trees of galaxy and MBH co-evolution. Orbital eccentricities in the LISA band will be clearly distinguishable from zero with $e \gtrsim 0.001-0.01$.

This paper studies the formation and evolution of binary supermassive black holes (SMBHs) in rotating galactic nuclei, focusing on the role of stellar dynamics. We present the first N-body simulations that follow the evolution of the SMBHs from kiloparsec separations all the way to their final relativistic coalescence, and that can robustly be scaled to real galaxies. The N-body code includes post-Newtonian (PN) corrections to the binary equations of motion up to order 2.5; we show that the evolution of the massive binary is only correctly reproduced if the conservative 1PN and 2PN terms are included. The orbital eccentricities of the massive binaries in our simulations are often found to remain large until shortly before coalescence. This directly affects not only their orbital evolution rates, but has important consequences as well for the gravitational waveforms emitted during the relativistic inspiral. We estimate gravitational wave amplitudes when the frequencies fall inside the band of the (planned) Laser Interferometer Space Antennae (LISA). We find significant contributions -- well above the LISA sensitivity curve -- from the higher-order harmonics.

We present preliminary results from self-consistent, high resolution direct {\it N}-body simulations of massive black hole binaries in mergers of galactic nuclei. The dynamics of the black hole binary includes the full Post-Newtonian corrections (up to 2.5PN) to its equations of motion. We show that massive black holes starting at separations of 100 pc can evolve down to gravitational-wave-induced coalescence in less than a Hubble time. The binaries, in our models, often form with very high eccentricity and, as a result, reach separations of 50 Schwarzschild radius with eccentricities which are clearly distinct from zero -- even though gravitational wave emission damps the eccentricity during the inspiral. These deviations from exact circular orbits, at such small separations, may have important consequences for LISA data analysis.

SUPERBOX is known as a very efficient particle-mesh code with highly-resolving sub-grids. Nevertheless, the height of a typical galactic disc is small compared to the size of the whole system. Consequently, the numerical resolution in z-direction, i.e. vertically with respect to the plane of the disc, remains poor. Here, we present a new version of SUPERBOX that allows for a considerably higher resolution along z. The improved code is applied to investigate disc heating by the infall of a galaxy satellite. We describe the improvement and communicate our results. As an important application we discuss the disruption of a dwarf galaxy within a disc-bulge-halo galaxy that consists of some 10^6 particles.

We examine the bar instability in galactic models with an exponential disk and a cuspy dark matter (DM) halo with a Navarro-Frenk-White (NFW) cosmological density profile. We construct equilibrium models from a 3-integral composite distribution function that are subject to the bar instability. We generate a sequence of models with a range of mass resolution from 1.8K to 18M particles in the disk and 10K to 100M particles in the halo along with a multi-mass model with an effective resolution of ~10^10 particles. We describe how mass resolution affects the bar instability, including its linear growth phase, the buckling instability, pattern speed decay through the resonant transfer of angular momentum to the DM halo, and the possible destruction of the halo cusp. Our higher resolution simulations show a converging spectrum of discrete resonance interactions between the bar and DM halo orbits. As the pattern speed decays, orbital resonances sweep through most of the DM halo phase space and widely distribute angular momentum among the halo particles. The halo does not develop a flat density core and preserves the cusp, except in the region dominated by gravitational softening. The formation of the bar increases the central stellar density and the DM is compressed adiabatically increasing the halo central density by 1.7X. Overall, the evolution of the bar displays a convergent behavior for halo particle numbers between 1M and 10M particles, when comparing bar growth, pattern speed evolution, the DM halo density profile and a nonlinear analysis of the orbital resonances. Higher resolution simulations clearly illustrate the importance of discrete resonances in transporting the angular momentum from the bar to the halo.

We present some preliminary results from recent numerical simulations that model the evolution of super-massive black hole (SMBH) binaries in galactic nuclei. Including the post-Newtonian terms for the binary system and adopting appropriate models for the galaxies allows us, for the first time, to follow the evolution of SMBH binaries from kpc scales down to the coalescence phase. We use our results to make predictions of the detectability of such events with the gravitational wave detector LISA.

We examine the bar instability in models with an exponential disk and a cuspy NFW-like dark matter (DM) halo inspired by cosmological simulations. Bar evolution is studied as a function of numerical resolution in a sequence of models spanning 10K to 100M DM particles - including a multi-mass model with an effective resolution of 10G. The goal is to find convergence in dynamical behaviour. We characterize the bar growth, the buckling instability, pattern speed decay through resonant transfer of angular momentum, and possible destruction of the DM halo cusp. Overall, most characteristics converge in behaviour in detail for halos containing more than 10M particles. Notably, the formation of the bar does not destroy the density cusp in this case. These higher resolution simulations clearly illustrate the importance of discrete resonances in transporting angular momentum from the bar to the halo.

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

We study the growth of galactic disks in live triaxial DM halos. The halos have been assembled through constrained realizations method and evolved from the linear regime using cosmological simulations. The `seed' disks have been inserted at redshift z=3 and increased in mass tenfold over various time periods, ~1-3 Gyr, with the halo responding quasi-adiabatically to this process. We follow the dynamical and secular evolution of the disk-halo system and analyze changes in the most important parameters, like 3-D DM shapes, stellar and DM radial density profiles, stellar bar development, etc. We find that a growing disk is responsible for washing out the halo prolateness and for diluting its flatness over a period of time comparable to the disk growth. Moreover, we find that a disk which contributes more to the overall rotation curve in the system is also more efficient in axisymmetrizing the halo, without accelerating the halo figure rotation. The observational corollary is that the maximal disks probably reside in nearly axisymmetric halos, while disks whose rotation is dominated by the halo at all radii are expected to reside in more prolate halos. The halo shape is sensitive to the final disk mass, but is independent of how the seed disk is introduced into the system. We also expect that the massive disks are subject to a bar instability, while light disks have this instability damped by the halo triaxiality. Implications to the cosmological evolution of disks embedded in asymmetric halos are discussed and so are the corollaries for the observed fraction of stellar bars. Finally, the halo responds to the stellar bar by developing a gravitational wake -- a `ghost' bar of its own which is almost in-phase with that in the disk.

We analyze the evolution of stellar bars in galactic disks in mildly triaxial flat-core and cuspy CDM halos. We use tailored simulations of rigid and live halos which include the feedback from disk/bar onto the halo in order to test the work by El-Zant & Shlosman (2002). The latter used the Liapunov exponents to analyze the fate of bars in analytical asymmetric halos. We find: (1) The bar growth is similar in all rigid axisymmetric and triaxial halos. (2) Bars in live models vertically buckle and form a pseudobulge with a boxy/peanut shape. (3) In live axisymmetric halos, the bar strength varies little during the secular evolution. The bar pattern speed anticorrelates with the halo core size. The bar strength is larger for smaller disk-to-halo mass ratios within disk radii, the bar size correlates with the halo core sizes, and the bar pattern speeds -- with the halo central mass concentration. Bars embedded in live triaxial halos have a starkly different fate: they dissolve on ~1.5-5 Gyr due to the onset of chaos over continuous zones, leaving behind a weak oval distortion. The onset of chaos is related to the halo triaxiality, the fast rotating bar and the halo cuspiness. Before the bar dissolves, the region outside it develops strong spiral structures, especially in the live triaxial halos. (4) More angular momentum is absorbed by the triaxial halos as compared to the axisymmetric models and its exchange is mediated by resonances. (5) Cuspy halos are more susceptible than flat-core halos to having their prolateness washed out by the bar. We analyze these results in terms of the stability of trajectories and development of chaos. We set constraints on the triaxiality of DM halos by comparing our predictions to recent observations of bars out to z~1.

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

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

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

We carry out a detailed orbit analysis of gravitational potentials selected at different times from an evolving self-consistent model galaxy consisting of a two-component disk (stars+gas) and a live halo. The results are compared with a pure stellar model, subject to nearly identical initial conditions, which are chosen as to make the models develop a large scale stellar bar. The bars are also subject to hose-pipe (buckling) instability which modifies the vertical structure of the disk. The diverging morphological evolution of both models is explained in terms of gas radial inflow, the resulting change in the gravitational potential at smaller radii, and the subsequent modification of the main families of orbits, both in and out of the disk plane. We find that dynamical instabilities become milder in the presence of the gas component, and that the stability of planar and 3D stellar orbits is strongly affected by the related changes in the potential -- both are destabilized with the gas accumulation at the center. This is reflected in the overall lower amplitude of the bar mode and in the substantial weakening of the bar, which appears to be a gradual process. The vertical buckling of the bar is much less pronounced and the characteristic peanut shape of the galactic bulge almost disappears when there is a substantial gas inflow towards the center. Milder instability results in a smaller bulge whose basic parameters are in agreement with observations. We also find that the overall evolution in the model with a gas component is accelerated due to the larger central mass concentration and resulting decrease in the characteristic dynamical time.