The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 April 2013. Here we present the genesis, setup and data reduction of RAVE as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in RAVE DR6. Furthermore, we present derived spectral classification and radial velocities for the RAVE targets, complemented by cross matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km/s, while 95% of the objects have radial velocities better than 4.0 km/s. Stellar atmospheric parameters, abundances and distances are presented in subsequent publication. The data can be accessed via the RAVE Web (http://rave-survey.org) or the Vizier database.

We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their completion on 4 April 2013. In the second of two publications, we present the data products derived from 518387 observations of 451783 unique stars using a suite of advanced reduction pipelines focussing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (Teff, log(g), and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [alpha/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The RAVE DR6 catalogs are cross matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the RAVE Web site (http://rave-survey.org) or the Vizier database.

We study the distribution of the Milky Way satellites stellar and dark matter debris. For the first time we address the question of the tidal disruption of satellites in simulations by utilising simultaneously a) a realistic set of orbits extracted from cosmological simulations, b) a three component host galaxy with live halo, disc and bulge components, and c) satellites from hydrodynamical simulations. We analyse the statistical properties of the satellite debris of all massive galaxies reaching the inner Milky Way on a timescale of 2 Gyr. Up to 80$\%$ of the dark matter is stripped from the satellites, while this happens for up to 30$\%$ of their stars. The stellar debris ends mostly in the inner Milky Way halo, whereas the dark matter debris shows a flat mass distribution over the full main halo. The dark matter debris follows a density profile with inner power law index $\alpha_{\rm DM}=-0.66$ and outer index $\beta_{\rm DM}=2.94$, while for stars $\alpha_{*}=-0.44$ and $\beta_{*}=6.17$. In the inner 25 kpc, the distribution of the stellar debris is flatter than that of the dark matter debris and the orientations of their short axes differ significantly. Changing the orientation of the stellar disc by 90$^{\rm{o}}$ has only a minor impact on the distribution of the satellite debris. Our results indicate that the dark matter is more easily stripped than stars from the Milky Way satellites. The structure of the debris is dominated by the satellite orbital properties. The radial profiles, the flattening and the orientation of the stellar and dark matter debris are significantly different, which prevents the prediction of the dark matter distribution from the observed stellar component.

A new so-called `gravitational loss-cone instability' in stellar systems has recently been investigated theoretically in the framework of linear perturbation theory and proved to be potentially important in understanding the physical processes in centres of galaxies, star clusters, and the Oort comet cloud. Using N-body simulations, we confirm previous findings and go beyond the linear theory. Unlike the well-known instabilities, the new one shows no notable change in spherical geometry of the cluster, but it significantly accelerates the speed of diffusion of particles in phase space leading to a repopulation of the loss cone and early instability saturation.

A new so-called `gravitational loss-cone instability' in stellar systems has recently been investigated theoretically in the framework of linear perturbation theory and proved to be potentially important in understanding the physical processes in centres of galaxies, star clusters, and the Oort comet cloud. Using N-body simulations of a toy model, we confirm previous findings for the harmonic dominating potential and go beyond the linear theory. Unlike the well-known instabilities, the new one shows no notable change in the spherical geometry of the cluster, but it significantly accelerates the speed of diffusion of particles in phase space leading to an early instability saturation.

We study the evolution of star clusters in the Galactic tidal field starting from their birth in molecular clumps. Our model clusters form according to the local-density-driven cluster formation model in which the stellar density profile is steeper than that of gas. As a result, clusters resist the gas expulsion better than predicted by earlier models. We vary the impact of the Galactic tidal field {\lambda}, considering different Galactocentric distances (3-18 kpc), as well as different cluster sizes. Our model clusters survive the gas expulsion independent of {\lambda}. We investigated the relation between the cluster mass at the onset of secular evolution and their dissolution time. The model clusters formed with a high star-formation efficiency (SFE) follow a tight mass-dependent dissolution relation, in agreement with previous theoretical studies. However, the low-SFE models present a shallower mass-dependent relation than high-SFE clusters, and most dissolve before reaching 1 Gyr (cluster teenage mortality).

We derive the circular velocity curves of the gaseous and stellar discs of UGC 4056, a giant barred galaxy with an active galactic nucleus (AGN). We analyse UGC 4056 using the 2D spectroscopy obtained within the framework of the Mapping Nearby Galaxies at APO (MaNGA) survey. Using images and the colour index g-r from the Sloan Digital Sky Survey (SDSS), we determined the tilt of the galaxy, which allows us to conclude that the galaxy rotates clockwise with trailing spiral arms. We found that the gas motion at the central part of the UGC 4056 shows peculiar features. The rotation velocity of the gaseous disc shows a bump within around three kiloparsecs while the rotation velocity of the stellar disc falls smoothly to zero with decreasing galactocentric distance. We demonstrate that the peculiar radial velocities in the central part of the galaxy may be caused by the inflow of the gas towards the nucleus of the galaxy. The unusual motion of the gas takes place at the region with the AGN-like radiation and can be explained by the gas response to the bar potential.

We report the discovery of tidal structures around the intermediate-aged ($\sim$ 700--800~Myr), nearby ($\sim85$~pc) star cluster Coma Berenices. The spatial and kinematic grouping of stars is determined with the {\it Gaia} DR2 parallax and proper motion data, by a clustering analysis tool, \textsc{StarGO}, to map 5D parameters ($X, Y, Z$, $\mu_\alpha \cos\delta, \mu_\delta$) onto a 2D neural network. A leading and a trailing tails, each with an extension of $\sim50$~pc are revealed for the first time around this disrupting star cluster. The cluster members, totaling to $\sim115^{+5}_{-3}\,\rm {M_\odot}$, are clearly mass segregated, and exhibit a flat mass function with $\alpha \sim 0.79\pm0.16$, in the sense of $dN/dm \propto m^{-\alpha}$, where $N$ is the number of member stars and $m$ is stellar mass, in the mass range of $m=0.25$--$2.51~{\rm M_\odot}$. Within the tidal radius of $\sim$6.9~pc, there are 77 member candidates with an average position, i.e., as the cluster center, of R.A.= 186.8110~deg, and decl.= 25.8112~deg, and an average distance of 85.8~pc. Additional 120 member candidates reside in the tidal structures, i.e., outnumbering those in the cluster core. The expansion of escaping members lead to an anisotropy in the velocity field of the tidal tails. Our analysis also serendipitously uncovers an adjacent stellar group, part of which has been cataloged in the literature. We identify 218 member candidates, 10 times more than previously known. This star group is some 65~pc away from, and $\sim400$~Myr younger than, Coma Ber, but is already at the final stage of disruption.

We study the impact of the tidal field on the survivability of star clusters following instantaneous gas expulsion. Our model clusters are formed with a centrally-peaked star-formation efficiency profile as a result of star-formation taking place with a constant efficiency per free-fall time. We define the impact of the tidal field as the ratio of the cluster half-mass radius to its Jacobi radius immediately after gas expulsion, $\lambda = r_{h}/R_{J}$. We vary $\lambda$ by varying either the Galactocentric distance, or the size (hence volume density) of star clusters. We propose a new method to measure the violent relaxation duration, in which we compare the total mass-loss rate of star clusters with their stellar evolutionary mass-loss rate. That way, we can robustly estimate the bound mass fraction of our model clusters at the end of violent relaxation. The duration of violent relaxation correlates linearly with the Jacobi radius, when considering identical clusters at different Galactocentric distances. In contrast, it is nearly constant for the solar neighbourhood clusters, slightly decreasing with $\lambda$. The violent relaxation does not last longer than 50 Myr in our simulations. Identical model clusters placed at different Galactocentric distances have the same final bound fraction, despite experiencing different impacts of the tidal field. The solar neighbourhood clusters with different densities experience only limited variations of their final bound fraction. In general, we conclude that the cluster survivability after instantaneous gas expulsion, as measured by their bound mass fraction at the end of violent relaxation, $F_{bound}$, is independent of the impact of the tidal field, $\lambda$.

We examine the possible dependence of the radial oxygen abundance distribution on non-axisymmetrical structures (bar/spirals) and other macroscopic parameters such as the mass, the optical radius R25, the color g-r, and the surface brightness of the galaxy. A sample of disk galaxies from the CALIFA DR3 is considered. We adopted the Fourier amplitude A2 of the surface brightness as a quantitative characteristic of the strength of non-axisymmetric structures in a galactic disk, in addition to the commonly used morphologic division for A, AB, and B types based on the Hubble classification. To distinguish changes in local oxygen abundance caused by the non-axisymmetrical structures, the multiparametric mass--metallicity relation was constructed as a function of parameters such as the bar/spiral pattern strength, the disk size, color index g-r in the SDSS bands, and central surface brightness of the disk. The gas-phase oxygen abundance gradient is determined by using the R calibration. We find that there is no significant impact of the non-axisymmetric structures such as a bar and/or spiral patterns on the local oxygen abundance and radial oxygen abundance gradient of disk galaxies. Galaxies with higher mass, however, exhibit flatter oxygen abundance gradients in units of dex/kpc, but this effect is significantly less prominent for the oxygen abundance gradients in units of dex/R25 and almost disappears when the inner parts are avoided. We show that the oxygen abundance in the central part of the galaxy depends neither on the optical radius R25 nor on the color g-r or the surface brightness of the galaxy. Instead, outside the central part of the galaxy, the oxygen abundance increases with g-r value and central surface brightness of the disk.

We study the dynamics and evolution of the Milky Way nuclear star cluster performing a high resolution direct one-million-body simulation. Focusing on the interactions between such stellar system and the central supermassive black hole, we find that different stellar components adapt their overall distribution differently. After 5 Gyr, stellar mass black holes are characterized by a spatial distribution with power-slope $-1.75$, fully consistent with the prediction of Bahcall-Wolf pioneering work. Using the vast amount of data available, we infer the rate for tidal disruption events, being $4 \times 10^{-6}$ per yr, and estimate the number of objects that emit gravitational waves during the phases preceding the accretion onto the super-massive black hole, $\sim 270$ per Gyr. We show that some of these sources could form extreme mass-ratio inspirals. We follow the evolution of binary stars population, showing that the initial binary fraction of $5\%$ drops down to $2.5\%$ inside the inner parsec. Also, we explored the possible formation of binary systems containing a compact object, discussing the implications for millisecond pulsars formation and the development of Ia Supernovae.

We test the performance of the semi-analytic self-consistent Just-Jahrei{\ss} disc model (JJ model) with the astrometric data from the Tycho-Gaia Astrometric Solution (TGAS) sub-catalogue of the first Gaia data release (Gaia DR1), as well as the radial velocities from the fifth data release of the Radial Velocity Experiment survey (RAVE DR5). We use a sample of 19,746 thin disc stars from the TGAS$\times$RAVE cross-match selected in the local solar cylinder of 300 pc radius and 1 kpc height below the Galactic plane and simulate this sample via the forward modelling technique. First, we convert the predicted vertical density laws of the thin disc populations into a mock sample. Then the obtained mock populations are reddened with a 3D dust map and are subjected to the selection criteria corresponding to the RAVE and TGAS observational limitations as well as to additional cuts applied to the data sample. We calculate the quantities of interest separately at different heights above the Galactic plane taking into account the distance error effects separately in horizontal and vertical directions. We investigate the simulated sample in terms of the vertical number density profiles, Hess diagrams and velocity distribution functions. Basing on a good agreement of our simulations with the data, we conclude that our fiducial disc model confidently reproduces the vertical trends in the thin disc stellar population properties. Thus, it can serve as a starting point for the future extension of the JJ model to other Galactocentric distances.

We have studied the long-term evolution of star clusters of the solar neighborhood, starting from their birth in gaseous clumps until their complete dissolution in the Galactic tidal field. We have combined the "local-density-driven cluster formation model" of Parmentier & Pfalzner (2013) with direct N-body simulations of clusters following instantaneous gas expulsion. We have studied the relation between cluster dissolution time, $t_{dis}$, and cluster "initial" mass, $M_{init}$, defined as the cluster {mass at the end of the dynamical response to gas expulsion (i.e. violent relaxation), when the cluster age is 20-30 Myr}. We consider the "initial" mass to be consistent with other works which neglect violent relaxation. The model clusters formed with a high star formation efficiency (SFE -- i.e. gas mass fraction converted into stars) follow a tight mass-dependent relation, in agreement with previous theoretical studies. However, the low-SFE models present a large scatter in both the "initial" mass and the dissolution time, and a shallower mass-dependent relation than high-SFE clusters. Both groups differ in their structural properties on the average. Combining two populations of clusters, high- and low-SFE ones, with domination of the latter, yields a cluster dissolution time for the solar neighborhood in agreement with that inferred from observations, without any additional destructive processes such as giant molecular cloud encounters. An apparent mass-independent relation may emerge for our low-SFE clusters when we neglect low-mass clusters (as expected for extra-galactic observations), although more simulations are needed to investigate this aspect.

The evolution of Supermassive Black Holes (SMBHs) initially embedded in the centres of merging galaxies realised with a stellar mass function (SMF) is studied from the onset of galaxy mergers till coalescence. We performed a large set of direct N-body simulations with three different slopes of the central stellar cusp and different random seeds. Post Newtonian terms up to order 3.5 are used to drive the SMBH binary evolution in the relativistic regime. The impact of a SMF on the hardening rate and the coalescence time is investigated. We find that SMBH binaries coalesce well within one billion years when our models are scaled to galaxies with a steep cusp at low redshift. Here higher central densities provide larger supply of stars to efficiently extract energy from the SMBH binary orbit and shrink it to the phase where gravitational wave (GW) emission becomes dominant leading to the coalescence of the SMBHs. Mergers of models with shallow cusps that are representative for giant elliptical galaxies having central cores result in less efficient extraction of binary orbital energy due to the lower stellar densities in the centre. However, high values of eccentricity witnessed for SMBH binaries in such galaxy mergers ensure that the GW emission dominated phase sets in earlier at larger values of the semi-major axis. This helps to compensate for the less efficient energy extraction during the phase dominated by stellar encounters resulting in mergers of SMBHs in about one Gyr after the formation of the binary. Additionally, we witness mass segregation in the merger remnant resulting in enhanced SMBH binary hardening rates. We show that at least the final phase of the merger in cuspy low mass galaxies would be observable with the GW detector eLISA.

We construct the rotation curve of the Milky Way in the extended solar neighbourhood using a sample of SEGUE (Sloan Extension for Galactic Understanding and Exploration) G-dwarfs. We investigate the rotation curve shape for the presence of any peculiarities just outside the solar radius as has been reported by some authors. We approach the problem in a framework of classical Jeans analysis. Using the most recent data from RAVE (RAdial Velocity Experiment), we determine the solar peculiar velocity and the radial scalelengths for the three populations of different metallicities representing the Galactic thin disc. Then with the same binning in metallicity for the SEGUE G-dwarfs, we construct the rotation curve in the range of Galactocentric distances 7-10 kpc. We derive the circular velocity by correcting the mean tangential velocity for the asymmetric drift in each distance bin. With SEGUE data we also calculate the radial scalelength of the thick disc taking as known the derived peculiar motion of the Sun and the slope of the rotation curve. The rotation curve constructed through SEGUE G-dwarfs appears to be smooth in the selected radial range. The local kinematics of the thin disc rotation as determined in the framework of our new careful analysis does not favour the presence of a massive overdensity ring just outside the solar radius.

The all-sky Milky Way Star Clusters (MWSC) survey provides uniform and precise ages and other parameters for a variety of clusters in the Solar Neighbourhood. We construct the cluster age distribution, investigate its spatial variations, and discuss constraints on cluster formation scenarios of the Galactic disk during the last 5 Gyrs. Due to the spatial extent of the MWSC, we consider spatial variations of the age distribution along galactocentric radius $R_G$, and along $Z$-axis. For the analysis of the age distribution we use 2242 clusters, which all lie within roughly 2.5 kpc of the Sun. To connect the observed age distribution to the cluster formation history we build an analytical model based on simple assumptions on the cluster initial mass function and on the cluster mass-lifetime relation, fit it to the observations, and determine the parameters of the cluster formation law. Comparison with the literature shows that earlier results strongly underestimated the number of evolved clusters with ages $t\gtrsim 100$ Myr. Recent studies based on all-sky catalogues agree better with our data, but still lack the oldest clusters with ages $t\gtrsim 1$ Gyr. We do not observe a strong variation in the age distribution along $R_G$, though we find an enhanced fraction of older clusters ($t>1$ Gyr) in the inner disk. In contrast, the distribution strongly varies along $Z$. The high altitude distribution practically does not contain clusters with $t<1$ Gyr. With simple assumptions on the cluster formation history, cluster initial mass function and cluster lifetime we can reproduce the observations. Cluster formation rate and cluster lifetime are strongly degenerate, which does not allow us to disentangle different formation scenarios. In all cases the cluster formation rate is strongly declining with time, and the cluster initial mass function is very shallow at the high mass end. (abridged)

We use data from the Radial Velocity Experiment (RAVE) and the Tycho-Gaia astrometric solution catalogue (TGAS) to compute the velocity fields yielded by the radial (VR), azimuthal (Vphi) and vertical (Vz) components of associated Galactocentric velocity. We search in particular for variation in all three velocity components with distance above and below the disc midplane, as well as how each component of Vz (line-of-sight and tangential velocity projections) modifies the obtained vertical structure. To study the dependence of velocity on proper motion and distance we use two main samples: a RAVE sample including proper motions from the Tycho-2, PPMXL and UCAC4 catalogues, and a RAVE-TGAS sample with inferred distances and proper motions from the TGAS and UCAC5 catalogues. In both samples, we identify asymmetries in VR and Vz. Below the plane we find the largest radial gradient to be dVR / dR = -7.01+- 0.61 km\s kpc, in agreement with recent studies. Above the plane we find a similar gradient with dVR / dR= -9.42+- 1.77 km\s kpc. By comparing our results with previous studies, we find that the structure in Vz is strongly dependent on the adopted proper motions. Using the Galaxia Milky Way model, we demonstrate that distance uncertainties can create artificial wave-like patterns. In contrast to previous suggestions of a breathing mode seen in RAVE data, our results support a combination of bending and breathing modes, likely generated by a combination of external or internal and external mechanisms.

We perform high resolution direct $N$-body simulations to study the effect of an accretion disc on stellar dynamics in an active galactic nucleus (AGN). We show that the interaction of the nuclear stellar cluster (NSC) with the gaseous disc (AD) leads to formation of a stellar disc in the central part of the NSC. The accretion of stars from the stellar disc onto the super-massive black hole is balanced by the capture of stars from the NSC into the stellar disc, yielding a stationary density profile. We derive the migration time through the AD to be 3\% of the half-mass relaxation time of the NSC. The mass and size of the stellar disc are 0.7\% of the mass and 5\% of the influence radius of the super-massive black hole. An AD lifetime shorter than the migration time would result in a less massive nuclear stellar disc. The detection of such a stellar disc could point to past activity of the hosting galactic nucleus.

We derive maps of the observed velocity of ionized gas, the oxygen abundance, and the extinction (Balmer decrement) across the area of the four spiral galaxies NGC36, NGC180, NGC6063, and NGC7653 from integral field spectroscopy obtained by the Calar Alto Legacy Integral Field Area (CALIFA) survey. We searched for spiral arms through Fourier analysis of the spatial distribution of three tracers (non-circular motion, enhancement of the oxygen abundance, and of the extinction) in the discs of our target galaxies. The spiral arms (two-armed logarithmic spirals in the deprojected map) are shown in each target galaxy for each tracer considered. The pitch angles of the spiral arms in a given galaxy obtained with the three different tracers are close to each other. The enhancement of the oxygen abundance in the spiral arms as compared to the abundance in the interarm regions at a given galactocentric distance is small; within a few per cent. We identified a metallicity gradient in our target galaxies. Both barred galaxies in our sample show flatter gradients than the two galaxies without bars. Galactic inclination, position angle of the major axis, and the rotation curve were also obtained for each target galaxy using the Fourier analysis of the two-dimensional velocity map.

We study the effect of the instantaneous gas expulsion on star clusters wherein the residual gas has a density profile shallower than that of the embedded cluster. This is expected if star formation proceeds with a given SFE per free-fall time in a centrally-concentrated molecular clump. We perform direct N-body simulations whose initial conditions are generated by the program "mkhalo" "falcON" adapted for our models. Our model clusters initially have a Plummer profile and are in virial equilibrium with the gravitational potential of the cluster-forming clump. The residual gas contribution is computed based on the model of Parmentier&Pfalzner(2013). Our simulations include mass loss by stellar evolution and the tidal field of the Galaxy. We find that a star cluster with a minimum global SFE of 15% is able to survive instantaneous gas expulsion and to produce a bound cluster. Its violent relaxation lasts no longer than 20 Myr, independently of its global SFE and initial stellar mass. At the end of violent relaxation the bound fractions of surviving clusters with the same global SFEs are similar regardless of their initial stellar mass. Their subsequent lifetime in the gravitational field of the Galaxy depends on their bound stellar masses. We therefore conclude that the critical SFE needed to produce a bound cluster is 15%, which is twice smaller than earlier estimates of 33%. Thus we have improved the survival likelihood of young clusters after instantaneous gas expulsion. Those can now survive instantaneous gas expulsion with global SFEs as low as those observed for embedded clusters of Solar Neighbourhood (15-30%). This is the consequence of the star cluster having a density profile steeper than that of the residual gas. However, in terms of the effective SFE, measured by the virial ratio of the cluster at gas expulsion, our results are in agreement with previous studies.

Elemental abundances of stars are the result of the complex enrichment history of their galaxy. Interpretation of observed abundances requires flexible modeling tools to explore and quantify the information about Galactic chemical evolution (GCE) stored in such data. Here we present Chempy, a newly developed code for GCE modeling, representing a parametrized open one-zone model within a Bayesian framework. A Chempy model is specified by a set of 5-10 parameters that describe the effective galaxy evolution along with the stellar and star-formation physics: e.g. the star-formation history, the feedback efficiency, the stellar initial mass function (IMF) and the incidence of supernova type Ia (SN Ia). Unlike established approaches, Chempy can sample the posterior probability distribution in the full model parameter space and test data-model matches for different nucleosynthetic yield sets. We extend Chempy to a multi-zone scheme. As an illustrative application, we show that interesting parameter constraints result from only the ages and elemental abundances of Sun, Arcturus and the present-day interstellar medium (ISM). For the first time, we use such information to infer IMF parameter via GCE modeling, where we properly marginalize over nuisance parameters and account for different yield sets. We find that of the IMF $11.6_{-1.6}^{+2.1}$ % explodes as core-collapse SN, compatible with Salpeter 1955. We also constrain the incidence of SN Ia per 10^3 Msun to 0.5-1.4. At the same time, this Chempy application shows persistent discrepancies between predicted and observed abundances for some elements, irrespective of the chosen yield set. These cannot be remedied by any variations of Chempy's parameters and could be an indication for missing nucleosynthetic channels. Chempy should be a powerful tool to confront predictions from stellar nucleosynthesis with far more complex abundance data sets.

Aims. We aim to show how encounters with low-mass satellite galaxies may alter the bar formation in a Milky Way-like disc galaxy. Methods. We use high-resolution N-body simulations of a disc galaxy prone to mild bar instability. For realistic initial conditions of satellites, we take advantage of cosmological simulations of Milky Way-like dark matter haloes. Results. The satellites may have a significant impact on the time of bar formation. Some runs with satellites demonstrate a delay, while others show an advancement in bar formation compared to the isolated run, with such time differences reaching $\sim$ 1 Gyr. Meanwhile, the final bar configuration, including its very appearance and the bar characteristics such as the pattern speed and the exponential growth rate of its amplitude are independent of the number of encounters and their orbits. The contribution of satellites with masses below $10^9 M_{\odot}$ is insignificant, unless their pericentre distances are small. We suggest that the encounters act indirectly via inducing perturbations across the disc that evolve to delayed waves in the central part and interfere with an emerging seed bar. The predicted effect for the present-day host galaxy is expected to be even more significant at redshifts $z \gtrsim 0.5$.

Many barred galaxies, possibly including the Milky Way, have cusps in the centres. There is a widespread belief, however, that usual bar instability taking place in bulgeless galaxy models is impossible for the cuspy models, because of the presence of the inner Lindblad resonance for any pattern speed. At the same time there are numerical evidences that the bar instability can form a bar. We analyse this discrepancy, by accurate and diverse N-body simulations and using the calculation of normal modes. We show that bar formation in cuspy galaxies can be explained by taking into account the disc thickness. The exponential growth time is moderate for typical current disc masses (about 250 Myr), but considerably increases (factor 2 or more) upon substitution of the live halo and bulge with a rigid halo/bulge potential; meanwhile pattern speeds remain almost the same. Normal mode analysis with different disc mass favours a young bar hypothesis, according to which the bar instability saturated only recently.

Supermassive black holes (SMBHs) are ubiquitous in galaxies with a sizable mass. It is expected that a pair of SMBHs originally in the nuclei of two merging galaxies would form a binary and eventually coalesce via a burst of gravitational waves. So far theoretical models and simulations have been unable to predict directly the SMBH merger timescale from ab-initio galaxy formation theory, focusing only on limited phases of the orbital decay of SMBHs under idealized conditions of the galaxy hosts. The predicted SMBH merger timescales are long, of order Gyrs, which could be problematic for future gravitational wave searches. Here we present the first multi-scale $\Lambda$CDM cosmological simulation that follows the orbital decay of a pair of SMBHs in a merger of two typical massive galaxies at $z\sim3$, all the way to the final coalescence driven by gravitational wave emission. The two SMBHs, with masses $\sim10^{8}$ M$_{\odot}$, settle quickly in the nucleus of the merger remnant. The remnant is triaxial and extremely dense due to the dissipative nature of the merger and the intrinsic compactness of galaxies at high redshift. Such properties naturally allow a very efficient hardening of the SMBH binary. The SMBH merger occurs in only $\sim10$ Myr after the galactic cores have merged, which is two orders of magnitude smaller than the Hubble time.

We show by numerical simulations that a purely stellar dynamical model composed of an exponential disc, a cuspy bulge, and an NFW halo with parameters relevant to the Milky Way Galaxy is subject to bar formation. Taking into account the finite disc thickness, the bar formation can be explained by the usual bar instability, in spite of the presence of an inner Lindblad resonance, that is believed to damp any global modes. The effect of replacing the live halo and bulge by a fixed external axisymmetric potential (rigid models) is studied. It is shown that while the e-folding time of bar instability increases significantly (from 250 to 500 Myr), the bar pattern speed remains almost the same. For the latter, our average value of 55 km/s/kpc agrees with the assumption that the Hercules stream in the solar neighbourhood is an imprint of the bar--disc interaction at the outer Lindblad resonance of the bar. Vertical averaging of the radial force in the central disc region comparable to the characteristic scale length allows us to reproduce the bar pattern speed and the growth rate of the rigid models, using normal mode analysis of linear perturbation theory in a razor thin disc. The strong increase of the e-folding time with decreasing disc mass predicted by the mode analysis suggests that bars in galaxies similar to the Milky Way have formed only recently.

We present evidence from the RAdial Velocity Experiment (RAVE) survey of chemically separated, kinematically distinct disc components in the solar neighbourhood. We apply probabilistic chemical selection criteria to separate our sample into $\alpha$-low (`thin disc') and $\alpha$-high (`thick disc') sequences. Using newly derived distances, which will be utilized in the upcoming RAVE DR5, we explore the kinematic trends as a function of metallicity for each of the disc components. For our $\alpha$-low disc, we find a negative trend in the mean rotational velocity ($V_{\mathrm{\phi}}$) as a function of iron abundance ([Fe/H]). We measure a positive gradient $\partial V_{\mathrm{\phi}}$/$\partial$[Fe/H] for the $\alpha$-high disc, consistent with results from high-resolution surveys. We also find differences between the $\alpha$-low and $\alpha$-high discs in all three components of velocity dispersion. We discuss the implications of an $\alpha$-low, metal-rich population originating from the inner Galaxy, where the orbits of these stars have been significantly altered by radial mixing mechanisms in order to bring them into the solar neighbourhood. The probabilistic separation we propose can be extended to other data sets for which the accuracy in [$\alpha$/Fe] is not sufficient to disentangle the chemical disc components a priori. For such datasets which will also have significant overlap with Gaia DR1, we can therefore make full use of the improved parallax and proper motion data as it becomes available to investigate kinematic trends in these chemical disc components.

We examine the effect of an accretion disc on the orbits of stars in the central star cluster surrounding a central massive black hole by performing a suite of 39 high-accuracy direct N-body simulations using state-of-the art software and accelerator hardware, with particle numbers up to 128k. The primary focus is on the accretion rate of stars by the black hole (equivalent to their tidal disruption rate for black holes in the small to medium mass range) and the eccentricity distribution of these stars. Our simulations vary not only the particle number, but disc model (two models examined), spatial resolution at the centre (characterised by the numerical accretion radius) and softening length. The large parameter range and physically realistic modelling allow us for the first time to confidently extrapolate these results to real galactic centres. While in a real galactic centre both particle number and accretion radius differ by a few orders of magnitude from our models, which are constrained by numerical capability, we find that the stellar accretion rate converges for models with N > 32k. The eccentricity distribution of accreted stars, however, does not converge. We find that there are two competing effects at work when improving the resolution: larger particle number leads to a smaller fraction of stars accreted on nearly-circular orbits, while higher spatial resolution increases this fraction. We scale our simulations to some nearby galaxies and find that the expected boost in stellar accretion (or tidal disruption, which could be observed as X-ray flares) in the presence of a gas disc is about a factor of 10. Even with this boost, the accretion of mass from stars is still a factor of ~ 100 slower than the accretion of gas from the disc. Thus, it seems accretion of stars is not a major contributor to black hole mass growth.

We present a high resolution study of the impact of realistic satellite galaxies, extracted from cosmological simulations of Milky Way haloes including 6 Aquarius suites and Via Lactea \rom{2}, on the dynamics of the galactic disc. The initial conditions for the multi-component Milky Way galaxy were generated using the GalIC code, to ensure a system in dynamical equilibrium state prior to addition of satellites. Candidate subhaloes that came closer than 25\,kpc to the centre of the host DM haloes with initial mass enclosed within the tidal radius, $M_\textrm{tid}$ $\ge$ 10$^{8} M_{\odot}$\,=\,0.003 $M_\textrm{disc}$, were identified, inserted into our high resolution N-body simulations and evolved for 2 Gyr. We quantified the vertical heating due to such impacts by measuring the disc thickness and squared vertical velocity dispersion $\sigma_{z}^{2}$ across the disc. According to our analysis the strength of heating is strongly dependent on the high mass end of the subhalo distribution from cosmological simulations. The mean increase of the vertical dispersion is $\sim$ 20\,km$^{2}$\,s$^{-2}$\,Gyr$^{-1}$ for R $>$ 4\,kpc with a flat radial profile while, excluding Aq-F2 results, the mean heating is $<$ 12\,km$^{2}$\,s$^{-2}$\,Gyr$^{-1}$, corresponding to 28\% and 17\% of the observed vertical heating rate in the solar neighbourhood. Taking into account the statistical dispersion around the mean we miss the observed heating rate by more than 3$\sigma$. We observed a general flaring of the disc height in the case of all 7 simulations in the outer disc.

Our detailed analytic local disc model (JJ-model) quantifies the interrelation between kinematic properties (e.g. velocity dispersions and asymmetric drift), spatial parameters (scale-lengths and vertical density profiles), and properties of stellar sub-populations (age and abundance distributions). Any consistent radial extension of the disc evolution model should predict specific features in the different distribution functions and in their correlations. Large spectroscopic surveys (SEGUE, RAVE, APOGEE, Gaia-ESO) allow significant constraints on the long-term evolution of the thin disc. We discuss the qualitative difference of correlations (like the alpha-enhancement as function of metallicity) and distribution functions (e.g. in [Mg/H] or [Fe/H]) for the construction of a disc model. In the framework of the JJ-model we build a local chemical enrichment model and show that significant vertical gradients for main sequence and red clump stars are expected in the thin disc. A Jeans analysis of the asymmetric drift provides a link to the radial structure of the disc. The derived metallicity-dependent radial scale-lengths can be combined in the future with the abundance distributions at different Galactocentric distances to construct full disc models. We expect to be able to constrain possible scenarios of inside-out growth of the thin disc and to characterise those populations, which require significant radial migration.

From direct N-body simulations we find that the dynamical evolution of star clusters is strongly influenced by the Roche volume filling factor. We present a parameter study of the dissolution of open star clusters with different Roche volume filling factors and different particle numbers. We study both Roche volume underfilling and overfilling models and compare with the Roche volume filling case. We find that in the Roche volume overfilling limit of our simulations two-body relaxation is no longer the dominant dissolution mechanism but the changing cluster potential. We call this mechnism "mass-loss driven dissolution" in contrast to "two-body relaxation driven dissolution" which occurs in the Roche volume underfilling regime. We have measured scaling exponents of the dissolution time with the two-body relaxation time. In this experimental study we find a decreasing scaling exponent with increasing Roche volume filling factor. The evolution of the escaper number in the Roche volume overfilling limit can be described by a log-logistic differential equation. We report the finding of a resonance condition which may play a role for the evolution of star clusters and may be calibrated by the main periodic orbit in the large island of retrograde quasiperiodic orbits in the Poincar\'e surfaces of section. We also report on the existence of a stability curve which may be of relevance with respect to the structure of star clusters.

We investigate the influence of stellar migration caused by minor mergers (mass ratio from 1:70 to 1:8) on the radial distribution of chemical abundances in the disks of Milky Way-like galaxies during the last four Gyr. A GPU-based pure N-body tree-code model without hydrodynamics and star formation was used. We computed a large set of mergers with different initial satellite masses, positions, and orbital velocities. We find that there is no significant metallicity change at any radius of the primary galaxy in the case of accretion of a low-mass satellite of 10$^9$ M$_{\odot}$ (mass ratio 1:70) except for the special case of prograde satellite motion in the disk plane of the host galaxy. The accretion of a satellite of a mass $\gtrsim3\times10^9$ M$_{\odot}$ (mass ratio 1:23) results in an appreciable increase of the chemical abundances at galactocentric distances larger than $\sim10$ kpc. The radial abundance gradient flattens in the range of galactocentric distances from 5 to 15 kpc in the case of a merger with a satellite with a mass $\gtrsim3\times10^9$ M$_{\odot}$. There is no significant change in the abundance gradient slope in the outer disk (from $\sim15$ kpc up to 25 kpc) in any merger while the scatter in metallicities at a given radius significantly increases for most of the satellite's initial masses/positions compared to the case of an isolated galaxy. This argues against attributing the break (flattening) of the abundance gradient near the optical radius observed in the extended disks of Milky Way-like galaxies only to merger-induced stellar migration.

A new sample of stars, representative of the solar neighbourhood luminosity function, is constructed from the Hipparcos catalogue and the Fifth Catalogue of Nearby Stars. We have cross-matched to sources in the 2MASS catalogue so that for all stars individually determined Near Infrared photometry (NIR) is available on a homogeneous system (typically K_s). The spatial completeness of the sample has been carefully determined by statistical methods, and the NIR luminosity function of the stars has been derived by direct star counts. We find a local volume luminosity of 0.121 +/- 0.004 L_K_sun/(pc**3), corresponding to a volumetric mass-to-light ratio of M/L_K = 0.31 +/- 0.02 M_sun/L_K_sun, where giants contribute 80 per cent to the light but less than 2 per cent to the stellar mass. We derive the surface brightness of the solar cylinder with the help of a vertical disc model. We find a surface brightness of 99 L_K_sun/(pc**2) with an uncertainty of approximately 10 %. This corresponds to a mass-to-light ratio for the solar cylinder of M/L_K = 0.34 M_sun/L_K_sun. The mass-to-light ratio for the solar cylinder is only 10% larger than the local value despite the fact that the local population has a much larger contribution of young stars. It turns out that the effective scale heights of the lower main sequence carrying most of the mass is similar to that of the giants, which are dominating the NIR light. The corresponding colour for the solar cylinder is V-K=2.89 mag compared to the local value of V-K = 2.46 mag. An extrapolation of the local surface brightness to the whole Milky Way yields a total luminosity of M_K = -24.2 mag. The Milky Way falls in the range of K band Tully-Fisher (TF) relations from the literature.

We use our vertical Milky Way disc model together with Galaxia to create mock observations of stellar samples in the solar neighbourhood. We compare these to the corresponding volume complete observational samples of dereddened and binary accounted data from Hipparcos and the Catalogue of Nearby Stars. Sampling the likelihood in the parameter space we determine a new fiducial IMF considering constraints from dwarf and giant stars. The resulting IMF observationally backed in the range from 0.5 to 10 Msun is a two slope broken power law with -1.49 +- 0.08 for the low mass slope, a break at 1.39 +- 0.05 Msun and a high mass slope of -3.02 +- 0.06. The Besancon group also converging to a similar IMF even though their observational sample being quite different to ours shows that the forward modelling technique is a powerful diagnostic to test theoretical concepts like the local field star IMF.

We test methods for the determination of unstable modes in stellar discs: a point collocation scheme in the action sub-space, a scheme based on expansion of the density and potential on the biorthonormal basis, and a finite element method. Using models of galaxies with low and high mass concentration to the center, the existence of two different kinds of spectra of unstable modes is demonstrated. Characteristic features of methods and obtained spectra are discussed. Despite ignoring any constraint on the continuity or differentiability of the perturbed DF, the collocation scheme is reliable for obtaining spectra and patterns in both kinds of models. The method based on the expansion of the perturbed potential and surface density over a biorthogonal basis was not applicable to a model with high mass concentration. The finite element method successfully used in various fields of science and engineering is currently sensitive to the presence of resonant orbits due to the choice of interpolation functions for the distribution function.

We provide new constraints on the chemo-dynamical models of the Milky Way by measuring the radial and vertical chemical gradients for the elements Mg, Al, Si, Ti, and Fe in the Galactic disc and the gradient variations as a function of the distance from the Galactic plane ($Z$). We selected a sample of giant stars from the RAVE database using the gravity criterium 1.7$<$log g$<$2.8. We created a RAVE mock sample with the Galaxia code based on the Besan\c con model and selected a corresponding mock sample to compare the model with the observed data. We measured the radial gradients and the vertical gradients as a function of the distance from the Galactic plane $Z$ to study their variation across the Galactic disc. The RAVE sample exhibits a negative radial gradient of $d[Fe/H]/dR=-0.054$ dex kpc$^{-1}$ close to the Galactic plane ($|Z|<0.4$ kpc) that becomes flatter for larger $|Z|$. Other elements follow the same trend although with some variations from element to element. The mock sample has radial gradients in fair agreement with the observed data. The variation of the gradients with $Z$ shows that the Fe radial gradient of the RAVE sample has little change in the range $|Z|\lesssim0.6$ kpc and then flattens. The iron vertical gradient of the RAVE sample is slightly negative close to the Galactic plane and steepens with $|Z|$. The mock sample exhibits an iron vertical gradient that is always steeper than the RAVE sample. The mock sample also shows an excess of metal-poor stars in the [Fe/H] distributions with respect to the observed data. These discrepancies can be reduced by decreasing the number of thick disc stars and increasing their average metallicity in the Besan\c con model.

Context: Using the recent observational data of R\"oser et al. we present $N$-body simulations of the Hyades open cluster. Aims: We make an attempt to determine initial conditions of the Hyades cluster at the time of its formation in order to reproduce the present-day cumulative mass profile, stellar mass and luminosity function (LF). Methods: We performed direct $N$-body simulations of the Hyades in an analytic Milky Way potential that account for stellar evolution and include primordial binaries in a few models. Furthermore, we applied a Kroupa (2001) IMF and used extensive ensemble-averaging. Results: We find that evolved single-star King initial models with King parameters $W_0 = 6-9$ and initial particle numbers $N_0 = 3000$ provide good fits to the observational present-day cumulative mass profile within the Jacobi radius. The best-fit King model has an initial mass of $1721\ M_\odot$ and an average mass loss rate of $-2.2 \ M_\odot/\mathrm{Myr}$. The K-band LFs of models and observations show a reasonable agreement. Mass segregation is detected in both observations and models. If 33% primordial binaries are included the initial particle number is reduced by 5% as compared to the model without primordial binaries. Conclusions: The present-day properties of the Hyades can be well reproduced by a standard King or Plummer initial model when choosing appropriate initial conditions. The degeneracy of good-fitting models can be quite high due to the large dimension of the parameter space. More simulations with different Roche-lobe filling factors and primordial binary fractions are required to explore this degeneracy in more detail.

We present the stellar atmospheric parameters (effective temperature, surface gravity, overall metallicity), radial velocities, individual abundances and distances determined for 425 561 stars, which constitute the fourth public data release of the RAdial Velocity Experiment (RAVE). The stellar atmospheric parameters are computed using a new pipeline, based on the algorithms of MATISSE and DEGAS. The spectral degeneracies and the 2MASS photometric information are now better taken into consideration, improving the parameter determination compared to the previous RAVE data releases. The individual abundances for six elements (magnesium, aluminum, silicon, titanium, iron and nickel) are also given, based on a special-purpose pipeline which is also improved compared to that available for the RAVE DR3 and Chemical DR1 data releases. Together with photometric information and proper motions, these data can be retrieved from the RAVE collaboration website and the Vizier database.

Aim: We aim at measuring the chemical gradients of the elements Mg, Al, Si, and Fe along the Galactic radius to provide new constraints on the chemical evolution models of the Galaxy and Galaxy models such as the Besancon model. Methods: We analysed three different samples selected from three independent datasets: a sample of 19,962 dwarf stars selected from the RAVE database, a sample of 10,616 dwarf stars selected from the Geneva-Copenhagen Survey (GCS) dataset, and a mock sample (equivalent to the RAVE sample) created by using the GALAXIA code, which is based on the Besancon model. We measured the chemical gradients as functions of the guiding radius (Rg) at different distances from the Galactic plane reached by the stars along their orbit (Zmax). Results: The chemical gradients of the RAVE and GCS samples are negative and show consistent trends, although they are not equal: at Zmax<0.4 kpc and 4.5<Rg(kpc)<9.5, the iron gradient for the RAVE sample is d[Fe/H]/dRg=-0.065 dex kpc^{-1}, whereas for the GCS sample it is d[Fe/H]/dRg=-0.043 dex kpc^{-1} with internal errors +-0.002 and +-0.004 dex kpc^{-1}, respectively. The gradients of the RAVE and GCS samples become flatter at larger Zmax. Conversely, the mock sample has a positive iron gradient of d[Fe/H]/dRg=+0.053+-0.003 dex kpc^{-1} at Zmax<0.4 kpc and remains positive at any Zmax. These positive and unrealistic values originate from the lack of correlation between metallicity and tangential velocity in the Besancon model. The discrepancies between the observational samples and the mock sample can be reduced by i) decreasing the density, ii) decreasing the vertical velocity, and iii) increasing the metallicity of the thick disc in the Besancon model.

