Ingo Thies

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.

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.

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, 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 on to 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.

Tidally induced brown dwarf and planet formation in circumstellar disks

Most stars are born in clusters and the resulting gravitational interactions between cluster members may significantly affect the evolution of circumstellar disks and therefore the formation of planets and brown dwarfs (BDs). Recent findings suggest that tidal perturbations of typical circumstellar disks due to close encounters may inhibit rather than trigger disk fragmentation and so would seem to rule out planet formation by external tidal stimuli. However, the disk models in these calculations were restricted to disk radii of 40 AU and disk masses below 0.1 M ☉. Here, we show that even modest encounters can trigger fragmentation around 100 AU in the sorts of massive (~0.5 M ☉), extended (≥100 AU) disks that are observed around young stars. Tidal perturbation alone can do this; no disk-disk collision is required. We also show that very low mass binary systems can form through the interaction of objects in the disk. In our computations, otherwise non-fragmenting massive disks, once perturbed, fragment into several objects between about 0.01 and 0.1 M ☉, i.e., over the whole BD 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 BDs and very massive planets which, interestingly, leads to a mass distribution consistent with the canonical substellar initial mass function. As a minor outcome, a possible explanation for the origin of misaligned extrasolar planetary systems is discussed.

A discontinuity in the low-mass IMF -the case of high multiplicity

ABSTRACT The empirical binary properties of brown dwarfs (BDs) differ from those of normalstars suggesting BDs form a separate population. Recent work by Thies and Krouparevealed a discontinuity of the initial mass function (IMF) in the very-low-mass starregime under the assumption of a low multiplicity of BDs of about 15 per cent. How-ever, previous observations had suggested that the multiplicity of BDs may be sig-nificantly higher, up to 45 per cent. This contribution investigates the implication ofa high BD multiplicity on the appearance of the IMF for the Orion Nebula Cluster,Taurus-Auriga, IC 348 and the Pleiades. We show that the discontinuity remains pro-nounced even if the observed MF appears to be continuous, even for a BD binaryfraction as high as 60%. We find no evidence for a variation of the BD IMF withstar-forming conditions. The BD IMF has a power-law index α BD ≈ +0.3 and about2 BDs form per 10 low-mass stars assuming equal-mass pairing of BDs.Key words: binaries: general — open clusters and associations: general — stars:low-mass, brown dwarfs — stars: luminosity function, mass function

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 GC 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, supermassive black holes 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 galaxys mass, half-mass radius and radial anisotropy. The fraction of eroded GCs is nearly 100 per cent in compact, M32-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.

CHARACTERIZING THE BROWN DWARF FORMATION CHANNELS FROM THE INITIAL MASS FUNCTION AND BINARY-STAR DYNAMICS

The stellar initial mass function (IMF) is a key property of stellar populations. There is growing evidence that the classical star-formation mechanism by the direct cloud fragmentation process has difficulties reproducing the observed abundance and binary properties of brown dwarfs and very-low-mass stars. In particular, recent analytical derivations of the stellar IMF exhibit a deficit of brown dwarfs compared to observational data. Here we derive the residual mass function of brown dwarfs as an empirical measure of the brown dwarf deficiency in recent star-formation models with respect to observations and show that it is compatible with the substellar part of the Thies-Kroupa IMF and the mass function obtained by numerical simulations. We conclude that the existing models may be further improved by including a substellar correction term that accounts for additional formation channels like disk or filament fragmentation. The term peripheral fragmentation is introduced here for such additional formation channels. In addition, we present an updated analytical model of stellar and substellar binarity. The resulting binary fraction and the dynamically evolved companion mass-ratio distribution are in good agreement with observational data on stellar and very-low-mass binaries in the Galactic field, in clusters, and in dynamically unprocessed groups of stars if all stars form as binaries with stellar companions. Cautionary notes are given on the proper analysis of mass functions and the companion mass-ratio distribution and the interpretation of the results. The existence of accretion disks around young brown dwarfs does not imply that these form just like stars in direct fragmentation.

Induced planet formation in stellar clusters: a parameter study of star–disc encounters

We present a parameter study of the possibility of tidally triggered disk instability. Using a restricted N-body model which allows for a survey of an extended parameter space, we show that a passing dwarf star with a mass between 0.1 and 1 M ⊙ can probably induce gravitational instabilities in the pre-planetary solar disk for prograde passages with minimum separations below 80–170 AU for isothermal or adiabatic disks. Inclined and retrograde encounters lead to similar results but require slightly closer passages. Such encounter distances are quite likely in young moderately massive star clusters (Scally & Clarke 2001; Bonnell et al. 2001). The induced gravitational instabilities may lead to enhanced planetesimal formation in the outer regions of the protoplanetary disk, and could therefore be relevant for the existence of Uranus and Neptune, whose formation timescale of about 100 Myr (Wuchterl et al. 2000) is inconsistent with the disk lifetimes of about a few Myr according to observational data by Haisch et al. (2001). The relatively small gas/solid ratio in Uranus and Neptune can be matched if the perturbing fly-by occurred after early gas depletion of the solar system, i.e. when the solar system was older than about 5 Myr. We also confirm earlier results by Heller (1993) that the observed 7 degree tilt of the solar equatorial plane relative to the ecliptic plane could be the consequence of such a close encounter.

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 “starlike” binaries with late-type primaries (≲ 1.5M⊙) 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.

MOND simulation suggests an origin for some peculiarities in the Local Group

The Milky Way (MW) and Andromeda (M31) galaxies possess rotating planes of satellites. Their formation has not been explained satisfactorily yet. It was suggested that the MW and M31 satellites are ancient tidal dwarf galaxies, which could explain their configuration. This suggestion gained support by an analytic backward calculation of the relative MW-M31 orbit in the MOND modified dynamics paradigm by Zhao et al. (2013) implying their close flyby 7-11 Gyr ago. Here we explore the Local Group history in MOND in more detail using a simplified first-ever self-consistent simulation. We note the features induced by the encounter in the simulation and identify their possible real counterparts. The simulation was set to approximately reproduce the observed MW and M31 masses, effective radii, separation, relative velocity and disk inclinations. We used the publicly available adaptive-mesh-refinement code Phantom of RAMSES. In the simulation, matter was transferred from the MW to M31 along a tidal tail. The encounter induced formation of several structures resembling the peculiarities of the Local Group. Most notably: 1) A rotating planar structure formed around M31 from the transferred material. It had a size similar to the observed satellite plane and was oriented edge-on to the simulated MW, just as the real one. 2) The same structure also resembled the tidal features observed around M31 by its size and morphology. 3) A warp in the MW developed with an amplitude and orientation similar to that observed. Moreover, a cloud of particles formed around the simulated MW, with the extent of the actual MW satellite system. The encounter did not end by merging in a Hubble time. The simulation thus demonstrated that MOND can possibly explain many peculiarities of the Local Group. Future more elaborate simulations should verify and expand our conclusions. (abridged)