Source author record

Pavel Kroupa

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

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.SR astro-ph gr-qc astro-ph.EP astro-ph.HE hep-ph astro-ph.IM hep-ex physics.comp-ph

Catalog footprint

What is connected

138works

11topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Intermediate-mass black hole incubators. Gas accretion onto stellar black hole clusters in galactic central molecular zones

The stellar dynamical evolution of massive star clusters formed during starburst periods leads to the segregation of $\gtrsim10^4 M_\odot$ stellar-mass black hole sub-clusters in their centres. In gas-rich environments, such as galactic central molecular zones, these black hole clusters are likely to accrete large amounts of the gas from their surroundings, which in turn affects their internal dynamics. In this Letter we estimated the corresponding accretion rate onto the black hole cluster and its radiative feedback. We assessed whether such an accretion flow can lead to the collapse of the black hole cluster into an intermediate-mass black hole. The estimates were obtained analytically, considering the astrophysical conditions and star formation history reported for the central molecular zone of our Galaxy. We find that a stellar black hole cluster with mass $\gtrsim10^4 M_\odot$ located in the twisted ring of molecular clouds with radius $\approx100$ pc that is observed in the central molecular zone of our Galaxy can accrete about the same mass in gas on a timescale of a few million years. We suggest that this is sufficient for its subsequent collapse into an intermediate-mass black hole. Based on an estimate of the dynamical friction inspiral time, we further argue that the locations of the intermediate-mass black hole candidates recently observed in the central molecular zone are compatible with their formation therein during the last starburst period reported to have occurred $\approx1$ Gyr ago.

preprint2025arXiv

Tidal Tails of Nearby Open Clusters. II. A Review of Simulated Properties and the Reliability of Observational Catalogues

Context: Recent studies using Gaia data have reported tidal tail detections for tens to hundreds of open clusters. However, a comprehensive assessment of the reliability & completeness of these detections is lacking. Aims: This work aims to summarise the expected properties of tidal tails based on N-body simulations, review the reliability of tidal tail detections in the literature, and grade them according to a set of diagnostic tests. Methods: We used a grid of 68-20000 Msun simulated clusters & analysed the formation & evolution of the tails. We compiled 122 catalogues (58 unique clusters) from literature, within 500 pc of the Sun. We employed tests based on photometric, morphological & dynamical signatures, & comparisons with simulations, to grade their tails. Results: Based on the simulations, we analysed the complex morphology of the tails & their properties (e.g., size, span, stellar types, number density, & mass function) at various cluster masses & ages. During the first 200 Myr of evolution, the tails typically form a characteristic S shape, with an amplitude that scales with cluster mass. The tail span increases at a rate of 4 times the initial velocity dispersion, & the near-tail is predominantly populated by recent escapees. Conclusions: In the 122 published tail catalogues, we found that 15 gold & 55 silver catalogues passed the majority of the tests. The remaining 51 were graded as bronze; care should be taken before using these for further analysis. The age, metallicity, binary fraction, & mass function of stars in the tails were generally consistent with those of their parent clusters. The gold/silver catalogues (69 catalogues of 40 clusters) represent reliable samples for detailed analyses of tails. Future data will be essential for further validation & for leveraging tidal tails as tracers of cluster dissolution & the Galactic potential.

preprint2022arXiv

3D hydrodynamic simulations for the formation of the Local Group satellite planes

The existence of mutually correlated thin and rotating planes of satellite galaxies around both the Milky Way (MW) and Andromeda (M31) calls for an explanation. Previous work in Milgromian dynamics (MOND) indicated that a past MW-M31 encounter might have led to the formation of these satellite planes. We perform the first-ever hydrodynamical MOND simulation of the Local Group using Phantom of RAMSES. We show that an MW-M31 encounter at $z \approx 1$, with a perigalactic distance of about 80 kpc, can yield two disc galaxies at $z=0$ oriented similarly to the observed galactic discs and separated similarly to the observed M31 distance. Importantly, the tidal debris are distributed in phase space similarly to the observed MW and M31 satellite planes, with the correct preferred orbital pole for both. The MW-M31 orbital geometry is consistent with the presently observed M31 proper motion despite this not being considered as a constraint when exploring the parameter space. The mass of the tidal debris around the MW and M31 at $z=0$ compare well with the mass observed in their satellite systems. The remnant discs of the two galaxies have realistic radial scale lengths and velocity dispersions, and the simulation naturally produces a much hotter stellar disc in M31 than in the MW. However, reconciling this scenario with the ages of stellar populations in satellite galaxies would require that a higher fraction of stars previously formed in the outskirts of the progenitors ended up within the tidal debris, or that the MW-M31 interaction occurred at $z>1$.

preprint2022arXiv

Do the majority of stars form as gravitationally unbound?

Some of the youngest stars (age $\lesssim 10$ Myr) are clustered, while many others are observed scattered throughout star forming regions or in complete isolation. It has been intensively debated whether the scattered or isolated stars originate in star clusters, or if they form truly isolated, which could help constrain the possibilities how massive stars are formed. We adopt the assumption that all stars form in gravitationally bound star clusters embedded in molecular cloud cores ($Γ$-$1\;$ model), which expel their natal gas, and compare the fraction of stars found in clusters with observational data. The star clusters are modelled by the code nbody6, which includes stellar and circumbinary evolution, gas expulsion, and the external gravitational field of their host galaxy. We find that small changes in the assumptions in the current theoretical model estimating the fraction, $Γ$, of stars forming in embedded clusters have a large influence on the results, and we present a counterexample as an illustration. This calls into question theoretical arguments about $Γ$ in embedded clusters, and it suggests that there is no firm theoretical ground for low $Γ$ in galaxies with lower star formation rates (SFRs). Instead, the assumption that all stars form in embedded clusters is in agreement with observational data for the youngest stars (age $\lesssim 10$ Myr). In the $Γ$-$1\;$ scenario, the observed fraction of the youngest stars in clusters increases with the SFR only weakly; the increase is caused by the presence of more massive clusters in galaxies with higher SFRs, which release fewer stars to the field in proportion to their mass. The $Γ$-$1\;$ model yields a higher fraction of stars in clusters for older stars (age between $10$ and $300$ Myr) than what is observed. This discrepancy can be caused by interactions with molecular clouds.

preprint2022arXiv

Far-ultraviolet investigation into the galactic globular cluster M30 (NGC 7099): I. Photometry and radial distributions

We present a far-ultraviolet (FUV) study of the globular cluster M30 (NGC 7099). The images were obtained using the Advanced Camera for Surveys (ACS/SBC, F150LP, FUV) and the Wide Field Planetary Camera 2 (WFPC2, F300W, UV) which were both on board the Hubble Space Telescope (HST). The FUV-UV colour-magnitude diagram (CMD) shows a main sequence (MS) turnoff at FUV $\approx$ 22 mag and FUV-UV $\approx$ 3 mag. The MS extends 4 mag below the turnoff, and a prominent horizontal branch (HB) and blue straggler (BS) sequence can be seen. A total of 1218 MS stars, 185 red giant branch stars, 47 BS stars and 41 HB stars are identified, along with 78 sources blueward of the MS which consist of white dwarfs (WDs) and objects in the gap between the WDs and the MS that include potential cataclysmic variable (CV) candidates. The radial distribution of the BS population is concentrated towards the cluster centre, indicating that mass segregation has occurred. The blue and red sub-populations of the double BS sequence appear mixed in the ultraviolet CMD, and no significant central concentration of CV candidates is seen in this cluster.

preprint2022arXiv

Far-ultraviolet investigation into the galactic globular cluster M30 (NGC 7099): II. Potential X-ray counterparts and variable sources

We present a far-ultraviolet (FUV) study of the globular cluster M30 (NGC 7099). The images were obtained using the Advanced Camera for Surveys (ACS/SBC, F150LP, FUV) and the Wide Field Planetary Camera 2 (WFPC2, F300W, UV) on board the Hubble Space Telescope (HST). We compare the catalogue of FUV objects to ten known X-ray sources and find six confident matches of two cataclysmic variables (CVs), one RS CVn, one red giant with strong FUV emission and two sources only detected in the FUV. We also searched for variable sources in our dataset and found a total of seven blue stragglers (BSs), four horizontal branch (HB) stars, five red giant branch stars, 28 main sequence stars and four gap objects that demonstrated variability. One BS star is a known W-UMa contact binary, one of the gap objects is a known CV identified in this work to be a dwarf nova, and the three other gap sources are weak variables. The periods and positions of two of the variable HB stars match them to two previously known RR Lyrae variables of types RRab and RRc.

preprint2022arXiv

Origin of the spectacular tidal shells of galaxy NGC474

The lenticular galaxy NGC474 hosts a rich system of tidal shells and streams, some of which are exceptionally bright. Two teams recently presented spectroscopic observations of the brightest shells. These were the first shell spectra ever observed in integrated starlight. The authors studied the stellar populations of the shell, of the center of the galaxy and of its globular clusters. The precise formation scenario for the tidal features of this prominent galaxy however still remained unclear. Here, we add further clues on their formation from the radii of the shells, and we present a scenario for the formation of the tidal features that seems to be unique and explaining all available data. Shell radii are analyzed with the shell identification method, and we run self-consistent simulations of the formation of the tidal features. We consider Newtonian as well as MOND gravity. Observations suggest that the tidal features originate from the accretion of a spiral galaxy. The shell identification method yields that the merging galaxies collided first 1.3Gyr ago and then again 0.9Gyr ago, thereby forming the shells in two generations. This would also explain the young ages of stellar populations in the center of the galaxy and the young age of the globular clusters. The analytic models of shell propagation, that underlie the shell identification method, are verified by a simulation. The simulations reproduce well the observed morphology of the tidal features. The accreted spiral likely reached NGC474 nearly radially, in the plane of the sky, from the south, its rotation axis pointing toward us. It should have had a stellar mass of around 1/6 of NGC474, i.e. $10^{9.8}\,M_\odot$. It seems that all tidal features in the galaxy originate from one merger.

preprint2022arXiv

Overestimated inclinations of Milgromian disc galaxies: the case of the ultradiffuse galaxy AGC 114905

We present hydrodynamical star-forming simulations in the Milgromian dynamics (MOND) framework of a gas-rich disc galaxy with properties similar to AGC 114905, which has recently been argued to have a rotation curve (RC) that is inconsistent with the MOND prediction. Our first model considers the galaxy in isolation, while our second model includes an external field of $0.05 \, a_{_0}$, the estimated gravitational field from large-scale structure. We show that isophotes in the face-on view can differ from circular at the 50% level. This could mislead observers into overestimating the inclination $i$ between disc and sky planes. Because RCs require a correction factor of $1/\sin i$, the actual RC could be much higher than reported by observers. This plausibly reconciles AGC 114905 with MOND expectations.

preprint2022arXiv

The formation of early-type galaxies through monolithic collapse of gas clouds in Milgromian gravity

Studies of stellar populations in early-type galaxies (ETGs) show that the more massive galaxies form earlier and have a shorter star formation history (SFH). In this study, we investigate the initial conditions of ETG formation. The study begins with the collapse of non-rotating post-Big-Bang gas clouds in Milgromian (MOND) gravitation. These produce ETGs with star-forming timescales (SFT) comparable to those observed in the real Universe. Comparing these collapse models with observations, we set constraints on the initial size and density of the post-Big-Bang gas clouds in order to form ETGs. The effective-radius-mass relation of the model galaxies falls short of the observed relation. Possible mechanisms for later radius expansion are discussed. Using hydrodynamic MOND simulations this work thus for the first time shows that the SFTs observed for ETGs may be a natural occurrence in the MOND paradigm. We show that different feedback algorithms change the evolution of the galaxies only to a very minor degree in MOND. The first stars have, however, formed more rapidly in the real Universe than possible just from the here studied gravitational collapse mechanism. Dark-matter-based cosmological structure formation simulations disagree with the observed SFTs at more than 5 sigma confidence.

preprint2022arXiv

The global stability of M33 in MOND

The dynamical stability of disk galaxies is sensitive to whether their anomalous rotation curves are caused by dark matter halos or Milgromian Dynamics (MOND). We investigate this by setting up a MOND model of M33. We first simulate it in isolation for 6 Gyr, starting from an initial good match to the rotation curve (RC). Too large a bar and bulge form when the gas is too hot, but this is avoided by reducing the gas temperature. A strong bar still forms in 1 Gyr, but rapidly weakens and becomes consistent with the observed weak bar. Previous work showed this to be challenging in Newtonian models with a live dark matter halo, which developed strong bars. The bar pattern speed implies a realistic corotation radius of 3 kpc. However, the RC still rises too steeply, and the central line of sight velocity dispersion (LOSVD) is too high. We then add a constant external acceleration field of $8.4 \times 10^{-12}$ m/s$^2$ at $30^\circ$ to the disk as a first order estimate for the gravity exerted by M31. This suppresses buildup of material at the centre, causing the RC to rise more slowly and reducing the central LOSVD. Overall, this simulation bears good resemblance to several global properties of M33, and highlights the importance of including even a weak external field on the stability and evolution of disk galaxies. Further simulations with a time-varying external field, modeling the full orbit of M33, will be needed to confirm its resemblance to observations.

preprint2022arXiv

The High Fraction of Thin Disk Galaxies Continues to Challenge ΛCDM Cosmology

Any viable cosmological framework has to match the observed proportion of early- and late-type galaxies. In this contribution, we focus on the distribution of galaxy morphological types in the standard model of cosmology (Lambda cold dark matter, $Λ$CDM). Using the latest state-of-the-art cosmological $Λ$CDM simulations known as Illustris, IllustrisTNG, and EAGLE, we calculate the intrinsic and sky-projected aspect ratio distribution of the stars in subhalos with stellar mass $M_* > 10^{10}\,M_\odot$ at redshift $z=0$. There is a significant deficit of intrinsically thin disk galaxies, which however comprise most of the locally observed galaxy population. Consequently, the sky-projected aspect ratio distribution produced by these $Λ$CDM simulations disagrees with the Galaxy And Mass Assembly (GAMA) survey and Sloan Digital Sky Survey at $\geq 12.52σ$ (TNG50-1) and $\geq 14.82σ$ (EAGLE50) confidence. The deficit of intrinsically thin galaxies could be due to a much less hierarchical merger-driven build-up of observed galaxies than is given by the $Λ$CDM framework. It might also arise from the implemented sub-grid models, or from the limited resolution of the above-mentioned hydrodynamical simulations. We estimate that an $8^5$ times better mass resolution realization than TNG50-1 would reduce the tension with GAMA to the $5.58σ$ level. Finally, we show that galaxies with fewer major mergers have a somewhat thinner aspect ratio distribution. Given also the high expected frequency of minor mergers in $Λ$CDM, the problem may be due to minor mergers. In this case, the angular momentum problem could be alleviated in Milgromian dynamics (MOND) because of a reduced merger frequency arising from the absence of dynamical friction between extended dark matter halos.

preprint2021arXiv

Estimating the ages of open star clusters from properties of their extended tidal tails

The most accurate current methods for determining the ages of open star clusters, stellar associations and stellar streams are based on isochrone fitting or the lithium depletion boundary. We propose another method for dating these objects based on the morphology of their extended tidal tails, which have been recently discovered around several open star clusters. Assuming that the early-appearing tidal tails, the so called tidal tails I, originate from the stars released from the cluster during early gas expulsion, or that they form in the same star forming region as the cluster (i.e. being coeval with the cluster), we derive the analytical formula for the tilt angle $β$ between the long axis of the tidal tail and the orbital direction for clusters or streams on circular trajectories. Since at a given Galactocentric radius, $β$ is only a function of age $t$ regardless of the initial properties of the cluster, we estimate the cluster age by inverting the analytical formula $β= β(t)$. We illustrate the method on a sample of $12$ objects, which we compiled from the literature, and we find a reasonable agreement with previous dating methods in $\approx 70$% of the cases. This can probably be improved by taking into account the eccentricity of the orbits and by revisiting the dating methods based on stellar evolution. The proposed morphological method is suitable for relatively young clusters (age $\lesssim 300$ Myr), where it provides a relative age error of the order of $10$ to $20$% for an error in the observed tilt angle of $5 ^\circ$.

preprint2021arXiv

HI mapping of the Leo Triplet: Morphologies and kinematics of tails and bridges

