Embedded Star Clusters Research Articles

Generation of inclined protoplanetary discs and misaligned planets through mass accretion – I. Coplanar secondary discs

We study the three-dimensional evolution of a viscous protoplanetary disc which accretes gas material from a second protoplanetary disc during a close encounter in an embedded star cluster. The aim is to investigate the capability of the mass accretion scenario to generate strongly inclined gaseous discs which could later form misaligned planets. We use smoothed particle hydrodynamics to study mass transfer and disc inclination for passing stars and circumstellar discs with different masses. We explore different orbital configurations to find the parameter space which allows significant disc inclination generation. \citet{Thi2011} suggested that significant disc inclination and disc or planetary system shrinkage can generally be produced by the accretion of external gas material with a different angular momentum. We found that this condition can be fullfilled for a large range of gas mass and angular momentum. For all encounters, mass accretion from the secondary disc increases with decreasing mass of the secondary proto-star. Thus, higher disc inclinations can be attained for lower secondary stellar masses. Variations of the secondary disc's orientation relative to the orbital plane can alter the disc evolution significantly. The results taken together show that mass accretion can change the three-dimensional disc orientation significantly resulting in strongly inclined discs. In combination with the gravitational interaction between the two star-disc systems, this scenario is relevant for explaining the formation of highly inclined discs which could later form misaligned planets.

Rotation in young massive star clusters

Hydrodynamical simulations of turbulent molecular clouds show that star clusters form from the hierarchical merger of several sub-clumps. We run smoothed-particle hydrodynamics simulations of turbulence-supported molecular clouds with mass ranging from 1700 to 43000 Msun. We study the kinematic evolution of the main cluster that forms in each cloud. We find that the parent gas acquires significant rotation, because of large-scale torques during the process of hierarchical assembly. The stellar component of the embedded star cluster inherits the rotation signature from the parent gas. Only star clusters with final mass < few X 100 Msun do not show any clear indication of rotation. Our simulated star clusters have high ellipticity (~0.4-0.5 at t=4 Myr) and are subvirial (Q_vir<~0.4). The signature of rotation is stronger than radial motions due to subvirial collapse. Our results suggest that rotation is common in embedded massive (>~1000 Msun) star clusters. This might provide a key observational test for the hierarchical assembly scenario.

Embedded star clusters as sources of high-energy cosmic rays

Massive stars are mainly found in stellar associations. These massive star clusters occur in the heart of giant molecular clouds. The strong stellar wind activity in these objects generates large bubbles and induces collective effects that could accelerate particles up to high energy and produce gamma rays. The best way to input an acceleration origin to the stellar wind interaction in massive stellar cluster is to observe young massive star clusters in which no supernova explosion has occurred yet. This work aims to constrain the part of stellar wind mechanical energy that is converted into energetic particles using the sensitivity of the ongoing Fermi/LAT instrument. This work further provides detailed predictions of expected gamma-ray fluxes in the view of the on-set of the next generation of imaging atmospheric Cherenkov telescopes. A one-zone model where energetic particles are accelerated by repeated interactions with strong supersonic shocks occurring in massive star clusters was developed. The particle escape from the star cluster and subsequent interaction with the surrounding dense material and magnetic fields of the HII region was computed. We applied this model to a selection of eight embedded star clusters constricted by existing observations. We evaluated the gamma-ray signal from each object, combining both leptonic and hadronic contributions. We searched for these emissions in the Fermi/LAT observations in the energy range from 3 to 300 GeV and compared them to the sensitivity of the Cherenkov Telescope Array. No significant gamma-ray emission from these star clusters has been found. Less than 10% of stellar wind luminosities are supplied to the relativistic particles. Some clusters even show acceleration efficiency of less than 1%. The CTA would be able to detect gamma-ray emission from several clusters in the case of an acceleration efficiency of close to 1%.

EMBEDDED CLUSTERS IN THE LARGE MAGELLANIC CLOUD USING THE VISTA MAGELLANIC CLOUDS SURVEY*

