BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andreas Just

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.SRastro-ph.IMastro-ph.HEastro-phastro-ph.EP

Catalog footprint

What is connected

23works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

23 published item(s)

preprint2024arXiv

Identification of new nearby white dwarfs using Gaia DR3

Based on the astrometry and photometry in Gaia DR3, we identified new nearby white dwarfs and validated those that had been missed from recent white dwarf catalogues despite being previously documented. To ensure the reliability of their astrometric solutions, we used a cut on just two parameters from Gaia DR3: the amplitude of the image parameter determination goodness-of-fit and the parallax-over-error ratio. In addition, we imposed photometric signal-to-noise requirements to ensure the reliable identification of white dwarfs when using the colour-magnitude diagram. We have identified nine previously unreported white dwarfs within the local population of 50 pc, and validated 21 previously reported white dwarfs missing from the GCWD21 (Gentile Fusillo et al. 2021) and other recent volume-limited white dwarf samples. A few of these objects belong to the rare class of ultra-cool white dwarfs. Four white dwarfs in our sample have an effective temperature of $T_{eff}\leq4000$ K within the $1σ$ interval, and two of them have an absolute magnitude of $M_G > 16.0$ mag. The identified white dwarfs are predominantly located in crowded fields, such as near the Galactic plane or in the foreground of the Large Magellanic Cloud. We also find that 19 of these white dwarfs have common proper motion companions with angular separations ranging from $1.1''$ to $7.1''$ and brightness differences between the components of up to 9.8 magnitudes. One of these systems is a triple system consisting of a white dwarf and two K dwarfs, while another is a double white dwarf system. We have identified 103 contaminants among the 2338 high-confidence white dwarfs in the 50 pc subsample of the GCWD21 and have found that their astrometric solutions in Gaia DR3 are spurious, improving the purity by 4.4%.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Addressing via N-body simulations the distribution of the satellite tidal debris in the Milky Way environment

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 $α_{\rm DM}=-0.66$ and outer index $β_{\rm DM}=2.94$, while for stars $α_{*}=-0.44$ and $β_{*}=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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- I: Survey Description, Spectra and Radial Velocities

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- II: Stellar Atmospheric Parameters, Chemical Abundances and Distances

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Star Clusters in the Galactic tidal field, from birth to dissolution

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 λ, considering different Galactocentric distances (3-18 kpc), as well as different cluster sizes. Our model clusters survive the gas expulsion independent of λ. 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).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Chemical separation of disc components using RAVE

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 $α$-low (`thin disc') and $α$-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 $α$-low disc, we find a negative trend in the mean rotational velocity ($V_{\mathrmϕ}$) as a function of iron abundance ([Fe/H]). We measure a positive gradient $\partial V_{\mathrmϕ}$/$\partial$[Fe/H] for the $α$-high disc, consistent with results from high-resolution surveys. We also find differences between the $α$-low and $α$-high discs in all three components of velocity dispersion. We discuss the implications of an $α$-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 [$α$/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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the bar formation mechanism in galaxies with cuspy bulges

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Star-disc interaction in galactic nuclei: orbits and rates of accreted stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Swift coalescence of supermassive black holes in cosmological mergers of massive galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamical and chemical evolution of the thin 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Influence of Minor Mergers on the Radial Abundance Gradient in Disks of Milky Way-like Galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Towards a fully consistent Milky Way disc model - III. Constraining the initial mass function

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bisymmetric normal modes in soft-centred and realistic galactic discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Detailed comparison of Milky Way models based on stellar population synthesis and SDSS star counts at the north Galactic pole

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Enhanced accretion rates of stars on Super-massive Black Holes by star-disk interactions in galactic nuclei

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]

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Formation and Hardening of Supermassive Black Hole Binaries in Minor Mergers of Disk Galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mergers of Unequal Mass Galaxies: Supermassive Black Hole Binary Evolution and Structure of Merger Remnants

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulating Sinking Satellites with SUPERBOX-10

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Efficient Merger of Binary Supermassive Black Holes in Merging Galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Formation and Evolution of the Disk System of the Milky Way: [alpha/Fe] Ratios and Kinematics of the SEGUE G-Dwarf Sample

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Towards a fully consistent Milky Way disc model - II. The local disc model and SDSS data of the NGP region

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the dissolution of star clusters in the Galactic centre. I. Circular orbits

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Escape from the vicinity of fractal basin boundaries of a star cluster

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

0 citations0 reviewsNo code linkedOpen access