A fully-sampled and hitherto highest resolution and sensitivity observation of neutral hydrogen (HI) in the Leo Triplet (NGC 3628, M 65/NGC 3623, and M 66/NGC 3627) reveals six HI structures beyond the three galaxies. We present detailed results of the morphologies and kinematics of these structures, which can be used for future simulations. In particular, we detect a two-arm structure in the plume of NGC 3628 for the first time, which can be explained by a tidal interaction model. The optical counterpart of the plume is mainly associated with the southern arm. The connecting part (base) of the plume (directed eastwards) with NGC 3628 is located at the blueshifted (western) side of NGC 3628. Two bases appear to be associated with the two arms of the plume. A clump with reversed velocity gradient (relative to the velocity gradient of M 66) and a newly detected tail, i.e. M 66SE, is found in the southeast of M 66. We suspect that M 66SE represents gas from NGC 3628 which was captured by M 66 in the recent interaction between the two galaxies. Meanwhile gas is falling toward M 66, resulting in features already previously observed in the southeastern part of M 66, e.g. large line widths and double peaks. An upside-down `Y'-shaped HI gas component (M 65S) is detected in the south of M 65 which suggests that M 65 may also have been involved in the interaction. We strongly encourage modern hydrodynamical simulations of this interacting group of galaxies to reveal the origin of the gaseous debris surrounding all three galaxies.

preprint2021arXiv

On the dynamical evolution of Cepheids in star clusters

We investigate the occurrence of classical (type-I) Cepheid variable stars (henceforth: Cepheids) in dynamically evolving star clusters from birth to an age of approximately 300 Myr. The clusters are modelled by the Aarseth code nbody6, and they feature a realistic stellar initial mass function and initial binary star population, single star and binary star evolution, expulsion of the primordial gas, and the tidal field of the galaxy. Our simulations provide the first detailed dynamical picture of how frequently Cepheids remain gravitationally bound to their birth clusters versus how frequently they occur in the field. They allow us to quantify the relevance of various cluster escape mechanisms and how they depend on stellar mass. Overall, the simulations agree with the empirical picture that a small fraction ($\approx 10\%$) of Cepheids reside in clusters, that cluster halo membership is relatively common, and that the majority of Cepheid hosting clusters have only a single Cepheid member. Additionally, the simulations predict that a) Cepheid progenitors are much more likely to escape from low-mass than higher-mass clusters; b) higher-mass (long-period) Cepheids are $\approx 30\%$ more likely to be found in clusters than low-mass (short-period) Cepheids; c) the clustered Cepheid fraction increases with galactocentric radius since cluster dispersal is less efficient at greater radii; d) a lower metallicity reduces the overall clustered Cepheid fraction; e) high-mass clusters are much more likely to have more than one Cepheid member at any given time, in particular at a lower metallicity. We interpret the results as outcomes of various aspects of star cluster dynamics. The comparison of predicted and observed clustered Cepheid fractions, $f_{\rm CC}$, highlights the need for additional cluster disruption mechanisms, most likely encounters with giant molecular clouds.

preprint2020arXiv

A correlation between the number of satellites and the baryonic bulge-to-total ratio extending beyond the Local Group

Recent observations of the fields surrounding a few Milky-Way-like galaxies in the local Universe have become deep enough to enable investigations of the predictions of the standard LCDM cosmological model down to small scales outside the Local Group. Motivated by an observed correlation between the number of dwarf satellites (N_sat) and the bulge-to-total baryonic mass ratios (B/T) of the three main galaxies in the Local Group, i.e. the Milky Way, Andromeda, and Triangulum (M33), we use published data of three well-studied galaxies outside the Local Group, namely M81, Centaurus A, and M101, and their confirmed satellites, and we find a strong and significant correlation between N_sat and B/T. This presents itself in contradiction with the hitherto published results from cosmological simulations reporting an absence of a correlation between N_sat and B/T in the LCDM model. We conclude that, based on the current data, the N_sat vs. B/T correlation is no longer a property confined to only the Local Group.

preprint2020arXiv

Chemical evolution of ultra-faint dwarf galaxies in the self-consistently calculated IGIMF theory

The galaxy-wide stellar initial mass function (gwIMF) of a galaxy in dependence of its metallicity and star formation rate (SFR) can be calculated by the integrated galactic IMF (IGIMF) theory. Lacchin et al. (2019) apply the IGIMF theory for the first time to study the chemical evolution of the ultra-faint dwarf (UFD) satellite galaxies and failed to reproduce the data. Here, we find that the IGIMF theory is naturally consistent with the data. We apply the time-evolving gwIMF calculated at each timestep. The number of type Ia supernova explosions per unit stellar mass formed is renormalised according to the gwIMF. The chemical evolution of Boötes I, one of the best observed UFD, is calculated. Our calculation suggests a mildly bottom-light and top-light gwIMF for Boötes I, and that this UFD has the same gas-consumption timescale as other dwarfs but was quenched about 0.1 Gyr after formation, being consistent with independent estimations and similar to Dragonfly 44. The recovered best fitting input parameters in this work are not covered in the work of Lacchin et al. (2019), creating the discrepancy between our conclusions. In addition, a detailed discussion of uncertainties is presented addressing how the results of chemical evolution models depend on applied assumptions. This study demonstrates the power of the IGIMF theory in understanding the star-formation in extreme environments and shows that UDFs are a promising pathway to constrain the variation of the low-mass stellar IMF.

preprint2020arXiv

Constraints on the star formation histories of galaxies in the Local Cosmological Volume

The majority of galaxies with current star-formation rates (SFRs), SFRo >= 10^-3 Msun/yr, in the Local Cosmological Volume where observations should be reliable, have the property that their observed SFRo is larger than their average star formation rate. This is in tension with the evolution of galaxies described by delayed-tau models, according to which the opposite would be expected. The tension is apparent in that local galaxies imply the star formation timescale tau approx 6.7 Gyr, much longer than the 3.5-4.5 Gyr obtained using an empirically determined main sequence at several redshifts. Using models where the SFR is a power law in time of the form propto (t - t1)^eta for t1 = 1.8 Gyr (with no stars forming prior to t1) implies that eta = 0.18 +- 0.03. This suggested near-constancy of a galaxy's SFR over time raises non-trivial problems for the evolution and formation time of galaxies, but is broadly consistent with the observed decreasing main sequence with increasing age of the Universe.

preprint2020arXiv

Scale-invariant dynamics in the Solar System

The covariant scale-invariant dynamics (SID) theory has recently been proposed as a possible explanation for the observed dynamical discrepancies in galaxies (Maeder & Gueorguiev 2020). SID implies that these discrepancies $-$ commonly attributed to dark matter $-$ arise instead from a non-standard velocity-dependent force that causes two-body near-Keplerian orbits to expand. We show that the predicted expansion of the Earth-Moon orbit is incompatible with lunar laser ranging data at $>200σ$. Moreover, SID predicts that the gravitating mass of any object was much smaller in the past. If true, a low-mass red giant star must be significantly older than in standard theory. This would make it much older than the conventional age of the Universe, which however is expected to be similarly old in SID. Moreover, it is not completely clear whether SID truly contains new physics beyond General Relativity, with several previous works arguing that the extra degree of freedom is purely mathematical. We conclude that the SID model is falsified at high significance by observations across a range of scales, even if it is theoretically well formulated.

preprint2020arXiv

Solar System limits on gravitational dipoles

The gravitational dipole theory of Hadjukovic (2010) is based on the hypothesis that antimatter has a negative gravitational mass and thus falls upwards on Earth. Astrophysically, the model is similar to but more fundamental than Modified Newtonian Dynamics (MOND), with the Newtonian gravity $g_{_N}$ towards an isolated point mass boosted by the factor $ν= 1 + \left( α/x \right) \tanh \left( \sqrt{x}/α\right)$, where $x \equiv g_{_N}/a_{_0}$ and $a_{_0} = 1.2 \times 10^{-10}$ m/s$^2$ is the MOND acceleration constant. We show that $α$ must lie in the range ${0.4-1}$ to acceptably fit galaxy rotation curves. In the Solar System, this interpolating function implies an extra Sunwards acceleration of ${αa_{_0}}$. This would cause Saturn to deviate from Newtonian expectations by ${7000 \left( α/0.4 \right)}$ km over 15 years, starting from known initial position and velocity on a near-circular orbit. We demonstrate that this prediction should not be significantly altered by the postulated dipole haloes of other planets due to the rather small region in which each planet's gravity dominates over that of the Sun. The orbit of Saturn should similarly be little affected by a possible ninth planet in the outer Solar System and by the Galactic gravity causing a non-spherical distribution of gravitational dipoles several kAU from the Sun. Radio tracking of the Cassini spacecraft orbiting Saturn yields a ${5σ}$ upper limit of 160 metres on deviations from its conventionally calculated trajectory. These measurements imply a much more stringent upper limit on $α$ than the minimum required for consistency with rotation curve data. Therefore, no value of $α$ can simultaneously match all available constraints, falsifying the gravitational dipole theory in its current form at extremely high significance.

preprint2020arXiv

The failure of testing for cosmic opacity via the distance-duality relation

The distance-duality relation (DDR) between the luminosity distance $D_L$ and the angular diameter distance $D_A$ is viewed as a powerful tool for testing for the opacity of the Universe, being independent of any cosmological model. It was applied by many authors, who mostly confirm its validity and report a negligible opacity of the Universe. Nevertheless, a thorough analysis reveals that applying the DDR in cosmic opacity tests is tricky. Its applicability is strongly limited because of a non-unique interpretation of the $D_L$ data in terms of cosmic opacity and a rather low accuracy and deficient extent of currently available $D_A$ data. Moreover, authors usually assume that cosmic opacity is frequency independent and parametrize it in their tests by a prescribed phenomenological function. In this way, they only prove that cosmic opacity does not follow their assumptions. As a consequence, no convincing evidence of transparency of the universe using the DDR has so far been presented.

preprint2020arXiv

The formation of exponential disk galaxies in MOND

The formation and evolution of galaxies is highly dependent on the dynamics of stars and gas, which is governed by the underlying law of gravity. To investigate how the formation and evolution of galaxies takes place in Milgromian gravity (MOND), we present full hydrodynamical simulations with the Phantom of Ramses (POR) code. These are the first-ever galaxy formation simulations done in MOND with detailed hydrodynamics, including star formation, stellar feedback, radiative transfer and supernovae. These models start from simplified initial conditions, in the form of isolated, rotating gas spheres in the early Universe. These collapse and form late-type galaxies obeying several scaling relations, which was not a priori expected. The formed galaxies have a compact bulge and a disk with exponentially decreasing surface mass density profiles and scale lengths consistent with observed galaxies, and vertical stellar mass distributions with distinct exponential profiles (thin and thick disk). This work thus shows for the first time that disk galaxies with exponential profiles in both gas and stars are a generic outcome of collapsing gas clouds in MOND. These models have a slight lack of stellar angular momentum because of their somewhat compact stellar bulge, which is connected to the simple initial conditions and the negligible later gas accretion. We also analyse how the addition of more complex baryonic physics changes the main resulting properties of the models and find this to be negligibly so in the Milgromian framework.

preprint2020arXiv

Tidal tails of open star clusters as probes to early gas expulsion I: A semi-analytic model

Star clusters form out of the densest parts of infrared dark clouds. The emergence of massive stars expels the residual gas, which has not formed stars yet. Gas expulsion lowers the gravitational potential of the embedded cluster, unbinding many of the cluster stars. These stars then move on their own trajectories in the external gravitational field of the Galaxy, forming a tidal tail. We investigate the formation and evolution of the tidal tail forming due to expulsion of primordial gas under various scenarios of gas expulsion to provide predictions for tidal tails around dynamically evolved (age > 100 Myr) galactic star clusters, which can be possibly detected by the Gaia mission. We provide a semi-analytical model for the tail evolution. We find that tidal tails released during gas expulsion have different kinematic properties than the tails gradually forming due to evaporation. The gas expulsion tidal tail shows non-monotonic expansion with time, where longer epochs of expansion are interspersed with shorter epochs of contraction. The tail thickness and velocity dispersions strongly, but not exactly periodically, vary with time. The times of minima of tail thickness and velocity dispersions are given only by the properties of the galactic potential, and not by the properties of the cluster. The estimates provided by the (semi-)analytical model for the extent of the tail, the minima of tail thickness, and velocity dispersions are in a very good agreement with the nbody6 simulations. This implies that the semi-analytic model can be used for estimating the properties of the gas expulsion tidal tail for a cluster of a given age and orbital parameters without the necessity of performing numerical simulations. A study with a more extended parameter space of the initial conditions is performed in the follow up paper.

preprint2020arXiv

Tidal tails of open star clusters as probes to early gas expulsion II: Predictions for Gaia

We study the formation and evolution of the tidal tail released from a young Pleiades-like star cluster, due to expulsion of primordial gas in a realistic gravitational field of the Galaxy. The tidal tails (as well as clusters) are integrated from their embedded phase for 300 Myr. We vary star formation efficiencies (SFEs) from 33% to 100% and the timescales of gas expulsion as free parameters, and provide predictions for the morphology and kinematics of the evolved tail for each of the models. The resulting tail properties are intended for comparison with anticipated Gaia observations in order to constrain the poorly understood early conditions during the gas phase and gas expulsion. The simulations are performed with the code Nbody6 including a realistic external gravitational potential of the Galaxy, and an analytical approximation for the natal gaseous potential. Assuming that the Pleiades formed with rapid gas expulsion and an SFE of $\approx 30$%, the current Pleiades are surrounded by a rich tail extending from $\approx 150$ to $\approx 350$ pc from the cluster and containing 0.7 to 2.7 times the number of stars in the present-day cluster. If the Pleiades formed with an SFE close to 100%, then the tail is shorter ($\lesssim 90$ pc) and substantially poorer, containing only $\approx 2$% of the number of present-day cluster stars. If the Pleiades formed with an SFE of $\approx 30$%, but the gas expulsion was adiabatic, the tail signatures are indistinguishable from the case of the model with 100% SFE. The model takes into account the estimated contamination due to the field stars and the Hyades-Pleiades stream, which constitutes a more limiting factor than the accuracy of the Gaia measurements.

preprint2020arXiv

Very high redshift quasars and the rapid emergence of super-massive black holes

The observation of quasars at very high redshift such as Poniuaena is a challenge for models of super-massive black hole (SMBH) formation. This work presents a study of SMBH formation via known physical processes in star-burst clusters formed at the onset of the formation of their hosting galaxy. While at the early stages hyper-massive star-burst clusters reach the luminosities of quasars, once their massive stars die, the ensuing gas accretion from the still forming host galaxy compresses its stellar black hole (BH) component to a compact state overcoming heating from the BH--BH binaries such that the cluster collapses, forming a massive SMBH-seed within about a hundred Myr. Within this scenario the SMBH--spheroid correlation emerges near-to-exactly. The highest-redshift quasars may thus be hyper-massive star-burst clusters or young ultra-compact dwarf galaxies (UCDs), being the precursors of the SMBHs that form therein within about 200 Myr of the first stars. For spheroid masses <10^9.6 Msun a SMBH cannot form and instead only the accumulated nuclear cluster remains. The number evolution of the quasar phases with redshift is calculated and the possible problem of missing quasars at very high redshift is raised. SMBH-bearing UCDs and the formation of spheroids are discussed critically in view of the high redshift observations. A possible tension is found between the high star-formation rates (SFRs) implied by downsizing and the observed SFRs, which may be alleviated within the IGIMF theory and if the downsizing times are somewhat longer.

preprint2019arXiv

The Milky Way's Disk of Classical Satellite Galaxies in Light of Gaia DR2

We study the correlation of orbital poles of the 11 classical satellite galaxies of the Milky Way, comparing results from previous proper motions with the independent data by Gaia DR2. Previous results on the degree of correlation and its significance are confirmed by the new data. A majority of the satellites co-orbit along the Vast Polar Structure, the plane (or disk) of satellite galaxies defined by their positions. The orbital planes of eight satellites align to $<20^\circ$ with a common direction, seven even orbit in the same sense. Most also share similar specific angular momenta, though their wide distribution on the sky does not support a recent group infall or satellites-of-satellites origin. The orbital pole concentration has continuously increased as more precise proper motions were measured, as expected if the underlying distribution shows true correlation that is washed out by observational uncertainties. The orbital poles of the up to seven most correlated satellites are in fact almost as concentrated as expected for the best-possible orbital alignment achievable given the satellite positions. Combining the best-available proper motions substantially increases the tension with $Λ$CDM cosmological expectations: <0.1 per cent of simulated satellite systems in IllustrisTNG contain seven orbital poles as closely aligned as observed. Simulated systems that simultaneously reproduce the concentration of orbital poles and the flattening of the satellite distribution have a frequency of <0.1 per cent for any number of k > 3 combined orbital poles, indicating that these results are not affected by a look-elsewhere effect. This compounds the Planes of Satellite Galaxies Problem.

preprint2019arXiv

The systematically varying stellar IMF

Some ultra-compact dwarf galaxies have large dynamical mass to light (M/L) ratios and also appear to contain an overabundance of LMXB sources, and some Milky Way globular clusters have a low concentration and appear to have a deficit of low-mass stars. These observations can be explained if the stellar IMF becomes increasingly top-heavy with decreasing metallicity and increasing gas density of the forming object. The thus constrained stellar IMF then accounts for the observed trend of metallicity and M/L ratio found amongst M31 globular star clusters. It also accounts for the overall shift of the observationally deduced galaxy-wide IMF from top-light to top-heavy with increasing star formation rate amongst galaxies. If the IMF varies similarly to deduced here, then extremely young very massive star-burst clusters observed at a high redshift would appear quasar-like (Jerabkova et al. 2017) .

preprint2019arXiv

The ultra-diffuse dwarf galaxies NGC 1052-DF2 and 1052-DF4 are in conflict with standard cosmology

Recently van Dokkum et al. (2018b) reported that the galaxy NGC 1052-DF2 (DF2) lacks dark matter if located at $20$ Mpc from Earth. In contrast, DF2 is a dark-matter-dominated dwarf galaxy with a normal globular cluster population if it has a much shorter distance near $10$ Mpc. However, DF2 then has a high peculiar velocity wrt. the cosmic microwave background of $886$ $\rm{km\,s^{-1}}$, which differs from that of the Local Group (LG) velocity vector by $1298$ $\rm{km\,s^{-1}}$ with an angle of $117 \, ^{\circ}$. Taking into account the dynamical $M/L$ ratio, the stellar mass, half-light radius, peculiar velocity, motion relative to the LG, and the luminosities of the globular clusters, we show that the probability of finding DF2-like galaxies in the lambda cold dark matter ($Λ$CDM) TNG100-1 simulation is at most $1.0\times10^{-4}$ at $11.5$ Mpc and is $4.8\times10^{-7}$ at $20.0$ Mpc. At $11.5$ Mpc, the peculiar velocity is in significant tension with the TNG100-1, TNG300-1, and Millennium simulations, but occurs naturally in a Milgromian cosmology. At $20.0$ Mpc, the unusual globular cluster population would challenge any cosmological model. Estimating that precise measurements of the internal velocity dispersion, stellar mass, and distance exist for $100$ galaxies, DF2 is in $2.6σ$ ($11.5$ Mpc) and $4.1σ$ ($20.0$ Mpc) tension with standard cosmology. Adopting the former distance for DF2 and assuming that NGC 1052-DF4 is at $20.0$ Mpc, the existence of both is in tension at $\geq4.8σ$ with the $Λ$CDM model. If both galaxies are at $20.0$ Mpc the $Λ$CDM cosmology has to be rejected by $\geq5.8σ$.

preprint2016arXiv

Declining rotation curves of galaxies as a test of gravitational theory

Unlike Newtonian dynamics which is linear and obeys the strong equivalence principle, in any nonlinear gravitation such as Milgromian dynamics (MOND), the strong version of the equivalence principle is violated and the gravitational dynamics of a system is influenced by the external gravitational field in which it is embedded. This so called External Field Effect (EFE) is one of the important implications of MOND and provides a special context to test Milgromian dynamics. Here, we study the rotation curves (RCs) of 18 spiral galaxies and find that their shapes constrain the EFE. We show that the EFE can successfully remedy the overestimation of rotation velocities in 80\% of the sample galaxies in Milgromian dynamics fits by decreasing the velocity in the outer part of the RCs. We compare the implied external field with the gravitational field for non-negligible nearby sources of each individual galaxy and find that in many cases it is compatible with the EFE within the uncertainties. We therefore argue that in the framework of Milgromian dynamics, one can constrain the gravitational field induced from the environment of galaxies using their RCs. We finally show that taking into account the EFE yields more realistic values for the stellar mass-to-light ratio in terms of stellar population synthesis than the ones implied without the EFE.

preprint2016arXiv

Dynamical ejections of massive stars from young star clusters under diverse initial conditions

We study the effects of initial conditions of star clusters and their massive star population on dynamical ejections of massive stars from star clusters up to an age of 3 Myr. We use a large set of direct N-body calculations for moderately massive star clusters (Mecl=$10^{3.5}$ Msun). We vary the initial conditions of the calculations such as the initial half-mass radius of the clusters, initial binary populations for massive stars and initial mass segregation. We find that the initial density is the most influential parameter for the ejection fraction of the massive systems. The clusters with an initial half-mass radius of 0.1 (0.3) pc can eject up to 50% (30%) of their O-star systems on average. Most of the models show that the average ejection fraction decreases with decreasing stellar mass. For clusters efficient at ejecting O stars, the mass function of the ejected stars is top-heavy compared to the given initial mass function (IMF), while the mass function of stars that remain in the cluster becomes slightly steeper (top-light) than the IMF. The top-light mass functions of stars in 3 Myr old clusters in our N-body models agree well with the mean mass function of young intermediate-mass clusters in M31, as reported previously. We show that the multiplicity fraction of the ejected massive stars can be as high as 60%, that massive high-order multiple systems can be dynamically ejected, and that high-order multiples become common especially in the cluster. We also discuss binary populations of the ejected massive systems. When a large kinematic survey of massive field stars becomes available, for instance through Gaia, our results may be used to constrain the birth configuration of massive stars in star clusters. (Abridged)

preprint2016arXiv

On the origin of the Schechter-like mass function of young star clusters in disk galaxies

The mass function of freshly formed star clusters is empirically often described as a power law. However the cluster mass function of populations of young clusters over the scale of a galaxy has been found to be described by a Schechter-function. Here we address this apparent discrepancy. We assume that in an annulus of an isolated self- regulated radially-exponential axially-symmetric disk galaxy, the local mass function of very young (embedded) clusters is a power law with an upper mass limit which depends on the local star formation rate density. Radial integration of this mass function yields a galaxy-wide embedded cluster mass function. This integrated embedded cluster mass function has a Schechter-type form, which results from the addition of many low mass clusters forming at all galactocentric distances and rarer massive clusters only forming close to the center of the galaxy.

preprint2016arXiv

Simulating disk galaxies and interactions in Milgromian dynamics

Since its publication 1983, Milgromian dynamics (aka MOND) has been very successful in modeling the gravitational potential of galaxies from baryonic matter alone. However, the dynamical modeling has long been an unsolved issue. In particular, the setup of a stable galaxy for Milgromian N-body calculations has been a major challenge. Here, we show a way to set up disc galaxies in MOND for calculations in the PHANTOM OF RAMSES (PoR) code by Lüghausen (2015) and Teyssier (2002). The method is done by solving the QUMOND Poisson equations based on a baryonic and a phantom dark matter component. The resulting galaxy models are stable after a brief settling period for a large mass and size range. Simulations of single galaxies as well as colliding galaxies are shown.

preprint2016arXiv

Star formation triggered by galaxy interactions in modified gravity

Together with interstellar turbulence, gravitation is one key player in star formation. It acts both at galactic scales in the assembly of gas into dense clouds, and inside those structures for their collapse and the formation of pre-stellar cores. To understand to what extent the large scale dynamics govern the star formation activity of galaxies, we present hydrodynamical simulations in which we generalise the behaviour of gravity to make it differ from Newtonian dynamics in the low acceleration regime. We focus on the extreme cases of interacting galaxies, and compare the evolution of galaxy pairs in the dark matter paradigm to that in the Milgromian Dynamics (MOND) framework. Following up on the seminal work by Tiret & Combes, this paper documents the first simulations of galaxy encounters in MOND with a detailed Eulerian hydrodynamical treatment of baryonic physics, including star formation and stellar feedback. We show that similar morphologies of the interacting systems can be produced by both the dark matter and MOND formalisms, but require a much slower orbital velocity in the MOND case. Furthermore, we find that the star formation activity and history are significantly more extended in space and time in MOND interactions, in particular in the tidal debris. Such differences could be used as observational diagnostics and make interacting galaxies prime objects in the study of the nature of gravitation at galactic scales.

preprint2016arXiv

The observed spatial distribution of matter on scales ranging from 100kpc to 1Gpc is inconsistent with the standard dark-matter-based cosmological models

The spatial arrangement of galaxies (of satellites on a scale of 100kpc) as well as their three-dimensional distribution in galaxy groups such as the Local Group (on a scale of 1Mpc), the distribution of galaxies in the nearby volume of galaxies (on a scale of 8Mpc) and in the nearby Universe (on a scale of 1Gpc) is considered. There is further evidence that the CMB shows irregularities and for anisotropic cosmic expansion. The overall impression one obtains, given the best data we have, is matter to be arranged as not expected in the dark-matter based standard model of cosmology (SMoC). There appears to be too much structure, regularity and organisation. Dynamical friction on the dark matter halos is a strong direct test for the presence of dark matter particles, but this process does not appear to be operative in the real Universe. This evidence suggests strongly that dynamically relevant dark matter does not exist and therefore cosmology remains largely not understood theoretically. More-accepted awareness of this case would by itself constitute a major advance in research providing fabulous opportunities for bright minds, and the observational data strongly suggest that gravitation must be effectively Milgromian, corresponding to a generalized Poisson equation in the classical limit. Thus, physical cosmology offers a significant historically relevant opportunity for ground-breaking work, at least for those daring to do so.

preprint2015arXiv

Comment on "Evidence for dark matter in the inner Milky Way"

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

preprint2015arXiv

Dependency of dynamical ejections of O stars on the masses of very young star clusters

Massive stars can be efficiently ejected from their birth clusters through encounters with other massive stars. We study how the dynamical ejection fraction of O star systems varies with the masses of very young star clusters, Mecl, by means of direct N -body calculations. We include diverse initial conditions by varying the half-mass radius, initial mass-segregation, initial binary fraction and orbital parameters of the massive binaries. The results show robustly that the ejection fraction of O star systems exhibits a maximum at a cluster mass of $10^{3.5}$ Msun for all models, even though the number of the ejected systems increases with cluster mass. We show that lower mass clusters (Mecl ~ 400 Msun ) are the dominant sources for populating the Galactic field with O stars by dynamical ejections, considering the mass function of embedded clusters. About 15 per cent (up to 38 per cent, depending on the cluster models) of O stars of which a significant fraction are binaries, and which would have formed in a 10 Myr epoch of star formation in a distribution of embedded clusters, will be dynamically ejected to the field. Individual clusters may eject 100 per cent of their original O star content. A large fraction of such O stars have velocities up to only 10 km/s. Synthesising a young star cluster mass function it follows, given the stellar-dynamical results presented here, that the observed fractions of field and runaway O stars, and the binary fractions among them can be well understood theoretically if all O stars form in embedded clusters.

preprint2015arXiv

Dynamical Evolution of Outer-Halo Globular Clusters

Outer-halo globular clusters show large half-light radii and flat stellar mass functions, depleted in low-mass stars. Using N-body simulations of globular clusters on eccentric orbits within a Milky Way-like potential, we show how a cluster's half-mass radius and its mass function develop over time. The slope of the central mass function flattens proportionally to the amount of mass a cluster has lost, and the half-mass radius grows to a size proportional to the average strength of the tidal field. The main driver of these processes is mass segregation of dark remnants. We conclude that the extended, depleted clusters observed in the Milky Way must have had small half-mass radii in the past, and that they expanded due to the weak tidal field they spend most of their lifetime in. Moreover, their mass functions must have been steeper in the past but flattened significantly as a cause of mass segregation and tidal mass loss.

preprint2015arXiv

Galaxies as simple dynamical systems: observational data disfavor dark matter and stochastic star formation

(Abridged) Galaxies are observed to be simple systems but the standard model of cosmology (SMoC) implies a haphazard merger history driven by dynamical friction on dark matter halos. The SMoC is tested here with the results that the dual dwarf galaxy theorem is in contradiction to the data if the SMoC is true, and the action of dynamical friction is not evident in the galaxy population. A consistency test for this conclusion comes from the significantly anisotropic distributions of satellite galaxies. Independently, the long history of failures of the SMoC have the likelihood that it describes the observed Universe to less than 10^-4 per cent. The implication for fundamental physics is that exotic dark matter particles do not exist and that consequently effective gravitational physics on the scales of galaxies and beyond ought to be non-Newtonian/non-Einsteinian. The data imply that scale-invariant dynamics, as discovered by Milgrom, is an excellent description of galaxies in the weak-field regime. Observations of stellar populations in galaxies suggest the galaxy-wide IMF, the IGIMF, to vary with star formation rate and that stochastic descriptions of star formation are inconsistent with the data. A consequence of this understanding of galactic astrophysics is that most dwarf satellite galaxies are formed as tidal dwarf galaxies in galaxy-galaxy encounters, that they follow the mass-metallicity relation, that galactic mergers are rare, that galaxies immersed in external potentials are physically larger than isolated galaxies and that star-forming galaxies form a main sequence. Eight predictions are offered which will allow the concepts raised here to be tested. A very conservative, cold- and warm-dark-matter-free cosmological model may be emerging from these considerations.

preprint2015arXiv

Globular Cluster Streams as Galactic High-Precision Scales - The Poster Child Palomar 5

Using the example of the tidal stream of the Milky Way globular cluster Palomar 5 (Pal 5), we demonstrate how observational data on streams can be efficiently reduced in dimensionality and modeled in a Bayesian framework. Our approach combines detection of stream overdensities by a Difference-of-Gaussians process with fast streakline models, a continuous likelihood function built from these models, and inference with MCMC. By generating $\approx10^7$ model streams, we show that the geometry of the Pal 5 debris yields powerful constraints on the solar position and motion, the Milky Way and Pal 5 itself. All 10 model parameters were allowed to vary over large ranges without additional prior information. Using only SDSS data and a few radial velocities from the literature, we find that the distance of the Sun from the Galactic Center is $8.30\pm0.25$ kpc, and the transverse velocity is $253\pm16$ km/s. Both estimates are in excellent agreement with independent measurements of these quantities. Assuming a standard disk and bulge model, we determine the Galactic mass within Pal 5's apogalactic radius of 19 kpc to be $(2.1\pm0.4)\times10^{11}$ M$_\odot$. Moreover, we find the potential of the dark halo with a flattening of $q_z = 0.95^{+0.16}_{-0.12}$ to be essentially spherical within the radial range that is effectively probed by Pal 5. We also determine Pal 5's mass, distance and proper motion independently from other methods, which enables us to perform vital cross-checks. We conclude that with more observational data and by using additional prior information, the precision of this method can be significantly increased.

preprint2015arXiv

How can young massive clusters reach their present-day sizes?

The classic question that how young massive star clusters attain their shapes and sizes, as we find them today, remains to be a challenge. Both observational and computational studies of star-forming massive molecular gas clouds infer that massive cluster formation is primarily triggered along the small-scale ($\lesssim0.3$ pc) filamentary substructures within the clouds. The present study is intended to investigate the possible ways in which a filament-like-compact, massive star cluster (effective radius 0.1-0.3 pc) can expand $\gtrsim10$ times, still remaining massive enough ($\gtrsim10^4 M_\odot$), to become a young massive star cluster, as we observe today. To that end, model massive clusters (of initially $10^4 M_\odot-10^5 M_\odot$) are evolved using Sverre Aarseth's state-of-the-art N-body code NBODY7. All the computed clusters expand with time, whose sizes (effective radii) are compared with those observed for young massive clusters, of age $\lesssim100$ Myr, in the Milky Way and other nearby galaxies. It is found that beginning from the above compact sizes, a star cluster cannot expand by its own, i.e., due to two-body relaxation, stellar-evolutionary mass loss, dynamical heating by primordial binaries and stellar-mass black holes, up to the observed sizes of young massive clusters; they always remain much more compact compared to the observed ones. This calls for additional mechanisms that can boost the expansion of a massive cluster after its assembly. Using further N-body calculations, it is shown that a substantial residual gas expulsion, with $\approx30$% star formation efficiency, can indeed swell the newborn embedded cluster adequately. The limitations of the present calculations and their consequences are discussed.

preprint2015arXiv

M-dwarf binaries as tracers of star and brown dwarf formation

The separation distribution for M-dwarf binaries in the ASTRALUX survey is narrower and peaking at smaller separations than the distribution for solar-type binaries. This is often interpreted to mean that M-dwarfs constitute a continuous transition from brown dwarfs (BDs) to stars. Here a prediction for the M-dwarf separation distribution is presented, using a dynamical population synthesis (DPS) model in which "star-like" binaries with late-type primaries ($\lesssim1.5 M_{\rm sun}$) follow universal initial distribution functions and are dynamically processed in their birth embedded clusters. A separate "BD-like" population has both its own distribution functions for binaries and initial mass function (IMF), which overlaps in mass with the IMF for stars. Combining these two formation modes results in a peak on top of a wider separation distribution for late M-dwarfs consistent with the late ASTRALUX sample. The DPS separation distribution for early M-dwarfs shows no such peak and is in agreement with the M-dwarfs in Multiples (MinMS) data. We note that the latter survey is potentially in tension with the early ASTRALUX data. Concluding, the ASTRALUX and MinMS data are unable to unambiguously distinguish whether or not BDs are a continuous extension of the stellar IMF. Future observational efforts are needed to fully answer this interesting question. The DPS model predicts that binaries outside the sensitivity range of the ASTRALUX survey remain to be detected. For application to future data, we present a means to observationally measure the overlap of the putative BD-like branch and the stellar branch. We discuss the meaning of universal star formation and distribution functions.

preprint2015arXiv

Mass distributions of star clusters for different star formation histories in a galaxy cluster environment

Clusters of galaxies usually contain rich populations of globular clusters (GCs). We investigate how different star formation histories (SFHs) shape the final mass distribution of star clusters. We assume that every star cluster population forms during a formation epoch of length dt at a constant star-formation rate (SFR). The mass distribution of such a population is described by the embedded cluster mass function (ECMF), which is a pure power law extending to an upper limit M_max. Since the SFR determines M_max, the ECMF implicitly depends on the SFR. Starting with different SFHs, each SFH is divided into formation epochs of length dt at different SFRs. The requested mass function arises from the superposition of the star clusters of all formation epochs. An improved optimal sampling technique is introduced that allows generating number and mass distributions, both of which accurately agree with the ECMF. Moreover, for each SFH the distribution function of all involved SFRs, F(SFR), is computed. For monotonically decreasing SFHs, F(SFR) always follows a power law. With F(SFR), we develope the theory of the integrated galactic embedded cluster mass function (IGECMF). It describes the distribution function of birth stellar masses of star clusters that accumulated over a formation episode much longer than dt. The IGECMF indeed reproduces the mass distribution of star clusters created according to the superposition principle. Interestingly, all considered SFHs lead to a turn-down with increasing star cluster mass in their respective IGECMFs. In the past, a turn-down at the high-mass end has been observed for GC systems in different galaxy clusters and in the cluster initial mass function. This behavior can be explained naturally if the observed star cluster ensembles are superpositions of several individual star cluster populations that formed at different times at different SFRs.

preprint2015arXiv

On the persistence of two small-scale problems in ΛCDM

We investigate the degree to which the inclusion of baryonic physics can overcome two long-standing problems of the standard cosmological model on galaxy scales: (i) the problem of satellite planes around Local Group galaxies, and (ii) the "too big to fail" problem. By comparing dissipational and dissipationless simulations, we find no indication that the addition of baryonic physics results in more flattened satellite distributions around Milky-Way-like systems. Recent claims to the contrary are shown to derive in part from a non-standard metric for the degree of flattening, which ignores the satellites' radial positions. If the full 3D positions of the satellite galaxies are considered, none of the simulations we analyse reproduce the observed flattening nor the observed degree of kinematic coherence of the Milky Way satellite system. Our results are consistent with the expectation that baryonic physics should have little or no influence on the structure of satellite systems on scales of hundreds of kiloparsecs. Claims that the "too big to fail" problem can be resolved by the addition of baryonic physics are also shown to be problematic.

preprint2015arXiv

Possible smoking-gun evidence for initial mass segregation in re-virialized post-gas expulsion globular clusters

We perform a series of direct $N$-body calculations to investigate the effect of residual gas expulsion from the gas-embedded progenitors of present-day globular clusters (GCs) on the stellar mass function (MF). Our models start either tidally filling or underfilling, and either with or without primordial mass segregation. We cover 100 Myr of the evolution of modeled clusters and show that the expulsion of residual gas from initially mass-segregated clusters leads to a significantly shallower slope of the stellar MF in the low- ($m\leq 0.50 M_\odot$) and intermediate-mass ($\simeq 0.50-0.85 M_\odot$) regime. Therefore, the imprint of residual gas expulsion and primordial mass segregation might be visible in the present-day MF. We find that the strength of the external tidal field, as an essential parameter, influences the degree of flattening, such that a primordially mass-segregated tidally-filling cluster with $r_h/r_t\geq 0.1$ shows a strongly depleted MF in the intermediate stellar mass range. Therefore, the shape of the present-day stellar MF in this mass range probes the birth place of clusters in the Galactic environment. We furthermore find that this flattening agrees with the observed correlation between the concentration of a cluster and its MF slope, as found by de Marchi et al.. We show that if the expansion through the residual gas expulsion in primordial mass segregated clusters is the reason for this correlation then GCs most probably formed in strongly fluctuating local tidal fields in the early proto-Milky Way potential, supporting the recent conclusion by Marks \& Kroupa.

preprint2015arXiv

Probing the isotropy of cosmic acceleration traced by Type Ia supernovae

We present a method to test the isotropy of the magnitude-redshift relation of Type Ia Supernovae (SNe Ia) and single out the most discrepant direction (in terms of the signal-to-noise ratio) with respect to the all-sky data. Our technique accounts for possible directional variations of the corrections for SNe Ia and yields all-sky maps of the best-fit cosmological parameters with arbitrary angular resolution. To show its potential, we apply our method to the recent Union2.1 compilation, building maps with three different angular resolutions. We use a Monte Carlo method to estimate the statistical significance with which we could reject the null hypothesis that the magnitude-redshift relation is isotropic based on the properties of the observed most discrepant directions. We find that, based on pure signal-to-noise arguments, the null hypothesis cannot be rejected at any meaningful confidence level. However, if we also consider that the strongest deviations in the Union2.1 sample closely align with the dipole temperature anisotropy of the cosmic microwave background, we find that the null hypothesis should be rejected at the $95-99$ per cent confidence level, slightly depending on the angular resolution of the study. If this result is not due to a statistical fluke, it might either indicate that the SN data have not been cleaned from all possible systematics or even point towards new physics. We finally discuss future perspectives in the field for achieving larger and more uniform data sets that will vastly improve the quality of the results and optimally exploit our method.

preprint2015arXiv

The number of tidal dwarf satellite galaxies in dependence of bulge index

We show that a significant correlation (up to 5sigma) emerges between the bulge index, defined to be larger for larger bulge/disk ratio, in spiral galaxies with similar luminosities in the Galaxy Zoo 2 of SDSS and the number of tidal-dwarf galaxies in the catalogue by Kaviraj et al. (2012). In the standard cold or warm dark-matter cosmological models the number of satellite galaxies correlates with the circular velocity of the dark matter host halo. In generalized-gravity models without cold or warm dark matter such a correlation does not exist, because host galaxies cannot capture in-falling dwarf galaxies due to the absence of dark-matter-induced dynamical friction. However, in such models a correlation is expected to exist between the bulge mass and the number of satellite galaxies, because bulges and tidal-dwarf satellite galaxies form in encounters between host galaxies. This is not predicted by dark matter models in which bulge mass and the number of satellites are a priori uncorrelated because higher bulge/disk ratios do not imply higher dark/luminous ratios. Hence, our correlation reproduces the prediction of scenarios without dark matter, whereas an explanation is not found readily from the a priori predictions of the standard scenario with dark matter. Further research is needed to explore whether some application of the standard theory may explain this correlation.

preprint2014arXiv

A census of the expected properties of classical Milky Way dwarfs in Milgromian dynamics

Prompted by the recent successful predictions of the internal dynamics of Andromeda's satellite galaxies (McGaugh & Milgrom 2013a,b), we revisit the classical Milky Way dwarf spheroidal satellites Draco, Sculptor, Sextans, Carina, and Fornax in the framework of Milgromian dynamics (MOND). We use for the first time a Poisson solver with adaptive mesh refinement (AMR) in order to account simultaneously for the gravitational influence of the Milky Way and its satellites. This allows to rigorously model the important external field effect (EFE) of Milgromian dynamics, which can reduce the effective acceleration significantly. We make predictions on the dynamical mass-to-light ratio (Mdyn/L) expected to be measured by an observer who assumes Newtonian dynamics to be valid. We show that Milgromian dynamics predicts typical Mdyn/L ~ 10-50 Msun/Lsun. The results for the most luminous ones, Fornax and Sculptor, agree well with available velocity dispersion data. Moreover, the central power-law slopes of the dynamical masses agree exceedingly well with values inferred observationally from velocity dispersion measurements. The results for Sextans, Carina and Draco are low compared to usually quoted observational estimates, as already pointed out by Angus (2008). For Milgromian dynamics to survive further observational tests in these objects, one would thus need that either (a) previous observational findings based on velocity dispersion measurements have overestimated the dynamical mass due to, e.g., binaries and contaminant outliers, (b) the satellites are not in virial equilibrium due to the Milky Way tidal field, or (c) the specific theory used here does not describe the EFE correctly (e.g., the EFE could be practically negligible in some other theories), or a combination of (a)-(c).

preprint2014arXiv

A perfect starburst cluster made in one go: the NGC 3603 young cluster

Understanding how distinct, near-spherical gas-free clusters of very young, massive stars shape out of vast, complex clouds of molecular hydrogen is one of the biggest challenges in astrophysics. A popular thought dictates that a single gas cloud fragments into many new-born stars which, in turn, energize and rapidly expel the residual gas to form a gas-free cluster. This study demonstrates that the above classical paradigm remarkably reproduces the well-observed central, young cluster (HD 97950) of the Galactic NGC 3603 star-forming region, in particular, its shape, internal motion and the mass distribution of stars, naturally and consistently follow from a single model calculation. Remarkably, the same parameters (star formation efficiency, gas expulsion time scale and delay) reproduce HD 97950 as were found to reproduce the Orion Nebula Cluster, Pleiades and R136. The present results thereby provide intriguing evidences of formation of star clusters through single-starburst events followed by significant residual gas expulsion.

preprint2014arXiv

Chemo-dynamical evolution of tidal dwarf galaxies. II. The long-term evolution and influence of a tidal field

In a series of papers, we present detailed chemo-dynamical simulations of tidal dwarf galaxies (TDGs). After the first paper, where we focused on the very early evolution, we present in this work simulations on the long-term evolution of TDGs, ranging from their formation to an age of 3 Gyr. Dark-matter free TDGs may constitute a significant component of the dwarf galaxy (DG) population. But it remains to be demonstrated that TDGs can survive their formation phase given stellar feedback processes, the time-variable tidal field of the post-encounter host galaxy and its dark matter halo and ram-pressure wind from the gaseous halo of the host. For robust results the maximally damaging feedback by a fully populated invariant stellar IMF in each star cluster is assumed, such that fractions of massive stars contribute during phases of low star-formation rates. The model galaxies are studied in terms of their star-formation history, chemical enrichment and rotational curves. All models evolve into a self-regulated long-term equilibrium star-formation phase lasting for the full simulation time, whereby the TDGs become significantly more compact and sustain significantly higher SFRs through compressive tides than the isolated model. None of the models is disrupted despite the unphysical extreme feedback, and none of the rotation curves achieves the high values observed in real TDGs, despite non-virial gas accretion phases.

preprint2014arXiv

Co-orbiting satellite galaxy structures are still in conflict with the distribution of primordial dwarf galaxies

Both major galaxies in the Local Group host planar distributions of co-orbiting satellite galaxies, the Vast Polar Structure (VPOS) of the Milky Way and the Great Plane of Andromeda (GPoA). The $Λ$CDM cosmological model did not predict these features. However, according to three recent studies the properties of the GPoA and the flattening of the VPOS are common features among sub-halo based $Λ$CDM satellite systems, and the GPoA can be naturally explained by satellites being acquired along cold gas streams. We point out some methodological issues in these studies: either the selection of model satellites is different from that of the observed ones, or an incomplete set of observational constraints has been considered, or the observed satellite distribution is inconsistent with basic assumptions. Once these issues have been addressed, the conclusions are different: features like the VPOS and GPoA are very rare (each with probability $\lesssim 10^{-3}$, and combined probability $< 10^{-5}$) if satellites are selected from a $Λ$CDM simulation combined with semi-analytic modelling, and accretion along cold streams is no natural explanation of the GPoA. The origin of planar dwarf galaxy structures remains unexplained in the standard paradigm of galaxy formation.

preprint2014arXiv

Direct $N$-body simulations of globular clusters - II. Palomar 4

We use direct $N$-body calculations to study the evolution of the unusually extended outer halo globular cluster Palomar 4 (Pal~4) over its entire lifetime in order to reproduce its observed mass, half-light radius, velocity dispersion and mass function slope at different radii. We find that models evolving on circular orbits, and starting from a non-mass segregated, canonical initial mass function (IMF) can reproduce neither Pal 4's overall mass function slope nor the observed amount of mass segregation. Including either primordial mass segregation or initially flattened IMFs does not reproduce the observed amount of mass segregation and mass function flattening simultaneously. Unresolved binaries cannot reconcile this discrepancy either. We find that only models with both a flattened IMF and primordial segregation are able to fit the observations. The initial (i.e. after gas expulsion) mass and half-mass radius of Pal~4 in this case are about 57000 M${\odot}$ and 10 pc, respectively. This configuration is more extended than most globular clusters we observe, showing that the conditions under which Pal~4 formed must have been significantly different from that of the majority of globular clusters. We discuss possible scenarios for such an unusual configuration of Pal~4 in its early years.

preprint2014arXiv

Erosion of Globular Cluster Systems: The Influence of Radial Anisotropy, Central Black Holes and Dynamical Friction

We present the adaptable MUESLI code for investigating dynamics and erosion processes of globular clusters (GCs) in galaxies. MUESLI follows the orbits of individual clusters and applies internal and external dissolution processes to them. Orbit integration is based on the self-consistent field method in combination with a time-transformed leapfrog scheme, allowing us to handle velocity-dependent forces like triaxial dynamical friction. In a first application, the erosion of globular cluster systems (GCSs) in elliptical galaxies is investigated. Observations show that massive ellipticals have rich, radially extended GCSs, while some compact dwarf ellipticals contain no GCs at all. For several representative examples, spanning the full mass scale of observed elliptical galaxies, we quantify the influence of radial anisotropy, galactic density profiles, SMBHs, and dynamical friction on the GC erosion rate. We find that GC number density profiles are centrally flattened in less than a Hubble time, naturally explaining observed cored GC distributions. The erosion rate depends primarily on a galaxy's mass, half-mass radius and radial anisotropy. The fraction of eroded GCs is nearly 100% in compact, M 32 like galaxies and lowest in extended and massive galaxies. Finally, we uncover the existence of a violent tidal disruption dominated phase which is important for the rapid build-up of halo stars.

preprint2014arXiv

Galactic rotation curves, the baryon-to-dark-halo-mass relation and space-time scale invariance

Low-acceleration space-time scale invariant dynamics (SID, Milgrom 2009a) predicts two fundamental correlations known from observational galactic dynamics: the baryonic Tully-Fisher relation (BTFR) and a correlation between the observed mass discrepancy and acceleration (MDA) in the low acceleration regime for disc galaxies. SID corresponds to the deep MOdified Newtonian Dynamics (MOND) limit. The MDA data emerging in cold/warm dark matter (C/WDM) cosmological simulations disagree significantly with the tight MDA correlation of the observed galaxies. Therefore, the most modern simulated disc galaxies, which are delicately selected to have a quiet merging history in a standard dark-matter-cosmological model, still do not represent the correct rotation curves. Also, the observed tight correlation contradicts the postulated stochastic formation of galaxies in low-mass DM halos. Moreover, we find that SID predicts a baryonic to apparent virial halo (dark matter) mass relation which agrees well with the correlation deduced observationally assuming Newtonian dynamics to be valid, while the baryonic to halo mass relation predicted from CDM models does not. The distribution of the observed ratios of dark-matter halo masses to baryonic masses may be empirical evidence for the external field effect, which is predicted in SID as a consequence of the forces acting between two galaxies depending on the position and mass of a third galaxy. Applying the external field effect, we predict the masses of galaxies in the proximity of the dwarf galaxies in the Miller et al. sample. Classical non-relativistic gravitational dynamics is thus best described as being Milgromian, rather than Newtonian.

preprint2014arXiv

Phantom of RAMSES (POR): A new Milgromian dynamics N-body code

Since its first formulation in 1983, Milgromian dynamics (MOND) has been very successful in predicting the gravitational potential of galaxies from the distribution of baryons alone, including general scaling relations and detailed rotation curves of large statistical samples of individual galaxies covering a large range of masses and sizes. Most predictions however rely on static models, and only a handful of N-body codes have been developed over the years to investigate the consequences of the Milgromian framework for the dynamics of complex evolving dynamical systems. In this work, we present a new Milgromian N-body code, which is a customized version of the RAMSES code (Teyssier 2002) and thus comes with all its features: it includes particles and gas dynamics, and importantly allows for high spatial resolution of complex systems due to the adaptive mesh refinement (AMR) technique. It further allows the direct comparison between Milgromian simulations and standard Newtonian simulations with dark matter particles. We provide basic tests of this customized code and demonstrate its performance by presenting N-body computations of dark-matter-free spherical equilibrium models as well as dark-matter-free disk galaxies in Milgromian dynamics.

preprint2014arXiv

Reproducing properties of MW dSphs as descendants of DM-free TDGs

The Milky Way (MW) dwarf spheroidal (dSph) satellites are known to be the most dark-matter (DM) dominated galaxies with estimates of dark to baryonic matter reaching even above one hundred. It comes from the assumption that dwarfs are dynamically supported by their observed velocity dispersions. However their spatial distributions around the MW is not at random and this could challenge their origin, previously assumed to be residues of primordial galaxies accreted by the MW potential. Here we show that alternatively, dSphs could be the residue of tidal dwarf galaxies (TDGs), which would have interacted with the Galactic hot gaseous halo and disk. TDGs are gas-rich and have been formed in a tidal tail produced during an ancient merger event at the M31 location, and expelled towards the MW. Our simulations show that low-mass TDGs are fragile to an interaction with the MW disk and halo hot gas. During the interaction, their stellar content is progressively driven out of equilibrium and strongly expands, leading to low surface brightness feature and mimicking high dynamical M/L ratios. Our modeling can reproduce the properties, including the kinematics, of classical MW dwarfs within the mass range of the Magellanic Clouds to Draco. An ancient gas-rich merger at the M31 location could then challenge the currently assumed high content of dark matter in dwarf galaxies. We propose a simple observational test with the coming GAIA mission, to follow their expected stellar expansion, which should not be observed within the current theoretical framework.

preprint2014arXiv

Sampling methods for stellar masses and the mmax-Mecl relation in the starburst dwarf galaxy NGC 4214

It has been claimed in the recent literature that a non-trivial relation between the mass of the most-massive star, mmax, in a star cluster and its embedded star cluster mass (the mmax-Mecl relation) is falsified by observations of the most-massive stars and the Halpha luminosity of young star clusters in the starburst dwarf galaxy NGC 4214. Here it is shown by comparing the NGC 4214 results with observations from the Milky Way that NGC 4214 agrees very well with the predictions of the the mmax-Mecl relation and the integrated galactic stellar initial mass function (IGIMF) theory and that this difference in conclusions is based on a high degree of degeneracy between expectations from random sampling and those from the mmax-Mecl relation, but are also due to interpreting mmax as a truncation mass in a randomly sampled IMF. Additional analysis of galaxies with lower SFRs than those currently presented in the literature will be required to break this degeneracy.

preprint2014arXiv

Simulations on the survivability of Tidal Dwarf Galaxies

We present detailed numerical simulations of the evolution of Tidal Dwarf Galaxies (TDGs) after they kinematically decouple from the rest of the tidal arm to investigate their survivability. Both the short-term (500 Myr) response of TDGs to the stellar feedback of different underlying stellar populations as well as the long-term evolution that is dominated by a time dependent tidal field is examined. All simulated TDGs survive until the end of the simulation time of up to 3 Gyr, despite their lack of a stabilising dark matter component.

preprint2014arXiv

The formation of NGC 3603 young starburst cluster: "prompt" hierarchical assembly or monolithic starburst?

The formation of very young massive clusters or "starburst" clusters is currently one of the most widely debated topic in astronomy. The classical notion dictates that a star cluster is formed in-situ in a dense molecular gas clump followed by a substantial residual gas expulsion. On the other hand, based on the observed morphologies of many young stellar associations, a hierarchical formation scenario is alternatively suggested. A very young (age $\approx$ 1 Myr), massive ($>10^4M_\odot$) star cluster like the Galactic NGC 3603 young cluster (HD 97950) is an appropriate testbed for distinguishing between such "monolithic" and "hierarchical" formation scenarios. A recent study by Banerjee and Kroupa (2014) demonstrates that the monolithic scenario remarkably reproduces the HD 97950 cluster. In the present work, we explore the possibility of the formation of the above cluster via hierarchical assembly of subclusters. These subclusters are initially distributed over a wide range of spatial volumes and have various modes of sub-clustering. Unlike the monolithic model that reproduces HD 97950 very well, the same is found to be prohibitive with hierarchical assembly alone. Only those systems which assemble promptly into a single cluster (in $\lesssim1$ Myr) from a close separation (all within $\lesssim2$ pc) could match the observed density profile of HD 97950 after $\approx70$% gas dispersal. These results therefore suggest that NGC 3603 young cluster has formed essentially monolithically, i.e., either in-situ or via a prompt assembly, followed by a substantial residual gas expulsion. Both scenarios are consistent with the inferred young age and the small age spread of this cluster. Future observations of molecular cloud filaments with ALMA and proper motion measurements of young clusters with Gaia will provide more direct tests of such birth environments. [Abridged]

preprint2014arXiv

The Puzzling Negative Orbit-Period Derivative of the Low-Mass X-Ray Binary 4U 1820-30 in NGC 6624

4U 1820-30 is a low-mass X-ray binary near the center of the globular cluster NGC 6624 consisting of, at least, one neutron star and one helium white dwarf. Analyzing 16 years of data from the Rossi X-ray Timing Explorer (RXTE) allows us to measure its orbital period and its time derivative with unprecedented accuracy to be P = 685.01197 +- 0.00003 s and dP/dt /P = -5.3 +- 0.3x10^-8 yr^-1. Hence, we confirm that the period derivative is significantly negative at the >17 sigma level, contrary to theoretical expectations for an isolated X-ray binary. We discuss possible scenarios that could explain this discrepancy, and conclude that the center of NGC 6624 most likely contains large amounts of non-luminous matter such as dark remnants. We also discuss the possibility of an IMBH inside NGC 6624, or that a dark remnant close to 4U 1820-30 causes the observed shift.

preprint2014arXiv

The state of globular clusters at birth II: primordial binaries

(abridged) In this paper, we constrain the properties of primordial binary populations in Galactic globular clusters using the MOCCA Monte Carlo code for cluster evolution. Our results are compared to the observations of Milone et al. (2012) using the photometric binary populations as proxies for the true underlying distributions, in order to test the hypothesis that the data are consistent with an universal initial binary fraction near unity and the binary orbital parameter distributions of Kroupa (1995). With the exception of a few possible outliers, we find that the data are to first-order consistent with the universality hypothesis. Specifically, the present-day binary fractions inside the half-mass radius r$_{\rm h}$ can be reproduced assuming either high initial binary fractions near unity with a dominant soft binary component as in the Kroupa distribution combined with high initial densities (10$^4$-10$^6$ M$_{\odot}$ pc$^{-3}$), or low initial binary fractions ($\sim$ 5-10%) with a dominant hard binary component combined with moderate initial densities near their present-day values (10$^2$-10$^3$ M$_{\odot}$ pc$^{-3}$). This apparent degeneracy can be broken using the binary fractions outside r$_{\rm h}$- only high initial binary fractions with a significant soft component combined with high initial densities can contribute to reproducing the observed anti-correlation between the binary fractions outside r$_{\rm h}$ and the total cluster mass. We further illustrate using the simulated present-day binary orbital parameter distributions and the technique introduced in Leigh et al. (2012) that the relative fractions of hard and soft binaries can be used to further constrain the initial cluster density and mass-density relation. Our results favour an initial mass-density relation of the form r$_{\rm h} \propto$ M$_{\rm clus}^α$ with $α<$ 1/3.

preprint2014arXiv

The state of globular clusters at birth: emergence from the gas-embedded phase

In this paper, we discuss the origin of the observed correlation between cluster concentration c and present-day mass function (PDMF) slope α reported by De Marchi, Paresce & Pulone. This relation can either be reproduced from universal initial conditions combined with some dynamical mechanism(s) that alter(s) the cluster structure and mass function over time, or it must arise early on in the cluster lifetime, such as during the gas-embedded phase of cluster formation. Using a combination of Monte Carlo and N-body models for globular cluster evolution performed with the MOCCA and NBODY6 codes, respectively, we explore a number of dynamical mechanisms that could affect the observed relation. For the range of initial conditions considered here, our results are consistent with an universal initial binary fraction \approx 10% (which does not, however, preclude 100%) and an universal initial stellar mass function resembling the standard Kroupa distribution. Most of the dispersion observed in the c-α relation can be attributed to two-body relaxation and Galactic tides. However, dynamical processes alone could not have reproduced the dispersion in concentration, and we require at least some correlation between the initial concentration and the total cluster mass. We argue that the origin of this trend could be connected to the gas-embedded phase of cluster evolution.

preprint2014arXiv

The Vast Polar Structure of the Milky Way Attains New Members

The satellite galaxies of the Milky Way (MW) align with and preferentially orbit in a vast polar structure (VPOS), which also contains globular clusters, stellar and gaseous streams. Similar alignments have been discovered around several other host galaxies. We test whether recently discovered objects in the MW halo, the satellite galaxy/globular cluster transition object named PSO J174.0675-10.8774 or Crater and three stellar streams, are part of the VPOS. Crater is situated close to the VPOS. Incorporating the new object in the VPOS-plane fit slightly improves the alignment of the plane with other features such as the Magellanic stream and the average orbital plane of the satellites co-orbiting in the VPOS. We predict Crater's proper motion by assuming that it, too, orbits in the VPOS. One of the three streams aligns well with the VPOS. Surprisingly, it appears to lie in the exact same orbital plane as the Palomar 5 stream and shares its distance, suggesting a direct connection between the two. The stream also crosses close to the Fornax dwarf galaxy and is oriented approximately along the galaxy's direction of motion. The two other streams cannot align closely with the VPOS because they were discovered in the direction of M31/M33, which is outside of the satellite structure. The VPOS thus attains two new members. This further emphasizes that the highly anisotropic and correlated distribution of satellite objects requires an explanation beyond the suggestion that the MW satellite system is an extreme statistical outlier of a $Λ$CDM sub-halo system.

preprint2013arXiv

Did the infant R136 and NGC 3603 clusters undergo residual gas expulsion?

Based on kinematic data observed for very young, massive clusters that appear to be in dynamical equilibrium, it has recently been argued that such young systems set examples where the early residual gas-expulsion did not happen or had no dynamical effect. The intriguing scenario of a star cluster forming through a single starburst has thereby been challenged. Choosing the case of the R136 cluster of the Large Magellanic Cloud, the most cited one in this context, we perform direct N-body computations that mimic the early evolution of this cluster including the gas-removal phase (on a thermal timescale). Our calculations show that under plausible initial conditions as consistent from observational data, a large fraction (> 60%) of a gas-expelled, expanding R136-like cluster is bound to regain dynamical equilibrium by its current age. Therefore, the recent measurements of velocity dispersion in the inner regions of R136, that indicate that the cluster is in dynamical equilibrium, are consistent with an earlier substantial gas expulsion of R136 followed by a rapid re-virialization (in approx. 1 Myr). Additionally, we find that the less massive Galactic NGC 3603 Young Cluster (NYC), with a substantially longer re-virialization time, is likely to be found deviated from dynamical equilibrium at its present age (approx. 1 Myr). The recently obtained stellar proper motions in the central part of the NYC indeed suggest this and are consistent with the computed models. This work significantly extends previous models of the Orion Nebula Cluster which already demonstrated that the re-virialization time of young post-gas-expulsion clusters decreases with increasing pre-expulsion density.

preprint2013arXiv

Dwarf Galaxy Planes: the discovery of symmetric structures in the Local Group

Both major galaxies in the Local Group (LG) are surrounded by thin planes of mostly co-orbiting satellite galaxies, the vast polar structure (VPOS) around the Milky Way (MW) and the Great Plane of Andromeda (GPoA) around M31. We summarize the current knowledge concerning these structures and compare their relative orientations and properties in a common coordinate system. The existence of coherent satellite structures motivates an investigation of the distribution of the more distant non-satellite galaxies in the LG. This results in the discovery of two planes (diameters of 1-2 Mpc) which contain almost all nearby non-satellite galaxies. The two LG planes are surprisingly symmetric, inclined by only 20 deg relative to the galactic disc of M31, similarly thin (heights of ~60 kpc) and have near-to-identical offsets from the MW and M31. They are inclined relative to each other by 35 deg. Comparing the plane orientations with each other and with additional features indicates an intimate connection between the VPOS and the GPoA. They are both polar to the MW, have similar orbital directions and are inclined by ~45 deg relative to each other. The Magellanic Stream approximately aligns with the VPOS and the GPoA, but also shares its projected position and line-of-sight velocity trend with the dominating structure of non-satellite dwarf galaxies. The proper motion of M31 indicates a prograde orbit of the MW-M31 system, the VPOS and the GPoA. The alignment with other features such as the Supergalactic Plane and the over-density in hypervelocity stars are discussed as well. We end with a summary of proposed scenarios trying to explain the LG galaxy structures as either originating from cosmological structures or from tidal debris of a past galaxy encounter. There currently exists no full detailed model which satisfactorily explains the existence of the thin symmetric LG planes. (abridged)

preprint2013arXiv

Gas surface density, star formation rate surface density, and the maximum mass of young star clusters in a disk galaxy. II. The grand-design galaxy M51

We analyze the relationship between maximum cluster mass, and surface densities of total gas (Sigma_gas), molecular gas (Sigma_H_2), neutral gas (Sigma_HI) and star formation rate (Sigma_SFR) in the grand design galaxy M51, using published gas data and a catalog of masses, ages, and reddenings of more than 1800 star clusters in its disk, of which 223 are above the cluster mass distribution function completeness limit. We find for clusters older than 25 Myr that M_3rd, the median of the 5 most massive clusters, is proportional to Sigma_HI^0.4. There is no correlation with Sigma_gas, Sigma_H2, or Sigma_SFR. For clusters younger than 10 Myr, M_3rd is proportional to Sigma_HI^0.6, M_3rd is proportional to Sigma_gas^0.5; there is no correlation with either Sigma_H_2 or Sigma_SFR. The results could hardly be more different than those found for clusters younger than 25 Myr in M33. For the flocculent galaxy M33, there is no correlation between maximum cluster mass and neutral gas, but M_3rd is proportional to Sigma_gas^3.8; M_3rd is proportional to Sigma_H_2^1.2; M_3rd proportional to Sigma_SFR^0.9. For the older sample in M51, the lack of tight correlations is probably due to the combination of the strong azimuthal variations in the surface densities of gas and star formation rate, and the cluster ages. These two facts mean that neither the azimuthal average of the surface densities at a given radius, nor the surface densities at the present-day location of a stellar cluster represent the true surface densities at the place and time of cluster formation. In the case of the younger sample, even if the clusters have not yet traveled too far from their birthsites, the poor resolution of the radio data compared to the physical sizes of the clusters results in measured Sigmas that are likely quite diluted compared to the actual densities relevant for the formation of the clusters.

preprint2013arXiv

Local Group timing in Milgromian dynamics. A past Milky Way-Andromeda encounter at z>0.8

The Local Group timing has been one of the first historical probes of the missing mass problem. Whilst modern cosmological probes indicate that pure baryonic dynamics is not sufficient on the largest scales, nearby galaxies and small galaxy groups persistently obey Milgrom's MOND law, which implies that dynamics at small scales is possibly entirely predicted by the baryons. Here, we investigate the Local Group timing in this context of Milgromian dynamics. Making use of the latest measured proper motions and radial velocities for Andromeda and the Magellanic clouds, we integrate backwards their orbits by making use of the Milgromian two-body equation of motion. We find that, with the currently measured proper motions and radial velocity of M31, MOND would imply that the Milky Way and M31 first moved apart via Hubble expansion after birth, but then necessarily got attracted again by the Milgromian gravitational attraction, and had a past fly-by encounter before coming to their present positions. This encounter would most probably have happened 7 to 11 Gyr ago (0.8<z<3). The absence of a dark matter halo and its associated dynamical friction is necessary for such a close encounter not to trigger a merger. Observational arguments which could exclude or favour such a past encounter would thus be very important in view of falsifying or vindicating Milgromian dynamics on the scale of the Local Group. Interestingly, the closest approach of the encounter is small enough (<55 kpc) to have severe consequences on the disk dynamics, including perhaps thick disk formation, and on the satellite systems of both galaxies. The ages of the satellite galaxies and of the young halo globular clusters, all of which form the vast polar structure around the Milky Way, are consistent with these objects having been born in this encounter.

preprint2013arXiv

R144 : a very massive binary likely ejected from R136 through a binary-binary encounter

R144 is a recently confirmed very massive, spectroscopic binary which appears isolated from the core of the massive young star cluster R136. The dynamical ejection hypothesis as an origin for its location is claimed improbable by Sana et al. due to its binary nature and high mass. We demonstrate here by means of direct N-body calculations that a very massive binary system can be readily dynamically ejected from a R136-like cluster, through a close encounter with a very massive system. One out of four N-body cluster models produces a dynamically ejected very massive binary system with a mass comparable to R144. The system has a system mass of $\approx$ 355 Msun and is located at 36.8 pc from the centre of its parent cluster, moving away from the cluster with a velocity of 57 km/s at 2 Myr as a result of a binary-binary interaction. This implies that R144 could have been ejected from R136 through a strong encounter with an other massive binary or single star. In addition, we discuss all massive binaries and single stars which are ejected dynamically from their parent cluster in the N-body models.

preprint2013arXiv

Recent advances on IMF research

Here I discuss recent work on brown dwarfs, massive stars and the IMF in general. The stellar IMF can be well described by an invariant two-part power law in present-day star-formation events within the Local Group of galaxies. It is nearly identical in shape to the pre-stellar core mass function. The majority of brown dwarfs follow a separate IMF. Evidence from globular clusters and ultra-compact dwarf galaxies has emerged that IMFs may have been top heavy depending on the star-formation rate density. The IGIMF then ranges from bottom heavy at low galaxy-wide star formation rates to being top-heavy in galaxy-scale star bursts.

preprint2013arXiv

The dynamical phase transitions of stellar systems and the corresponding kinetimacs

External fields in Migromian dynamics (MD or MOND, Milgrom 1983) break the Strong Equivalence Principle (SEP) and change the dynamics of self-bound stellar systems moving in space-varying background gravitational fields. We study two kinds of re-virialisation of the stellar systems: the violent phase transition and the adiabatic phase transition for systems moving on radial orbits, where the external field evolves from strong to weak and whose corresponding dynamics changes from Newtonian to Milgromian. We find that the time scale for the phase tranformation from Newtonian to Milgromian gravity lies only within one to a few crossing times for low density globular clusters with masses ranging from $10^4\msun$ to $10^6\msun$. Thus a globular cluster can appear frozen in the Newtonian regime despite being in the Milgromian regime for not longer than a few crossing times. We also study the kinematics and anisotropy profiles of the systems. The velocity dispersions of the systems are larger after the phase transitions, especially for the outer regions of the stellar systems. Moreover, the isotropic systems become radially anisotropic, especially for the outer parts, after the process caused by the dynamical phase transition. Deeper Milgromian systems have more radially anisotropic velocity dispersion functions. We also find that the final profiles of density, velocity dispersion and anisotropy do not depend on the details of the phase transition. I.e., the mass distribution and kinematics of the end-states of the globular clusters do not depend on the rapidity of the transition from Newtonian to Milgromian gravity. Thus, the transition from the Newtonian to the Milgromian regime naturally induces a significant radially anisotropic velocity distribution in a globular cluster.

preprint2013arXiv

The failures of the standard model of cosmology require a new paradigm

Cosmological models that invoke warm or cold dark matter can not explain observed regularities in the properties of dwarf galaxies, their highly anisotropic spatial distributions, nor the correlation between observed mass discrepancies and acceleration. These problems with the standard model of cosmology have deep implications, in particular in combination with the observation that the data are excellently described by Modified Newtonian Dynamics (MOND). MOND is a classical dynamics theory which explains the mass discrepancies in galactic systems, and in the universe at large, without invoking dark entities. MOND introduces a new universal constant of nature with the dimensions of acceleration, a0, such that the pre-MONDian dynamics is valid for accelerations a >> a0, and the deep MONDian regime is obtained for a << a0, where space-time scale invariance is invoked. Remaining challenges for MOND are (i) explaining fully the observed mass discrepancies in galaxy clusters, and (ii) the development of a relativistic theory of MOND that will satisfactorily account for cosmology. The universal constant a0 turns out to have an intriguing connection with cosmology: \bar a0 == 2 pi a0 \approx c H_0 \approx c^2(Lambda/3)^{1/2}. This may point to a deep connection between cosmology and internal dynamics of local systems.

preprint2013arXiv

The galactocentric radius dependent upper mass limit of young star clusters: stochastic star formation ruled out

It is widely accepted that the distribution function of the masses of young star clusters is universal and can be purely interpreted as a probability density distribution function with a constant upper mass limit. As a result of this picture the masses of the most-massive objects are exclusively determined by the size of the sample. Here we show, with very high confidence, that the masses of the most-massive young star clusters in M33 decrease with increasing galactocentric radius in contradiction to the expectations from a model of a randomly sampled constant cluster mass function with a constant upper mass limit. Pure stochastic star formation is thereby ruled out. We use this example to elucidate how naive analysis of data can lead to unphysical conclusions.

preprint2013arXiv

The galaxy-wide IMF of dwarf late-type to massive early-type galaxies

Observational studies are showing that the galaxy-wide stellar initial mass function are top-heavy in galaxies with high star-formation rates (SFRs). Calculating the integrated galactic stellar initial mass function (IGIMF) as a function of the SFR of a galaxy, it follows that galaxies which have or which formed with SFRs > 10 Msol yr^-1 would have a top-heavy IGIMF in excellent consistency with the observations. Consequently and in agreement with observations, elliptical galaxies would have higher M/L ratios as a result of the overabundance of stellar remnants compared to a stellar population that formed with an invariant canonical stellar initial mass function (IMF). For the Milky Way, the IGIMF yields very good agreement with the disk- and the bulge-IMF determinations. Our conclusions are that purely stochastic descriptions of star formation on the scales of a pc and above are falsified. Instead, star formation follows the laws, stated here as axioms, which define the IGIMF theory. We also find evidence that the power-law index beta of the embedded cluster mass function decreases with increasing SFR. We propose further tests of the IGIMF theory through counting massive stars in dwarf galaxies.

preprint2013arXiv

The mmax-Mecl relation, the IMF and IGIMF: probabilistically sampled functions?

We introduce a new method to measure the dispersion of mmax values of star clusters and show that the observed sample of mmax is inconsistent with random sampling from an universal stellar initial mass function (IMF) at a 99.9% confidence level. The scatter seen in the mmax-Mecl data can be mainly (76%) understood as being the result of observational uncertainties only. The scatter of mmax values at a given Mecl are consistent with mostly measurement uncertainties such that the true (physical) scatter may be very small. Additionally, new data on the local star-formation regions Taurus-Auriga and L1641 in Orion make stochastically formed stellar populations rather unlikely. The data are however consistent with the local IGIMF (integrated galactic stellar initial mass function) theory according to which a stellar population is a sum of individual star-forming events each of which is described by well defined physical laws. Randomly sampled IMFs and henceforth scale-free star formation seems to be in contradiction to observed reality.

preprint2013arXiv

The rotationally stabilized VPOS and predicted proper motions of the Milky Way satellite galaxies

The satellite galaxies of the Milky Way (MW) define a vast polar structure (VPOS), a thin plane perpendicular to the MW disc. Proper motion (PM) measurements are now available for all of the 11 brightest, `classical' satellites and allow an updated analysis of the alignment of their orbital poles with this spatial structure. The coherent orbital alignment of 7 to 9 out of 11 satellites demonstrates that the VPOS is a rotationally stabilized structure and not only a pressure-supported, flattened ellipsoid. This allows us to empirically and model independently predict the PMs of almost all satellite galaxies by assuming that the MW satellite galaxies orbit within the VPOS. As a test of our method, the predictions are best met by satellites whose PMs are already well constrained, as expected because more uncertain measurements tend to deviate more from the true values. Improved and new PM measurements will further test these predictions. A strong alignment of the satellite galaxy orbital poles is not expected in dark matter based simulations of galaxy formation. Coherent orbital directions of satellite galaxies are, however, a natural consequence of tidal dwarf galaxies formed together in the debris of a galaxy collision. The orbital poles of the MW satellite galaxies therefore lend further support to tidal scenarios for the origin of the VPOS and are a very significant challenge for the standard LCDM model of cosmology. We also note that the dependence of the MW satellite speeds on Galactocentric distance appear to map an effective potential with a constant velocity of approximately 240 km/s to about 250 kpc. The individual satellite velocities are only mildly radial.

