Source author record

P. Kroupa

P. Kroupa 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

49works
6topics
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

49 published item(s)

preprint2022arXiv

First detection of a magnetic field in low-luminosity B[e] stars: New scenarios for the nature and evolutionary stages of FS CMa stars

We report the first detection of the magnetic field in a star of FS CMa type, a subgroup of objects characterized by the B[e] phenomenon. The split of magnetically sensitive lines in IRAS 17449+2320 determines the magnetic field modulus of 6.2+/-0.2 kG. Spectral lines and their variability reveal the presence of a B-type spectrum and a hot continuum source in the visible. The hot source confirms GALEX UV photometry. Because there is a lack of spectral lines for the hot source in the visible, the spectral fitting gives only the lower temperature limit of the hot source, which is 50 000 K, and the upper limit for the B-type star of 11 100 K. The V/R ratio of the H alpha line shows quasiperiodic behavior on timescale of 800 days. We detected a strong red-shifted absorption in the wings of Balmer and OI lines in some of the spectra. The absorption lines of helium and other metals show no, or very small, variations, indicating unusually stable photospheric regions for FS CMa stars. We detected two events of material infall, which were revealed to be discrete absorption components of resonance lines. The discovery of the strong magnetic field together with the Gaia measurements of the proper motion show that the most probable nature of this star is that of a post-merger object created after the leaving the binary of the birth cluster. Another possible scenario is a magnetic Ap star around Terminal-Age Main Sequence (TAMS). On the other hand, the strong magnetic field defies the hypothesis that IRAS 17449+2320 is an extreme classical Be star. Thus, IRAS 17449+2320 provides a pretext for exploring a new explanation of the nature of FS CMa stars or, at least, a group of stars with very similar spectral properties.

preprint2021arXiv

A possible solution to the Milky Way's binary-deficient retrograde stellar population. Evidence that $ω$ Centauri has formed in an extreme starburst

Context. The fraction of field binaries on retrograde orbits about the Milky Way is significantly lower compared to its prograde counterpart. Chemical and dynamical evidence suggests that the retrograde stellar population originates from $ω$ Centauri, which is either the most massive globular cluster (GC) of the Milky Way or the putative core of a former dwarf galaxy. Aims. Star formation conditions required to produce the retrograde binary population are constrained assuming that the retrograde stellar population originates from $ω$ Centauri's progenitor. Methods. We match the observed low binary fraction with dynamical population synthesis models, including a universal initial binary population and dynamical processing in star clusters, making use of the publicly available binary population synthesis tool BiPoS1. Results. It is found that either the GC progenitor of $ω$ Centauri must have formed with a stellar density of $\approx 10^8 \; M_{sun} \; pc^{-3}$ or the $ω$ Centauri dwarf galaxy's progenitor star cluster population must have formed in an extreme starburst with a star formation rate exceeding $1000 \; M_{sun} \; yr^{-1}$ and probably a top-heavy embedded-cluster mass function with suppressed low-mass cluster formation. The separation and mass ratio distribution for retrograde field binaries are predicted for comparison with future observations. Conclusions. A viable solution for the deficiency of binaries on retrograde orbits is presented, and star formation conditions for $ω$ Centauri as well as orbital parameter distributions for the Milky Way's retrograde binary population are predicted. The dwarf galaxy origin for $ω$ Centauri is tentatively preferred within the present context.

preprint2016arXiv

DGSAT: Dwarf Galaxy Survey with Amateur Telescopes I. Discovery of low surface brightness systems around nearby spiral galaxies

Context: We introduce the Dwarf Galaxy Survey with Amateur Telescopes (DGSAT) project and report the discovery of eleven low surface brightness (LSB) galaxies in the fields of the nearby galaxies NGC 2683, NGC 3628, NGC 4594 (M104), NGC 4631, NGC 5457 (M101), and NGC 7814. Aims: The DGSAT project aims to use the potential of small-sized telescopes to probe LSB features around large galaxies and to increase the sample size of the dwarf satellite galaxies in the Local Volume. Methods: Using long exposure images, fields of the target spiral galaxies are explored for extended low surface brightness objects. After identifying dwarf galaxy candidates, their observed properties are extracted by fitting models to their light profiles. Results: We find three, one, three, one, one, and two new LSB galaxies in the fields of NGC 2683, 3628, 4594, 4631, 5457, and 7814, respectively. In addition to the newly found galaxies, we analyse the structural properties of nine already known galaxies. All of these 20 dwarf galaxy candidates have effective surface brightnesses in the range $25.3\lesssimμ_{e}\lesssim28.8$ mag.arcsec$^{-2}$ and are fit with Sersic profiles with indices $n\lesssim 1$. Assuming that they are in the vicinity of the above mentioned massive galaxies, their $r$-band absolute magnitudes, their effective radii, and their luminosities are in the ranges $-15.6 \lesssim M_r \lesssim -7.8$, 160 pc $\lesssim R_e \lesssim$ 4.1 kpc, and $0.1\times 10^6 \lesssim\left(\frac{L}{L_{\odot}}\right)_r\lesssim127 \times 10^6$, respectively. To determine whether these LSB galaxies are indeed satellites of the above mentioned massive galaxies, their distances need to be determined via further observations. Conclusions: Using small telescopes, we are readily able to detect LSB galaxies with similar properties to the known dwarf galaxies of the Local Group.

preprint2016arXiv

Distribution of star formation rates during the rapid assembly of NGC 1399 as deduced from its globular cluster system

Ultra-compact dwarf galaxies (UCDs) share many properties with globular clusters (GCs) and are found in similar environments. A large sample of UCDs and GCs around NGC 1399, the central giant elliptical of the Fornax galaxy cluster, is used to infer their formation history and also that of NGC 1399. We assumed that all GCs and UCDs in our sample are star clusters (SCs) and used them as tracers of past star formation activities. After correcting our GC/UCD sample for mass loss, we interpreted their overall mass function to be a superposition of SC populations that formed coevally during different times. The SC masses of each population were distributed according to the embedded cluster mass function (ECMF), a pure power law with the slope $-β$ and a stellar upper mass limit, $M_{\mathrm{max}}$, which depended on the star formation rate (SFR). We decomposed the observed GC/UCD mass function into individual SC populations and converted $M_{\mathrm{max}}$ of each SC population to an SFR. The overall distribution of SFRs reveals how the GC/UCD sample formed. Considering the age of the GCs/UCDs and the present stellar mass of NGC 1399, we found that the formation of the GCs/UCDs can be well explained for $β<2.3$. This agrees very well with the observation in young SCs where $β\approx2.0$ is usually found. Even if taking into account that some of the most massive objects might not be genuine SCs and applying different corrections for the mass loss, the outcome is not influenced much. We found peak SFRs between approximately 300 and 3000 $M_{\odot}\mathrm{yr}^{-1}$, which matches the SFRs observed in massive high-$z$ sub-mm galaxies and an SFR estimate inferred from NGC 1399 based on "downsizing", i.e. more massive galaxies must have formed over shorter times. Our results suggest that NGC 1399 and its GC/UCD system formed in an early, short, and intense starburst.

preprint2016arXiv

Understanding the internal dynamics of elliptical galaxies without non-baryonic dark matter

Assuming virial equilibrium and Newtonian dynamics, low-mass early-type galaxies have larger velocity dispersions than expected from the amount of baryons they contain. The conventional interpretation of this finding is that their dynamics is dominated by non-baryonic matter. However, there is also strong evidence that many low-mass early-type galaxies formed as tidal dwarf galaxies, which would contain almost no dark matter. Using an extensive catalogue of early-type galaxies, we therefore discuss how the internal dynamics of early-type galaxies in general can be understood by replacing the assumption of non-baryonic dark matter with two alternative assumptions. The first assumption is that Milgromian dynamics (i.e., MOND) is valid, which changes the effective gravitational force in the weak-field limit. The second assumption is that binary stars affect the observed line-of-sight velocity dispersions. Some moderate discrepancies between observed and predicted velocity dispersions remain also when these effects are implemented. Nevertheless, the observed velocity dispersions in early-type galaxies can then easily be explained without invoking the presence of non-baryonic dark matter in them, but with already documented variations of the galaxy-wide stellar initial mass function and non-equilibrium dynamics in some of the low-mass early-type galaxies.

