Source author record

A. Just

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

19 published item(s)

preprint2021arXiv

Towards a fully consistent Milky Way disk model -- IV. The impact of Gaia DR2 and APOGEE

We present an updated version of the semi-analytic Just-Jahreiß (JJ) model of the Galactic disk and constrain its parameters in the Solar neighbourhood. The new features of the JJ model include a simple two-component gaseous disk, a star-formation rate (SFR) function of the thick disk that has been extended in time, and a correlation between the kinematics of molecular gas and thin-disk populations. Here, we study the vertical number density profiles and W-velocity distributions determined from ~2 million local stars of the Gaia DR2. We also investigate an apparent Hess diagram of the Gaia stars selected in a conic volume towards the Galactic poles. Using a stellar evolution library, we synthesise stellar populations with a four-slope broken power-law initial mass function (IMF), the SFR, and an age-metallicity relation (AMR). The latter is consistently derived with the observed metallicity distribution of the local Red Clump (RC) giants from the APOGEE. Working within a Bayesian approach, we sample the posterior probability distribution in a multidimensional parameter space using the MCMC method. We find that the spatial distribution and motion of the Gaia stars imply two recent SF bursts centered at ages of ~0.5 Gyr and ~3 Gyr and characterised by a ~30% and ~55% SF enhancement, respectively, relative to a monotonously declining SFR continuum. The stellar populations associated with this SF excess are found to be dynamically hot for their age: they have W-velocity dispersions of ~12.5 km/s and ~26 km/s. The new JJ model is able to reproduce the local star counts with an accuracy of ~5 %. Using Gaia DR2 data, we self-consistently constrained 22 parameters of the updated JJ model. Our optimised model predicts two SF bursts within the last ~4 Gyr, which may point to recent episodes of gas infall.

preprint2019arXiv

Simulation of the loss-cone instability in spherical systems. I. Dominating harmonic potential

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.

preprint2019arXiv

Simulation of the loss-cone instability in spherical systems. II. Dominating Keplerian potential

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.

preprint2016arXiv

Impact of cosmological satellites on the vertical heating of the Milky Way disc

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 $σ_{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$σ$. We observed a general flaring of the disc height in the case of all 7 simulations in the outer disc.

preprint2015arXiv

Roche volume filling and the dissolution of open star clusters

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

preprint2015arXiv

The local stellar luminosity function and mass-to-light ratio in the NIR

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.

preprint2014arXiv

Chemical gradients in the Milky Way from the RAVE data. II. Giant stars

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.

preprint2013arXiv

Chemical gradients in the Milky Way from the RAVE data. I. Dwarf stars

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.

preprint2013arXiv

Sextans' cold substructures as a dynamical judge: Core, Cusp or MOND?

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.

preprint2013arXiv

Simulations of the Hyades

Context: Using the recent observational data of Röser 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.

preprint2013arXiv

The asymmetric drift, the local standard of rest, and implications from RAVE data

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ömberg 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ömberg 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önrich value of V_\sun=12.24 km/s for the LSR, the new Strömberg equation yields much larger individual radial scalelengths of the RAVE subpopulations, which seem unphysical in part. The new Strömberg 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.

preprint2013arXiv

The RAdial Velocity Experiment (RAVE): Fourth data release

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.

preprint2012arXiv

Dark Matter Subhalos in the Ursa Minor Dwarf Galaxy

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 Λ-cold dark matter scenario at dwarf galaxy scale.

preprint2012arXiv

The properties of the local spiral arms from RAVE data: two-dimensional density wave approach

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 Ω_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.

preprint2010arXiv

Calibration of radii and masses of open clusters with a simulation

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.

preprint2010arXiv

Dynamical friction of massive objects in galactic centres

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.

preprint2010arXiv

The local star formation history of the thin disc derived from kinematic data

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.

preprint2009arXiv

Towards a fully consistent Milky Way disc model: Part 1 The local model based on kinematic and photometric data

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.

preprint2008arXiv

Quantitative analysis of clumps in the tidal tails of star clusters

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.