preprint2013arXiv

The specific frequency and the globular cluster formation efficiency in Milgromian dynamics

Previous studies of globular cluster (GC) systems show that there appears to be a universal specific GC formation efficiency $η$ which relates the total mass of GCs to the virial mass of host dark matter halos, $M_{vir}$ (Georgiev et al 2010, Spitler & Forbes2009). In this paper, the specific frequency, $S_N$, and specific GC formation efficiency, $η$, are derived as functions of $M_{vir}$ in Milgromian dynamics, i.e., in modified Newtonian dynamics (MOND). In Milgromian dynamics, for the galaxies with GCs, the mass of the GC system, $M_{GC}$, is a two-component function of $M_{vir}$ instead of a simple linear relation. An observer in a Milgromian universe, who interprets this universe as being Newtonian/Einsteinian, will incorrectly infer a universal constant fraction between the mass of the GC system and a (false) dark matter halo of the baryonic galaxy. In contrast to a universal constant of $η$, in a Milgromian universe, for galaxies with $M_{vir} <= 10^{12}\msun$, $η$ decreases with the increase of $M_{vir}$, while for massive galaxies with $M_{vir}>10^{12}\msun$, $η$ increases with the increase of $M_{vir}$.

preprint2013arXiv

The vast thin plane of M31 co-rotating dwarfs: an additional fossil signature of the M31 merger and of its considerable impact in the whole Local Group