preprint2015arXiv

Globular Cluster Streams as Galactic High-Precision Scales

Tidal streams of globular clusters are ideal tracers of the Galactic gravitational potential. Compared to the few known, complex and diffuse dwarf-galaxy streams, they are kinematically cold, have thin morphologies and are abundant in the halo of the Milky Way. Their coldness and thinness in combination with potential epicyclic substructure in the vicinity of the stream progenitor turns them into high-precision scales. With the example of Palomar 5, we demonstrate how modeling of a globular cluster stream allows us to simultaneously measure the properties of the disrupting globular cluster, its orbital motion, and the gravitational potential of the Milky Way.

preprint2015arXiv

The mass-metallicity relation of tidal dwarf galaxies

Dwarf galaxies generally follow a mass-metallicity (MZ) relation, where more massive objects retain a larger fraction of heavy elements. Young tidal dwarf galaxies (TDGs), born in the tidal tails produced by interacting gas-rich galaxies, have been thought to not follow the MZ relation, because they inherit the metallicity of the more massive parent galaxies. We present chemical evolution models to investigate if TDGs that formed at very high redshifts, where the metallicity of their parent galaxy was very low, can produce the observed MZ relation. Assuming that galaxy interactions were more frequent in the denser high-redshift universe, TDGs could constitute an important contribution to the dwarf galaxy population. The survey of chemical evolution models of TDGs presented here captures for the first time an initial mass function (IMF) of stars that is dependent on both the star formation rate and the gas metallicity via the integrated galactic IMF (IGIMF) theory. As TDGs form in the tidal debris of interacting galaxies, the pre-enrichment of the gas, an underlying pre-existing stellar population, infall, and mass dependent outflows are considered. The models of young TDGs that are created in strongly pre-enriched tidal arms with a pre-existing stellar population can explain the measured abundance ratios of observed TDGs. The same chemical evolution models for TDGs, that form out of gas with initially very low metallicity, naturally build up the observed MZ relation. The modelled chemical composition of ancient TDGs is therefore consistent with the observed MZ relation of satellite galaxies.

preprint2015arXiv

Young tidal dwarf galaxies cannot be used to probe dark matter in galaxies

The location of dark-matter free, tidal dwarf galaxies (TDGs) in the baryonic Tully Fisher (bTF) diagram has been used to test cosmological scenarios, leading to various and controversial results. Using new high-resolution 3D spectroscopic data, we re-investigate the morpho-kinematics of these galaxies to verify whether or not they can be used for such a purpose. We find that the three observed TDGs are kinematically not virialized and show complex morphologies and kinematics, leading to considerable uncertainties about their intrinsic rotation velocities and their locations on the bTF. Only one TDG can be identify as a (perturbed) rotation disk that it is indeed a sub-component of NGC5291N and that lies at $<$1$σ$ from the local bTF relation. It results that the presently studied TDGs are young, dynamically forming objects, which are not enough virialized to robustly challenge cosmological scenarios.

preprint2014arXiv

The chemical evolution of galaxies with a variable IGIMF

Standard analytical chemical evolution modelling of galaxies has been assuming the stellar initial mass function (IMF) to be invariant and fully sampled allowing fractions of massive stars to contribute even in dwarf galaxies with very low star formation rates (SFRs). Recent observations show the integrated galactic initial mass function (IGIMF) of stars, i.e. the galaxy-wide IMF, to become systematically top-heavy with increasing SFR. This has been predicted by the IGIMF theory, which is here used to develop the analytical theory of the chemical evolution of galaxies. This theory is non-linear and requires the iterative solution of implicit integral equations due to the dependence of the IGIMF on the metallicity and on the SFR. It is shown that the mass-metallicity relation of galaxies emerges naturally, although at low masses the theoretical predictions overestimate the observations by 0.3--0.4 dex. A good agreement with the observation can be obtained only if gas flows are taken into account. In particular, we are able to reproduce the mass--metallicity relation observed by Lee et al. (2006) with modest amounts of infall and with an outflow rate which decreases as a function of the galactic mass. The outflow rates required to fit the data are considerably smaller than required in models with invariant IMFs.

preprint2013arXiv

Chemo-dynamical evolution of tidal dwarf galaxies. I. Method and IMF dependence

We present high-resolution simulations of tidal dwarf galaxies (TDG) to investigate their early chemo-dynamical evolution and test their survivability. In this work the simulation setup is introduced and the response of TDGs to self-consistent star formation (SF) and an external tidal field is examined. Throughout the simulation star cluster particles with variable masses down to $5\,M_{\odot}$ form, depending on the local gas reservoir. For low cluster masses $M_{\mathrm{cl}}$, the stellar initial mass function (IMF) is considered to be either filled or truncated at a maximal star mass $m_\mathrm{max}$ to represent the observed $m_{\mathrm{max}} - M_{\mathrm{cl}}$ relation (IGIMF theory). The evolution of TDGs with fully-populated and truncated IMFs are compared to study the impact of stellar energy feedback on their survivability. Both TDGs experience an initial starburst but after a dynamical time they evolve into dwarf galaxies with self-regulated and continuous SF. At this stage the truncated-IMF model contains about 6 times more stellar mass than the invariant IMF models, but the final bound gas mass is comparable in both models. In spite of their significantly different SF histories, both TDG models are not disrupted within the first 500 Myr. We conclude that TDGs can survive an early starburst, independent of the underlying IMF description, even though they do not harbor a stabilizing dark matter halo.

preprint2013arXiv

Main Sequence Star Populations in the Virgo Overdensity Region

We present deep CMDs for two Subaru Suprime-Cam fields in the Virgo Stellar Stream(VSS)/Virgo Overdensity(VOD) and compare them to a field centred on the highest concentration of Sagittarius (Sgr) Tidal Stream stars in the leading arm, Branch A of the bifurcation. A prominent population of MS stars is detected in all three fields and can be traced as faint as g~24 mag. Using theoretical isochrone fitting we derive an age of 9.1(+1.0;-1.1)Gyr, a median abundance of [Fe/H]=-0.70 (+0.15; -0.20)dex and d_helio of 30.9+-3.0kpc for the MS of the Sgr Stream Branch A. The dominant main sequence populations in the two VSS/VOD fields (Lsun ~265 deg, Bsun ~13 deg) are located at a mean distance of 23.3+-1.6kpc and have an age ~8.2Gyr and an abundance [Fe/H]=-0.67(+0.16;-0.12)dex similar to the Sgr Stream stars. These parameters are also in good agreement with the age of the main population in the Sagittarius dwarf (8.0+-1.5Gyr). They also agree with the peak in the metallicity distribution of 2-3Gyr old M-giants, [Fe/H] ~ -0.6dex, in the Sgr leading arm north. We then compare the results from the VSS/VOD fields with the Sgr Tidal Stream model by Law & Majewski based on a triaxial Galactic halo that is calibrated with SDSS Sgr A-branch and 2MASS M-giant stars. We find that the most prominent feature in the CMDs, the MS population at 23kpc, is not explained by the model. Instead the model predicts in these directions a low density filamentary structure of Sgr debris stars at ~9kpc and a slightly higher concentration of Sgr stars spread from 42-53kpc. At best there is only marginal evidence for the presence of these populations in our data. Our findings then suggest that while there are probably some Sgr debris stars present, the dominant stellar population in the VOD originates from a different halo structure that has almost identical age and metallicity as some sections of the Sgr tidal stream.

preprint2013arXiv

Phase mixing due to the Galactic potential: steps in the position and velocity distributions of popped star clusters

