We present $GalevNB$ (Galev for $N$-body simulations), an utility that converts fundamental stellar properties of $N$-body simulations into observational properties using the $GALEV$ (GAlaxy EVolutionary synthesis models) package, and thus allowing direct comparisons between observations and $N$-body simulations. It works by converting fundamental stellar properties, such as stellar mass, temperature, luminosity and metallicity into observational magnitudes for a variety of filters of mainstream instruments/telescopes, such as HST, ESO, SDSS, 2MASS, etc.), and into spectra that spans from far-UV (90 $\rm \AA$) to near-IR (160 $\rm \mu$m). As an application, we use $GalevNB$ to investigate the secular evolution of spectral energy distribution (SED) and color-magnitude diagram (CMD) of a simulated star cluster over a few hundred million years. With the results given by $GalevNB$ we discover an UV-excess in the SED of the cluster over the whole simulation time. We also identify four candidates that contribute to the FUV peak, core helium burning stars, thermal pulsing asymptotic giant branch (TPAGB) stars, white dwarfs and naked helium stars.

We report on GalevNB (Galev for N-Body simulations), an integrated software solution that provides N-body users direct access to the software package GALEV (GALaxy EVolutionary synthesis models). GalevNB is developed for the purpose of a direct comparison between N-body simulations and observations. It converts the fundamental stellar properties of N-body simulations, i.e., stellar mass, temperature, stellar luminosity and metallicity, into observational magnitudes for a variety of filters of widely used instruments/telescopes (HST, ESO, SDSS, 2MASS), and into spectra that span from far-UV (90 $\rm \AA$) to near-IR (160 $\rm \mu$m).

Tidal interactions of galaxies in galaxy clusters have been proposed as one potential explanation of the morphology-density relation at low masses. The efficiency of tidal transformation is expected to depend strongly on the orbit of a galaxy within the cluster halo. The orbit determines both the strength of the cluster's global tidal field and the probability of encounters with other cluster members. Additionally, we investigate the effect of an inclination between the disk of the infalling galaxy and its orbital plane. We also compare our results to observational data. We find that galaxies that entered a cluster from the outskirts are unlikely to be significantly transformed. Closer to the cluster centre, tidal interactions are a more efficient mechanism for producing harassed galaxies. But the inclination of the disk can reduce the mass loss significantly, yet it amplifies the thickening.

We investigate the circumstellar disc fraction as determined from L-band excess observations of the young, massive Arches and Quintuplet clusters residing in the central molecular zone of the Milky Way. The Quintuplet cluster was searched for L-band excess sources for the first time. We find a total of 26 excess sources in the Quintuplet cluster and 21 in the Arches cluster, of which 13 are new detections. With the aid of proper motion membership samples, the disc fraction of the Quintuplet cluster was derived for the first time to be 4.0 +/- 0.7%. There is no evidence for a radially varying disc fraction in this cluster. In the case of the Arches cluster, a disc fraction of 9.2 +/- 1.2% approximately out to the cluster's predicted tidal radius, r < 1.5 pc, is observed. This excess fraction is consistent with our previously found disc fraction in the cluster in the radial range 0.3 < r < 0.8 pc. In both clusters, the host star mass range covers late A- to early B-type stars, 2 < M < 15 Msun, as derived from J-band photospheric magnitudes. We discuss the unexpected finding of dusty circumstellar discs in these UV intense environments in the context of primordial disc survival and formation scenarios of secondary discs. We consider the possibility that the L-band excess sources in the Arches and Quintuplet clusters could be the high-mass counterparts to T Tauri pre-transitional discs. As such a scenario requires a long pre-transitional disc lifetime in a UV intense environment, we suggest that mass transfer discs in binary systems are a likely formation mechanism for the B-star discs observed in these starburst clusters.

