We conducted infrared spectroscopic observations of bright stars in the direction of the molecular clouds W33 and GMC G23.3 − 0.3. We compared stellar spectro-photometric distances with parallactic distances to these regions, and we were able to assess the association of the detected massive stars with these molecular complexes. The spatial and temporal distributions of the detected stars enabled us to locate sources of ionizing radiation and to gather precise information on the star formation history of these clouds. The studied clouds present different distributions of massive stars.

Here, we report on the first results from the e-MERLIN Cyg OB2 Radio Survey (COBRaS), which is designed to exploit e-MERLIN’s enhanced capabilities to conduct deep-field mapping of the tremendously rich Cyg OB2 association. The project aims to deliver the most detailed radio census of the most massive OB association in the northern hemisphere. There exists considerable evidence for clumping in the winds of hot stars, which has hugely important consequences for mass-loss determinations. The amount of mass lost from a massive star is a crucial parameter required for stellar and cluster evolution models that are paramount to our understanding of the formation and evolution of massive star clusters. Presenting some of the first 20cm (L band) detections of massive stars in Cyg OB2, both accurate mass-loss determinations and constraints upon clumping within their winds are made. These data substantially increase the observational detections of the outer wind of massive stars and in combination with other observations at different wavelengths, COBRaS will greatly advance our knowledge of clumping as a function of radial distance around massive stars.

We investigate the competitive role of the different dissipative phenomena acting on the onset of star formation of gravitationally bound systems in an external environment. Ram pressure, Kelvin-Helmholtz and Rayleigh-Taylor instabilities, and tidal forces are accounted for separately in an analytical framework and compared in their role in influencing the star forming regions. We present an analytical criterion to elucidate the dependence of star formation in a spherical stellar system on its surrounding environment. We consider the different signatures of these phenomena in synthetically realized colour-magnitude diagrams (CMDs) of the orbiting system thus investigating the detectability limits of these different effects for future observational projects and their relevance. The developed theoretical framework has direct applications to the cases of massive star clusters, dwarf galaxies in galaxy clusters and dwarf galaxies orbiting our Milky Way system, as well as any primordial gas-rich cluster of stars orbiting within its host galaxy.

New CO J=1–0 observations with NANTEN and NANTEN2 reveal that extensive collisions between two molecular clouds at relative velocity of 15 km s−1 triggered the O star formation in the Galactic mini-starbursts NGC 6357 and NGC 6334. Correlated/anti-correlated gas distributions and intermediate velocity features between the two clouds lend support for the cloud-cloud collision scenario. The timescale of the collision and high-mass star formation is as short as less than 0.5 Myrs, suggesting rapid O star formation.

N131 is an infrared dust bubble residing in a molecular filament. We aim to study the formation and fragmentation of this bubble with multi-wavelength dust and gas observations. Towards the bubble N131, we analyzed archival multi-wavelength observations including 3.6, 4.5, 5.8, 8.0, 24, 70, 160, 250, 350, 500 μm, 1.1 mm, and 21 cm. In addition, we performed new observations of CO (2-1), CO (1-0), and 13CO (1-0) with the IRAM 30-m telescope. Multi-wavelength dust and gas observations reveal a ringlike shell with compact fragments, two filamentary structures, and a secondary bubble N131-A. The bubble N131 is a rare object with a large hole at 24 μm and 21 cm in the direction of its center. The dust and gas clumps are compact and might have been compressed at the inner edge of the ringlike shell, while they are extended and might be pre-existing at the outer edge. The column density, excitation temperature, and velocity show a potentially hierarchical distribution from the inner to outer edge of the ringlike shell. We also detected the front and back sides of the secondary bubble N131-A in the direction of its center. The derived Lyman-continuum ionizing photon flux within N131-A is equivalent to an O9.5 star. Based on the above, we suggest that the bubble N131 might be triggered by the strong stellar winds from a group of massive stars inside the bubble. We propose a scenario in which the bubble N131 forms from the disruption of a gas filament by expansion of H II region, strong stellar winds, and fragments under self-gravity.

Massive star clusters are of fundamental importance both observationally, since they are visible at such great distances, and theoretically, because of their influence on the large–scale ISM. Understanding stellar feedback is a prerequisite for making sense of their formation and early evolution, since feedback influences cluster structure, star formation efficiency, and sets the timescales on which clusters emerge from their parent clouds to become optically visible. I review the progress made in understanding these issues from a numerical perspective.