As star clusters are expected to form with low star formation efficiencies, the gas in the cluster is expelled quickly and early in their development: the star cluster "pops." This leads to an unbound stellar system, evolving in the Galactic potential. Previous N-body simulations have demonstrated the existence of a stepped number density distribution of cluster stars after popping, both in vertical position and vertical velocity, with a passing resemblance to a Christmas tree. Using numerical and analytical methods, we investigate the source of this structure, which arises due to the phase mixing of the out-of-equilibrium stellar system, determined entirely by the background analytic potential. Considering only the vertical motions, we construct a theoretical model to describe the time evolution of the phase space distribution of stars in a Miyamoto-Nagai disk potential and a full Milky-Way type potential comprising bulge, halo and disk components, which is then compared with N-body simulations. Using our theoretical model, we investigate the possible observational signatures and the feasibility of detection.

preprint2013arXiv

Polar ring galaxies as tests of gravity

Polar ring galaxies are ideal objects with which to study the three-dimensional shapes of galactic gravitational potentials since two rotation curves can be measured in two perpendicular planes. Observational studies have uncovered systematically larger rotation velocities in the extended polar rings than in the associated host galaxies. In the dark matter context, this can only be explained through dark halos that are systematically flattened along the polar rings. Here, we point out that these objects can also be used as very effective tests of gravity theories, such as those based on Milgromian dynamics (MOND). We run a set of polar ring models using both Milgromian and Newtonian dynamics to predict the expected shapes of the rotation curves in both planes, varying the total mass of the system, the mass of the ring with respect to the host, as well as the size of the hole at the center of the ring. We find that Milgromian dynamics not only naturally leads to rotation velocities being typically higher in the extended polar rings than in the hosts, as would be the case in Newtonian dynamics without dark matter, but that it also gets the shape and amplitude of velocities correct. Milgromian dynamics thus adequately explains this particular property of polar ring galaxies.

preprint2013arXiv

The [alpha/Fe] ratios of very metal-poor stars within the IGIMF theory

The aim of this paper is to quantify the amplitude of the predicted plateau in [alpha/Fe] ratios associated with the most metal-poor stars of a galaxy. We assume that the initial mass function in galaxies is steeper if the star formation rate (SFR) is low -- as per the integrated galactic initial mass function (IGIMF) theory. A variant of the theory, in which the IGIMF depends upon the metallicity of the parent galaxy, is also considered. The IGIMF theory predicts low [alpha/Fe] plateaus in dwarf galaxies, characterised by small SFRs. The [alpha/Fe] plateau is up to 0.7dex lower than the corresponding plateau of the Milky Way. For a universal IMF one should expect instead that the [alpha/Fe] plateau is the same for all the galaxies, irrespective of their masses or SFRs. Assuming a strong dependence of the IMF on the metallicity of the parent galaxy, dwarf galaxies can show values of the [alpha/Fe] plateau similar to those of the Milky Way, and almost independent on the SFR. The [Mg/Fe] ratios of the most metal-poor stars in dwarf galaxies satellites of the Milky Way can be reproduced either if we consider metallicity-dependent IMFs or if the early SFRs of these galaxies were larger than we presently think. Present and future observations of dwarf galaxies can help disentangle between these different IGIMF formulations.

preprint2013arXiv

The Curious Case of Palomar 13: The Influence of the Orbital Phase on the Appearance of Galactic Satellites

We investigate the dynamical status of the low-mass globular cluster Palomar 13 by means of N-body computations to test whether its unusually high mass-to-light ratio of about 40 and its peculiarly shallow surface density profile can be caused by tidal shocking. Alternatively, we test - by varying the assumed proper motion - if the orbital phase of Palomar 13 within its orbit about the Milky Way can influence its appearance and thus may be the origin of these peculiarities, as has been suggested by Kuepper et al. (2010). We find that, of these two scenarios, only the latter can explain the observed mass-to-light ratio and surface density profile. We note, however, that the particular orbit that best reproduces those observed parameters has a proper motion inconsistent with the available literature value. We discuss this discrepancy and suggest that it may be caused by an underestimation of the observational uncertainties in the proper motion determination. We demonstrate that Palomar 13 is most likely near apogalacticon, which makes the cluster appear supervirial and blown-up due to orbital compression of its tidal debris. Since the satellites of the Milky Way are on average closer to apo- than perigalacticon, their internal dynamics may be influenced by the same effect, and we advocate that this needs to be taken into account when interpreting their kinematical data. Moreover, we briefly discuss the influence of a possible binary population on such measurements.

preprint2012arXiv

A Catalog of Extended Clusters and Ultra-Compact Dwarf Galaxies - An Analysis of their Parameters in Early- and Late-Type Galaxies

Extended stellar clusters with effective radii larger than 10 pc have been found in various environments. Objects with masses comparable to globular clusters (GCs) are called extended clusters (ECs), while objects with masses in the dwarf galaxy regime are called ultra-compact dwarf galaxies (UCDs). The paper analyses the observational parameters luminosity, effective radius, and projected distance to the host galaxy, of all known ECs and UCDs and the dependence of these parameters on the type and the luminosity of their host galaxy. We searched the available literature to compile a catalog of star clusters larger than 10 pc. As there is no clear distinction between ECs and UCDs, both types of objects will be called extended stellar objects (EOs). In total, we found 813 EOs of which 171 are associated with late-type and 642 with early-type galaxies. EOs cover a luminosity range from about MV = -4 to -14 mag. However, the vast majority of EOs brighter than -10 mag are associated with elliptical galaxies. At each magnitude EOs are found with effective radii between 10 pc and an upper size limit, which shows a clear trend: the more luminous the object the larger is the upper size limit. For EOs associated with early-and late-type galaxies, the luminosity functions peak at -6.40 and -6.47 mag, respectively, which is about one magnitude fainter than the peak of the GC luminosity function. EOs and GCs form a coherent structure in the reff vs. MV parameter space, while there is a clear gap between EOs and early type dwarf galaxies. However, there is a small potential overlap at the high-mass end, where the most extended EOs are close to the parameters of compact elliptical galaxies. We compare the EO sample with numerical models and conclude that the parameters of the EO sample as a whole can be very well explained by a star cluster origin, where EOs are the results of merged cluster complexes.

preprint2012arXiv

Can filamentary accretion explain the orbital poles of the Milky Way satellites?

Several scenarios have been suggested to explain the phase-space distribution of the Milky Way (MW) satellite galaxies in a disc of satellites (DoS). To quantitatively compare these different possibilities, a new method analysing angular momentum directions in modelled data is presented. It determines how likely it is to find sets of angular momenta as concentrated and as close to a polar orientation as is observed for the MW satellite orbital poles. The method can be easily applied to orbital pole data from different models. The observed distribution of satellite orbital poles is compared to published angular momentum directions of subhalos derived from six cosmological state-of-the-art simulations in the Aquarius project. This tests the possibility that filamentary accretion might be able to naturally explain the satellite orbits within the DoS. For the most likely alignment of main halo and MW disc spin, the probability to reproduce the MW satellite orbital pole properties turns out to be less than 0.5 per cent in Aquarius models. Even an isotropic distribution of angular momenta has a higher likelihood to produce the observed distribution. The two Via Lactea cosmological simulations give results similar to the Aquarius simulations. Comparing instead with numerical models of galaxy-interactions gives a probability of up to 90 per cent for some models to draw the observed distribution of orbital poles from the angular momenta of tidal debris. This indicates that the formation as tidal dwarf galaxies in a single encounter is a viable, if not the only, process to explain the phase-space distribution of the MW satellite galaxies.

preprint2012arXiv

Dwarf spheroidal satellites of the Milky Way from dark matter free tidal dwarf galaxy progenitors: maps of orbits