Context: Using the recent observational data of R\"oser et al. we present $N$-body simulations of the Hyades open cluster. Aims: We make an attempt to determine initial conditions of the Hyades cluster at the time of its formation in order to reproduce the present-day cumulative mass profile, stellar mass and luminosity function (LF). Methods: We performed direct $N$-body simulations of the Hyades in an analytic Milky Way potential that account for stellar evolution and include primordial binaries in a few models. Furthermore, we applied a Kroupa (2001) IMF and used extensive ensemble-averaging. Results: We find that evolved single-star King initial models with King parameters $W_0 = 6-9$ and initial particle numbers $N_0 = 3000$ provide good fits to the observational present-day cumulative mass profile within the Jacobi radius. The best-fit King model has an initial mass of $1721\ M_\odot$ and an average mass loss rate of $-2.2 \ M_\odot/\mathrm{Myr}$. The K-band LFs of models and observations show a reasonable agreement. Mass segregation is detected in both observations and models. If 33% primordial binaries are included the initial particle number is reduced by 5% as compared to the model without primordial binaries. Conclusions: The present-day properties of the Hyades can be well reproduced by a standard King or Plummer initial model when choosing appropriate initial conditions. The degeneracy of good-fitting models can be quite high due to the large dimension of the parameter space. More simulations with different Roche-lobe filling factors and primordial binary fractions are required to explore this degeneracy in more detail.

We perform coagulation & fragmentation simulations to understand grain growth in T Tauri & brown dwarf discs. We present a physically-motivated approach using a probability distribution function for the collision velocities and separating the deterministic & stochastic velocities. We find growth to larger sizes compared to other models. Furthermore, if brown dwarf discs are scaled-down versions of T Tauri discs (in terms of stellar & disc mass, and disc radius), growth at the same location with respect to the outer edge occurs to similar sizes in both discs.

The Hyades cluster is an ideal target to study the dynamical evolution of a star cluster over the entire mass range due to its intermediate age and proximity to the Sun. We wanted to extend the Hyades mass function towards lower masses down to 0.1 Msol and to use the full three-dimensional spatial information to characterize the dynamical evolution of the cluster. We performed a kinematic and photometric selection using the PPMXL and Pan-STARRS1 sky surveys, to search for cluster members up to 30 pc from the cluster centre. We determined our detection efficiency and field star contamination rate to derive the cluster luminosity and mass functions down to masses of 0.1 Msol. The thorough astrometric and photometric constraints minimized the contamination. A minimum spanning tree algorithm was used to quantify the mass segregation. We discovered 43 new Hyades member candidates with velocity perpendicular to the Hyades motion up to 2 km/s. They have mass estimates between 0.43 and 0.09 Msol, for a total mass of 10 Msol. This doubles the number of Hyades candidates with masses smaller than 0.15 Msol. We provide an additional list of 11 possible candidates with velocity perpendicular to the Hyades motion up to 4 km/s. The cluster is significantly mass segregated. The extension of the mass function towards lower masses provided an even clearer signature than estimated in the past. We also identified as likely Hyades member an L0 dwarf previously assumed to be a field dwarf. Finally we question the membership of a number of previously published candidates, including a L2.5-type dwarf.

We perform coagulation and fragmentation simulations using the new physically-motivated model by Garaud et al. (2013) to determine growth locally in brown dwarf disks. We show that large grains can grow and that if brown dwarf disks are scaled down versions of T Tauri disks (in terms of stellar mass, disk mass and disk radius) growth at an equivalent location with respect to the disk truncation radius can occur to the same size in both disks. We show that similar growth occurs because the collisional timescales in the two disks are comparable. Our model may therefore potentially explain the recent observations of grain growth to millimetre sizes in brown dwarf disks, as seen in T Tauri disks.