ABSTRACT We present initial results of the first large-scale survey of embedded star clusters in molecular clouds in the Large Magellanic Cloud (LMC) using near-infrared imaging from the Visible and Infrared Survey Telescope for Astronomy Magellanic Clouds Survey. We explored a ∼1.65 deg2 area of the LMC, which contains the well-known star-forming region 30 Doradus as well as ∼14% of the galaxy’s CO clouds, and identified 67 embedded cluster candidates, 45 of which are newly discovered as clusters. We have determined the sizes, luminosities, and masses for these embedded clusters, examined the star formation rates (SFRs) of their corresponding molecular clouds, and made a comparison between the LMC and the Milky Way. Our preliminary results indicate that embedded clusters in the LMC are generally larger, more luminous, and more massive than those in the local Milky Way. We also find that the surface densities of both embedded clusters and molecular clouds is ∼3 times higher than in our local environment, the embedded cluster mass surface density is ∼40 times higher, the SFR is ∼20 times higher, and the star formation efficiency is ∼10 times higher. Despite these differences, the SFRs of the LMC molecular clouds are consistent with the SFR scaling law presented in Lada et al. This consistency indicates that while the conditions of embedded cluster formation may vary between environments, the overall process within molecular clouds may be universal.

Observational constraints on star cluster formation theory

Stars form predominantly in groups usually denoted as clusters or associations. The observed stellar groups display a broad spectrum of masses, sizes and other properties, so it is often assumed that there is no underlying structure in this diversity. Here we show that the assumption of an unstructured multitude of cluster or association types might be misleading. Current data compilations of clusters show correlations between cluster mass, size, age, maximum stellar mass etc. In this first paper we take a closer look at the correlation of cluster mass and radius. We use literature data to explore relations in cluster and molecular core properties in the solar neighborhood. We show that for embedded clusters in the solar neighborhood there exists a clear correlation between cluster mass and half-mass radius of the form $M_c = C R_c^{\gamma}$ with gamma = 1.7 +/-0.2. This correlation holds for infra red K band data as well as X-ray sources and for clusters containing a hundred stars up to those consisting of a few tens of thousands of stars. The correlation is difficult to verify for clusters containing <30 stars due to low-number statistics. Dense clumps of gas are the progenitors of the embedded clusters. We find a similar slope for the mass-size relation of dense, massive clumps as for the embedded star clusters. This might point at a direct translation from gas to stellar mass: however, it is difficult to relate size measurements for clusters (stars) to those for gas profiles. Taking into account multiple paths for clump mass into cluster mass, we obtain an average star-formation efficiency of 18%{+9.3}{-5.7} for the embedded clusters in the solar neighborhood. The derived mass-radius relation gives constraints for the theory of clustered star formation. Analytical models and simulations of clustered star formation have to reproduce this relation in order to be realistic (abridged)

Fundamental Properties of a Large, Unbiased Sample of Massive, Young, Embedded Star Clusters in the Milky Way

AbstractIn this project, we focus on the analysis of infrared observations of the clumps defined with the Galactic Census of High- and Medium-mass Protostars (ChaMP) (Barneset al.2011). We derive line of sight infrared extinction values, star counts and protostar candidates around the molecular gas emission obtained with the Mopra telescope. Then, we examine the correlation between radio and infrared properties of the clumps. For this stage of the project, we use the Vela-Carina and 2MASS catalogs to obtain a preliminary understanding of the final results. For the later stages, we will extract infrared photometry from our deep AAT near-IR andSpitzer3.6 and 4.5 μm images. With the final deep photometry results, we will compile the properties of individual clusters.

The difficult early stages of embedded star clusters and the importance of the pre-gas expulsion virial ratio