The long term time evolution of tidal dwarf satellite galaxies with two different initial densities orbiting a host galaxy that resembles the Milky Way has been studied using a large set of Newtonian N-Body simulations. From the simulations two maps of the orbital conditions that lead to quasi-equilibrium objects were constructed. It has been found that several orbits of the satellites allow for the existence, for about 1 Gyr or more, of out-of-equilibrium bodies with high apparent mass-to-light ratios. Within this framework the satellites in the quasi-stable phase reproduce the observed satellite properties for about 16% of the orbit for high density progenitors, and for about 66% for progenitors with lower densities An additional simulation for a single satellite with initial mass of 10^7 Msun and Plummer radius of 0.15 kpc leads to remnants in the quasi- equilibrium phase that simultaneously reproduce remarkably well the observational quantities of the UFDGs of the Milky Way. This satellite in the quasi-stable phase reproduces the observed satellite properties for about 42% of the orbit. The results suggest that a fraction of the observed satellites could plausibly be galaxies without dark matter that have true M/L ratios much lower than those measured. The inflated M/L ratios arise because they are observed at the right time, along the right orbit and during the quasi-equilibrium phase of their evolution. This is a viable explanation for the high M/L ratios observed in all satellites as long as the satellites are preferentially on certain orbits and are observed at certain times. This could arise within the TDG scenario if all satellites are created at the same time along a few specific orbits that are particularly susceptible to the quasi-equilibrium phase.

preprint2012arXiv

Field O stars: formed in situ or as runaways?

A significant fraction of massive stars in the Milky Way and other galaxies are located far from star clusters. It is known that some of these stars are runaways and therefore most likely were formed in embedded clusters and then ejected into the field because of dynamical few-body interactions or binary-supernova explosions. However, there exists a group of field O stars whose runaway status is difficult to prove via direct proper motion measurements or whose low space velocities and/or young ages appear to be incompatible with their large separation from known star clusters. The existence of this group led some authors to believe that field O stars can form in situ. In this paper, we examine the runaway status of the best candidates for isolated formation of massive stars in the Milky Way and the Magellanic Clouds by searching for bow shocks around them, by using the new reduction of the Hipparcos data, and by searching for stellar systems from which they could originate within their lifetimes. We show that most of the known O stars thought to have formed in isolation are instead very likely runaways. We show also that the field must contain a population of O stars whose low space velocities and/or young ages are in apparent contradiction with the large separation of these stars from their parent clusters and/or the ages of these clusters. These stars (the descendants of runaway massive binaries) cannot be traced back to their parent clusters and therefore can be mistakenly considered as having formed in situ. We argue also that some field O stars could be detected in optical wavelengths only because they are runaways, while their cousins residing in the deeply embedded parent clusters might still remain totally obscured. The main conclusion of our study is that there is no significant evidence whatsoever in support of the in situ proposal on the origin of massive stars.

preprint2012arXiv

The VPOS: a vast polar structure of satellite galaxies, globular clusters and streams around the Milky Way

It has been known for a long time that the satellite galaxies of the Milky Way (MW) show a significant amount of phase-space correlation, they are distributed in a highly inclined Disc of Satellites (DoS). We have extended the previous studies on the DoS by analysing for the first time the orientations of streams of stars and gas, and the distributions of globular clusters within the halo of the MW. It is shown that the spatial distribution of MW globular clusters classified as young halo clusters (YH GC) is very similar to the DoS, while 7 of the 14 analysed streams align with the DoS. The probability to find the observed clustering of streams is only 0.3 per cent when assuming isotropy. The MW thus is surrounded by a vast polar structure (VPOS) of subsystems (satellite galaxies, globular clusters and streams), spreading from Galactocentric distances as small as 10 kpc out to 250 kpc. These findings demonstrate that a near-isotropic infall of cosmological sub-structure components onto the MW is essentially ruled out because a large number of infalling objects would have had to be highly correlated, to a degree not natural for dark matter sub-structures. The majority of satellites, streams and YH GCs had to be formed as a correlated population. This is possible in tidal tails consisting of material expelled from interacting galaxies. We discuss the tidal scenario for the formation of the VPOS, including successes and possible challenges. The potential consequences of the MW satellites being tidal dwarf galaxies are severe. If all the satellite galaxies and YH GCs have been formed in an encounter between the young MW and another gas-rich galaxy about 10-11 Gyr ago, then the MW does not have any luminous dark-matter substructures and the missing satellites problem becomes a catastrophic failure of the standard cosmological model.

preprint2011arXiv

A New Formation Scenario for the Milky Way Cluster NGC2419

We present a new formation scenario for NGC2419, which is one of the most luminous,one of the most distant, and as well one of the most extended globular clusters of the Milky Way. We propose that NGC2419 is the remnant of a merged star cluster complex, which was possibly formed during an interaction between a gas-rich galaxy and the Milky Way. To test this hypothesis, we performed numerical simulations of 27 different models of star cluster complexes (CCs) moving on a highly eccentric orbit in the Galactic halo. We vary the CC mass, the CC size, and the initial distribution of star clusters in the CC to analyze the influence of these parameters on the resulting objects. In all cases, the vast majority of star clusters merged into a stable object. The derived parameters mass, absolute V-band magnitude, effective radius, velocity dispersion and the surface brightness profile are, for a number of models, in good agreement with those observed for NGC2419. Despite the large range of CC sizes, the effective radii of the merger objects are constrained to a relatively small interval. A turnover in the r_eff vs. M_encl space leads to degenerate states, i.e. relatively compact CCs can produce an object with the same structural parameters as a more massive and larger CC. In consequence, a range of initial conditions can form a merger object comparable to NGC2419 preventing us to pinpoint the exact parameters of the original CC, which formed NGC2419. We conclude that NGC2419 can be well explained by the merged cluster complex scenario. Some of the recently discovered stellar streams in the Galactic halo might be related to the parent galaxy, which produced the cluster complex in our scenario. Measurements of the proper motion of NGC2419 are necessary to prove an association with one of the stellar streams.

preprint2011arXiv

A parametric study on the formation of extended star clusters and ultra-compact dwarf galaxies

In the last decade, extended stellar clusters with masses in the range from a few 10^4 to 10^8 M_sun have been found in various types of galaxies in different environments. Objects with masses comparable to normal globular clusters (GCs) are called extended clusters (ECs), while objects with masses in the dwarf galaxy regime are called ultra-compact dwarf galaxies (UCDs). In heavily interacting galaxies star clusters tend to form in larger conglomerations called star cluster complexes (CCs). In this work we systematically scan a suitable parameter space for CCs and perform numerical simulations to study their further fate. The varied sizes and masses of the CCs cover a matrix of 5x6 values with CC Plummer radii between 10 - 160 pc and CC masses between 10^5.5 - 10^8 M_sun, which are consistent with observed CC parameters. The CCs of the parametric study are on orbits with distances between 20 kpc and 60 kpc. In addition, we studied also the evolution of CCs on a circular orbit at a distance of 60 kpc to verify that also extremely extended ECs and UCDs can be explained by our formation scenario. All 54 simulations end up with stable merger objects, wherein 26 to 97% of the initial CC mass is bound. The objects show a general trend of increasing effective radii with increasing mass. Despite the large range of input Plummer radii of the CCs (10 to 160 pc) the effective radii of the merger objects are constrained to values between 10 and 20 pc at the low mass end and to values between 15 and 55 pc at the high mass end. The structural parameters of the models are comparable to those of the observed ECs and UCDs. The results of the circular orbits demonstrate that even very extended objects like the M31 ECs found by Huxor in 2005 and the very extended (r_eff > 80 pc), high-mass UCDs can be explained by merged cluster complexes in regions with low gravitational fields at large distances.

preprint2011arXiv

Making Counter-Orbiting Tidal Debris: The Origin of the Milky Way Disc of Satellites

Using stellar-dynamical calculations it is shown for the first time that counter-orbiting material emerges naturally in tidal interactions of disc galaxies. Model particles on both pro- and retrograde orbits can be formed as tidal debris in single encounters with disc galaxies of 1-to-1 and 4-to-1 mass ratios. A total of 74 model calculations are performed for a range of different initial parameters. Interactions include fly-by and merger cases. The fraction of counter-orbiting material produced varies over a wide range (from a few up to 50 percent). All fly-by models show a similar two-phase behaviour, with retrograde material forming first. Properties of the prograde and retrograde populations are extracted to make an observational discrimination possible. During such encounters the tidal debris occupies a certain region in phase space. In this material, tidal-dwarf galaxies may form. The modelling therefore can explain why galaxies may have dwarf galaxies orbiting counter to the bulk of their dwarf galaxies. An example is the Sculptor dwarf of the Milky Way, which orbits counter to the bulk of the disc of satellites. The modelling thus supports the scenario of the MW satellites being ancient tidal-dwarf galaxies formed from gaseous material stripped from another galaxy during an encounter with the young MW. A possible candidate for this galaxy is identified as the Magellanic Cloud progenitor galaxy. Its angular motion fits the angular motion of the MW disc of satellites objects. This scenario is in agreement with Lynden-Bell's original suggestion for the origin of the dSph satellites and the near-unbound orbit of the LMC.