Planet formation occurs within the gas and dust rich environments of protoplanetary disks. Observations of these objects show that the growth of primordial sub micron sized particles into larger aggregates occurs at the earliest stages of the disks. However, theoretical models of particle growth that use the Smoluchowski equation to describe collisional coagulation and fragmentation have so far failed to produce large particles while maintaining a significant populations of small grains. This has been generally attributed to the existence of two barriers impeding growth due to bouncing and fragmentation of colliding particles. In this paper, we demonstrate that the importance of these barriers has been artificially inflated through the use of simplified models that do not take into account the stochastic nature of the particle motions within the gas disk. We present a new approach in which the relative velocities between two particles is described by a probability distribution function that models both deterministic motion and stochastic motion. Taking both into account can give quite different results to what has been considered recently in other studies. We demonstrate the vital effect of two "ingredients" for particle growth: the proper implementation of a velocity distribution function that overcomes the bouncing barrier and, in combination with mass transfer in high-mass-ratio collisions, boosts the growth of larger particles beyond the fragmentation barrier. A robust result of our simulations is the emergence of two particle populations (small and large), potentially explaining simultaneously a number of long-standing problems in protoplanetary disks, including planetesimal formation close to the central star, the presence of mm to cm size particles far out in the disk, and the persistence of micron-size grains for millions of years.

Most stars form in a cluster environment. These stars are initially surrounded by discs from which potentially planetary systems form. Of all cluster environments starburst clusters are probably the most hostile for planetary systems in our Galaxy. The intense stellar radiation and extreme density favour rapid destruction of circumstellar discs via photoevaporation and stellar encounters. Evolving a virialized model of the Arches cluster in the Galactic tidal field we investigate the effect of stellar encounters on circumstellar discs in a prototypical starburst cluster. Despite its proximity to the deep gravitational potential of the Galactic centre only a moderate fraction of members escapes to form an extended pair of tidal tails. Our simulations show that encounters destroy one third of the circumstellar discs in the cluster core within the first 2.5 Myr of evolution, preferentially affecting the least and most massive stars. A small fraction of these events causes rapid ejection and the formation of a weaker second pair of tidal tails that is overpopulated by disc-poor stars. Two predictions arise from our study: (i) If not destroyed by photoevaporation protoplanetary discs of massive late B- and early O-type stars represent the most likely hosts of planet formation in starburst clusters. (ii) Multi-epoch K- and L-band photometry of the Arches cluster would provide the kinematically selected membership sample required to detect the additional pair of disc-poor tidal tails.

The realization that most stars form in clusters, raises the question of whether star/planet formation are influenced by the cluster environment. The stellar density in the most prevalent clusters is the key factor here. Whether dominant modes of clustered star formation exist is a fundamental question. Using near-neighbour searches in young clusters Bressert et al. (2010) claim this not to be the case and conclude that star formation is continuous from isolated to densely clustered. We investigate under which conditions near-neighbour searches can distinguish between different modes of clustered star formation. Near-neighbour searches are performed for model star clusters investigating the influence of the combination of different cluster modes, observational biases, and types of diagnostic and find that the cluster density profile, the relative sample sizes, limitations in observations and the choice of diagnostic method decides whether modelled modes of clustered star formation are detected. For centrally concentrated density distributions spanning a wide density range (King profiles) separate cluster modes are only detectable if the mean density of the individual clusters differs by at least a factor of ~65. Introducing a central cut-off can lead to underestimating the mean density by more than a factor of ten. The environmental effect on star and planet formation is underestimated for half of the population in dense systems. A analysis of a sample of cluster environments involves effects of superposition that suppress characteristic features and promotes erroneous conclusions. While multiple peaks in the distribution of the local surface density imply the existence of different modes, the reverse conclusion is not possible. Equally, a smooth distribution is not a proof of continuous star formation, because such a shape can easily hide modes of clustered star formation (abridged)

Using HST/WFPC2 data, we have performed a comparative study of the Blue Straggler Star (BSS) populations in the central regions of the globular clusters AM 1, Eridanus, Palomar 3, and Palomar 4. Located at distances RGC > 50 kpc from the Galactic Centre, these are (together with Palomar 14 and NGC 2419) the most distant clusters in the Halo. We determine their colour-magnitude diagrams and centres of gravity. The four clusters turn out to have similar ages (10.5-11 Gyr), significantly smaller than those of the inner-Halo globulars, and similar metallicities. By exploiting wide field ground based data, we build the most extended radial density profiles from resolved star counts ever published for these systems. These are well reproduced by isotropic King models of relatively low concentration. BSSs appear to be significantly more centrally segregated than red giants in all globular clusters, in agreement with the estimated core and half-mass relaxation times which are smaller than the cluster ages. Assuming that this is a signature of mass segregation, we conclude that AM 1 and Eridanus are slightly dynamically more evolved than Pal 3 and Pal 4.