The recent discovery by Ibata et al. (2013) of a vast thin disk of satellites (VTDS) around M31 offers a new challenge for the understanding of the Local Group properties. This comes in addition to the unexpected proximity of the Magellanic Clouds (MCs) to the Milky Way (MW), and to another vast polar structure (VPOS), which is almost perpendicular to our Galaxy disk. We find that the VTDS plane is coinciding with several stellar, tidally-induced streams in the outskirts of M31, and, that its velocity distribution is consistent with that of the Giant Stream (GS). This is suggestive of a common physical mechanism, likely linked to merger tidal interactions, knowing that a similar argument may apply to the VPOS at the MW location. Furthermore, the VTDS is pointing towards the MW, being almost perpendicular to the MW disk, as the VPOS is. We compare these properties to the modelling of M31 as an ancient, gas-rich major merger, which has been successfully used to predict the M31 substructures and the GS origin. We find that without fine tuning, the induced tidal tails are lying in the VTDS plane, providing a single and common origin for many stellar streams and for the vast stellar structures surrounding both the MW and M31. The model also reproduces quite accurately positions and velocities of the VTDS dSphs. Our conjecture leads to a novel interpretation of the Local Group past history, as a gigantic tidal tail due to the M31 ancient merger is expected to send material towards the MW, including the MCs. Such a link between M31 and the MW is expected to be quite exceptional, though it may be in qualitative agreement with the reported rareness of MW-MCs systems in nearby galaxies.