We investigate the formation of stellar clusters from a Galactic scale SPH simulation. The simulation traces star formation over a 5.6 Myr timescale, with local gravitational instabilities resulting in ~ 105 solar masses of star formation in the form of sink particles. We investigate the time evolution of the physical properties of the forming clusters including their half-mass radii, their energies and the depletion time of the gas. Star formation is driven by the large scale flows which compress the gas to higher densities where self gravity takes over and collapse occurs. We show that the more massive clusters (up to ~ 2 × 104 solar masses) gather their material from of order 10 pc due to these large scale motions associated with the spiral arm passage and shock. The bulk of the gas becomes gravitationally bound near 1-2 Myr before sink formation, and in the absence of feedback, significant accretion ongoing on longer timescales. We trace the hierarchical merging process of cluster formation which naturally results in age spreads of order the crossing time of the original region which provides the gas reservoir for the cluster.

Stellar feedback from massive stars can unbind and disperse large amount of molecular gas, affecting the star formation efficiency. Based on ALMA and VLT observations in the Antennae galaxies we study a massive (~ 107 M⊙) and young (~ 3 Myr) SSC, B1, associated with compact molecular and ionized emission, which suggests that it is embedded in its parent cloud. However, we found contradictories and puzzling results on the structure and dynamics of the matter around the cluster, indicating that SSC B1 is not embedded in its parent cloud after all. We propose that radiation pressure was highly enhanced at the early stages of the SSC formation, disrupting the parent cloud in < 3 Myr. We show evidences of outflowing gas from the parent cloud in the more extended CO gas. Higher angular resolution observations are needed to validate this interpretation and to understand the origin and fate of the component seen to be associated with SSC B1.

We discuss the mechanism of cluster formation in hierarchically collapsing molecular clouds. Recent evidence, both observational and numerical, suggests that molecular clouds (MCs) may be undergoing global, hierarchical gravitational collapse. The “hierarchical” regime consists of small-scale collapses within larger-scale ones. The latter implies that the star formation rate increases systematically during the early stages of evolution, and occurs via filamentary flows onto “hubs” of higher density, mass, and velocity dispersion, and culminates a few Myr after than the small-scale collapses have started to form stars. In turn, the small-scale collapses occur in clumps embedded in the filaments, and are themselves falling into the larger potential well of the still-ongoing large-scale collapse. The stars formed in the early, small-scale collapses share the infall motion of their parent clumps towards the larger potential trough, so that the filaments feed both gaseous and stellar material to the hubs. This leads to the presence of older stars in a region where new protostars are still forming, to a scale-free or fractal structure of the clusters, in which each unit is composed of smaller-scale ones, and to the eventual merging of the subunits, explaining the observed structural features of open clusters.

There is now evidence that some aspects of compact star cluster formation and destruction are quasi-universal in nature, and some aspects depend on environment. But what do we mean by these terms, environmental and universal? Is one the dominant influence? How can things be both universal and environmentally dependent? In this contribution we first provide a brief historical overview, then examine evidence for both universality and environmental dependences, and finish by examining a new approach that both demonstrates the degree to which cluster mass functions are universal (i.e., to a level of roughly 0.2 in the Log over three orders of magnitude when normalizing by the star formation rate), and enables a method for quantifying 2nd-order environmental effects.

RCW 38 is the youngest super star cluster in the Galaxy and is located at a distance of 1.7 kpc. Molecular observations revealed that the cluster is associated with two molecular clouds having velocity difference of 12 km s−1. We interpret that the two clouds are colliding with each other and the collision triggered the cluster formation. The natal molecular gas still survives within ~ 0.5 pc of the central O stars which have an age of 0.1 Myrs as inferred from the collision morphology. We suggest that the high column density of one of the clouds 1023 cm−2 enabled formation of ~ 20 O stars in the cluster center and discuss the implications on massive cluster formation.