preprint2011arXiv

Mass segregation and fractal substructure in young massive clusters: (I) the McLuster code and method calibration

By analysing models of the young massive cluster R136 in 30 Doradus, set-up using the herewith introduced and publicly made available code McLuster, we investigate and compare different methods for detecting and quantifying mass segregation and substructure in non-seeing limited N-body data. For this purpose we generate star cluster models with different degrees of mass segregation and fractal substructure and analyse them. We quantify mass segregation by measuring, from the projected 2d model data, the mass function slope in radial annuli, by looking for colour gradients in radial colour profiles, by measuring Allison's Lambda parameter, and by determining the local stellar surface density around each star. We find that these methods for quantifying mass segregation often produce ambiguous results. Most reliable for detecting mass segregation is the mass function slope method, whereas the colour gradient method is the least practical in an R136-like configuration. The other two methods are more sensitive to low degrees of mass segregation but are computationally much more demanding. We also discuss the effect of binaries on these measures. Moreover, we quantify substructure by looking at the projected radial stellar density profile, by comparing projected azimuthal stellar density profiles, and by determining Cartwright & Whitworth's Q parameter. We find that only high degrees of substructure affect the projected radial density profile, whereas the projected azimuthal density profile is very sensitive to substructure. The Q parameter is also sensitive to substructure but its absolute value shows a dependence on the radial density gradient of the cluster and is strongly influenced by binaries. (abridged)

preprint2011arXiv

Popping star clusters as building blocks of the Milky Way Thick Disc

It is widely believed that star clusters form with low star formation efficiencies. With the onset of stellar winds by massive stars or finally when the first super nova blows off, the residual gas is driven out of the embedded star cluster. Due to this fact a large amount, if not all, of the stars become unbound and disperse in the gravitational potential of the galaxy. In this context, Kroupa (2002) suggested a new mechanism for the emergence of thickened Galactic discs. Massive star clusters add kinematically hot components to the galactic field populations, building up in this way, the Galactic thick disc as well. In this work we perform, for the first time, numerical simulations to investigate this scenario for the formation of the galactic discs of the Milky Way. We find that a significant kinematically hot population of stars may be injected into the disk of a galaxy such that a thick disk emerges. For the MW the star clusters that formed the thick disk must have had masses of about 10^6 Msol.

preprint2011arXiv

Search for OB stars running away from young star clusters. II. The NGC 6357 star-forming region

Dynamical few-body encounters in the dense cores of young massive star clusters are responsible for the loss of a significant fraction of their massive stellar content. Some of the escaping (runaway) stars move through the ambient medium supersonically and can be revealed via detection of their bow shocks (visible in the infrared, optical or radio). In this paper, which is the second of a series of papers devoted to the search for OB stars running away from young (several Myr) Galactic clusters and OB associations, we present the results of the search for bow shocks around the star-forming region NGC 6357. Using the archival data of the Midcourse Space Experiment (MSX) satellite and the Spitzer Space Telescope, and the preliminary data release of the Wide-Field Infrared Survey Explorer (WISE), we discovered seven bow shocks, whose geometry is consistent with the possibility that they are generated by stars expelled from the young star clusters, Pismis 24 and AH03 J1725-34.4, associated with NGC 6357. Two of the seven bow shocks are driven by the already known O stars. Follow-up spectroscopy of three other bow shock-producing stars showed that they are O-type stars as well, while the 2MASS photometry of the remaining two stars suggests that they could be B0 V stars, provided that both are located at the same distance as NGC 6357. Detection of numerous massive stars ejected from the very young clusters is consistent with the theoretical expectation that star clusters can effectively lose massive stars at the very beginning of their dynamical evolution and lends strong support to the idea that probably all field OB stars have been dynamically ejected from their birth clusters. A by-product of our search for bow shocks around NGC 6357 is the detection of three circular shells typical of luminous blue variable and late WN-type Wolf-Rayet stars.

preprint2011arXiv

The evolution of the surface brightness of a star cluster as a result of residual star-forming gas expulsion

Direct N-body calculations are presented of the early evolution of exposed clusters to quantify the influence of gas expulsion on the time-varying surface brightness. By assuming that the embedded OB stars drive out most of the gas after a given time delay, the change of the surface brightness of expanding star clusters is studied. The influence of stellar dynamics and stellar evolution is discussed. The growth of the core radii of such models shows a remarkable core re-virialisation. The decrease of the surface mass density during gas expulsion is large and is only truncated by this re-virialisation process. However, the surface brightness within a certain radius does not increase noticeably. Thus, an embedded star cluster cannot reappear in observational surveys after re-virialisation. This finding has a bearing on the observed infant mortality fraction.

preprint2011arXiv

The origin of very wide binary stars

A large population of fragile, wide (> 1000 AU) binary systems exists in the Galactic field and halo. These wide binary stars cannot be primordial because of the high stellar density in star forming regions, while formation by capture in the Galactic field is highly improbable. We propose that these binary systems were formed during the dissolution phase of star clusters (see Kouwenhoven et al. 2010, for details). Stars escaping from a dissolving star cluster can have very similar velocities, which can lead to the formation of a wide binary systems. We carry out N-body simulations to test this hypothesis. The results indicate that this mechanism explains the origin of wide binary systems in the Galaxy. The resulting wide binary fraction and semi-major axis distribution depend on the initial conditions of the dissolving star cluster, while the distributions in eccentricity and mass ratio are universal. Finally, since most stars are formed in (relatively tight) primordial binaries, we predict that a large fraction of the wide "binary stars" are in fact higher-order multiple systems.

preprint2011arXiv

Tidal disruption rate of stars by supermassive black holes obtained by direct N-body simulations

The disruption rate of stars by supermassive black holes (SMBHs) is calculated numerically with a modified version of Aarseth's NBODY6 code. The initial stellar distribution around the SMBH follows a Sérsic n=4 profile representing bulges and early type galaxies. In order to infer relaxation driven effects and to increase the statistical significance, a very large set of N-body integrations with different particle numbers N, ranging from 10^{3} to 0.5 \cdot 10^{6} particles, is performed. Three different black hole capture radii are taken into account, enabling us to scale these results to a broad range of astrophysical systems with relaxation times shorter than one Hubble time, i.e. for SMBHs up to M_bh \approx 10^{7} M_sun. The computed number of disrupted stars are driven by diffusion in angular momentum space into the loss cone of the black hole and the rate scales with the total number of particles as dN/dt \propto N^{b}, where b is as large as 0.83. This is significantly steeper than the expected scaling dN/dt \propto ln(N) derived from simplest energy relaxation arguments. Only a relatively modest dependence of the tidal disruption rate on the mass of the SMBH is found and we discuss our results in the context of the M_bh/sigma relation. The number of disrupted stars contribute a significant part to the mass growth of black holes in the lower mass range as long as a significant part of the stellar mass becomes swallowed by the SMBH. This also bears direct consequences for the search and existence of IMBHs in globular clusters. For SMBHs similar to the galactic center black hole SgrA*, a tidal disruption rate of 55 \pm 27 events per Myr is deduced. Finally relaxation driven stellar feeding can not account for the masses of massive black holes M_bh \geq 10^{7} M_sun. (abridged)

preprint2010arXiv

A top-heavy stellar initial mass function in starbursts as an explanation for the high mass-to-light ratios of ultra compact dwarf galaxies