preprint2012arXiv

A natural formation scenario for misaligned and short-period eccentric extrasolar planets

Recent discoveries of strongly misaligned transiting exoplanets pose a challenge to the established planet formation theory which assumes planetary systems to form and evolve in isolation. However, the fact that the majority of stars actually do form in star clusters raises the question how isolated forming planetary systems really are. Besides radiative and tidal forces the presence of dense gas aggregates in star-forming regions are potential sources for perturbations to protoplanetary discs or systems. Here we show that subsequent capture of gas from large extended accretion envelopes onto a passing star with a typical circumstellar disc can tilt the disc plane to retrograde orientation, naturally explaining the formation of strongly inclined planetary systems. Furthermore, the inner disc regions may become denser, and thus more prone to speedy coagulation and planet formation. Pre-existing planetary systems are compressed by gas inflows leading to a natural occurrence of close-in misaligned hot Jupiters and short-period eccentric planets. The likelihood of such events mainly depends on the gas content of the cluster and is thus expected to be highest in the youngest star clusters.

preprint2012arXiv

Catch me if you can: is there a runaway-mass black hole in the Orion Nebula Cluster?

We investigate the dynamical evolution of the Orion Nebula Cluster (ONC) by means of direct N-body integrations. A large fraction of residual gas was probably expelled when the ONC formed, so we assume that the ONC was much more compact when it formed compared to its current size, in agreement with the embedded cluster radius-mass relation from Marks & Kroupa (2012). Hence, we assume that few-body relaxation played an important role during the initial phase of evolution of the ONC. In particular, three body interactions among OB stars likely led to their ejection from the cluster and, at the same time, to the formation of a massive object via runaway physical stellar collisions. The resulting depletion of the high mass end of the stellar mass function in the cluster is one of the important points where our models fit the observational data. We speculate that the runaway-mass star may have collapsed directly into a massive black hole (Mbh > 100Msun). Such a dark object could explain the large velocity dispersion of the four Trapezium stars observed in the ONC core. We further show that the putative massive black hole is likely to be a member of a binary system with appr. 70 per cent probability. In such a case, it could be detected either due to short periods of enhanced accretion of stellar winds from the secondary star during pericentre passages, or through a measurement of the motion of the secondary whose velocity would exceed 10 km/s along the whole orbit.

preprint2012arXiv

Dwarf elliptical galaxies as ancient tidal dwarf galaxies

The formation of tidal dwarf galaxies (TDGs) is triggered by the encounters of already existing galaxies. Their existence is predicted from numerical calculations of encountering galaxies and is also well documented with observations. The numerical calculations on the formation of TDGs furthermore predict that TDGs cannot contain significant amounts of non-baryonic dark matter. In this paper, the first exhaustive sample of TDG-candidates from observations and numerical calculations is gathered from the literature. These stellar systems are gas-rich at the present, but they will probably evolve into gas-poor objects that are indistinguishable from old dwarf elliptical galaxies (dEs) based on their masses and radii. Indeed, known gas-poor TDGs appear as normal dEs. According to the currently prevailing cosmological paradigm, there should also be a population of primordial galaxies that formed within haloes of dark matter in the same mass range. Due to their different composition and origin, it would be expected that objects belonging to that population would have a different structure than TDGs and would thus be distinguishable from them, but such a population cannot be identified from their masses and radii. Moreover, long-lived TDGs could indeed be numerous enough to account for all dEs in the Universe. Downsizing, i.e. that less massive galaxies tend to be younger, would then be a natural consequence of the nature of the dEs. If these claims can be kept up in the light of future observations, the presently prevailing understanding of galaxy formation would need to be revised.

preprint2012arXiv

Evidence for top-heavy stellar initial mass functions with increasing density and decreasing metallicity

Residual-gas expulsion after cluster formation has recently been shown to leave an imprint in the low-mass present-day stellar mass function (PDMF) which allowed the estimation of birth conditions of some Galactic globular clusters (GCs) such as mass, radius and star formation efficiency. We show that in order to explain their characteristics (masses, radii, metallicity, PDMF) their stellar initial mass function (IMF) must have been top-heavy. It is found that the IMF is required to become more top-heavy the lower the cluster metallicity and the larger the pre-GC cloud-core density are. The deduced trends are in qualitative agreement with theoretical expectation. The results are consistent with estimates of the shape of the high-mass end of the IMF in the Arches cluster, Westerlund 1, R136 and NGC 3603, as well as with the IMF independently constrained for ultra-compact dwarf galaxies (UCDs). The latter suggests that GCs and UCDs might have formed along the same channel or that UCDs formed via mergers of GCs. A fundamental plane is found which describes the variation of the IMF with density and metallicity of the pre-GC cloud-cores simultaneously. The implications for the evolution of galaxies and chemical enrichment over cosmological times are expected to be major.

preprint2012arXiv

Gas surface density, star formation rate surface density, and the maximum mass of young star clusters in a disk galaxy. I. The flocculent galaxy M33

We analyze the relationship between maximum cluster mass, M_max, and surface densities of total gas (Sigma_gas), molecular gas (Sigma_H2) and star formation rate (Sigma_SFR) in the flocculent galaxy M33, using published gas data and a catalog of more than 600 young star clusters in its disk. By comparing the radial distributions of gas and most massive cluster masses, we find that M_max is proportional to Sigma_gas^4.7, M_max is proportional Sigma_H2^1.3, and M_max is proportional to Sigma_SFR^1.0. We rule out that these correlations result from the size of sample; hence, the change of the maximum cluster mass must be due to physical causes.

preprint2012arXiv

Inverse dynamical population synthesis: Constraining the initial conditions of young stellar clusters by studying their binary populations

Binary populations in young star clusters show multiplicity fractions both lower and up to twice as high as those observed in the Galactic field. We follow the evolution of a population of binary stars in dense and loose star clusters starting with an invariant initial binary population and a formal multiplicity fraction of unity, and demonstrate that these models can explain the observed binary properties in Taurus, Rho-Ophiuchus, Chamaeleon, Orion, IC 348, Upper Scorpius A, Praesepe, and the Pleiades. The model needs to consider solely different birth densities for these regions. The evolved theoretical orbital-parameter distributions are highly probable parent distributions for the observed ones. We constrain the birth conditions (stellar mass, M_ecl, and half-mass radius, r_h) for the derived progenitors of the star clusters and the overall present-day binary fractions allowed by the present model. The results compare very well with properties of molecular cloud clumps on the verge of star formation. Combining these with previously and independently obtained constraints on the birth densities of globular clusters, we identify a weak stellar mass -- half-mass radius correlation for cluster-forming cloud clumps, r_h / pc ~ (M_ecl / M_sun)^(0.13+-0.04). The ability of the model to reproduce the binary properties in all the investigated young objects, covering present-day densities from 1-10 stars pc^-3 (Taurus) to 2x10^4 stars pc^-3 (Orion), suggests that environment-dependent dynamical evolution plays an important role in shaping the present-day properties of binary populations in star clusters, and that the initial binary properties may not vary dramatically between different environments.

preprint2012arXiv

Low-mass X-ray binaries indicate a top-heavy stellar initial mass function in ultra compact dwarf galaxies

It has been shown before that the high mass-to-light ratios of ultra compact dwarf galaxies (UCDs) can be explained if their stellar initial mass function (IMF) was top-heavy, i.e. that the IMF was skewed towards high mass stars. In this case, neutron stars and black holes would provide unseen mass in the UCDs. In order to test this scenario with an independent method, we use data on which fraction of UCDs has a bright X-ray source. These X-ray sources are interpreted as low-mass X-ray binaries (LMXBs), i.e. binaries where a neutron star accretes matter from an evolving low-mass star. We find that LMXBs are indeed up to 10 times more frequent in UCDs than expected if the IMF was invariant. The top-heavy IMF required to account for this overabundance is the same as needed to explain the unusually high mass-to-light ratios of UCDs and a top-heavy IMF appears to be the only simultaneous explanation for both findings. Furthermore, we show that the high rate of type II supernovae (SNII) in the star-burst galaxy Arp 220 suggests a top-heavy IMF in that system. This finding is consistent with the notion that star-burst galaxies are sites where UCDs are likely to be formed and that the IMF of UCDs is top-heavy. It is estimated that the IMF becomes top-heavy whenever the star formation rate per volume surpasses 0.1 Solar masses per year and cubic parsec in pc-scale regions.

preprint2012arXiv

On the true shape of the upper end of the stellar initial mass function: The case of R136

Context: The shape of the stellar initial mass function (IMF) of a star cluster near its upper mass limit is a focal topic of investigation as it determines the high mass stellar content and hence the dynamics of the cluster at its embedded phase as well as during its young gas-free phase. The massive stellar content of a young cluster, however, can be substantially modified due to the dynamical ejections of the massive stars so that the present-day high-mass stellar mass function (hereafter MF) can be different than that with which the cluster is born. Aims: In the present study, we provide a preliminary estimate of this evolution of the high-mass IMF of a young cluster due to early ejections of massive stars, using the Large Magellanic Cloud massive, young cluster R136 as an example. Methods: To that end, we utilize the results of the state-of-the-art calculations by Banerjee, Kroupa & Oh (2012) comprising direct N-body computations of realistic, binary-rich, mass-segregated models of R136. In particular, these calculations provide the ejection fraction of stars as a function of stellar mass. Results: We find that if the measured IMF of R136 is granted to be canonical, as observations indicate, then the "true" high-mass IMF of R136 at its birth must be at least moderately top-heavy when corrected for the dynamical escape of massive stars. Conclusions: The top-heaviness of the true high-mass IMF over the observationally determined one is a general feature of massive, young clusters where the dynamical ejection of massive stars is efficient. We discuss its implications and possible improvements over our current estimate.

preprint2012arXiv

The dark matter crisis: falsification of the current standard model of cosmology

The current standard model of cosmology (SMoC) requires The Dual Dwarf Galaxy Theorem to be true according to which two types of dwarf galaxies must exist: primordial dark-matter (DM) dominated (type A) dwarf galaxies, and tidal-dwarf and ram-pressure-dwarf (type B) galaxies void of DM. Type A dwarfs surround the host approximately spherically, while type B dwarfs are typically correlated in phase-space. Type B dwarfs must exist in any cosmological theory in which galaxies interact. Only one type of dwarf galaxy is observed to exist on the baryonic Tully-Fisher plot and in the radius-mass plane. The Milky Way satellite system forms a vast phase-space-correlated structure that includes globular clusters and stellar and gaseous streams. Other galaxies also have phase-space correlated satellite systems. Therefore, The Dual Dwarf Galaxy Theorem is falsified by observation and dynamically relevant cold or warm DM cannot exist. It is shown that the SMoC is incompatible with a large set of other extragalactic observations. Other theoretical solutions to cosmological observations exist. In particular, alone the empirical mass-discrepancy--acceleration correlation constitutes convincing evidence that galactic-scale dynamics must be Milgromian. Major problems with inflationary big bang cosmologies remain unresolved.

preprint2012arXiv

The emergence of super-canonical stars in R136-type star-burst clusters

[abridged] Among the most remarkable features of the stellar population of R136, the central, young, massive star cluster in the 30 Doradus complex of the Large Magellanic Cloud, are the single stars whose masses substantially exceed the canonical stellar upper mass limit of 150 M_sun. A recent study by us, viz., that of Banerjee, Kroupa & Oh (2012; Paper I), which involves realistic N-body computations of star clusters mimicking R136, indicates that such "super-canonical" (SC) stars can be formed out of a dense stellar population with a canonical initial mass function (IMF) through dynamically induced mergers of the most massive binaries. Here we study the formation of SC stars in the R136 models of Paper I in detail. To avoid forming extraneous SC stars from initially highly eccentric primordial binaries as in Paper I, we compute additional models with only initially circular primordial binaries. We also take into account the mass-evolution of the SC stars using detailed stellar evolutionary models that incorporate updated treatments of stellar winds. We find that SC stars begin to form via dynamical mergers of massive binaries from approx. 1 Myr cluster age. We obtain SC stars with initial masses up to approx. 250 M_sun from these computations. Multiple SC stars are found to remain bound to the cluster simultaneously within a SC-lifetime. These properties of the dynamically formed SC stars are consistent with those observed in R136. In fact, the stellar evolutionary models of SC stars imply that had they formed primordially along with the rest of the R136 cluster, i.e., violating the canonical upper limit, they would have evolved below the canonical 150 M_sun limit by approx. 3 Myr, the likely age of R136, and would not have been observable as SC stars at the present time in R136. This strongly supports the dynamical formation scenario of the observed SC stars in R136.

preprint2012arXiv

The influence of stellar-dynamical ejections and collisions on the relation between the maximum-star and star-cluster-mass