The cold dark matter model predicts cuspy dark matter halos. However, it has been found that, in some low-mass galaxies, cored dark halos provide a better description of their internal dynamics. Here we give constraints on the dark halo profile in the Sextans dwarf spheroidal galaxy by studying the longevity of two cold kinematic substructures detected in this galaxy. We perform N-body simulations of a stellar clump in the Sextans dwarf galaxy, including a live dark matter halo and the main stellar component. We find that, if the dark halo is cuspy, stellar clumps orbiting with semi-major axis ~400 pc are disrupted in ~5 Gyr, even if the clump is initially as compact stellar cluster with a radius of r_c=5 pc. Stellar clusters in an initial orbit with semi-major axis \leq 250 pc may survive to dissolution but their orbits decay towards the center by dynamical friction. In contrast, the stellar clumps can persist for a Hubble time within a cored dark matter halo, even if the initial clump's radius is as extended as r_c=80 pc. We also study the evolution of the clump in the MONDian context. In this scenario, we find that even an extended stellar clump with radius r_c=80 pc survives for a Hubble time, but an unrealistic value for the stellar mass-to-light ratio of 9.2 is needed.

The determination of the LSR is still a matter of debate. The classical value of the tangential peculiar motion of the Sun with respect to the LSR was challenged in recent years, claiming a significantly larger value. We show that the RAdial Velocity Experiment (RAVE) sample of dwarf stars is an excellent data set to derive tighter boundary conditions to chemodynamical evolution models of the extended solar neighbourhood. We present an improved Jeans analysis, which allows a better interpretation of the measured kinematics of stellar populations in the Milky Way disc. We propose an improved version of the Str\"omberg relation with the radial scalelengths as the only unknown. Binning RAVE stars in metallicity reveals a bigger asymmetric drift (corresponding to a smaller radial scalelength) for more metal-rich populations. With the standard assumption of velocity-dispersion independent radial scalelengths in each metallicity bin, we redetermine the LSR. The new Str\"omberg equation yields a joint LSR value of V_\sun=3.06 \pm 0.68 km/s, which is even smaller than the classical value based on Hipparcos data. The corresponding radial scalelength increases from 1.6 kpc for the metal-rich bin to 2.9 kpc for the metal-poor bin, with a trend of an even larger scalelength for young metal-poor stars. When adopting the recent Sch\"onrich value of V_\sun=12.24 km/s for the LSR, the new Str\"omberg equation yields much larger individual radial scalelengths of the RAVE subpopulations, which seem unphysical in part. The new Str\"omberg equation allows a cleaner interpretation of the kinematic data of disc stars in terms of radial scalelengths. Lifting the LSR value by a few km/s compared to the classical value results in strongly increased radial scalelengths with a trend of smaller values for larger velocity dispersions.