Investigations of stellar encounters in cluster environments have demonstrated their potential influence on the mass and angular momentum of protoplanetary discs around young stars. In this study it is investigated in how far the initial surface density in the disc surrounding a young star influences the outcome of an encounter. Based on a power-law ansatz for the surface density, $\Sigma(r) \propto r^{-p}$, a parameter study of star-disc encounters with different initial disc-mass distributions has been performed using N-body simulations. It is demonstrated that the shape of the disc-mass distribution has a significant impact on the quantity of the disc-mass and angular momentum losses in star-disc encounters. Most sensitive are the results where the outer parts of the disc are perturbed by high-mass stars. By contrast, disc-penetrating encounters lead more or less independently of the disc-mass distribution always to large losses. However, maximum losses are generally obtained for initially flat distributed disc material. Based on the parameter study a fit formula is derived, describing the relative mass and angular momentum loss dependent on the initial disc-mass distribution index p. Generally encounters lead to a steepening of the density profile of the disc. The resulting profiles can have a r^{-2}-dependence or even steeper independent of the initial distribution of the disc material. From observations the initial density distribution in discs remains unconstrained, so the here demonstrated strong dependence on the initial density distribution might require a revision of the effect of encounters in young stellar clusters. The steep surface density distributions induced by some encounters might be the prerequisite to form planetary systems similar to our own solar system.

The young star clusters we observe today are the building blocks of a new generation of stars and planets in our Galaxy and beyond. Despite their fundamental role we still lack knowledge about the conditions under which star clusters form and the impact of these often harsh environments on the evolution of their stellar and substellar members. We demonstrate the vital role numerical simulations play to uncover both key issues. Using dynamical models of different star cluster environments we show the variety of effects stellar interactions potentially have. Moreover, our significantly improved measure of mass segregation reveals that it can occur rapidly even for star clusters without substructure. This finding is a critical step to resolve the controversial debate on mass segregation in young star clusters and provides strong constraints on their initial conditions.

Investigations of mass segregation are of vital interest for the understanding of the formation and dynamical evolution of stellar systems on a wide range of spatial scales. Our method is based on the minimum spanning tree (MST) that serves as a geometry-independent measure of concentration. Compared to previous such approaches we obtain a significant refinement by using the geometrical mean as an intermediate-pass. It allows the detection of mass segregation with much higher confidence and for much lower degrees of mass segregation than other approaches. The method shows in particular very clear signatures even when applied to small subsets of the entire population. We confirm with high significance strong mass segregation of the five most massive stars in the Orion Nebula Cluster (ONC). Our method is the most sensitive general measure of mass segregation so far and provides robust results for both data from simulations and observations. As such it is ideally suited for tracking mass segregation in young star clusters and to investigate the long standing paradigm of primordial mass segregation by comparison of simulations and observations.

Context: Field stars are not always single stars, but can often be found in bound double systems. Since binary frequencies in the birth places of stars, young embedded clusters, are sometimes even higher than on average the question arises of how binary stars form in young dense star clusters and how their properties evolve to those observed in the field population. Aims: We assess, the influence of stellar dynamical interactions on the primordial binary population in young dense cluster environments. Methods: We perform numerical N-body simulations of the Orion Nebula Cluster like star cluster models including primordial binary populations using the simulation code nbody6++. Results: We find two remarkable results that have yet not been reported: The first is that the evolution of the binary frequency in young dense star clusters is independent predictably of its initial value. The time evolution of the normalized number of binary systems has a fundamental shape. The second main result is that the mass of the primary star is of vital importance to the evolution of the binary. The more massive a primary star, the lower the probability that the binary is destroyed by gravitational interactions. This results in a higher binary frequency for stars more massive than $2\Msun$ compared to the binary frequency of lower mass stars. The observed increase in the binary frequency with primary mass is therefore most likely not due to differences in the formation process but can be entirely explained as a dynamical effect. Conclusions: Our results allow us to draw conclusions about the past and the future number of binary systems in young dense star clusters and demonstrate that the present field stellar population has been influenced significantly by its natal environments.