We perform the largest currently available set of direct N-body calculations of young star cluster models to study the dynamical influence, especially through the ejections of the most massive star in the cluster, on the current relation between the maximum-stellar-mass and the star-cluster-mass. We vary several initial parameters such as the initial half-mass radius of the cluster, the initial binary fraction, and the degree of initial mass segregation. Two different pairing methods are used to construct massive binaries for more realistic initial conditions of massive binaries. We find that lower mass clusters (<= 10^2.5 Msun) do not shoot out their heaviest star. In the case of massive clusters (>= 1000 Msun), no most-massive star escapes the cluster within 3 Myr regardless of the initial conditions if clusters have initial half-mass radii, r_0.5, >= 0.8 pc. However, a few of the initially smaller sized clusters (r_0.5 = 0.3 pc), which have a higher density, eject their most massive star within 3 Myr. If clusters form with a compact size and their massive stars are born in a binary system with a mass-ratio biased towards unity, the probability that the mass of the most massive star in the cluster changes due to the ejection of the initially most massive star can be as large as 20 per cent. Stellar collisions increase the maximum-stellar-mass in a large number of clusters when clusters are relatively dense (M_ecl >= 10^3 Msun and r_0.5 = 0.3 pc) and binary-rich. Overall, we conclude that dynamical effects hardly influence the observational maximum-stellar-mass -- cluster mass relation.

preprint2012arXiv

The stellar and sub-stellar IMF of simple and composite populations

The current knowledge on the stellar IMF is documented. It appears to become top-heavy when the star-formation rate density surpasses about 0.1Msun/(yr pc^3) on a pc scale and it may become increasingly bottom-heavy with increasing metallicity and in increasingly massive early-type galaxies. It declines quite steeply below about 0.07Msun with brown dwarfs (BDs) and very low mass stars having their own IMF. The most massive star of mass mmax formed in an embedded cluster with stellar mass Mecl correlates strongly with Mecl being a result of gravitation-driven but resource-limited growth and fragmentation induced starvation. There is no convincing evidence whatsoever that massive stars do form in isolation. Various methods of discretising a stellar population are introduced: optimal sampling leads to a mass distribution that perfectly represents the exact form of the desired IMF and the mmax-to-Mecl relation, while random sampling results in statistical variations of the shape of the IMF. The observed mmax-to-Mecl correlation and the small spread of IMF power-law indices together suggest that optimally sampling the IMF may be the more realistic description of star formation than random sampling from a universal IMF with a constant upper mass limit. Composite populations on galaxy scales, which are formed from many pc scale star formation events, need to be described by the integrated galactic IMF. This IGIMF varies systematically from top-light to top-heavy in dependence of galaxy type and star formation rate, with dramatic implications for theories of galaxy formation and evolution.

preprint2012arXiv

Top-heavy IMFs in ultra compact dwarf galaxies?

Ultra compact dwarf galaxies (UCDs) are dense stellar systems at the border between massive star-clusters and small galaxies. Their on average high optical mass-to-light (M/L) ratio cannot be explained by stellar populations with the canonical stellar initial mass function (IMF), while it is doubtful that non-baryonic dark matter can accumulate enough on the scales of UCDs for influencing their dynamics significantly. UCDs in the Virgo galaxy cluster apparently also have an over-abundance of neutron stars, strongly suggesting a top-heavy IMF, which would explain both findings. This is because a top-heavy IMF can provide the unseen mass through an abundance of stellar remnants. The suggested variation of the IMF can be understood if UCDs represent a case of rapid star-formation in an extremely dense environment. While top-heavy IMFs imply a much heavier mass-loss shortly after the formation of a stellar system, this process does not necessarily dissolve the UCDs. Their formation with a top-heavy IMF would therefore not contradict their existence.

preprint2011arXiv

A new type of compact stellar population: dark star clusters

Among the most explored directions in the study of dense stellar systems is the investigation of the effects of the retention of supernova remnants, especially that of the massive stellar remnant black holes (BHs), in star clusters. By virtue of their eventual high central concentration, these stellar mass BHs potentially invoke a wide variety of physical phenomena, the most important ones being emission of gravitational waves (GWs), formation of X-ray binaries, and modification of the dynamical evolution of the cluster. Here we propose, for the first time, that rapid removal of stars from the outer parts of a cluster by the strong tidal field in the inner region of our Galaxy can unveil its BH sub-cluster, which appears as a star cluster that is gravitationally bound by an invisible mass. We study the formation and properties of such systems through direct N-body computations and estimate that they can be present in significant numbers in the inner region of the Milky Way. We call such objects "dark star clusters" (DSCs) as they appear dimmer than normal star clusters of similar mass and they comprise a predicted, new class of entities. The finding of DSCs will robustly cross-check BH retention; they will not only constrain the uncertain natal kicks of BHs, thereby the widely debated theoretical models of BH formation, but will also pinpoint star clusters as potential sites for GW emission for forthcoming ground-based detectors such as the Advanced LIGO. Finally, we also discuss the relevance of DSCs for the nature of IRS 13E.

preprint2011arXiv

An analytical description of the evolution of binary orbital-parameter distributions in N-body computations of star clusters

A new method is presented to describe the evolution of the orbital-parameter distributions for an initially universal binary population in star clusters by means of the currently largest existing library of N-body models. It is demonstrated that a stellar-dynamical operator exists, which uniquely transforms an initial orbital parameter distribution function for binaries into a new distribution depending on the initial cluster mass and half-mass radius, after some time of dynamical evolution. For the initial distribution the distribution functions derived by Kroupa (1995a,b) are used, which are consistent with constraints for pre-main sequence and Class I binary populations. Binaries with a lower energy and a higher reduced-mass are dissolved preferentially. The stellar-dynamical operator can be used to efficiently calculate and predict binary properties in clusters and whole galaxies without the need for further N-body computations. For the present set of N-body models it is found that the binary populations change their properties on a crossing time-scale such that the stellar dynamical operator can be well parametrized as a function of the initial cluster density. Furthermore it is shown that the binary-fraction in clusters with similar initial velocity dispersions follows the same evolutionary tracks as a function of the passed number of relaxation-times. Present-day observed binary populations in star clusters put constraints on their initial stellar densities which are found to be in the range 10^2 - 2x10^5 M_sun pc^-3 for open clusters and a few x 10^3 - 10^8 M_sun pc^-3 for globular clusters, respectively.

preprint2011arXiv

Disruption of Dwarf Satellite Galaxies without Dark Matter

The evolution of a satellite galaxy of a Milky Way like galaxy has been studied using N-Body simulations. The initial satellites, containing one million particles, have been simulated by a Plummer sphere, while the potential of the host galaxy is a three component rigid potential: disc, bulge and dark matter halo. It has been found that several orbits of the satellites allow for the existence, for about 1 Gyr or more, of an out-of-equilibrium body that could be interpreted as a dSph satellite galaxy of the Milky Way. In addition, from the study of the evolution of the mass-to-light ratios of satellites that show a disrupted phase it has been found that it is possible that some dSph galaxies of the Milky Way with large M/L ratios might not be dark matter dominated and that their high mass to light ratios are observed because they are out of equilibrium objects.

preprint2011arXiv

Distant star clusters of the Milky Way in MOND

We determine the mean velocity dispersion of six Galactic outer halo globular clusters, AM 1, Eridanus, Pal 3, Pal 4, Pal 15, and Arp 2 in the weak acceleration regime to test classical vs. modified Newtonian dynamics (MOND). Owing to the non-linearity of MOND's Poisson equation, beyond tidal effects, the internal dynamics of clusters is affected by the external field in which they are immersed. For the studied clusters, particle accelerations are much lower than the critical acceleration a_0 of MOND, but the motion of stars is neither dominated by internal accelerations (a_i >> a_e) nor external accelerations (a_e >> a_i). We use the N-body code N-MODY in our analysis, which is a particle-mesh-based code with a numerical MOND potential solver developed by Ciotti, Londrillo, and Nipoti (2006) to derive the line-of-sight velocity dispersion by adding the external field effect. We show that Newtonian dynamics predicts a low-velocity dispersion for each cluster, while in modified Newtonian dynamics the velocity dispersion is much higher. We calculate the minimum number of measured stars necessary to distinguish between Newtonian gravity and MOND with the Kolmogorov-Smirnov test. We also show that for most clusters it is necessary to measure the velocities of between 30 to 80 stars to distinguish between both cases. Therefore the observational measurement of the line-of-sight velocity dispersion of these clusters will provide a test for MOND.

preprint2011arXiv

Do all O stars form in star clusters?

The question whether or not massive stars can form in isolation or only in star clusters is of great importance for the theory of (massive) star-formation as well as for the stellar initial mass function of whole galaxies (IGIMF-theory). While a seemingly easy question it is rather difficult to answer. Several physical processes (e.g. star-loss due to stellar dynamics or gas expulsion) and observational limitations (e.g. dust obscuration of young clusters, resolution) pose severe challenges to answer this question. In this contribution we will present the current arguments in favour and against the idea that all O stars form in clusters.

preprint2011arXiv

Dynamical population synthesis: Constructing the stellar single and binary contents of galactic field populations

[abridged] The galactic field's late-type stellar single and binary population is calculated on the supposition that all stars form as binaries in embedded star clusters. A recently developed tool (Marks, Kroupa & Oh) is used to evolve the binary star distributions in star clusters for a few Myr so that a particular mixture of single and binary stars is achieved. On cluster dissolution the population enters the galactic field with these characteristics. The different contributions of single stars and binaries from individual star clusters which are selected from a power-law embedded star cluster mass function are then added up. This gives rise to integrated galactic field binary distribution functions (IGBDFs) resembling a galactic field's stellar content (Dynamical Population Synthesis). It is found that the binary proportion in the galactic field of a galaxy is larger the lower the minimum cluster mass, the lower the star formation rate, the steeper the embedded star cluster mass function and the larger the typical size of forming star clusters in the considered galaxy. In particular, period-, mass-ratio- and eccentricity IGBDFs for the Milky Way are modelled. The afore mentioned theoretical IGBDFs agree with independently observed distributions. Of all late-type binaries, 50% stem from M<300M_sun clusters, while 50% of all single stars were born in M>10^4M_sun clusters. Comparison of the G-dwarf and M-dwarf binary population indicates that the stars formed in mass-segregated clusters. In particular it is pointed out that although in the present model all M-dwarfs are born in binary systems, in the Milky Way's Galactic field the majority ends up being single stars. This work predicts that today's binary frequency in elliptical galaxies is lower than in spiral and in dwarf-galaxies. The period and mass-ratio distributions in these galaxies are explicitly predicted.

preprint2011arXiv

Runaway massive stars from R136: VFTS 682 is very likely a "slow runaway"

We conduct a theoretical study on the ejection of runaway massive stars from R136 --- the central massive, star-burst cluster in the 30 Doradus complex of the Large Magellanic Cloud. Specifically, we investigate the possibility of the very massive star (VMS) VFTS 682 being a runaway member of R136. Recent observations of the above VMS, by virtue of its isolated location and its moderate peculiar motion, have raised the fundamental question whether isolated massive star formation is indeed possible. We perform the first realistic N-body computations of fully mass-segregated R136-type star clusters in which all the massive stars are in primordial binary systems. These calculations confirm that the dynamical ejection of a VMS from a R136-like cluster, with kinematic properties similar to those of VFTS 682, is common. Hence the conjecture of isolated massive star formation is unnecessary to account for this VMS. Our results are also quite consistent with the ejection of 30 Dor 016, another suspected runaway VMS from R136. We further note that during the clusters' evolution, mergers of massive binaries produce a few single stars per cluster with masses significantly exceeding the canonical upper-limit of 150 solar mass. The observations of such single super-canonical stars in R136, therefore, do not imply an IMF with an upper limit greatly exceeding the accepted canonical 150 solar mass limit, as has been suggested recently, and they are consistent with the canonical upper limit.

preprint2011arXiv

Star cluster formation and some implications for GAIA

Stars form in spatially and temporarily correlated star formation events (CSFEs) and the dynamical processes within these "embedded clusters" leave imprints in the stellar populations in galactic fields. Such imprints are correlations in phase space (e.g. gravitationally bound star clusters, tidal streams), in the binary properties of stars and in the present-day stellar mass functions in the surviving clusters. The dynamical processes include expulsion of massive stars from cluster cores, disruption of CSFEs due to residual gas expulsion and energy-equipartition driven evaporation of stars from clusters leading to dark star clusters and cold kinematical streams with epicyclic overdensities. The properties of such phase-space structures in the Milky Way (MW) field depend on the effective gravitational potential of the MW. GAIA data will significantly constrain all of these aspects, and will in particular impact on gravitational dynamics via the properties of cold streams and on star-formation via the constraint on the gas expulsion process through the expanding unbound populations that must be associated with every CSFE.

preprint2011arXiv

The dynamical fingerprint of gas-expulsion: Insights into the assembly of the Milky Ways' old GC system

Since the oldest globular clusters (GCs) are early residuals from the formation of the Milky Way (MW), GCs were exposed to the likely evolving potential of our Galaxy. The expulsion of the residual-gas from the GC's embedded progenitors is sensitive to the conditions in the pre-MW gas cloud. By means of N-body computations it is shown that gas throw-out from initially mass-segregated GCs affect the shape of the low-mass stellar mass function (MF) and that its imprint might still be visible in the present-day MF (PDMF). The strength of the tidal-field at birth influences the degree of gas-expulsion driven low-mass-star depletion and therefore the PDMF probes the MW potential at the time of GC formation. It is argued that among the old GC population in the MW, younger GCs show stronger low-mass-star loss than older GCs. This is shown to be consistent with a contracting and self-gravitating cloud in which fluctuations in the pre-MW potential grow with time. An initially relatively smooth tidal field evolved into a grainy potential within a dynamical time-scale of the collapsing cloud (based on Marks & Kroupa 2010, MNRAS, 406, 2000).

preprint2011arXiv

The initial period function of late-type binary stars and its variation

The variation of the period distribution function of late-type binaries is studied. It is shown that the Taurus--Auriga pre-main sequence population and the main sequence G dwarf sample do not stem from the same parent period distribution with better than 95 per cent confidence probability. The Lupus, Upper Scorpius A and Taurus--Auriga populations are shown to be compatible with being drawn from the same initial period function (IPF), which is inconsistent with the main sequence data. Two possible IPF forms are used to find parent distributions to various permutations of the available data which include Upper Scorpius B (UScB), Chameleon and Orion Nebula Cluster pre-main sequence samples. All the pre-main sequence samples studied here are consistent with the hypothesis that there exists a universal IPF which is modified through binary-star disruption if it forms in an embedded star cluster leading to a general decline of the observed period function with increasing period. The pre-main sequence data admit a log-normal IPF similar to that arrived at by Duquennoy & Mayor (1991) for main sequence stars, provided the binary fraction among pre-main sequence stars is significantly higher. But, for consistency with proto-stellar data, the possibly universal IPF ought to be flat in log-P or log-semi-major axis and must be similar to the K1 IPF form derived through inverse dynamical population synthesis, which has been shown to lead to the main sequence period function if most stars form in typical embedded clusters.

preprint2011arXiv

The Role of the IGIMF in the chemical evolution of the solar neighbourhood

The integrated galactic initial mass function (IGIMF) is computed from the combination of the stellar initial mass function (IMF) and the embedded cluster mass function, described by a power law with index beta. The result of the combination is a time-varying IMF which depends on the star formation rate. We applied the IGIMF formalism to a chemical evolution model for the solar neighbourhood and compared the results obtained by assuming three possible values for beta with the ones 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, higher Type Ia and II supernova rates, higher mass ejection rates and higher [alpha/Fe] values at a given metallicity. Our suggested fiducial value for beta is 2, since with this value we can account for most of the local observables. We discuss our results in a broader perspective, with some implications regarding the possible universality of the IMF and the importance of the star formation threshold.

preprint2011arXiv

Using dwarf satellite proper motions to determine their origin

The highly organised distribution of satellite galaxies surrounding the Milky Way is a serious challenge to the concordance cosmological model. Perhaps the only remaining solution, in this framework, is that the dwarf satellite galaxies fall into the Milky Way's potential along one or two filaments, which may or may not plausibly reproduce the observed distribution. Here we test this scenario by making use of the proper motions of the Fornax, Sculptor, Ursa Minor and Carina dwarf spheroidals, and trace their orbits back through several variations of the Milky Way's potential and account for dynamical friction. The key parameters are the proper motions and total masses of the dwarf galaxies. Using a simple model we find no tenable set of parameters that can allow Fornax to be consistent with filamentary infall, mainly because the 1 sigma error on its proper motion is relatively small. The other three must walk a tightrope between requiring a small pericentre (less than 20 kpc) to lose enough orbital energy to dynamical friction and avoiding being tidally disrupted. We then employed a more realistic model with host halo mass accretion, and found that the four dwarf galaxies must have fallen in at least 5 Gyrs ago. This time interval is longer than organised distribution is expected to last before being erased by the randomisation of the satellite orbits.

preprint2011arXiv

What is a Galaxy? Cast your vote here...

Although originally classified as galaxies, Ultra Compact Dwarfs (UCDs) share many properties in common with globular star clusters. The debate on the origin and nature of UCDs, and the recently discovered ultra-faint dwarf spheroidal (dSph) galaxies which contain very few stars, has motivated us to ask the question `what is a galaxy?' Our aim here is to promote further discussion of how to define a galaxy and, in particular, what separates it from a star cluster. Like most previous definitions, we adopt the requirement of a gravitationally bound stellar system as a minimum. In order to separate a dwarf galaxy from a globular cluster, we discuss other possible requirements, such as a minimum size, a long two-body relaxation time, a satellite system, the presence of complex stellar populations and non-baryonic dark matter. We briefly mention the implications of each of these definitions if they are adopted. Some special cases of objects with an ambiguous nature are also discussed. Finally, we give our favoured criteria, and in the spirit of a `collective wisdom', invite readers to vote on their preferred definition of a galaxy via a dedicated website.

preprint2010arXiv

Applications of the IGIMF-theory to the astrophysics of galaxies

The functional form of the galaxy-wide stellar initial mass function is of fundamental importance for understanding galaxies. So far this stellar initial mass function has been assumed to be identical to the IMF observed directly in star clusters. But because stars form predominantly in embedded groups rather than uniformly distributed over the whole galaxy, the galaxy-wide IMF needs to be calculated by adding all IMFs of all embedded groups. This integrated galactic stellar initial mass function (IGIMF) is steeper than the canonical IMF and steepens with decreasing SFR, leading to fundamental new insights and understanding of star forming properties of galaxies. This contribution reviews the existing applications of the IGIMF theory to galactic astrophysics, while the parallel contribution by Weidner, Pflamm-Altenburg & Kroupa (this volume) introduces the IGIMF theory.

preprint2010arXiv

Direct N-body simulations of globular clusters: (I) Palomar 14

We present the first ever direct $N$-body computations of an old Milky Way globular cluster over its entire life time on a star-by-star basis. Using recent GPU hardware at Bonn University, we have performed a comprehensive set of $N$-body calculations to model the distant outer halo globular cluster Palomar 14 (Pal 14). By varying the initial conditions we aim at finding an initial $N$-body model which reproduces the observational data best in terms of its basic parameters, i.e. half-light radius, mass and velocity dispersion. We furthermore focus on reproducing the stellar mass function slope of Pal 14 which was found to be significantly shallower than in most globular clusters. While some of our models can reproduce Pal 14's basic parameters reasonably well, we find that dynamical mass segregation alone cannot explain the mass function slope of Pal 14 when starting from the canonical Kroupa initial mass function (IMF). In order to seek for an explanation for this discrepancy, we compute additional initial models with varying degrees of primordial mass segregation as well as with a flattened IMF. The necessary degree of primordial mass segregation turns out to be very high. This modelling has shown that the initial conditions of Pal 14 after gas expulsion must have been a half-mass radius of about 20 pc, a mass of about 50000 M$_{\odot}$, and possibly some mass segregation or an already established non-canonical IMF depleted in low-mass stars. Such conditions might be obtained by a violent early gas-expulsion phase from an embedded cluster born with mass segregation. Only at large Galactocentric radii are clusters likely to survive as bound entities the destructive gas-expulsion process we seem to have uncovered for Pal 14. In addition we compute a model with a 5% primordial binary fraction to test if such a population has an effect on the cluster's evolution.

preprint2010arXiv

Initial conditions for globular clusters and assembly of the old globular cluster population of the Milky Way

Comparing N-body calculations that include primordial residual-gas expulsion with the observed properties of 20 Galactic globular clusters (GCs) for which the stellar mass function (MF) has been measured, we constrain the time-scale over which the gas of their embedded counterparts must have been removed, the star formation efficiency the progenitor cloud must have had and the strength of the tidal-field the clusters must have formed in. The three parameters determine the expansion and mass-loss during residual-gas expulsion. After applying corrections for stellar and dynamical evolution we find birth cluster masses, sizes and densities for the GC sample and the same quantities for the progenitor gas clouds. The pre-cluster cloud core masses were between 10^5-10^7 Msun and half-mass radii were typically below 1 pc and reach down to 0.2 pc. We show that the low-mass present day MF slope, initial half-mass radius and initial density of clusters correlates with cluster metallicity, unmasking metallicity as an important parameter driving cluster formation and the gas expulsion process. This work predicts that PD low-concentration clusters should have a higher binary fraction than PD high-concentration clusters. Since the oldest GCs are early residuals from the formation of the Milky Way (MW) and the derived initial conditions probe the environment in which the clusters formed, we use the results as a new tool to study the formation of the inner GC system of the Galaxy. We achieve time-resolved insight into the evolution of the pre-MW gas cloud on short time-scales (a few hundred Myr) via cluster metallicities. The results are shown to be consistent with a contracting and self-gravitating cloud in which fluctuations in the pre-MW potential grow with time. An initially relatively smooth tidal-field evolved into a grainy potential within a dynamical time-scale of the collapsing cloud.