We examine the ratios $r_h/r_J$ of projected half-mass and Jacobi radius as well as $r_t/r_J$ of tidal and Jacobi radius for open and globular clusters in the Milky Way using data of both observations and simulations. We applied an improved calculation of $r_J$ for eccentric orbits of globular clusters. A sample of 236 open clusters of Piskunov et al. within the nearest kiloparsec around the Sun has been used. For the Milky Way globular clusters, data are taken from the Harris catalogue. We particularly use the subsample of 38 Milky Way globular clusters for which orbits have been integrated by Dinescu et al. We aim to quantify the differences between open and globular clusters and to understand, why they form two intrinsically distinct populations. We find under certain assumptions, or, in other words, in certain approximations, (i) that globular clusters are presently Roche volume underfilling and (ii) with at least $3\sigma$ confidence that the ratio $r_h/r_J$ of half-mass and Jacobi radius is $3 - 5$ times larger at present for an average open cluster in our sample than for an average globular cluster in our sample and (iii) that a significant fraction of globular clusters may be Roche volume overfilling at pericentre with $r_t > r_J$. Another aim of this paper is to throw light on the underlying theoretical reason for the existence of the van den Bergh correlation between half-mass and galactocentric radius.

SUPERBOX-10 is the successor of SUPERBOX, a particle-mesh code where additional grids and sub-grids are applied to regions of high particle density. Previous limitations have been solved. For instance, the vertical resolution is improved considerably when flattened grids are used. Since the computationally most intensive part is the Fast Fourier Transform, we introduce a parallelised version using the library FFTW, resulting in a speed-up of a few. The new features are tested using a galaxy model consisting of an exponential disc, a bulge and a dark matter halo. We demonstrate that the use of flattened grids efficiently reduces numerical heating. We simulate the merging of disc-bulge-halo galaxies with small spherical satellites. As a result, satellites on orbits with both low eccentricity and inclination heat the disc most efficiently. Moreover, we find that most of the satellite's energy and angular momentum is transfered to the halo.