Stellar encounters potentially affect the evolution of the protoplanetary discs in the Orion Nebula Cluster (ONC). However, the role of encounters in other cluster environments is less known. We investigate the effect of the encounter-induced disc-mass loss in different cluster environments. Starting from an ONC-like cluster we vary the cluster size and density to determine the correlation of collision time scale and disc-mass loss. We use the NBODY6++ code to model the dynamics of these clusters and analyze the effect of star-disc encounters. We find that the disc-mass loss depends strongly on the cluster density but remains rather unaffected by the size of the stellar population. The essential outcome of the simulations are: i) Even in clusters four times sparser than the ONC the effect of encounters is still apparent. ii) The density of the ONC itself marks a threshold: in less dense and less massive clusters it is the massive stars that dominate the encounter-induced disc-mass loss whereas in denser and more massive clusters the low-mass stars play the major role for the disc mass removal. It seems that in the central regions of young dense star clusters -- the common sites of star formation -- stellar encounters do affect the evolution of the protoplanetary discs. With higher cluster density low-mass stars become more heavily involved in this process. This finding allows for the extrapolation towards extreme stellar systems: in case of the Arches cluster one would expect stellar encounters to destroy the discs of most of the low- and high-mass stars in several hundred thousand years, whereas intermediate mass stars are able to retain to some extant their discs even under these harsh environmental conditions.

The external destruction of protoplanetary discs in a clustered environment acts mainly due to two mechanisms: gravitational drag by stellar encounters and evaporation by strong stellar winds and radiation. If encounters play a role in disc destruction, one would expect that stars devoid of disc material would show unexpectedly high velocities as an outcome of close interactions. We want to quantify this effect by numerical simulations and compare it to observations. As a model cluster we chose the Orion Nebula Cluster (ONC). We found from the observational data that 8 to 18 stars leave the ONC with velocities several times the velocity dispersion. The majority of these high-velocity stars are young low-mass stars, among them several lacking infrared excess emission. Interestingly, the high-velocity stars are found only in two separate regions of the ONC. Our simulations give an explanation for the location of the high-velocity stars and provide evidence for a strong correlation between location and disc destruction. The high-velocity stars can be explained as the outcome of close three-body encounters; the partial lack of disc signatures is attributed to gravitational interaction. The spatial distribution of the high-velocity stars reflects the initial structure and dynamics of the ONC. Our approach can be generalized to study the evolution of other young dense star clusters, like the Arches cluster, back in time.

Observations indicate that in young stellar clusters the binary fraction for massive stars is higher than for solar mass stars. For the Orion Nebula Cluster (ONC) there is a binary frequency of ~ 50% for solar-mass stars compared to 70-100% for the massive O- and B-stars. We explore the reasons for this discrepancy and come up with two possible answers: a) a primordially higher binarity of massive stars could be inherent to the star formation process or b) the primordial binary rate might be the same for solar-mass and massive stars, but the higher capture cross section of the massive stars possibly leads to the formation of additional massive binaries in the early cluster development. Here we investigate the likelihood of the latter using the ONC as an example. N-body simulations are performed to track the capture events in an ONC-like cluster. We find that whereas low-mass stars rarely form bound systems through capture, the dynamics of the massive stars - especially in the first 0.5 Myrs - is dominated by a rapid succession of ``transient binary or multiple systems''. In observations the transient nature of these systems would not be apparent, so that they would be rated as binaries. At 1-2 Myrs, the supposed age of the ONC, the ``transient'' massive systems become increasingly stable, lasting on average several 10^6 yrs. Despite the ONC being so young, the observed binary frequency for massive stars -- unlike that of solar-mass stars -- is not identical to the primordial binary frequency but is increased by at least 10-15% through dynamical interaction processes. This value might be increased to at least 20-25% by taking disc effects into account. The primordial binary frequency could well be the same for massive and solar mass stars because the observed difference can be explained by capture processes alone.