It has been shown recently that the dynamical V-band mass-to-light ratios of compact stellar systems with masses from 10^6 to 10^8 Solar masses are not consistent with the predictions from simple stellar population (SSP) models. Top-heavy stellar initial mass functions (IMFs) in these so-called ultra compact dwarf galaxies (UCDs) offer an attractive explanation for this finding, the stellar remnants and retained stellar envelopes providing the unseen mass. We therefore construct a model which quantifies by how much the IMFs of UCDs would have to deviate in the intermediate-mass and high-mass range from the canonical IMF in order to account for the enhanced M/L_V ratio of the UCDs. The deduced high-mass IMF in the UCDs depends on the age of the UCDs and the number of faint products of stellar evolution retained by them. Assuming that the IMF in the UCDs is a three-part power-law equal to the canonical IMF in the low-mass range and taking 20% as a plausible choice for the fraction of the remnants of high-mass stars retained by UCDs, the model suggests the exponent of the high-mass IMF to be approximately 1.6 if the UCDs are 13 Gyr old (i.e. almost as old as the Universe) or approximately 1.0 if the UCDs are 7 Gyr old, in contrast to 2.3 for the Salpeter-Massey IMF. If the IMF was as top-heavy as suggested here, the stability of the UCDs might have been threatened by heavy mass loss induced by the radiation and evolution of massive stars. The central densities of UCDs must have been in the range 10^6 to 10^7 Solar masses per cubic parsec when they formed with star formation rates of 10 to 100 Solar masses per year.

preprint2010arXiv

Effects of the integrated galactic IMF on the chemical evolution of the solar neighbourhood

The initial mass function determines the fraction of stars of different intial mass born per stellar generation. In this paper, we test the effects of the integrated galactic initial mass function (IGIMF) on the chemical evolution of the solar neighbourhood. The IGIMF (Weidner & Kroupa 2005) is computed from the combination of the stellar intial mass function (IMF), i.e. the mass function of single star clusters, and the embedded cluster mass function, i.e. a power law with index beta. By taking into account also the fact that the maximum achievable stellar mass is a function of the total mass of the cluster, the IGIMF becomes a time-varying IMF which depends on the star formation rate. We applied this formalism to a chemical evolution model for the solar neighbourhood and compared the results obtained by assuming three possible values for beta with the results obtained by means of a standard, well-tested, constant IMF. In general, a lower absolute value of beta implies a flatter IGIMF, hence a larger number of massive stars and larger metal ejection rates. This translates into higher type Ia and II supernova rates, higher mass ejection rates from massive stars and a larger amount of gas available for star formation, coupled with lower present-day stellar mass densities. (abridged) We also discuss the importance of the present day stellar mass function (PDMF) in providing a way to disentangle among various assumptions for beta. Our results indicate that the model adopting the IGIMF computed with beta ~2 should be considered the best since it allows us to reproduce the observed PDMF and to account for most of the chemical evolution constraints considered in this work.

preprint2010arXiv

Evolution of two stellar populations in globular clusters II. Effects of primordial gas expulsion

We investigate the early evolution of two distinct populations of low-mass stars in globular clusters under the influence of primordial gas expulsion driven by supernovae to study if this process can increase the fraction of second generation stars at the level required by observations. We analyse N-body models that take into account the effect of primordial gas expulsion. We divide the stars into two populations which mimic the chemical and dynamical properties of stars in globular clusters so that second generation stars start with a more centrally concentrated distribution. The main effect of gas expulsion is to eject preferentially first generation stars while second generation stars remain bound to the cluster. In the most favourable cases second generation stars can account for 60% of the bound stars we see today. We also find that at the end of the gas expulsion phase, the radial distribution of the two populations is still different, so that long-term evolution will further increase the fraction of second generation stars. The large fraction of chemically anomalous stars is readily explainable as a second generation of stars formed out of the slow winds of rapidly rotating massive stars if globular clusters suffer explosive residual gas expulsion for a star formation efficiency of about 0.33.

preprint2010arXiv

Local-Group tests of dark-matter Concordance Cosmology: Towards a new paradigm for structure formation

(abridged) Predictions of the Concordance Cosmological Model (CCM) of the structures in the environment of large spiral galaxies are compared with observed properties of Local Group galaxies. Five new most probably irreconcilable problems are uncovered. However, the Local Group properties provide hints that may lead to a solution of the above problems The DoS and bulge--satellite correlation suggest that dissipational events forming bulges are related to the processes forming phase-space correlated satellite populations. Such events are well known to occur since in galaxy encounters energy and angular momentum are expelled in the form of tidal tails, which can fragment to form populations of tidal-dwarf galaxies (TDGs) and associated star clusters. If Local Group satellite galaxies are to be interpreted as TDGs then the sub-structure predictions of CCM are internally in conflict. All findings thus suggest that the CCM does not account for the Local Group observations and that therefore existing as well as new viable alternatives have to be further explored. These are discussed and natural solutions for the above problems emerge.

preprint2010arXiv

Mass-loss and expansion of ultra compact dwarf galaxies through gas expulsion and stellar evolution for top-heavy stellar initial mass functions

(abridged) The dynamical V-band mass-to-light ratios of ultra compact dwarf galaxies (UCDs) are higher than predicted by simple stellar population models with the canonical stellar initial mass function (IMF). One way to explain this finding is a top-heavy IMF, so that the unseen mass is provided by additional remnants of high-mass stars. A possible explanation for why the IMF in UCDs could be top-heavy while this is not the case in less massive stellar systems is that encounters between proto-stars and stars become probable in forming massive systems. However, the required number of additional stellar remnants proves to be rather high, which raises the question of how their progenitors would affect the early evolution of a UCD. We have therefore calculated the first 200 Myr of the evolution of the UCDs, using the particle-mesh code Superbox. It is assumed that the stellar populations of UCDs were created in an initial starburst, which implies heavy mass loss during the following approximately 40 Myr due to primordial gas expulsion and supernova explosions. We find at the end of the simulations for various initial conditions and (tabulated) mass-loss histories objects that roughly resemble UCDs. Thus, the existence of UCDs does not contradict the notion that their stellar populations formed rapidly and with a top-heavy IMF. We find tentative evidence that the UCDs may have had densities as high as 10^8 M_sun/pc^3 at birth.

preprint2010arXiv

Massive runaway stars in the Large Magellanic Cloud

The origin of massive field stars in the Large Magellanic Cloud (LMC) has long been an enigma. The recent measurements of large offsets (~100 km/s) between the heliocentric radial velocities of some very massive (O2-type) field stars and the systemic LMC velocity provides a possible explanation of this enigma and suggests that the field stars are runaway stars ejected from their birth places at the very beginning of their parent cluster's dynamical evolution. A straightforward way to prove this explanation is to measure the proper motions of the field stars and to show that they are moving away from one of the nearby star clusters or OB associations. This approach however is complicated by the large distance to the LMC, which makes accurate proper motion measurements difficult. We use an alternative approach for solving the problem, based on the search for bow shocks produced by runaway stars. The geometry of detected bow shocks would allow us to infer the direction of stellar motion and thereby to determine their possible parent clusters. In this paper we present the results of a search for bow shocks around six massive field stars which were suggested in the literature as candidate runaway stars. Using archival (Spitzer Space Telescope) data, we found a bow shock associated with one of our program stars, the O2 V((f*)) star BI 237, which is the first-ever detection of bow shocks in the LMC. Orientation of the bow shock suggests that BI 237 was ejected from the OB association LH 82 (located at ~120 pc in projection from the star). A by-product of our search is the detection of bow shocks generated by four OB stars in the field of the LMC and an arc-like structure attached to the candidate luminous blue variable R81 (HD 269128). The geometry of two of these bow shocks is consistent with the possibility that their associated stars were ejected from the 30 Doradus star forming complex.

preprint2010arXiv

Massive runaway stars in the Small Magellanic Cloud

Using archival Spitzer Space Telescope data, we identified for the first time a dozen runaway OB stars in the Small Magellanic Cloud (SMC) through the detection of their bow shocks. The geometry of detected bow shocks allows us to infer the direction of motion of the associated stars and to determine their possible parent clusters and associations. One of the identified runaway stars, AzV 471, was already known as a high-velocity star on the basis of its high peculiar radial velocity, which is offset by ~40 km/s from the local systemic velocity. We discuss implications of our findings for the problem of the origin of field OB stars. Several of the bow shock-producing stars are found in the confines of associations, suggesting that these may be "alien" stars contributing to the age spread observed for some young stellar systems. We also report the discovery of a kidney-shaped nebula attached to the early WN-type star SMC-WR3 (AzV 60a). We interpreted this nebula as an interstellar structure created owing to the interaction between the stellar wind and the ambient interstellar medium.