We investigate the dynamical interaction of a central star cluster surrounding a super-massive black hole and a central accretion disk. The dissipative force acting on stars in the disk leads to an enhanced mass flow towards the super-massive black hole and to an asymmetry in the phase space distribution due to the rotating accretion disk. The accretion disk is considered as a stationary Keplerian rotating disk, which is vertically extended in order to employ a fully self-consistent treatment of stellar dynamics including the dissipative force originating from star-gas ram pressure effects. The stellar system is treated with a direct high-accuracy N-body integration code. A star-by-star representation, desirable in N-body simulations, cannot be extended to real particle numbers yet. Hence, we carefully discuss the scaling behavior of our model with regard to particle number and tidal accretion radius. The main idea is to find a family of models for which the ratio of two-body relaxation time and dissipation time (for kinetic energy of stellar orbits) is constant, which then allows us to extrapolate our results to real parameters of galactic nuclei. Our model is derived from basic physical principles and as such it provides insight into the role of physical processes in galactic nuclei, but it should be regarded as a first step towards more realistic and more comprehensive simulations. Nevertheless, the following conclusions appear to be robust: the star accretion rate onto the accretion disk and subsequently onto the super-massive black hole is enhanced by a significant factor compared to purely stellar dynamical systems neglecting the disk. This process leads to enhanced fueling of central disks in active galactic nuclei and to an enhanced rate of tidal stellar disruptions. [Abridged]