preprint2010arXiv

Is the distant globular cluster Pal 14 in a deep-freeze?

We investigate the velocity dispersion of Pal 14, an outer Milky-Way globular cluster at Galactocentric distance of 71 kpc with a very low stellar density (central density 0.1-0.2 Msun/pc^3). Due to this low stellar density the binary population of Pal 14 is likely to be close to the primordial binary population. Artificial clusters are generated with the observed properties of Pal 14 and the velocity dispersion within these clusters is measured as Jordi et al. (2009) have done with 17 observed stars of Pal 14. We discuss the effect of the binary population on these measurements and find that the small velocity dispersion of 0.38 km/s which has been found by Jordi et al. (2009) would imply a binary fraction of less than 0.1, even though from the stellar density of Pal 14 we would expect a binary fraction of more than 0.5. We also discuss the effect of mass segregation on the velocity dispersion as possible explanation for this discrepancy, but find that it would increase the velocity dispersion further. Thus, either Pal 14 has a very unusual stellar population and its birth process was significantly different than we see in today's star forming regions, or the binary population is regular and we would have to correct the observed 0.38 km/s for binarity. In this case the true velocity dispersion of Pal 14 would be much smaller than this value and the cluster would have to be considered as "kinematically frigid", thereby possibly posing a challenge for Newtonian dynamics but in the opposite sense to MOND.

preprint2010arXiv

Simple stellar population models including blue stragglers

Observations show that nearly all star clusters and stellar populations contain blue straggler stars (BSs). BSs in a cluster can significantly enhance the integrated spectrum of the host population, preferentially at short wavelengths, and render it much bluer in photometric colours. Current theoretical simple stellar population (SSP) models constructed within the traditional framework of single and binary stellar evolution cannot fully account for the impact of these objects on the integrated spectral properties of stellar populations. Using conventional SSP models without taking into account BS contributions may significantly underestimate a cluster's age and/or metallicity, simply because one has to balance the observed bluer colours (or a bluer spectrum) with a younger age and/or a lower metallicity. Therefore, inclusion of BS contributions in SSP models is an important and necessary improvement for population synthesis and its applications. Here, we present a new set of SSP models, which include BS contributions based on our analysis of individual star clusters. The models cover the wavelength range from 91~Å to 160~$μ$m, ages from 0.1 to 20 Gyr and metallicities $Z=0.0004, 0.004, 0.008, 0.02$ (solar) and 0.05. We use the observed integrated spectra of several Magellanic Cloud star clusters to cross-check and validate our models. The results show that the age predictions from our models are closer to those from isochrone fitting in the clusters' colour-magnitude diagrams compared to age predictions based on standard SSP models.

preprint2010arXiv

The [alpha/Fe] Ratios in Dwarf Galaxies: Evidence for a Non-universal Stellar Initial Mass Function?

It is well established that the [alpha/Fe] ratios in elliptical galaxies increase with galaxy mass. This relation holds also for early-type dwarf galaxies, although it seems to steepen at low masses. The [alpha/Fe] vs. mass relation can be explained assuming that smaller galaxies form over longer timescales (downsizing), allowing a larger amount of Fe (mostly produced by long-living Type Ia Supernovae) to be released and incorporated into newly forming stars. Another way to obtain the same result is by using a flatter initial mass function (IMF) in large galaxies, increasing in this way the number of Type II Supernovae and therefore the production rate of alpha-elements. The integrated galactic initial mass function (IGIMF) theory predicts that the higher the star formation rate, the flatter the IMF. We have checked, by means of semi-analytical calculations, that the IGIMF theory, combined with the downsizing effect (i.e. the shorter duration of the star formation in larger galaxies), well reproduces the observed [alpha/Fe] vs. mass relation. In particular, we show a steepening of this relation in dwarf galaxies, in accordance with the available observations.

preprint2010arXiv

The coupling of a young stellar disc with the molecular torus in the Galactic centre

The Galactic centre hosts, according to observations, a number of early-type stars. About one half of those which are orbiting the central supermassive black hole on orbits with projected radii $\gtrsim$ 0.03 pc form a coherently rotating disc. Observations further reveal a massive gaseous torus and a significant population of late-type stars. In this paper, we investigate, by means of numerical N-body computations, the orbital evolution of the stellar disc, which we consider to be initially thin. We include the gravitational influence of both the torus and the late-type stars, as well as the self-gravity of the disc. Our results show that, for a significant set of system parameters, the evolution of the disc leads, within the lifetime of the early-type stars, to a configuration compatible with the observations. In particular, the disc naturally reaches a specific - perpendicular - orientation with respect to the torus, which is indeed the configuration observed in the Galactic centre. We, therefore, suggest that all the early-type stars may have been born within a single gaseous disc.

preprint2010arXiv

The galaxy-wide IMF - from star clusters to galaxies

Over the past years observations of young and populous star clusters have shown that the stellar initial mass function (IMF) can be conveniently described by a two-part power-law with an exponent alpha_2 = 2.3 for stars more massive than about 0.5 Msol and an exponent of alpha_1 = 1.3 for less massive stars. A consensus has also emerged that most, if not all, stars form in stellar groups and star clusters, and that the mass function of these can be described as a power-law (the embedded cluster mass function, ECMF) with an exponent beta ~2. These two results imply that the integrated galactic IMF (IGIMF) for early-type stars cannot be a Salpeter power-law, but that they must have a steeper exponent. An application to star-burst galaxies shows that the IGIMF can become top-heavy. This has important consequences for the distribution of stellar remnants and for the chemo-dynamical and photometric evolution of galaxies. In this contribution the IGIMF theory is described, and the accompanying contribution by Pflamm-Altenburg, Weidner & Kroupa (this volume) documents the applications of the IGIMF theory to galactic astrophysics.

preprint2010arXiv

The puzzle of the cluster-forming core mass-radius relation and why it matters

We highlight how the mass-radius relation of cluster-forming cores combined with an external tidal field can influence infant weight-loss and disruption likelihood of clusters after gas expulsion. Specifically, we study how the relation between the bound fraction of stars staying in clusters at the end of violent relaxation and the cluster-forming core mass is affected by the slope and normalization of the core mass-radius relation. Assuming mass-independent star formation efficiency and gas-expulsion time-scale $τ_{GExp}/τ_{cross}$ and a given external tidal field, it is found that constant surface density cores and constant radius cores have the potential to lead to the preferential removal of high- and low-mass clusters, respectively. In contrast, constant volume density cores result in mass-independent cluster infant weight-loss, as suggested by observations. Our modelling includes predictions about the evolution of high-mass cluster-forming cores, a regime not yet covered by the observations. An overview of various issues directly affected by the nature of the core mass-radius relation is presented (e.g. cluster mass function, galaxy star formation histories, globular cluster self-enrichment). Finally, we emphasize that observational mass-radius data-sets of dense gas regions must be handled with caution as they may be the imprint of the molecular tracer used to map them, rather than reflecting cluster formation conditions. [Abridged]

preprint2010arXiv

The star formation history of the Large Magellanic Cloud as seen by star clusters and stars

The aim of this work is to test to what extent the star cluster population of a galaxy can be utilised to constrain or estimate the star formation history, with the Large Magellanic Cloud as our testbed. We follow two methods to extract information about the star formation rate from star clusters, either using only the most massive clusters (Maschberger & Kroupa 2007) or using the whole cluster population, albeit this is only possible for a shorter age span. We compare these results with the star formation history derived from colour-magnitude diagrams and find good overall agreement for the most recent approximately 1 Gyr. For later ages, and especially during the "cluster age gap", there is a deficiency of star clusters in relation to the star formation rate derived from the colour-magnitude diagram. The star formation rates following from the whole cluster population lie a factor of approximately 10 lower than the star formation rates deduced from the most massive clusters or from the colour-magnitude diagram, suggesting that only approximately 10 per cent of all stars form in long-lived bound star clusters.

preprint2010arXiv

The universality hypothesis: binary and stellar populations in star clusters and galaxies

It is possible to extract, from the observations, distribution functions of the birth dynamical properties of a stellar population, and to also infer that these are quite invariant to the physical conditions of star formation. The most famous example is the stellar IMF, and the initial binary population (IBP) seems to follow suit. A compact mathematical formulation of the IBP can be derived from the data. It has three broad parts: the IBP of the dominant stellar population (0.08-2 M_sol), the IBP of the more-massive stars and the IBP of brown dwarfs. These three mass regimes correspond to different physical regimes of star formation but not to structure in the IMF. With this formulation of the IBP it becomes possible to synthesise the stellar-population of whole galaxies.

preprint2010arXiv

Tidally induced brown dwarf and planet formation in circumstellar discs

Most stars are born in clusters and the resulting gravitational interactions between cluster members may significantly affect the evolution of circumstellar discs and therefore the formation of planets and brown dwarfs. Recent findings suggest that tidal perturbations of typical circumstellar discs due to close encounters may inhibit rather than trigger disc fragmentation and so would seem to rule out planet formation by external tidal stimuli. However, the disc models in these calculations were restricted to disc radii of 40 AU and disc masses below 0.1 M_sun. Here we show that even modest encounters can trigger fragmentation around 100 AU in the sorts of massive (~0.5 M_sun), extended (>=100 AU) discs that are observed around young stars. Tidal perturbation alone can do this, no disc-disc collision is required. We also show that very-low-mass binary systems can form through the interaction of objects in the disc. In our computations, otherwise non-fragmenting massive discs, once perturbed, fragment into several objects between about 0.01 and 0.1 M_sun, i.e., over the whole brown dwarf mass range. Typically these orbit on highly eccentric orbits or are even ejected. While probably not suitable for the formation of Jupiter- or Neptune-type planets, our scenario provides a possible formation mechanism for brown dwarfs and very massive planets which, interestingly, leads to a mass distribution consistent with the canonical substellar IMF. As a minor outcome, a possible explanation for the origin of misaligned extrasolar planetary systems is discussed.

preprint2010arXiv

Top-heavy integrated galactic stellar initial mass functions (IGIMFs) in starbursts

Star formation rates (SFR) larger than 1000 Msun/ yr are observed in extreme star bursts. This leads to the formation of star clusters with masses > 10^6 Msun in which crowding of the pre-stellar cores may lead to a change of the stellar initial mass function (IMF). Indeed, the large mass-to-light ratios of ultra-compact dwarf galaxies and recent results on globular clusters suggest the IMF to become top-heavy with increasing star-forming density. We explore the implications of top-heavy IMFs in these very massive and compact systems for the integrated galactic initial mass function (IGIMF), which is the galaxy-wide IMF, in dependence of the star-formation rate of galaxies. The resulting IGIMFs can have slopes, alpha_3, for stars more massive than about 1 Msun between 1.5 and the Salpeter slope of 2.3 for an embedded cluster mass function (ECMF) slope (beta) of 2.0, but only if the ECMF has no low-mass clusters in galaxies with major starbursts. Alternatively, beta would have to decrease with increasing SFR >10 Msun/ yr such that galaxies with major starbursts have a top-heavy ECMF. The resulting IGIMFs are within the range of observationally deduced IMF variations with redshift.

preprint2010arXiv

Two step ejection of massive stars and the issue of their formation in isolation

In this paper we investigate the combined effect of massive binary ejection from star clusters and a second acceleration of a massive star during a subsequent supernova explosion. We call this the "two-step-ejection" scenario. The main results are: i) Massive field stars produced via the two-step-ejection process can not in the vast majority of cases be traced back to their parent star clusters. These stars can be mistakenly considered as massive stars formed in isolation. ii) The expected O star fraction produced via the two-step-ejection process is of the order of 1-4 per cent, in quantitative agreement with the observed fraction of candidates for isolated-O-star formation. iii) Stars ejected via the two-step-ejection process can get a higher final velocity (up to 1.5-2 times higher) than the pre-supernova velocity of the massive-star binary.

preprint2009arXiv

Clustered star formation as a natural explanation of the Halpha cutoff in disc galaxies

Star formation is mainly determined by the observation of H$α$ radiation which is related to the presence of short lived massive stars. Disc galaxies show a strong cutoff in H$α$ radiation at a certain galactocentric distance which has led to the conclusion that star formation is suppressed in the outer regions of disc galaxies. This is seemingly in contradiction to recent UV observations (Boissier et al., 2007) that imply disc galaxies to have star formation beyond the Halpha cutoff and that the star-formation-surface density is linearly related to the underlying gas surface density being shallower than derived from Halpha luminosities (Kennicutt, 1998). In a galaxy-wide formulation the clustered nature of star formation has recently led to the insight that the total galactic Halpha luminosity is non-linearly related to the galaxy-wide star formation rate (Pflamm-Altenburg et al., 2007d). Here we show that a local formulation of the concept of clustered star formation naturally leads to a steeper radial decrease of the Halpha surface luminosity than the star-formation-rate surface density in quantitative agreement with the observations, and that the observed Halpha cutoff arises naturally.

preprint2009arXiv

Constraining the initial mass function of stars in the Galactic Centre

(abridged) Here we discuss the question whether the extreme circumstances in the centre of the Milky Way may be the reason for a significant variation of the IMF. By means of stellar evolution models using different codes we show that the observed luminosity in the central parsec is too high to be explained by a long-standing top-heavy IMF, considering the limited amount of mass inferred from stellar kinematics in this region. In contrast, continuous star formation over the Galaxy's lifetime following a canonical IMF results in a mass-to-light ratio and a total mass of stellar black holes (SBHs) consistent with the observations. Furthermore, these SBHs migrate towards the centre due to dynamical friction, turning the cusp of visible stars into a core as observed in the Galactic Centre. For the first time here we explain the luminosity and dynamical mass of the central cluster and both the presence and extent of the observed core, since the number of SBHs expected from a canonical IMF is just enough to make up for the missing luminous mass. We conclude that the Galactic Centre is consistent with the canonical IMF and do not suggest a systematic variation as a result of the region's properties such as high density, metallicity, strong tidal field etc.

preprint2009arXiv

Current data on the globular cluster Palomar 14 are not inconsistent with MOND

Certain types of globular clusters have the very important property that the predictions for their kinematics in the Newtonian and modified Newtonian dynamics (MOND) contexts are divergent. Here, we caution the recent claim that the stellar kinematics data (using 17 stars) of the globular cluster Palomar 14 are inconsistent with MOND. We compare the observations to the theoretical predictions using a Kolmogorov-Smirnov test, which is appropriate for small samples. We find that, with the currently available data, the MOND prediction for the velocity distribution can only be excluded with a very low confidence level, clearly insufficient to claim that MOND is falsified.

preprint2009arXiv

Did the Milky Way dwarf satellites enter the halo as a group?

The dwarf satellite galaxies in the Local Group are generally considered to be hosted in dark matter subhalos that survived the disruptive processes during infall onto their host halos. It has recently been argued that if the majority of satellites entered the Milky Way halo in a group rather than individually, this could explain the spatial and dynamical peculiarities of its satellite distribution. Such groups were identified as dwarf galaxy associations that are found in the nearby Universe. In this paper we address the question whether galaxies in such associations can be the progenitors of the Milky Way satellite galaxies. We find that the dwarf associations are much more extended than would be required to explain the disk-like distribution of the Milky Way and Andromeda satellite galaxies. We further identify a possible minor filamentary structure, perpendicular to the supergalactic plane, in which the dwarf associations are located, that might be related to the direction of infall of a progenitor galaxy of the Milky Way satellites, if they are of tidal origin.

preprint2009arXiv

Do binaries in clusters form in the same way as in the field?

We examine the dynamical destruction of binary systems in star clusters of different densities. We find that at high densities (10^4 - 10^5 Msun pc^-3) almost all binaries with separations > 10^3 AU are destroyed after a few crossing times. At low densities (order(10^2) Msun pc^-3) many binaries with separations > 10^3 AU are destroyed, and no binaries with separations > 10^4 AU survive after a few crossing times. Therefore the binary separations in clusters can be used as a tracer of the dynamical age and past density of a cluster. We argue that the central region of the Orion Nebula Cluster was around 100 times denser in the past with a half-mass radius of only 0.1 - 0.2 pc as (a) it is expanding, (b) it has very few binaries with separations > 10^3 AU, and (c) it is well-mixed and therefore dynamically old. We also examine the origin of the field binary population. Binaries with separations < 10^2 AU are not significantly modified in any cluster, therefore at these separations the field reflects the sum of all star formation. Binaries with separations in the range 10^2 - 10^4 AU are progressively more and more heavily affected by dynamical disruption in increasingly dense clusters. If most star formation is clustered, these binaries must be over-produced relative to the field. Finally, no binary with a separation > 10^4 AU can survive in any cluster and so must be produced by isolated star formation, but only if all isolated star formation produces extremely wide binaries.

preprint2009arXiv

Estimators for the exponent and upper limit, and goodness-of-fit tests for (truncated) power-law distributions

Many objects studied in astronomy follow a power law distribution function, for example the masses of stars or star clusters. A still used method by which such data is analysed is to generate a histogram and fit a straight line to it. The parameters obtained in this way can be severely biased, and the properties of the underlying distribution function, such as its shape or a possible upper limit, are difficult to extract. In this work we review techniques available in the literature and present newly developed (effectively) bias-free estimators for the exponent and the upper limit. The software packages are made available as downloads. Furthermore we discuss various graphical representations of the data and powerful goodness-of-fit tests to assess the validity of a power law for describing the distribution of data. As an example, we apply the presented methods to the data set of massive stars in R136 and the young star clusters in the Large Magellanic Cloud. (abridged)

preprint2009arXiv

Influence of a stellar cusp on the dynamics of young stellar discs and the origin of the S-stars in the Galactic Centre

Observations of the Galactic Centre show evidence of one or two disc-like structures of very young stars orbiting the central super-massive black hole within a distance of a few 0.1 pc. A number of analyses have been carried out to investigate the dynamical behaviour and consequences of these discs, including disc thickness and eccentricity growth as well as mutual interaction and warping. However, most of these studies have neglected the influence of the stellar cusp surrounding the black hole, which is believed to be 1-2 orders of magnitude more massive than the disc(s). By means of N-body integrations using our bhint code, we study the impact of stellar cusps of different compositions. We find that although the presence of a cusp does have an important effect on the evolution of an otherwise isolated flat disc, its influence on the evolution of disc thickness and warping is rather mild in a two-disc configuration. However, we show that the creation of highly eccentric orbits strongly depends on the graininess of the cusp (i.e. the mean and maximum stellar masses): While Chang (2009) recently found that full cycles of Kozai resonance are prevented by the presence of an analytic cusp, we show that relaxation processes play an important role in such highly dense regions and support short-term resonances. We thus find that young disc stars on initially circular orbits can achieve high eccentricities by resonant effects also in the presence of a cusp of stellar remnants, yielding a mechanism to create S-stars and hyper-velocity stars. Furthermore, we discuss the underlying initial mass function (IMF) of the young stellar discs and find no definite evidence for a non-canonical IMF.

preprint2009arXiv

Recurrent gas accretion by massive star clusters, multiple stellar populations and mass thresholds for spherical stellar systems

We explore the gravitational influence of pressure supported stellar systems on the internal density distribution of a gaseous environment. We conclude that compact massive star clusters with masses >= 10^6 M_sun act as cloud condensation nuclei and are able to accrete gas recurrently from a warm interstellar medium which may cause further star formation events and account for multiple stellar populations in the most massive globular and nuclear star clusters. The same analytical arguments can be used to decide whether an arbitrary spherical stellar system is able to keep warm or hot interstellar material or not. These mass thresholds coincide with transition masses between pressure supported galaxies of different morphological types.

preprint2009arXiv

Stellar-mass black holes in star clusters: implications for gravitational wave radiation

We study the dynamics of stellar-mass black holes (BH) in star clusters with particular attention to the formation of BH-BH binaries, which are interesting as sources of gravitational waves (GW). We examine the properties of these BH-BH binaries through direct N-body simulations of star clusters using the GPU-enabled NBODY6 code. We perform simulations of N <= 10^5 Plummer clusters of low-mass stars with an initial population of BHs. Additionally, we do several calculations of star clusters confined within a reflective boundary mimicking only the core of a massive cluster. We find that stellar-mass BHs with masses ~ 10 solar mass segregate rapidly into the cluster core and form a sub-cluster of BHs within typically 0.2 - 0.5 pc radius, which is dense enough to form BH-BH binaries through 3-body encounters. While most BH binaries are ejected from the cluster by recoils received during super-elastic encounters with the single BHs, few of them harden sufficiently so that they can merge via GW emission within the cluster. We find that for clusters with $N \ga 5\times 10^4$, typically 1 - 2 BH-BH mergers occur within them during the first ~ 4 Gyr of evolution. Also for each of these clusters, there are a few escaping BH binaries that can merge within a Hubble time, most of the merger times being within a few Gyr. These results indicate that intermediate-age massive clusters constitute the most important class of candidates for producing dynamical BH-BH mergers. Old globular clusters cannot contribute significantly to the present-day BH-BH merger rate since most of the mergers from them would have occurred earlier. In contrast, young massive clusters are too young to produce significant number of BH-BH mergers. Our results imply significant BH-BH merger detection rates for the proposed "Advanced LIGO" GW detector. (Abridged)