We examine the effects of gas-expulsion on initially substructured distributions of stars. We perform N-body simulations of the evolution of these distributions in a static background potential to mimic the gas. We remove the static potential instantaneously to model gas-expulsion. We find that the exact dynamical state of the cluster plays a very strong role in affecting a cluster's survival, especially at early times: they may be entirely destroyed or only weakly affected. We show that knowing both detailed dynamics and relative star-gas distributions can provide a good estimate of the post-gas expulsion state of the cluster, but even knowing these is not an absolute way of determining the survival or otherwise of the cluster.

Highly efficient star formation in NGC 5253 possibly from stream-fed accretion.

Gas clouds in present-day galaxies are inefficient at forming stars. Low star-formation efficiency is a critical parameter in galaxy evolution: it is why stars are still forming nearly 14billion years after the Big Bang and why star clusters generally do not survive their births, instead dispersing to form galactic disks or bulges. Yet the existence of ancient massive bound star clusters (globular clusters) in the Milky Way suggests that efficiencies were higher when they formed ten billion years ago. A local dwarf galaxy, NGC5253, has a young star cluster that provides an example of highly efficient star formation. Here we report the detection of the J = 3→2 rotational transition of CO at the location of the massive cluster. The gas cloud is hot, dense, quiescent and extremely dusty. Its gas-to-dust ratio is lower than the Galactic value, which we attribute to dust enrichment by the embedded star cluster. Its star-formation efficiency exceeds 50 per cent, tenfold that of clouds in the Milky Way. We suggest that high efficiency results from the force-feeding of star formation by a streamer of gas falling into the galaxy.

Erratum: 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 Hα 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 mmaxMecl relationand with the integrated galactic stellar initial mass function (IGIMF) theory. The difference in conclusions is based on a high degree of degeneracy between expectations from random sampling and those from 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.

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.

Evolution of star clusters in a cosmological tidal field

We present a method to couple N-body star cluster simulations to a cosmological tidal field, using the Astrophysical Multipurpose Software Environment. We apply this method to star clusters embedded in the CosmoGrid dark matter-only LambdaCDM simulation. Our star clusters are born at z = 10 (corresponding to an age of the Universe of about 500 Myr) by selecting a dark matter particle and initializing a star cluster with 32,000 stars on its location. We then follow the dynamical evolution of the star cluster within the cosmological environment. We compare the evolution of star clusters in two Milky-Way size haloes with a different accretion history. The mass loss of the star clusters is continuous irrespective of the tidal history of the host halo, but major merger events tend to increase the rate of mass loss. From the selected two dark matter haloes, the halo that experienced the larger number of mergers tends to drive a smaller mass loss rate from the embedded star clusters, even though the final masses of both haloes are similar. We identify two families of star clusters: native clusters, which become part of the main halo before its final major merger event, and the immigrant clusters, which are accreted upon or after this event; native clusters tend to evaporate more quickly than immigrant clusters. Accounting for the evolution of the dark matter halo causes immigrant star clusters to retain more mass than when the z=0 tidal field is taken as a static potential. The reason for this is the weaker tidal field experienced by immigrant star clusters before merging with the larger dark matter halo.

Probing the super star cluster environment of NGC 1569 using FISICA

We present near-IR JH spectra of the central regions of the dwarf starburst galaxy NGC 1569 using the Florida Image Slicer for Infrared Cosmology and Astrophysics (FISICA). The dust-penetrating properties and available spectral features of the near-IR, combined with the integral field unit (IFU) capability to take spectra of a field, make FISICA an ideal tool for this work. We use the prominent [He I] (1.083\mu m) and Pa\beta (1.282 \mu m) lines to probe the dense star forming regions as well as characterize the general star forming environment around the super star clusters (SSCs) in NGC 1569. We find [He I] coincident with CO clouds to the north and west of the SSCs, which provides the first, conclusive evidence for embedded star clusters here.

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

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

The embedded cluster or association Trumpler 37 in IC 1396: a search for evolutionary constraints