Through numerical simulations, we study the dissolution timescale of the Ursa Minor cold stellar clump, due to the combination of phase-mixing and gravitational encounters with compact dark substructures in the halo of Ursa Minor. We compare two scenarios; one where the dark halo is made up by a smooth mass distribution of light particles and one where the halo contains 10% of its mass in the form of substructures (subhalos). In a smooth halo, the stellar clump survives for a Hubble time provided that the dark matter halo has a big core. In contrast, when the point-mass dark substructures are added, the clump survives barely for \sim 1.5 Gyr. These results suggest a strong test to the \Lambda-cold dark matter scenario at dwarf galaxy scale.

We test the ability of the TRILEGAL and Besancon models to reproduce the CMD of SDSS data at the north Galactic pole (NGP). We show that a Hess diagram analysis of colour-magnitude diagrams is much more powerful than luminosity functions (LFs) in determining the Milky Way structure. We derive a best-fitting TRILEGAL model to simulate the NGP field in the (g-r, g) CMD of SDSS filters via Hess diagrams. For the Besancon model, we simulate the LFs and Hess diagrams in all SDSS filters. We use a chi2 analysis and determine the median of the relative deviations in the Hess diagrams to quantify the quality of the fits by the TRILEGAL models and the Besancon model in comparison and compare this with the Just-Jahreiss model. The input isochrones in the colour-absolute magnitude diagrams of the thick disc and halo are tested via the observed fiducial isochrones of globular clusters (GCs). We find that the default parameter set lacking a thick disc component gives the best representation of the LF in TRILEGAL. The Hess diagram reveals that a metal-poor thick disc is needed. In the Hess diagram, the median relative deviation of the TRILEGAL model and the SDSS data amounts to 25 percent, whereas for the Just-Jahreiss model the deviation is only 5.6 percent. The isochrone analysis shows that the representation of the MS of (at least metal-poor) stellar populations in the SDSS system is reliable. In contrast, the RGBs fail to match the observed fiducial sequences of GCs. The Besancon model shows a similar median relative deviation of 26 percent in (g-r, g). In the u band, the deviations are larger. There are significant offsets between the isochrone set used in the Besancon model and the observed fiducial isochrones. In contrast to Hess diagrams, LFs are insensitive to the detailed structure of the Milky Way components due to the extended spatial distribution along the line of sight.