preprint2009arXiv

Testing fundamental physics with distant star clusters: theoretical models for pressure-supported stellar systems

We investigate the mean velocity dispersion and the velocity dispersion profile of stellar systems in MOND, using the N-body code N-MODY, which is a particle-mesh based code with a numerical MOND potential solver developed by Ciotti, Londrillo and Nipoti (2006). We have calculated mean velocity dispersions for stellar systems following Plummer density distributions with masses in the range of $10^4 M_\odot$ to $10^9 M_\odot$ and which are either isolated or immersed in an external field. Our integrations reproduce previous analytic estimates for stellar velocities in systems in the deep MOND regime ($a_i, a_e \ll a_0$), where the motion of stars is either dominated by internal accelerations ($a_i \gg a_e$) or constant external accelerations ($a_e \gg a_i$). In addition, we derive for the first time analytic formulae for the line-of-sight velocity dispersion in the intermediate regime ($a_i \sim a_e \sim a_0$). This allows for a much improved comparison of MOND with observed velocity dispersions of stellar systems. We finally derive the velocity dispersion of the globular cluster Pal 14 as one of the outer Milky Way halo globular clusters that have recently been proposed as a differentiator between Newtonian and MONDian dynamics.

preprint2009arXiv

The chemical evolution of galaxies within the IGIMF theory: the [alpha/Fe] ratios and downsizing

The chemical evolution of galaxies is investigated within the framework of the star formation rate (SFR) dependent integrated galactic initial mass function (IGIMF). We study how the global chemical evolution of a galaxy and in particular how [alpha/Fe] abundance ratios are affected by the predicted steepening of the IGIMF with decreasing SFR. We use analytical and semi-analytical calculations to evaluate the mass-weighted and luminosity-weighted [alpha/Fe] ratios in early-type galaxies of different masses. The models with the variable IGIMF produce a [alpha/Fe] vs. velocity dispersion relation which has the same slope as the observations of massive galaxies, irrespective of the model parameters, provided that the star formation duration inversely correlates with the mass of the galaxy (downsizing). These models also produce steeper [alpha/Fe] vs. sigma relations in low-mass early-type galaxies and this trend is consistent with the observations. Constant IMF models are able to reproduce the [alpha/Fe] ratios in large elliptical galaxies as well, but they do not predict this change of slope for small galaxies. In order to obtain the best fit between our results and the observations, the downsizing effect (i.e. the shorter duration of the star formation in larger galaxies) must be milder than previously thought.

preprint2009arXiv

The fundamental gas depletion and stellar-mass buildup times of star forming galaxies

Stars do not form continuously distributed over star forming galaxies. They form in star clusters of different masses. This nature of clustered star formation is taken into account in the theory of the integrated galactic stellar initial mass function (IGIMF) in which the galaxy-wide IMF (the IGIMF) is calculated by adding all IMFs of young star clusters. For massive stars the IGIMF is steeper than the universal IMF in star clusters and steepens with decreasing SFR which is called the IGIMF-effect. The current SFR and the total Halpha luminosity of galaxies therefore scale non-linearly in the IGIMF theory compared to the classical case in which the galaxy-wide IMF is assumed to be constant and identical to the IMF in star clusters. We here apply for the first time the revised SFR-L_Halpha relation on a sample of local volume star forming galaxies with measured Halpha luminosities. The fundamental results are: i) the SFRs of galaxies scale linearly with the total galaxy neutral gas mass, ii) the gas depletion time scales of dwarf irregular and large disk galaxies are about 3 Gyr implying that dwarf galaxies do not have lower star formation efficiencies than large disk galaxies, and iii) the stellar mass buildup times of dwarf and large galaxies are only in agreement with downsizing in the IGIMF context, but contradict downsizing within the traditional framework that assumes a constant galaxy-wide IMF.

preprint2008arXiv

The influence of star clusters on galactic disks: new insights on star-formation in galaxies

Stars form in embedded star clusters which play a key role in determining the properties of a galaxy's stellar population. Physical mechanisms discussed in this paper are runaway stars shot out from young clusters, binary-star disruption in clusters, gas blow-out from clusters and the origin of thick galactic disks. I emphasise that the SNIa rate per low-mass star depends on the star-clusters formed in a galaxy and I discuss the IGIMF theory. Based on the IGIMF theory, the re-calibrated Halpha-luminosity--SFR relation implies dwarf irregular galaxies to have the same gas-depletion time-scale as major disk galaxies, suggesting a major change in our understanding of dwarf-galaxy evolution. The IGIMF-theory also naturally leads to the observed radial Halpha cutoff in disk galaxies without a radial star-formation cutoff. It emerges that the thorough understanding of the physics and distribution of star clusters may be leading to a major paradigm shift in our understanding of galaxy evolution.

preprint2007arXiv

A discontinuity in the low-mass initial mass function

The origin of brown dwarfs (BDs) is still an unsolved mystery. While the standard model describes the formation of BDs and stars in a similar way recent data on the multiplicity properties of stars and BDs show them to have different binary distribution functions. Here we show that proper treatment of these uncovers a discontinuity of the multiplicity-corrected mass distribution in the very-low-mass star (VLMS) and BD mass regime. A continuous IMF can be discarded with extremely high confidence. This suggests that VLMSs and BDs on the one hand, and stars on the other, are two correlated but disjoint populations with different dynamical histories. The analysis presented here suggests that about one BD forms per five stars and that the BD-star binary fraction is about 2%-3% among stellar systems.

preprint2005arXiv

Limits on the primordial stellar multiplicity

Most stars - especially young stars - are observed to be in multiple systems. Dynamical evolution is unable to pair stars efficiently, which leads to the conclusion that star-forming cores must usually fragment into \geq 2 stars. However, the dynamical decay of systems with \geq 3 or 4 stars would result in a large single-star population that is not seen in the young stellar population. Additionally, ejections would produce a significant population of hard binaries that are not observed. This leads to a strong constraint on star formation theories that cores must typically produce only 2 or 3 stars. This conclusion is in sharp disagreement with the results of currently available numerical simulations that follow the fragmentation of molecular cores and typically predict the formation of 5--10 seeds per core. In addition, open cluster remnants may account for the majority of observed highly hierarchical higher-order multiple systems in the field.

preprint2004arXiv

Monte-Carlo experiments on star-cluster induced integrated-galaxy IMF variations

As most if not all stars are born in stellar clusters the shape of the mass function of the field stars is not only determined by the initial mass function of stars (IMF) but also by the cluster mass function (CMF). In order to quantify this Monte-Carlo simulations were carried out by taking cluster masses randomly from a CMF and then populating these clusters with stars randomly taken from an IMF. Two cases were studied. Firstly the star masses were added randomly until the cluster mass was reached. Secondly a number of stars, given by the cluster mass divided by an estimate of the mean stellar mass and sorted by mass, were added until the desired cluster mass was reached. Both experiments verified the analytical results of Kroupa & Weidner (2003) that the resulting integrated stellar initial mass function is a folding of the IMF with the CMF and therefore steeper than the input IMF above 1 Msol.

preprint2004arXiv

The great disk of Milky-Way satellites and cosmological sub-structures

We show that the shape of the observed distribution of Milky Way (MW) satellites is inconsistent with being drawn from a cosmological sub-structure population with a confidence of 99.5 per cent. Most of the MW satellites therefore cannot be related to dark-matter dominated satellites.

preprint2004arXiv

Variation of the IMF

(abridged) The {stellar IMF} has been found to be essentially invariant. While some apparent differences are seen, the uncertainties inherent to this game do not allow a firm conclusion to be made that the IMF varies systematically with conditions. The IMF integrated over entire galaxies, however, is another matter. Chemical and photometric properties of various galaxies do hint at {galaxial IMFs} being steeper than the stellar IMF, as is also deduced from direct star-count analysis in the MW. These results are sensitive to the modelling of stellar populations and to corrections for stellar evolution, and are thus also uncertain. However, by realising that galaxies are made from dissolving star clusters, star clusters being viewed as {the fundamental building blocks of galaxies}, the result is found that galaxial IMFs must be significantly steeper than the stellar IMF, because the former results from a folding of the latter with the star-cluster mass function. Furthermore, this notion leads to the important insight that galaxial IMFs must vary with galaxy mass, and that the galaxial IMF is a strongly varying function of the star-formation history for galaxies that have assembled only a small mass in stars. Cosmological implications of this are discussed.

preprint2001arXiv

On the variation of the Initial Mass Function

(shortened) In this contribution an average or Galactic-field IMF is defined, stressing that there is evidence for a change in the power-law index at only two masses: near 0.5 Msun and 0.08 Msun. Using this supposed universal IMF, the uncertainty inherent to any observational estimate of the IMF is investigated, by studying the scatter introduced by Poisson noise and the dynamical evolution of star clusters. It is found that this apparent scatter reproduces quite well the observed scatter in power-law index determinations, thus defining the fundamental limit within which any true variation becomes undetectable. Determinations of the power-law indices alpha are subject to systematic errors arising mostly from unresolved binaries. The systematic bias is quantified here, with the result that the single-star IMFs for young star-clusters are systematically steeper by d_alpha=0.5 between 0.1 and 1 Msun than the Galactic-field IMF, which is populated by, on average, about 5 Gyr old stars. The MFs in globular clusters appear to be, on average, systematically flatter than the Galactic-field IMF, and the recent detection of ancient white-dwarf candidates in the Galactic halo and absence of associated low-mass stars suggests a radically different IMF for this ancient population. Star-formation in higher-metallicity environments thus appears to produce relatively more low-mass stars.

preprint1997arXiv

The Theoretical Mass--Magnitude Relation of Low-Mass Stars and its Metallicity Dependence

We investigate the dependence of theoretically generated mass - (absolute magnitude) relations on stellar models. Using up to date physics we compute models in the mass range 0.1 < m < 1M_sun. We compare the solar-metallicity models with our older models, with recent models computed by others, and also with an empirical mass - (absolute magnitude) relation that best fits the observed data. At a given mass below 0.6M_sun the effective temperatures differ substantially from model to model. However taken individually each set of models is in good agreement with observations in the mass - luminosity plane. A minimum in the derivative dm/dM_V at M_V = 11.5, which is due to H_2 formation and establishment of a fully convective stellar interior, is present in all photometric bands, for all models. This minimum leads to a maximum in the stellar luminosity function for Galactic disk stars at M_V = 11.5, M_bol = 9.8. Stellar models should locate this maximum in the stellar luminosity function at the same magnitude as observations. Models which incorporate the most realistic theoretical atmospheres and the most recent equation of state and opacities can satisfy this constraint. These models are also in best agreement with the most recent luminosity - (effective temperature) and mass-luminosity data. Each set of our models of a given metallicity (with 0.2 > [Fe/H] > -2.3) shows a maximum in -dm/dM_bol, which moves to brighter bolometric magnitudes with decreasing metallicity. The change in location of the maximum, as a function of [Fe/H], follows the location of structure in luminosity functions for stellar populations with different metal abundances. This structure seen in all observed stellar populations can be accounted for by the mass--luminosity relation.

preprint1995arXiv

Inverse Dynamical Population Synthesis and Star Formation

We investigate the hypothesis that stars form in aggregates of binary systems and that the dynamical evolution of these aggregates leads to the observed properties of binary stars in the Galactic field. We assume that the initial distribution of periods is flat in logP, where P is the orbital period in days, and 3<logP<7.5. We distribute 200 binaries in aggregates with half mass radii R corresponding to the range from tightly clustered to isolated star formation, and follow the subsequent evolution of the stellar systems by direct N-body integration. Hardening and softening of binary systems do not significantly increase the number of orbits with logP<3 and logP>7.5. After the cluster with R=0.8pc disintegrates we obtain a population which consists of about 60 per cent binary systems with a period distribution for logP>4 as is observed and in which the G-dwarf binaries have a mass ratio distribution which agrees with the observed distribution. This result indicates that the majority of Galactic field stars may originate from a clustered star formation mode. We invert the orbit depletion function and obtain an approximation to the initial binary star period distribution for star formation in the dominant mode cluster. Comparison with the measured distribution of orbits for pre-main sequence stars suggests that the initial distribution may not depend on the star formation environment. Inverse dynamical population synthesis suggests that the Galactic field stellar mass function may be related to the stellar density at birth in the most common, or dominant, mode of star formation.

preprint1995arXiv

Star Cluster Evolution, Dynamical Age Estimation and the Kinematical Signature of Star Formation

We use N-body integration to follow the evolution of clusters of 200 binary systems with different initial half mass radii $R_{0.5}$. We also simulate single-star clusters. All clusters evolve according to the same $n(t)$ curve, where $n(t)$ is the number density of stars in the central 2~pc sphere at time $t$. $n(t)$ and the lifetime are independent of (i) the inital proportion of binaries and (ii) the initial $R_{0.5}$. Mass segregation measures the dynamical age of the cluster. The proportion of binaries in the central cluster region is a sensitive indicator of the initial cluster concentration. If most stars form in binaries in a typical embedded cluster which is located at the edge of a giant molecular cloud, then we estimate that at most about 10~per cent of all pre-main sequence stars achieve near escape velocities from the molecular cloud. The large ejection velocities resulting from close encounters between binary systems imply a `halo' distribution of young stars over large areas surrounding star forming sites which is expected to have a significantly reduced binary proportion and a significantly increased proportion of stars with depleted circumstellar disks. We compare the time dependent model single star and system luminosity function in the central cluster region with the observational Hyades and Pleiades luminosity functions and find no evidence for different dynamical properties of stellar systems at birth in the Hyades, Pleiades and Galactic field stellar samples. The observed proportion of binary stars in the very young Trapezium Cluster is consistent with the early dynamical evolution of a cluster with a very high initial stellar number density.

preprint1995arXiv

The Dynamical Properties of Stellar Systems in the Galactic Disc

We postulate that most stars are born in aggregates of binary systems which are dynamically equivalent to the `dominant mode cluster'. The initial binary orbits are consitent with pre-main sequence data. Stellar masses are paired at random from the KTG(1.3) mass function. We develop a simple model which establishes the observed correlations between eccentricity, mass ratio and period. Stimulated evolution of orbital parameters within the cluster places a few per cent of all orbits into the eigenevolution region of the eccentricity- period diagram. Observations should reveal the odd binary with eccentric orbits at periods smaller than the circularisation cutoff period. We show that the birth eccentricity distribution must be approximately dynamically relaxed because stimulated evolution in the dominant mode cluster cannot sufficiently thermalise a significantly different distribution. After disintegration of the dominant mode cluster we have a proportion of binaries which compares favourably with the observed proportion. The shape of the period distribution agrees well with the observed distribution. The mass-ratio distribution of G dwarf binaries is depleted at small mass ratios and has the shape of the main sequence distribution. We predict and tabulate the mass ratio distribution for main sequence binaries with a primary star less massive than $1.1\,M_\odot$. Our model Galactic field stellar population has a binary proportion among G-, K- and M-dwarfs in good agreement with the observational data. Too few triple and quadruple systems form by capture to account for the number of observed systems. We compare the specific angular momentum distribution of our initial binary star population with the observed distribution of specific angular momenta of molecular cloud cores.

Pavel Kroupa

What is connected

Connect this record

See the researcher in context

Building this map preview

138 published item(s)

Intermediate-mass black hole incubators. Gas accretion onto stellar black hole clusters in galactic central molecular zones

Tidal Tails of Nearby Open Clusters. II. A Review of Simulated Properties and the Reliability of Observational Catalogues

3D hydrodynamic simulations for the formation of the Local Group satellite planes

Do the majority of stars form as gravitationally unbound?

Far-ultraviolet investigation into the galactic globular cluster M30 (NGC 7099): I. Photometry and radial distributions

Far-ultraviolet investigation into the galactic globular cluster M30 (NGC 7099): II. Potential X-ray counterparts and variable sources

Origin of the spectacular tidal shells of galaxy NGC474

Overestimated inclinations of Milgromian disc galaxies: the case of the ultradiffuse galaxy AGC 114905

The formation of early-type galaxies through monolithic collapse of gas clouds in Milgromian gravity

The global stability of M33 in MOND

The High Fraction of Thin Disk Galaxies Continues to Challenge ΛCDM Cosmology

Estimating the ages of open star clusters from properties of their extended tidal tails

HI mapping of the Leo Triplet: Morphologies and kinematics of tails and bridges

On the dynamical evolution of Cepheids in star clusters

A correlation between the number of satellites and the baryonic bulge-to-total ratio extending beyond the Local Group

Chemical evolution of ultra-faint dwarf galaxies in the self-consistently calculated IGIMF theory

Constraints on the star formation histories of galaxies in the Local Cosmological Volume

Scale-invariant dynamics in the Solar System

Solar System limits on gravitational dipoles

The failure of testing for cosmic opacity via the distance-duality relation

The formation of exponential disk galaxies in MOND

Tidal tails of open star clusters as probes to early gas expulsion I: A semi-analytic model

Tidal tails of open star clusters as probes to early gas expulsion II: Predictions for Gaia

Very high redshift quasars and the rapid emergence of super-massive black holes

The Milky Way's Disk of Classical Satellite Galaxies in Light of Gaia DR2

The systematically varying stellar IMF

The ultra-diffuse dwarf galaxies NGC 1052-DF2 and 1052-DF4 are in conflict with standard cosmology

Declining rotation curves of galaxies as a test of gravitational theory

Dynamical ejections of massive stars from young star clusters under diverse initial conditions

On the origin of the Schechter-like mass function of young star clusters in disk galaxies

Simulating disk galaxies and interactions in Milgromian dynamics

Star formation triggered by galaxy interactions in modified gravity

The observed spatial distribution of matter on scales ranging from 100kpc to 1Gpc is inconsistent with the standard dark-matter-based cosmological models

Comment on "Evidence for dark matter in the inner Milky Way"

Dependency of dynamical ejections of O stars on the masses of very young star clusters

Dynamical Evolution of Outer-Halo Globular Clusters

Galaxies as simple dynamical systems: observational data disfavor dark matter and stochastic star formation

Globular Cluster Streams as Galactic High-Precision Scales - The Poster Child Palomar 5

How can young massive clusters reach their present-day sizes?

M-dwarf binaries as tracers of star and brown dwarf formation

Mass distributions of star clusters for different star formation histories in a galaxy cluster environment

On the persistence of two small-scale problems in ΛCDM

Possible smoking-gun evidence for initial mass segregation in re-virialized post-gas expulsion globular clusters

Probing the isotropy of cosmic acceleration traced by Type Ia supernovae

The number of tidal dwarf satellite galaxies in dependence of bulge index

A census of the expected properties of classical Milky Way dwarfs in Milgromian dynamics

A perfect starburst cluster made in one go: the NGC 3603 young cluster

Chemo-dynamical evolution of tidal dwarf galaxies. II. The long-term evolution and influence of a tidal field

Co-orbiting satellite galaxy structures are still in conflict with the distribution of primordial dwarf galaxies

Direct $N$-body simulations of globular clusters - II. Palomar 4

Erosion of Globular Cluster Systems: The Influence of Radial Anisotropy, Central Black Holes and Dynamical Friction

Galactic rotation curves, the baryon-to-dark-halo-mass relation and space-time scale invariance

Phantom of RAMSES (POR): A new Milgromian dynamics N-body code

Reproducing properties of MW dSphs as descendants of DM-free TDGs

Sampling methods for stellar masses and the mmax-Mecl relation in the starburst dwarf galaxy NGC 4214

Simulations on the survivability of Tidal Dwarf Galaxies

The formation of NGC 3603 young starburst cluster: "prompt" hierarchical assembly or monolithic starburst?

The Puzzling Negative Orbit-Period Derivative of the Low-Mass X-Ray Binary 4U 1820-30 in NGC 6624

The state of globular clusters at birth II: primordial binaries

The state of globular clusters at birth: emergence from the gas-embedded phase

The Vast Polar Structure of the Milky Way Attains New Members

Did the infant R136 and NGC 3603 clusters undergo residual gas expulsion?

Dwarf Galaxy Planes: the discovery of symmetric structures in the Local Group

Gas surface density, star formation rate surface density, and the maximum mass of young star clusters in a disk galaxy. II. The grand-design galaxy M51

Local Group timing in Milgromian dynamics. A past Milky Way-Andromeda encounter at z>0.8

R144 : a very massive binary likely ejected from R136 through a binary-binary encounter

Recent advances on IMF research

The dynamical phase transitions of stellar systems and the corresponding kinetimacs

The failures of the standard model of cosmology require a new paradigm

The galactocentric radius dependent upper mass limit of young star clusters: stochastic star formation ruled out

The galaxy-wide IMF of dwarf late-type to massive early-type galaxies

The mmax-Mecl relation, the IMF and IGIMF: probabilistically sampled functions?

The rotationally stabilized VPOS and predicted proper motions of the Milky Way satellite galaxies

The specific frequency and the globular cluster formation efficiency in Milgromian dynamics