Scaling relations for globular clusters (GC) differ from the scaling relations for pressure supported (elliptical) galaxies. In this contribution we discuss the relative importance of nature and nurture in the establishment of the scaling between cluster density (or radius), mass and Galactocentric distance for the Milky Way GCs. We show that energy diffusion by stellar encounters (i.e. two-body relaxation) is the dominant mechanism in shaping the bivariate dependence of density on mass and Galactocentric distance for GCs with masses ≲ 106M⊙, and it can not be excluded that GCs formed with similar scaling relations as the more massive ultra-compact dwarf galaxies (UCDs). To explore the initial properties that give rise to the distributions of these quantities, we developed a fast cluster evolution model (Evolve Me A Cluster of StarS, emacss) and use it in a hierarchical Bayesian framework to fit a parameterised model for the initial properties of Milky Way GCs to the observed present-day properties. The best-fit cluster initial mass function is substantially flatter (power-law with index − 0.6 ± 0.2) than what is observed for young massive clusters (YMCs) forming in the nearby Universe (power-law with index − 2). This result is driven by the metal-poor GCs, a slightly steeper CIMF is allowed when considering the metal-rich GCs separately (α ≃ −1.2 ± 0.4). If stellar mass loss and two-body relaxation in the Milky Way tidal field are the dominant disruption mechanisms, then GCs formed differently from YMCs.

Characterization of the binary fractions in star clusters is of fundamental importance for many fields in astrophysics. Observations indicate that the majority of stars are found in binary systems, while most stars with masses greater than 0.5M⊙ are formed in star clusters. In addition, since binaries are on average more massive than single stars, in resolved star clusters these systems are thought to be good tracers of (dynamical) mass segregation. Over time, dynamical evolution through two-body relaxation will cause the most massive objects to migrate to the cluster center, while the relatively lower-mass objects remain in or migrate to orbits at greater radii. This process will globally dominate a cluster’s stellar distribution. However, close encounters involving binary systems may disrupt ‘soft’ binaries. This process will occur more frequently in a cluster’s central, dense region than in its periphery, which may mask the effects of mass segregation. Using high resolution Hubble Space Telescope observations, combined with sophisticated N-body simulations, we investigate the radial distributions of the main-sequence binary fractions in massive young Large Magellanic Cloud star clusters. We show that binary disruption may play an important role on very short timescales, depending on the environmental conditions in the cluster cores. This may lead to radial binary fractions that initially decline in the cluster centers, which is contrary to the effects expected from dynamical mass segregation.

How starburst clusters form out of molecular clouds is still an open question. In this article, I highlight some of the key constraints in this regard, that one can get from the dynamical evolutionary properties of dense stellar systems. I particularly focus on secular expansion of massive star clusters and hierarchical merging of sub-clusters, and discuss their implications vis-á-vis the observed properties of young massive clusters. The analysis suggests that residual gas expulsion is necessary for shaping these clusters as we see them today, irrespective of their monolithic or hierarchical mode of formation.

Recent observations have revealed the existence of stellar mass black hole (BH) candidates in some globular clusters (GC) in the Milky Way and in other galaxies. Given that the detection of BHs is challenging, these detections likely indicate the existence of large populations of BHs in these clusters. This is in direct contrast to the past understanding that at most a handful of BHs may remain in old GCs due to quick mass segregation and rapid mutual dynamical ejection. Modern realistic star-by-star numerical simulations suggest that the retention fraction of BHs is typically much higher than what was previously thought. The BH dynamics near the cluster center leads to dynamical formation of new binaries and dynamical ejections, and acts as a persistent and significant energy source for these clusters. We have started exploring effects of BHs on the global evolution and survival of star clusters. We find that the evolution as well as survival of massive star clusters can critically depend on the details of the initial assumptions related to BH formation physics, such as natal kick distribution, and the initial stellar mass function (IMF). In this article we will present our latest results.

The study of intermediate-mass black holes (IMBHs) is a young and promising field of research. If IMBH exist, they could explain the rapid growth of supermassive black holes by acting as seeds in the early stage of galaxy formation. Formed by runaway collisions of massive stars in young and dense stellar clusters, intermediate-mass black holes could still be present in the centers of globular clusters, today. We measured the inner kinematic profiles with integral-field spectroscopy for 10 Galactic globular cluster and determined masses or upper limits of central black holes. In combination with literature data we further studied the positions of our results on known black-hole scaling relations (such as M• − σ) and found a similar but flatter correlation for IMBHs. Applying cluster evolution codes, the change in the slope could be explained with the stellar mass loss occurring in clusters in a tidal field over its life time. Furthermore, we present results from several numerical simulations on the topic of IMBHs and integral field units (IFUs). N-body simulations were used to simulate IFU data cubes. For the specific case of NGC 6388 we simulated two different IFU techniques and found that velocity dispersion measurements from individual velocities are strongly biased towards lower values due to blends of neighbouring stars and background light. In addition, we use the Astrophysical Multipurpose Software Environment (AMUSE) to combine gravitational physics, stellar evolution and hydrodynamics to simulate the accretion of stellar winds onto a black hole. We find that the S-stars need to provide very strong winds in order to explain the accretion rate in the galactic center.