Using the RAVE survey, we recently brought to light a gradient in the mean galactocentric radial velocity of stars in the extended solar neighbourhood. This gradient likely originates from non-axisymmetric perturbations of the potential, among which a perturbation by spiral arms is a possible explanation. Here, we apply the traditional density wave theory and analytically model the radial component of the two-dimensional velocity field. Provided that the radial velocity gradient is caused by relatively long-lived spiral arms that can affect stars substantially above the plane, this analytic model provides new independent estimates for the parameters of the Milky Way spiral structure. Our analysis favours a two-armed perturbation with the Sun close to the inner ultra-harmonic 4:1 resonance, with a pattern speed \Omega_p=18.6^{+0.3}_{-0.2} km/s/kpc and a small amplitude A=0.55 \pm 0.02% of the background potential (14% of the background density). This model can serve as a basis for numerical simulations in three dimensions, additionally including a possible influence of the galactic bar and/or other non-axisymmetric modes.

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.

We employ measurements of the [alpha/Fe] ratio derived from low-resolution (R~2000) spectra of 17,277 G-type dwarfs from the SEGUE survey to separate them into likely thin- and thick-disk subsamples. Both subsamples exhibit strong gradients of orbital rotational velocity with metallicity, of opposite signs, -20 to -30 km/s/dex for the thin-disk and +40 to +50 km/s/dex for the thick-disk population. The rotational velocity is uncorrelated with Galactocentric distance for the thin-disk subsample, and exhibits a small trend for the thick-disk subsample. The rotational velocity decreases with distance from the plane for both disk components, with similar slopes (-9.0 {\pm} 1.0 km/s/kpc). Thick-disk stars exhibit a strong trend of orbital eccentricity with metallicity (about -0.2/dex), while the eccentricity does not change with metallicity for the thin-disk subsample. The eccentricity is almost independent of Galactocentric radius for the thin-disk population, while a marginal gradient of the eccentricity with radius exists for the thick-disk population. Both subsamples possess similar positive gradients of eccentricity with distance from the Galactic plane. The shapes of the eccentricity distributions for the thin- and thick-disk populations are independent of distance from the plane, and include no significant numbers of stars with eccentricity above 0.6. Among several contemporary models of disk evolution we consider, radial migration appears to have played an important role in the evolution of the thin-disk population, but possibly less so for the thick disk, relative to the gas-rich merger or disk heating scenarios. We emphasize that more physically realistic models and simulations need to be constructed in order to carry out the detailed quantitative comparisons that our new data enable.