preprint2010arXiv

Peculiarities in Velocity Dispersion and Surface Density Profiles of Star Clusters

Based on our recent work on tidal tails of star clusters (Kuepper et al. 2009) we investigate star clusters of a few 10^4 Msun by means of velocity dispersion profiles and surface density profiles. We use a comprehensive set of $N$-body computations of star clusters on various orbits within a realistic tidal field to study the evolution of these profiles with time, and ongoing cluster dissolution From the velocity dispersion profiles we find that the population of potential escapers, i.e. energetically unbound stars inside the Jacobi radius, dominates clusters at radii above about 50% of the Jacobi radius. Beyond 70% of the Jacobi radius nearly all stars are energetically unbound. The velocity dispersion therefore significantly deviates from the predictions of simple equilibrium models in this regime. We furthermore argue that for this reason this part of a cluster cannot be used to detect a dark matter halo or deviations from Newtonian gravity. By fitting templates to the about 10^4 computed surface density profiles we estimate the accuracy which can be achieved in reconstructing the Jacobi radius of a cluster in this way. We find that the template of King (1962) works well for extended clusters on nearly circular orbits, but shows significant flaws in the case of eccentric cluster orbits. This we fix by extending this template with 3 more free parameters. Our template can reconstruct the tidal radius over all fitted ranges with an accuracy of about 10%, and is especially useful in the case of cluster data with a wide radial coverage and for clusters showing significant extra-tidal stellar populations. No other template that we have tried can yield comparable results over this range of cluster conditions. All templates fail to reconstruct tidal parameters of concentrated clusters, however. (abridged)

preprint2010arXiv

Properties of hierarchically forming star clusters

We undertake a systematic analysis of the early (< 0.5 Myr) evolution of clustering and the stellar initial mass function in turbulent fragmentation simulations. These large scale simulations for the first time offer the opportunity for a statistical analysis of IMF variations and correlations between stellar properties and cluster richness. The typical evolutionary scenario involves star formation in small-n clusters which then progressively merge; the first stars to form are seeds of massive stars and achieve a headstart in mass acquisition. These massive seeds end up in the cores of clusters and a large fraction of new stars of lower mass is formed in the outer parts of the clusters. The resulting clusters are therefore mass segregated at an age of 0.5 Myr, although the signature of mass segregation is weakened during mergers. We find that the resulting IMF has a smaller exponent (alpha=1.8-2.2) than the Salpeter value (alpha=2.35). The IMFs in subclusters are truncated at masses only somewhat larger than the most massive stars (which depends on the richness of the cluster) and an universal upper mass limit of 150 Msun is ruled out. We also find that the simulations show signs of the IGIMF effect proposed by Weidner & Kroupa, where the frequency of massive stars is suppressed in the integrated IMF compared to the IMF in individual clusters. We identify clusters through the use of a minimum spanning tree algorithm which allows easy comparison between observational survey data and the predictions of turbulent fragmentation models. In particular we present quantitative predictions regarding properties such as cluster morphology, degree of mass segregation, upper slope of the IMF and the relation between cluster richness and maximum stellar mass. [abridged]

preprint2010arXiv

The formation of very wide binaries during the star cluster dissolution phase

Over the past few decades, numerous wide (>1000 au) binaries in the Galactic field and halo have been discovered. Their existence cannot be explained by the process of star formation or by dynamical interactions in the field, and their origin has long been a mystery. We explain the origin of these wide binaries by formation during the dissolution phase of young star clusters: an initially unbound pair of stars may form a binary when their distance in phase-space is small. Using N-body simulations, we find that the resulting wide binary fraction in the semi-major axis range 1000 au - 0.1 pc for individual clusters is 1-30%, depending on the initial conditions. The existence of numerous wide binaries in the field is consistent with observational evidence that most clusters start out with a large degree of substructure. The wide binary fraction decreases strongly with increasing cluster mass, and the semi-major axis of the newly formed binaries is determined by the initial cluster size. The resulting eccentricity distribution is thermal, and the mass ratio distribution is consistent with gravitationally-focused random pairing. As a large fraction of the stars form in primordial binaries, we predict that a large number of the observed 'wide binaries' are in fact triple or quadruple systems. By integrating over the initial cluster mass distribution, we predict a binary fraction of a few per cent in the semi-major axis range 1000 au - 0.1 pc in the Galactic field, which is smaller than the observed wide binary fraction. However, this discrepancy may be solved when we consider a broad range of cluster morphologies.

preprint2009arXiv

Angular momentum transfer and the size - mass relation in early - type galaxies

Early - type galaxies (ETGs) define a narrow strip in the size - mass plane because of the observed correlation between the effective radius $R_{eff}$ and the total stellar mass $M_{\star}$. When expressed in logarithmic units, a linear relation, $\log{R_{eff}} \propto γ\log{M_{\star}}$, is indeed observationally found, but the slope $γ$ deviates from the canonical $γ= 1/2$ value which can be naively predicted for a spherically symmetric isotropic system. We propose here that a transfer of angular momentum to the stellar component induces an anisotropy in the velocity space thus leading to a modified distribution function (DF). Assuming an Osipkov - Merritt like anisotropic DF, we derive an analytic relation between the slope $γ$ of the size - mass relation and the slope $α$ of the angular momentum term in the DF. With this simple model, we are then able to recover the observed $γ$ value provided $α$ is suitably set. It turns out that an anisotropy profile which is tangential inside $\sim 0.6 r_a$ and radial outside, with $r_a$ the anisotropy radius, is able to reproduce the observed size - mass relation observed for massive ($M_{\star} \ge 2 \times 10^{10} \ h^{-1} {\rm M}_{\odot}$) elliptical galaxies.

preprint2009arXiv

Density profiles of dark matter haloes on Galactic and Cluster scales

In the present paper, we improve the "Extended Secondary Infall Model" (ESIM) of Williams et al. (2004) to obtain further insights on the cusp/core problem. The model takes into account the effect of ordered and random angular momentum, dynamical friction and baryon adiabatic contraction in order to obtain a secondary infall model more close to the collapse reality. The model is applied to structures on galactic scales (normal and dwarf spiral galaxies) and on cluster of galaxies scales. The results obtained suggest that angular momentum and dynamical friction are able, on galactic scales, to overcome the competing effect of adiabatic contraction eliminating the cusp. The NFW profile can be reobtained, in our model only if the system is constituted just by dark matter and the magnitude of angular momentum and dynamical friction are reduced with respect to the values predicted by the model itself. The rotation curves of four LSB galaxies from de Blok & Bosma (2002) are compared to the rotation curves obtained by the model in the present paper obtaining a good fit to the observational data. On scales smaller than $\simeq 10^{11} h^{-1} M_{\odot}$ the slope $α\simeq 0$ and on cluster scales we observe a similar evolution of the dark matter density profile but in this case the density profile slope flattens to $α\simeq 0.6$ for a cluster of $\simeq 10^{14} h^{-1} M_{\odot}$. The total mass profile, differently from that of dark matter, shows a central cusp well fitted by a NFW model.

preprint2009arXiv

Faint Fuzzy Star Clusters in NGC1023 as Remnants of Merged Star Cluster Complexes

In the lenticular galaxy NGC1023 a third population of globular clusters (GCs), called faint fuzzies (FFs), was discovered next to the blue and red GC populations by Larsen & Brodie. While these FFs have colors comparable to the red population, the new population is fainter, larger (R_eff > 7 pc) and, most importantly, shows clear signs of co-rotation with the galactic disk of NGC1023. We present N-body simulations verifying the hypothesis that these disk-associated FFs are related to the young massive cluster complexes (CCs) observed by Bastian et. al in M51, who discovered a mass-radius relation for these CCs. Our models have an initial configuration based on the observations from M51 and are placed on various orbits in a galactic potential derived for NGC1023. All computations end up with a stable object containing 10 to 60% of the initial CC mass after an integration time of 5 Gyr. A conversion to visual magnitudes demonstrates that the resulting objects cover exactly the observed range for FFs. Moreover, the simulated objects show projected half-mass radii between 3.6 and 13.4 pc, in good agreement with the observed FF sizes. We conclude that objects like the young massive CCs in M51 are likely progenitors of the FFs observed in NGC1023.