It is currently widely accepted that open star clusters and stellar associations result from the evolution of embedded star clusters. Parameters such star formation efficiency, time-scale of gas removal and velocity dispersion can be determinants of their future as bound or unbound systems. Finding objects at an intermediate evolution state can provide constraints to model the embedded cluster evolution. In the HII region IC1396, Trumpler 37 is an extended young cluster that presents characteristics of an association. We employed the Two Micron All Sky Survey (2MASS) photometry to analysing its structure and stellar content, and determining its astrophysical parameters. We also analysed 11 bright-rimmed clouds in IC1396 in order to search for young infrared star clusters, and the background open star cluster Teutsch 74, to verify whether it has any contribution to the observed stellar density profile of Trumpler 37. The derived parameters and comparison with template objects from other studies lead us to conclude that Trumpler 37, rather than as a star cluster, will probably emerge from its molecular cloud as an OB association.

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.

Popping star clusters as building blocks of the Milky Way’s thick disc

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

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.

The role of stellar collisions for the formation of massive stars

We use direct N-body simulations of gas embedded star clusters to study the importance of stellar collisions for the formation and mass accretion history of high-mass stars. Our clusters start in virial equilibrium as a mix of gas and proto-stars. Proto-stars then accrete matter using different mass accretion rates and the amount of gas is reduced in the same way as the mass of stars increases. During the simulations we check for stellar collisions and we investigate the role of these collisions for the build-up of high-mass stars and the formation of runaway stars. We find that a significant number of collisions only occur in clusters with initial half-mass radii r_h < 0.1 pc. After emerging from their parental gas clouds, such clusters end up too compact compared to observed young, massive open clusters. In addition, collisions lead mainly to the formation of a single runaway star instead of the formation of many high mass stars with a broad mass spectrum. We therefore conclude that massive stars form mainly by gas accretion, with stellar collisions only playing a minor role if any at all. Collisions of stars in the pre-main sequence phase might however contribute to the formation of the most massive stars in the densest star clusters and possibly to the formation of intermediate-mass black holes with masses up to a few 100 Msun.

An Investigation of the Loss of Planet-Forming Potential in Intermediate-Sized Young Embedded Star Clusters

ABSTRACTA large fraction of stars forming in our galaxy are born within clusters embedded in giant molecular clouds. In these environments, the background UV radiation fields impinging upon circumstellar disks can often dominate over the radiation fields produced by each disk’s central star. As a result, this background radiation can drive the evaporation of circumstellar disks and lead to the loss of planet-forming potential within a cluster. This article presents a detailed analysis of this process for clusters whose stellar membership falls within the range of 100 ≤ N ≤ 1000. For these intermediate-sized clusters, the background UV field is often dominated by the most massive stellar member. Due to the steep slope of the initial mass function, the amount of background UV light that bathes clusters of similar size displays significant variance. As a result, we perform a statistical analysis of this problem by calculating distributions of FUV flux values impinging upon star/disk systems for several cluster scenarios. We find that in the absence of dust attenuation, giant planet formation would likely be inhibited in approximately half of the systems forming within intermediate-sized clusters, regardless of stellar membership. In contrast, the presence of dust can significantly lower this value, with the effect considerably more pronounced in more populated clusters.

INDICATORS FOR CLUSTER SURVIVABILITY IN A DISPERSING CLOUD

We use N-body simulations to survey the response of embedded star clusters to the dispersal of their parent molecular cloud. The final stages of the clusters can be divided into three classes: the cluster (i) is destroyed, (ii) has a loose structure, and (iii) has a compact core. We are interested in three of the governing parameters of the parent cloud: (i) the mass, (ii) the size, and (iii) the dispersing rate. It is known that the final stage of the cluster is well correlated with the star formation efficiency (SFE) for systems with the same cluster and cloud profile. We deem that the SFE alone is not enough to address systems with clouds of different sizes. Our result shows that the initial cluster-cloud mass ratio at a certain Lagrangian radius, and the initial kinetic energy are better indicators for the survivability of embedded clusters.