In spherical galaxies, binary supermassive black holes (SMBHs) have difficulty reaching sub-parsec separations due to depletion of stars on orbits that intersect the massive binary - the final-parsec problem. Galaxies that form via major mergers are substantially nonspherical, and it has been argued that the centrophilic orbits in triaxial galaxies might provide stars to the massive binary at a high enough rate to avoid stalling. Here we test that idea by carrying out fully self-consistent merger simulations of galaxies containing central SMBHs. We find hardening rates of the massive binaries that are indeed much higher than in spherical models, and essentially independent of the number of particles used in the simulations. Binary eccentricities remain high throughout the simulations. Our results constitute a fully stellar-dynamical solution to the final-parsec problem and imply a potentially high rate of events for low-frequency gravitational wave detectors like LISA.

We present an evolutionary disc model for the thin disc in the solar cylinder based on a continuous star formation history and a continuous dynamical heating of the stellar subpopulations. The vertical distribution of the stellar subpopulations are calculated self-consistently in dynamical equilibrium. The SFR and AVR of the stellar subpopulations are determined by fitting the velocity distribution functions of main sequence stars. The SFR shows a maximum 10 Gyr ago declining by a factor of 10 until present time. A constant SFR can be ruled out. The velocity dispersion of the stellar subpopulations increase with age according to a power law with index 0.375. The new scale heights lead to a best fit IMF with power-law indices of 1.5 below and 4.0 above 1.6 Msun, which has no kink around 1 Msun. Including a thick disc component results in slight variations of the thin disc properties, but has a negligible influence on the SFR. A variety of predictions are made concerning the number density, age and metallicity distributions of stellar subpopulations as a function of z above the galactic plane.

We have used the self-consistent vertical disc models of the solar neighbourhood presented in Just & Jahreiss (2010), which are based on different star formation histories (SFR) and fit the local kinematics of main sequence stars equally well, to predict star counts towards the North Galactic Pole (NGP). We combined these four different models with the local main sequence in the filter system of the SDSS and predicted the star counts in the NGP field with b>80deg. All models fit the Hess diagrams in the F-K dwarf regime better than 20 percent and the star number densities in the solar neighbourhood are consistent with the observed values. The chi^2 analysis shows that model A is clearly preferred with systematic deviations of a few percent only. The SFR of model A is characterised by a maximum at an age of 10Gyr and a decline by a factor of four to the present day value of 1.4Msun/pc^2/Gyr. The thick disc can be modelled very well by an old isothermal simple stellar population. The density profile can be approximated by a sech^(alpha_t) function. We found a power law index alpha_t=1.16 and a scale height of 800pc corresponding to a vertical velocity dispersion of 45.3km/s. About 6 percent of the stars in the solar neighbourhood are thick disc stars.

Dynamical friction leads to an orbital decay of massive objects like young compact star clusters or Massive Black Holes in central regions of galaxies. The dynamical friction force can be well approximated by Chandrasekhar's standard formula, but recent investigations show, that corrections to the Coulomb logarithm are necessary. With a large set of N-body simulations we show that the improved formula for the Coulomb logarithm fits the orbital decay very well for circular and eccentric orbits. The local scale-length of the background density distribution serves as the maximum impact parameter for a wide range of power-law indices of -1 ... -5. For each type of code the numerical resolution must be compared to the effective minimum impact parameter in order to determine the Coulomb logarithm. We also quantify the correction factors by using self-consistent velocity distribution functions instead of the standard Maxwellian often used. These factors enter directly the decay timescale and cover a range of 0.5 ... 3 for typical orbits. The new Coulomb logarithm combined with self-consistent velocity distribution functions in the Chandrasekhar formula provides a significant improvement of orbital decay times with correction up to one order of magnitude compared to the standard case. We suggest the general use of the improved formula in parameter studies as well as in special applications.

Context: A recent new approach to apply a simple dynamical mass estimate of tidally limited star clusters is based on the identification of the tidal radius in a King profile with the dynamical Jacobi radius. The application to an unbiased open cluster catalogue yields significantly higher cluster masses compared to the classical methods. Aims: We quantify the bias in the mass determination as function of projection direction and cluster age by analysing a simulated star cluster. Methods: We use direct $N$-body simulations of a star cluster including stellar evolution in an analytic Milky Way potential and apply a best fit to the projected number density of cluster stars. Results: We obtain significantly overestimated star cluster masses which depend strongly on the viewing direction. The overestimation is typically in the range of 10-50 percent and reaches a factor of 3.5 for young clusters. Mass segregation reduces the derived limiting radii systematically.

We present a fully consistent evolutionary disc model of the solar cylinder. The model is based on a sequence of stellar sub-populations described by the star formation history (SFR) and the dynamical heating law (given by the age-velocity dispersion relation AVR). The combination of kinematic data from Hipparcos and the finite lifetimes of main sequence (MS) stars enables us to determine the detailed vertical disc structure independent of individual stellar ages and only weakly dependent on the IMF. The disc parameters are determined by applying a sophisticated best fit algorithm to the MS star velocity distribution functions in magnitude bins. We find that the AVR is well constrained by the local kinematics, whereas for the SFR the allowed range is larger. A simple chemical enrichment model is included in order to fit the local metallicity distribution of G dwarfs. In our favoured model A the power law index of the AVR is 0.375 with a minimum and maximum velocity dispersion of 5.1 km/s and 25.0 km/s, respectively. The SFR shows a maximum 10 Gyr ago and declines by a factor of four to the present day value of 1.5 M_sun/pc^2/Gyr. A best fit of the IMF leads to power-law indices of -1.46 below and -4.16 above 1.72 M_sun avoiding a kink at 1 M_sun. An isothermal thick disc component with local density of ~6% of the stellar density is included. A thick disc containing more than 10% of local stellar mass is inconsistent with the local kinematics of K and M dwarfs.

We present N-body simulations of dissolving star clusters close to galactic centres. For this purpose, we developed a new N-body program called nbody6gc based on Aarseth's series of N-body codes. We describe the algorithm in detail. We report about the density wave phenomenon in the tidal arms which has been recently explained by Kuepper et al. (2008). Standing waves develop in the tidal arms. The wave knots or clumps develop at the position, where the emerging tidal arm hits the potential wall of the effective potential and is reflected. The escaping stars move through the wave knots further into the tidal arms. We show the consistency of the positions of the wave knots with the theory in Just et al. (2009). We also demonstrate a simple method to study the properties of tidal arms. By solving many eigenvalue problems along the tidal arms, we construct numerically a 1D coordinate system whose direction is always along a principal axis of the local tensor of inertia. Along this coordinate system, physical quantities can be evaluated. The half-mass or dissolution times of our models are almost independent of the particle number which indicates that two-body relaxation is not the dominant mechanism leading to the dissolution. This may be a typical situation for many young star clusters. We propose a classification scheme which sheds light on the dissolution mechanism.

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.

Tidal tails of star clusters are not homogeneous but show well defined clumps in observations as well as in numerical simulations. Recently an epicyclic theory for the formation of these clumps was presented. A quantitative analysis was still missing. We present a quantitative derivation of the angular momentum and energy distribution of escaping stars from a star cluster in the tidal field of the Milky Way and derive the connection to the position and width of the clumps. For the numerical realization we use star-by-star $N$-body simulations. We find a very good agreement of theory and models. We show that the radial offset of the tidal arms scales with the tidal radius, which is a function of cluster mass and the rotation curve at the cluster orbit. The mean radial offset is 2.77 times the tidal radius in the outer disc. Near the Galactic centre the circumstances are more complicated, but to lowest order the theory still applies. We have also measured the Jacobi energy distribution of bound stars and showed that there is a large fraction of stars (about 35%) above the critical Jacobi energy at all times, which can potentially leave the cluster. This is a hint that the mass loss is dominated by a self-regulating process of increasing Jacobi energy due to the weakening of the potential well of the star cluster, which is induced by the mass loss itself.

For an understanding of Galactic stellar populations in the SDSS filter system well defined stellar samples are needed. The nearby stars provide a complete stellar sample representative for the thin disc population. We compare the filter transformations of different authors applied to the main sequence stars from F to K dwarfs to SDSS filter system and discuss the properties of the main sequence. The location of the mean main sequence in colour-magnitude diagrams is very sensitive to systematic differences in the filter transformation. A comparison with fiducial sequences of star clusters observed in g',r',i' show good agreement. Theoretical isochrones from Padua and from Dartmouth have still some problems especially in (r-i)-colour.