Stellar clusters form with a high level of substructure, inherited from the molecular cloud and the star formation process. Evidence from observations and simulations also indicate the stars in such young clusters form a subvirial system. The subsequent dynamical evolution can cause important mass loss, ejecting a large part of the birth population in the field. It can also imprint the stellar population and still be inferred from observations of evolved clusters. Nbody simulations allow a better understanding of these early twists and turns, given realistic initial conditions. Nowadays, substructured, clumpy young clusters are usually obtained through pseudo-fractal growth and velocity inheritance. We introduce a new way to create clumpy initial conditions through a ”Hubble expansion” which naturally produces self consistent clumps, velocity-wise. In depth analysis of the resulting clumps shows consistency with hydrodynamical simulations of young star clusters. We use these initial conditions to investigate the dynamical evolution of young subvirial clusters. We find the collapse to be soft, with hierarchical merging leading to a high level of mass segregation. The subsequent evolution is less pronounced than the equilibrium achieved from a cold collapse formation scenario.

Young and massive star clusters above a critical mass form thermally unstable clumps reducing locally the temperature and pressure of the hot 107 K cluster wind. The matter reinserted by stars, and mass loaded in interactions with pristine gas and from evaporating circumstellar disks, accumulate on clumps that are ionized with photons produced by massive stars. We discuss if they may become self-shielded when they reach the central part of the cluster, or even before it, during their free fall to the cluster center. Here we explore the importance of heating efficiency of stellar winds.

UCDs are super massive star clusters found largely in dense regions but have also been found around individual galaxies and in smaller groups. Their origin is still under debate but consensus is that they formed either during major galaxy mergers as mergers of super massive star clusters, are simply the high mass end of the globular cluster luminosity function and formed in the same way as globular clusters, or that they formed from the threshing of galaxies and are remnant nuclear star clusters, which themselves may have formed from the mergers of globular star clusters within galaxies. We are attempting to disentangle these competing formation scenarios with a large survey of UCDs in the Coma cluster. Using ACS two-passband imaging from the HST/ACS Coma Cluster Treasury Survey, we are using colors and sizes to identify the UCD cluster members. With a large sample within the core region of the Coma cluster, we will use the population size, properties, and spatial distribution, and comparison with the Coma globular cluster and nuclear star cluster populations to discriminate between the threshing and globular cluster scenarios. In particular, previously we have found a possible correlation of UCD colors with host galaxy and a possible excess of UCDs around a non-central giant galaxy with an unusually large globular cluster population, both suggestive of a globular cluster origin. With a larger sample size, we are investigating whether the color correlation with host persists and whether the UCD population is consistent with, or in excess of, the bright end of the GCLF. We present initial results from the survey.

Studying the internal dynamics of stellar clusters is conducted primarily through N-Body simulations. One of the major inputs into N-Body simulations is the binary star frequency and mass distribution, which is currently constrained by relations derived from field binary stars. However to truly understand how clustered environments evolve, binary data from within star clusters is needed including masses. Detailed information on binaries masses, primary and secondary, in star clusters has been limited to date. The primary technique currently available has been radial velocity surveys that are limited in depth. Using previous two-band photometry-based studies that may cover different mass ranges produce potentially discrepant interpretations of the observed binary population. We introduce a new binary detection method, Binary INformation from Open Clusters Using SEDs (BINOCS) that covers the wide mass range needed to improve cluster N-body simulation inputs and comparisons. Using newly-observed multi-wavelength photometric catalogs (0.3 - 8 microns) of the key open clusters with a range of ages, we can show that the BINOCS method determines accurate binary component masses for unresolved cluster binaries through comparison to available RV-based studies. Using this method, we present results on the dynamical evolution of binaries from 0.4 - 2.5 solar masses within five prototypical clusters, spaning 30 Myr to 3.5 Gyr, and how the binary populations evolve as a function of mass.