preprint2009arXiv

Merging time-scales of stellar sub-clumps in young star-forming regions

Recent observations and hydrodynamical simulations of star formation inside a giant molecular cloud have revealed that, within a star forming region, stars do not form evenly distributed throughout this region, but rather in small sub-clumps. It is generally believed that these sub-clumps merge and form a young star cluster. The time-scale of this merging process is crucial for the evolution and the possible survival of the final star cluster. The key issue is whether this merging process happens faster than the time needed to remove the residual gas of the cloud. A merging time-scale shorter than the gas-removal time would enhance the survival chances of the resulting star cluster. In this paper we show by means of numerical simulations that the time-scale of the merging is indeed very fast. Depending on the details of the initial sub-clump distribution, the merging may occur before the gas is expelled from the newly-formed cluster either via supernovae or the winds from massive stars. Our simulations further show that the resulting merger-objects have a higher effective star formation efficiency than the overall star forming region and confirm the results that mass-segregated sub-clumps form mass-segregated merger-objects.

preprint2009arXiv

The origin of very wide binary systems

The majority of stars in the Galactic field and halo are part of binary or multiple systems. A significant fraction of these systems have orbital separations in excess of thousands of astronomical units, and systems wider than a parsec have been identified in the Galactic halo. These binary systems cannot have formed through the 'normal' star-formation process, nor by capture processes in the Galactic field. We propose that these wide systems were formed during the dissolution phase of young star clusters. We test this hypothesis using N-body simulations of evolving star clusters and find wide binary fractions of 1-30%, depending on initial conditions. Moreover, given that most stars form as part of a binary system, our theory predicts that a large fraction of the known wide 'binaries' are, in fact, multiple systems.

preprint2009arXiv

The relation between the most-massive star and its parental star cluster mass

We present a thorough literature study of the most-massive star, m_max, in several young star clusters in order to assess whether or not star clusters are populated from the stellar initial mass function (IMF) by random sampling over the mass range 0.01 < m < 150 M_sol without being constrained by the cluster mass, M_ecl. The data reveal a partition of the sample into lowest mass objects (M_ecl < 10^2 M_sol), moderate mass clusters (10^2 M_sol < M_ecl < 10^3 M_sol) and rich clusters above 10^3 M_sol. Additionally, there is a plateau of a constant maximal star mass (m_max ~ 25 M_sol) for clusters with masses between 10^3 M_sol and 4 10^3 M_sol. Statistical tests of this data set reveal that the hypothesis of random sampling from the IMF between 0.01 and 150 M_sol is highly unlikely for star clusters more massive than 10^2 M_sol with a probability of p ~ 2 10^-7 for the objects with M_ecl between 10^2 M_sol and 10^3 M_sol and p ~ 3 10^-9 for the more massive star clusters. Also, the spread of m_max values at a given M_ecl is smaller than expected from random sampling. We suggest that the basic physical process able to explain this dependence of stellar inventory of a star cluster on its mass may be the interplay between stellar feedback and the binding energy of the cluster-forming molecular cloud core. Given these results, it would follow that an integrated galactic initial mass function (IGIMF) sampled from such clusters would automatically be steeper in comparison to the IMF within individual star clusters.

preprint2009arXiv

Tidal Tails of Star Clusters

Based on recent findings of a formation mechanism of substructure in tidal tails by Kuepper, Macleod & Heggie (2008) we investigate a more comprehensive set of N-body models of star clusters on orbits about a Milky-Way-like potential. We find that the predicted epicyclic overdensities arise in any tidal tail no matter which orbit the cluster follows as long as the cluster lives long enough for the overdensities to build up. The distance of the overdensities along the tidal tail from the cluster centre depends for circular orbits only on the mass of the cluster and the strength of the tidal field, and therefore decreases monotonically with time, while for eccentric orbits the orbital motion influences the distance, causing a periodic compression and stretching of the tails and making the distance oscillate with time. We provide an approximation for estimating the distance of the overdensities in this case. We describe an additional type of overdensity which arises in extended tidal tails of clusters on eccentric orbits, when the acceleration of the tidal field on the stellar stream is no longer homogeneous. Moreover, we conclude that a pericentre passage or a disk shock is not the direct origin of an overdensity within a tidal tail. Escape due to such tidal perturbations does not take place immediately after the perturbation but is rather delayed and spread over the orbit of the cluster. All observable overdensities are therefore of the mentioned two types. In particular, we note that substructured tidal tails do not imply the existence of dark-matter sub-structures in the haloes of galaxies.

preprint2008arXiv

The ages of Galactic globular clusters in the context of self-enrichment

A significant fraction of stars in globular clusters (about 70%-85%) exhibit peculiar chemical patterns with strong abundance variations in light elements along with constant abundances in heavy elements. These abundance anomalies can be created in the H-burning core of a first generation of fast rotating massive stars and the corresponding elements are convoyed to the stellar surface thanks to rotational induced mixing. If the rotation of the stars is fast enough this matter is ejected at low velocity through a mechanical wind at the equator. It then pollutes the ISM from which a second generation of chemically anomalous stars can be formed. The proportion of anomalous to normal star observed today depends on at least two quantities : (1) the number of polluter stars; (2) the dynamical history of the cluster which may lose during its lifetime first and second generation stars in different proportions. Here we estimate these proportions based on dynamical models for globular clusters. When internal dynamical evolution and dissolution due to tidal forces are accounted for, starting from an initial fraction of anomalous stars of 10% produces a present day fraction of about 25%, still too small with respect to the observed 70-85%. In case gas expulsion by supernovae is accounted for, much higher fraction is expected to be produced. In this paper we also address the question of the evolution of the second generation stars that are He-rich, and deduce consequences for the age determination of globular clusters.

preprint2006arXiv

A highly abnormal massive star mass function in the Orion Nebula cluster and the dynamical decay of trapezia systems

The ONC appears to be unusual on two grounds: The observed constellation of the OB stars of the entire Orion Nebula cluster and its Trapezium at its centre implies a time-scale problem given the age of the Trapezium, and an IMF problem for the whole OB star population in the ONC. Given the estimated crossing time of the Trapezium, it ought to have totally dynamically decayed by now. Furthermore, by combining the lower limit of the ONC mass with a standard IMF it emerges that the ONC should have formed at least about 40 stars heavier than 5 M_sun while only ten are observed. Using N-body experiments we (i) confirm the expected instability of the trapezium and (ii) show that beginning with a compact OB-star configuration of about 40 stars the number of observed OB stars after 1 Myr within 1 pc radius and a compact trapezium configuration can both be reproduced. These two empirical constraints thus support our estimate of 40 initial OB stars in the cluster. Interestingly, a more-evolved version of the ONC resembles the Upper Scorpius OB association. The N-body experiments are performed with the new C-code CATENA by integrating the equations of motion using the chain-multiple-regularisation method. In addition, we present a new numerical formulation of the initial mass function.

preprint2002arXiv

The impact of mass loss on star cluster formation. I. Analytic results

We study analytically the disruptive effect of instantaneous gas removal from a cluster containing O stars. We setup an iterative calculation based on the stellar velocity distribution function to compute the fraction of stars that remain bound once the cluster has ejected the gas and is out of equilibrium. We show that the stellar bound fraction is a function of the initial cluster distribution function as well as the star formation efficiency, $ε$, taken constant throughout the cluster. The case of the Plummer sphere is dealt with in greater details. We find for this case that up to ~ 50% of the stars may remain bound when $ε$ assumes values < 1/2, contrary to expectations derived from the virial theorem. The fraction of bound stars is expressed algebraically for polytropic distribution functions.