The dissolution process of star clusters is rather intricate for theory. We investigate it in the context of chaotic dynamics. We use the simple Plummer model for the gravitational field of a star cluster and treat the tidal field of the Galaxy within the tidal approximation. That is, a linear approximation of tidal forces from the Galaxy based on epicyclic theory in a rotating reference frame. The Poincar\'e surfaces of section reveal the effect of a Coriolis asymmetry. The system is non-hyperbolic which has important consequences for the dynamics. We calculated the basins of escape with respect to the Lagrangian points $L_1$ and $L_2$. The longest escape times have been measured for initial conditions in the vicinity of the fractal basin boundaries. Furthermore, we computed the chaotic saddle for the system and its stable and unstable manifolds. The chaotic saddle is a fractal structure in phase space which has the form of a Cantor set and introduces chaos into the system.

We present N-body simulations of unstable spiral modes in a dynamically cool collisionless disc. We show that spiral modes grow in a thin collisionless disk in accordance with the analytical perturbation theory. We use the particle-mesh code SUPERBOX with nested grids to follow the evolution of unstable spirals that emerge from an unstable equilibrium state. We use a large number of particles (up to 40 million particles) and high-resolution spatial grids in our simulations (128^3 cells). These allow us to trace the dynamics of the unstable spiral modes until their wave amplitudes are saturated due to nonlinear effects. In general, the results of our simulations are in agreement with the analytical predictions. The growth rate and the pattern speed of the most unstable bar-mode measured in N-body simulations agree with the linear analysis. However the parameters of secondary unstable modes are in lesser agreement because of the still limited resolution of our simulations.

We have studied the dynamical evolution of rotating globular clusters with direct $N$-body models. Our initial models are rotating King models; we obtained results for both equal-mass systems and systems composed out of two mass components. Previous investigations using a Fokker-Planck solver have revealed that rotation has a noticeable influence on stellar systems like globular clusters, which evolve by two-body relaxation. In particular, it accelerates their dynamical evolution through the gravogyro instability. We have validated the occurence of the gravogyro instability with direct $N$-body models. In the case of systems composed out of two mass components, mass segregation takes place, which competes with the rotation in the acceleration of the core collapse. The "accelerating" effect of rotation has not been detected in our isolated two-mass $N$-body models. Last, but not least, we have looked at rotating $N$-body models in a tidal field within the tidal approximation. It turns out that rotation increases the escape rate significantly. A difference between retrograde and prograde rotating star clusters occurs with respect to the orbit of the star cluster around the Galaxy, which is due to the presence of a ``third integral'' and chaotic scattering, respectively.

We present a physical model that explains the two disparate observational facts: 1) the exponential vertical disc structure in the optical and NIR of the non-obscured part of the stellar disc and 2) the enhanced FIR/submm luminosity by about a factor of four near the obscured mid-plane, which requires additional dust and also stellar light to heat the dust component. We use multi-band photometry in U, B, V, R, and I- band combined with radiative transfer through a dust component to fit simultaneously the vertical surface-brightness and colour index profiles in all bands adopting a reasonable star formation history and dynamical heating function. The final disc model reproduces the surface-brightness profiles in all bands with a moderately declining star formation rate and a slowly starting heating function for young stars. The total dust mass is 57 million solar masses as required from the FIR/submm measurements. Without a recent star burst we find in the midplane an excess of 5.2-, 4.0-, and 3.0-times more stellar light in the U-, B-, and V-band, respectively. The corresponding stellar mass-to-light ratios are 0.91 in V- and 1.0 in R-band. The central face-on optical depth in V-band is 0.81 and the radial scale length of the dust is 40% larger than that of the stellar disc. Evolutionary disc models are a powerful method to understand the vertical structure of edge-on galaxies. Insights to the star formation history and the dynamical evolution of stellar discs can be gained. FIR/submm observations are necessary to restrict the parameter space for the models.

We investigate the effect of a density gradient on Chandrasekhar's dynamical friction formula based on the method of 2-body encounters in the local approximation. We apply these generalizations to the orbit evolution of satellite galaxies in Dark Matter haloes. We find from the analysis that the main influence occurs through a position-dependent maximum impact parameter in the Coulomb logarithm, which is determined by the local scale-length of the density distribution. We also show that for eccentric orbits the explicit dependence of the Coulomb logarithm on position yields significant differences for the standard homogeneous force. Including the velocity dependence of the Coulomb logarithm yields ambigous results. The orbital fits in the first few periods are further improved, but the deviations at later times are much larger. The additional force induced by the density gradient, the inhomogeneous force, is not antiparallel to the satellite motion and can exceed 10% of the homogeneous friction force in magnitude. However, due to the symmetry properties of the inhomogeneous force, there is a deformation and no secular effect on the orbit at the first order. Therefore the inhomogeneous force can be safely neglected for the orbital evolution of satellite galaxies. For the homogeneous force we compare numerical N-body calculations with semi-analytical orbits to determine quantitatively the accuracy of the generalized formulae of the Coulomb logarithm in the Chandresekhar approach. With the local scale-length as the maximum impact parameter we find a significant improvement of the orbital fits and a better interpretation of the quantitative value of the Coulomb logarithm.

First we present a recently developed 3D chemodynamical code for galaxy evolution from the K**2 collaboration. It follows the evolution of all components of a galaxy such as dark matter, stars, molecular clouds and diffuse interstellar matter (ISM). Dark matter and stars are treated as collisionless N-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse (hot) gas and a sticky particle scheme for the (cool) molecular clouds. Physical processs such as star formation, stellar death or condensation and evaporation processes of clouds interacting with the ISM are described locally. An example application of the model to a star forming dwarf galaxy will be shown for comparison with other codes. Secondly we will discuss new kinds of exotic chemodynamical processes, as they occur in dense gas-star systems in galactic nuclei, such as non-standard ``drag''-force interactions, destructive and gas producing stellar collisions. Their implementation in 1D dynamical models of galactic nuclei is presented. Future prospects to generalize these to 3D are work in progress and will be discussed.

We carry out a set of self-consistent N-body calculations to investigate how important the velocity anisotropy in non-spherical dark-matter haloes is for dynamical friction. For this purpose we allow satellite galaxies to orbit within flattened and live dark-matter haloes (DMHs) with aspect-ratio q_h=0.6 and compare the resulting orbit evolution with a semi-analytic code. This code solves the equation of motion of satellites orbiting within the same DMH and implements, according to the dynamical friction model, either Chandrasekhar's formula, which does not incorporate the velocity anisotropy, or Binney's description of dynamical friction in anisotropic systems. In the numerical and the two semi-analytic models the satellites are given different initial orbital inclinations and orbital eccentricities. We find that Binney's approach successfully describes the overall satellite decay and orbital inclination decrease for the whole set of orbits, with an averaged discrepancy of less than 4 per cent in orbital radius during the first 3 orbits. If Chandrasekhar's expression is used instead, the discrepancy increases to 20 per cent and no inclination decrease can be observed. Binney's treatment therefore appears to provide a significantly improved treatment of dynamical friction in anisotropic systems.

Highly condensed gaseous objects with masses larger than 5x10^4 M_sun are called super-massive stars. In the quasistationary contraction phase, the hydrostatic equilibrium is determined by radiation pressure and gravitation. The global structure is that of an n=3 polytrope at the stability limit. Small relativistic corrections for example can initiate a free fall collapse due to the 'post Newtonian' instability. Since the outcome of the final collapse -A super-massive black hole or hypernova- depends sensitively on the structure and the size of the object, when the instability sets in, it is important to investigate in more detail the contraction phase of the SMS. If the gaseous object is embedded in a dense stellar system, the central star cluster, the interaction and coupling of both components due to dynamical friction changes the energy balance and evolution of the SMS dramatically. Dynamical friction between stars and gas, which can be estimated semi-analytically (see Just et al. 1986), has three different effects on the two-component system. We discuss in which evolutionary stages and parameter range these interaction processes are relevant and how they can influence the stability and evolution of the SMS.

The concept of central super-massive stars (${\cal M} \ge 5 \times 10^4 M_{\odot}$, where ${\cal M}$ is the mass of the super-massive star) embedded in dense stellar systems was suggested as a possible explanation for high- energy emissions phenomena occurring in active galactic nuclei and quasars (Vilkoviski 1976, Hara 1978), such as X-ray emissions (Bahcall and Ostriker, 1975). SMSs and super-massive black holes are two possibilities to explain the nature of super-massive central objects, and super-massive stars may be an intermediate step towards the formation of super-massive black holes (Rees 1984). Therefore it is important to study such a dense gas-star system in detail. We address here the implementation of radiative transfer in a model which was presented in former work (Amaro-Seoane and Spurzem 2001, Amaro-Seoane et al. 2002). In this sense, we extend here and improve the work done by Langbein et al. (1990) by describing the radiative transfer in super-massive stars using previous work on this subject (Castor 1972).

We use a gaseous model and a semi-analytical approach to study the evolution of a super-massive central gaseous object (a super-massive star, (SMS from now on) in an AGN and its evolution by interactions with the surrounding stellar system. Our future work in this field is outlined, which aims at a more detailed study of energy flows in the interstellar medium, stellar evolution and the relation between QSOs and galaxy formation.