AIM:The aim of this work is to understand to what extend gravitational interactions between the stars in high-density young stellar clusters, like the Orion Nebula Cluster (ONC), change the angular momentum in their protoplanetary discs. METHOD:Two types of simulations were combined -- N-body simulations of the dynamics of the stars in the ONC, and angular momentum loss results from simulations of star-disc encounters. RESULTS:It is shown that in a star-disc encounter the angular momentum loss is usually larger than the mass loss, so that the disc remnant has a lower specific angular momentum. Assuming an age of 1-2 Myr for the ONC, the disc angular momentum in the higher density region of the Trapezium is reduced by 15-20% on average. Encounters therefore play an important part in the angular momentum transport in these central regions but are not the dominant process. More importantly, even in the outer cluster regions the angular momentum loss is on average 3-5%. Here it is shown that a 3-5% loss in angular momentum might be enough to trigger gravitational instabilities even in low-mass discs - a possible prerequisite for the formation of planetary systems.

AIMS: The aim of this work is to understand whether there is a difference in the dispersion of discs around stars in high-density young stellar clusters like the Orion Nebula Cluster (ONC) according to the mass of the star. METHODS: Two types of simulations were combined -- N-body simulations of the dynamics of the stars in the ONC and mass loss results from simulations of star-disc encounters, where the disc mass loss of all stars is determined as a function of time. RESULTS: We find that in the Trapezium, the discs around high-mass stars are dispersed much more quickly and to a larger degree by their gravitational interaction than for intermediate-mass stars. This is consistent with the very recent observations of IC 348, where a higher disc frequency was found around solar mass stars than for more massive stars, suggesting that this might be a general trend in large young stellar clusters.

The relevance of encounters on the destruction of protoplanetary discs in the Orion Nebula Cluster (ONC) is investigated by combining two different types of numerical simulation. First, star-cluster simulations are performed to model the stellar dynamics of the ONC, the results of which are used to investigate the frequency of encounters, the mass ratio and separation of the stars involved, and the eccentricity of the encounter orbits. The results show that interactions that could influence the star-surrounding disc are more frequent than previously assumed in the core of the ONC, the so-called Trapezium cluster. Second, a parameter study of star-disc encounters is performed to determine the upper limits of the mass loss of the discs in encounters. For simulation times of $\sim$ 1-2 Myr (the likely age of the ONC) the results show that gravitational interaction might account for a significant disc mass loss in dense clusters. Disc destruction is dominated by encounters with high-mass stars, especially in the Trapezium cluster, where the fraction of discs destroyed due to stellar encounters can reach 10-15%. These estimates are in accord with observations of (Lada et al. 2000) who determined a stellar disc fraction of 80-85%. Thus, it is shown that in the ONC - a typical star-forming region - stellar encounters do have a significant effect on the mass of protoplanetary discs and thus affect the formation of planetary systems.

Interactions between disc-surrounded stars might play a vital role in the formation of planetary systems. Here a first parameter study of the effects of encounters on low-mass discs is presented. The dependence of the mass and angular momentum transport on the periastron distance, the relative mass of the encountering stars and eccentricity of the encounter is investigated in detail. This is done for prograde and retrograde coplanar encounters as well as non-coplanar encounters. For distant coplanar encounters our simulation results agree with the analytical approximation of the angular momentum loss by Ostriker(1994). However, for close or high-mass encounters, significant differences to this approximation are found. This is especially so in the case of retrograde encounters, where the analytical result predict no angular momentum loss regardless of the periastron distance whereas the simulations find up to ~ 20% loss for close encounters. For the non-coplanar case a more complex dependency on the inclination between orbital path and disc plane is found than for distant encounters. For the coplanar prograde case new fitting formulae for the mass and angular momentum loss are obtained, which cover the whole range from grazing to distant encounters. In addition, the final disc size and the mass exchange between discs is examined, demonstrating that for equal mass stars in encounters as close as 1.5 the disc radius, the disc size only is reduced by approximately 10%.