© ESO, 2017

1. Introduction

Super star clusters (SSCs), potentially the progenitors of globular clusters (e.g., Portegies Zwart et al. 2010), are one of the most extreme forms of star formation. They are compact (a few parsec size) and massive (>10⁵M_⊙) star clusters. SSCs are found in galaxy interactions (e.g., Whitmore & Schweizer 1995) and starburst dwarf galaxies (e.g., Turner et al. 1998; Beck et al. 2002), as well as in the Milky Way (Clark et al. 2005). The number of such objects greatly increases in galaxy interactions and mergers, common phenomena in the Universe. Therefore, SSCs must be ubiquitous in the Universe.

Understanding the formation and early evolution of SSCs is observationally challenging, thus it is mainly based on theoretical work. Theoretical studies suggest that high external pressures (~ 10⁷k_B−10⁸k_B cm^-3 K) are needed to form such compact and massive clusters, conditions typically found in galaxy interactions (Elmegreen & Efremov 1997; Ashman & Zepf 2001). The absence of rotational shear in interacting and irregular galaxies may be an additional contribution to SSC formation, which, in spirals, may prevent giant molecular clouds (GMCs) from collapsing into dense clusters (Weidner et al. 2010). High external pressures can accelerate the collapse of molecular clouds, enhancing the star-formation efficiency (SFE, the fraction of gas which turns into stars). If the local (≲1 pc) SFE is high (≳30–50%), bound cluster formation will occur (Elmegreen & Efremov 1997; Ashman & Zepf 2001). Recent hydrodynamical models indicate that in dense star-forming regions, even if the global SFE in a GMC is below 30%, the local SFE can be higher than 50%, allowing the formation of bound clusters (Kruijssen et al. 2012). However, observational evidence constraining the SFE in young massive clusters is difficult to obtain. For instance, the SFE of a young (~1 Myr) SSC in the dwarf galaxy NGC 5253 is estimated to be 50% (Turner et al. 2015) or ~20% (Smith et al. 2016), depending on the assumed initial mass function.

High local SFE will occur if the local dynamical timescale of the parent molecular cloud (cluster formation timescale) is shorter than the timescale of disruptive processes (i.e., Portegies Zwart et al. 2010). The most invoked disruptive process is the stellar feedback, that is, the interaction of stars with the interstellar medium. The energy and momentum injected by the newly formed stars will unbind and disperse the matter around them. How quickly this occurs is key to the SFE. Studying feedback in SSCs is essential for understanding how such clusters form and evolve.

The specific role of each stellar feedback mechanism on impeding star formation and driving cloud disruption is a debated topic (see review paper Krumholz et al. 2014). Several stellar feedback mechanisms have been proposed and quantified analytically. Matzner & McKee (2000) discuss the ability of winds from low-mass stars to drive turbulence and balance dissipation. Matzner (2002) discusses the disruption of clouds by photoionization and expansion of the H ii region by the thermal gas pressure. More recent studies have included the impact of the radiation pressure on the expansion of the H ii region and the disruption of the parent molecular clouds (Murray et al. 2010; Fall et al. 2010). The mechanical energy associated with supernova explosions and winds from massive stars are also potentially able to disrupt clouds. The mechanical energy would be thermalized in shocks producing high pressure hot plasma, which can drive the expansion of the surrounding matter, provided it remains within the cloud. Observations and theoretical work indicate that this may not be as efficient as originally proposed (Castor et al. 1975), because the plasma is able to flow out through holes in the bubbles very early in the expansion (Harper-Clark & Murray 2009).

Analytical models and numerical simulations have attempted to explain how parent molecular clouds of massive (~10⁵M_⊙) star clusters are disrupted (e.g., Murray et al. 2010; Dale & Bonnell 2011; Dale et al. 2012; Skinner & Ostriker 2015). Analytical models suggest that radiation pressure is the dominant stellar feedback mechanism close to the clusters where the escape velocity is larger than the sound speed of the H ii gas (Krumholz & Matzner 2009; Murray et al. 2010). But, observational evidence of feedback in individual young (≲3 Myr) and massive (≳10⁵M_⊙) clusters is still scarce.

SSCs in nearby galaxies are ideal sites to investigate feedback mechanisms of massive clusters. Lopez et al. (2011) have used observations of the 30 Doradus star-forming region in the Large Magellanic Cloud, powered by one of the closest SSCs (≳5 × 10⁴M_⊙, Andersen et al. 2009), to compare the relative values of the pressures associated with radiation, ionized gas, and hot plasma, as a function of position in the nebula. They claim that the embedded massive cluster has broken apart its parent cloud within a time scale shorter than the main sequence life time of their most massive stars, thus before SN explosions become significant. Within the shells close to the core of the central R136 SSC, radiation pressure dominates. In this paper, we carry out a similar study in the Antennae galaxy merger for the most massive and youngest SSC across the overlap region.

The Antennae galaxy (NGC 4038/39) is a nearby (22 Mpc, Schweizer et al. 2008) and well studied merger between two gas-rich spiral galaxies. The region where the two spiral disks interact is known as the overlap region. The molecular gas content is principally constrained to both nuclei and to the overlap region, where it is observed to be fragmented into Super-Giant Molecular Complexes (SGMCs, Wilson et al. 2000). Interferometric observations showed that SGMCs are clumpy and present different velocity components (SMA observations by Ueda et al. 2012, ALMA observations by Herrera et al. 2012 and Whitmore et al. 2014). The Antennae galaxies host a large number of SSCs (Whitmore & Schweizer 1995; Mengel et al. 2005; Whitmore et al. 2010), which have masses of up to almost 10⁷M_⊙. We refer to Adamo & Bastian (2015) for a recent discussion of the fraction of star formation that happens in clusters. In this study we focus on the overlap region, where the most massive (≥10⁵M_⊙), and also the youngest (<10 Myr) SSCs are located (Mengel et al. 2005), which are spatially associated with SGMCs.

The paper is organized as follows. In Sect. 2, we introduce the datasets used in this paper. In Sect. 3, we present the search in our fields-of-view (FOVs) for SSCs associated with molecular and ionized gas on scales of their parent GMCs. We focus on a unique SSC with both molecular and ionized gas emission, suggesting that it may still be embedded in its parent cloud. This SSC, B1, is the most massive SSC in the overlap region. Section 4 presents the physical characteristics of SSC B1 and its parent cloud, and the measurement of their gas emission lines. In Sect. 5, we describe the physical structure of the matter surrounding SSC B1. In Sect. 6, we discuss the role of radiation pressure in disrupting the parent molecular cloud and argue that the bound natal cloud is already disrupted and outflowing gas from this cloud can still be observed. In Sect. 7, we discuss our interpretation of the data. Finally, Sect. 8 gives the conclusions of this paper.

2. Observations

Our investigation combines two datasets on the Antennae overlap region. We used observations obtained in February 2011 with the SINFONI imager spectrometer facility on the ESO Very Large Telescope, in the near-IR K-band which includes Brγ and several H₂ rovibrational lines. The spectral resolving power is R = 4000 at λ = 2.2 μm. The SINFONI FOV is 8′′ × 8′′ in size. In Herrera et al. (2011) and Herrera et al. (2012), we introduced our four SINFONI pointings across the overlap region, each of them covers a single SGMC. Figure 1 shows the K-band emission from the Antennae overlap region, where we mark our four SINFONI FOVs with dashed boxes. In this paper, we use the SINFONI observations of SGMC 4/5. The angular resolution is 0.̋7 × 0.̋6 (full width at half maximum, FWHM). We computed the velocity resolution at different wavelengths by fitting Gaussian curves to the sky-lines (see Table 3). Details of the data reduction can be found in Herrera et al. (2012).

We combined the previous dataset with archive observations done with the ALMA interferometer in the Cycle 0 early science process between July and November, 2012, which covered the entire Antennae overlap region (PI: B. Whitmore, Whitmore et al. 2014). The integration time on source was 3 h and the number of antennas varied between 14 and 24. These observations, done in Band 7 (345 GHz), include the CO(3−2) line emission as well as dust continuum emission. The data was retrieved from the ALMA science archive fully reduced by the ALMA team. We imaged the data using the CLEAN algorithm in CASA, by defining cleaning boxes enclosing the gas emission at each channel. The native spectral resolution of the observations corresponds to a channel width of ~0.5 MHz (~0.5 km s^-1). To increase the signal-to-noise ratio, we smoothed the channels to a velocity resolution of 10 km s^-1. The synthesized beam FWHM is 0.̋55 × 0.̋43, which corresponds to a lineal scale of 59 pc × 46 pc at the distance of the Antennae. We finally corrected the image by the primary beam response. We also used archive CO(2−1) (230 GHz) observations obtained with ALMA in the Band 6 during the science verification process. These observations were presented by Espada et al. (2012). Data reduction was performed by the ALMA team. We corrected the data cube by the primary beam response. The synthesized beam is 0.̋86 × 1.̋67 (PA = 77°), and the velocity resolution is 20 km s^-1.

3. Search for embedded clusters

In this section, we show how we selected a unique SSC across the Antennae overlap region to study the disruption of its parent cloud. We first introduce the SSCs population over the overlap region and then we explain how we used the SINFONI data, which cover four FOVs in the overlap region, to search for compact line emission associated with SSCs.

All clusters younger than ~10 Myr have significant ionizing flux. For SSCs embedded in their parent cloud, it is expected to detect Brγ line emission from its H ii region, which should be compact (~100 pc) and barely resolved by the arcsecond spatial resolution of near-IR ground-based observations. We also expect to detect H₂ line emission from photo dissociated regions (PDRs; e.g., Habart et al. 2011). For older SSCs, the Brγ emission would be more extended because the ionizing photons travel into the SGMCs and the more diffuse ISM, over distances larger than the size of their parent cloud. The detection of line emission associated with SSCs in the SINFONI data-cube of the Antennae overlap region is not straightforward because the SGMCs, where SSCs are found, show complex velocity structures and some of them are associated with several SSCs (Herrera et al. 2012). We face two main difficulties. The first one is to isolate the emission associated with the SSC from that of SGMCs, and the second one is to isolate different sources within SGMCs. We focus on SSCs, which are isolated. We selected a few SSCs previously studied by Gilbert & Graham (2007), based on single-slit spectroscopic Brγ observations. Gilbert & Graham (2007) estimated ages and masses from the observed Brγ equivalent width and the cluster stellar magnitude, employing the stellar population synthesis model Starburst99¹ and adopting a Kroupa initial mass function (IMF). These SSCs are young (≲6 Myr) and massive (>8 × 10⁵M_⊙), and lie within the four SINFONI fields presented in Herrera et al. (2012). These are the clusters D, D1, and D2 in SGMC 1, C in SGMC 2 and B1 in SGMC 4/5, which are highlighted in Fig. 1 and listed in Table 1. We also list in Table 1 the properties of these clusters estimated by Whitmore et al. (2010) using integrated photometry in the filters UBVIH_α complemented with population synthesis models. For each of these SSCs, we searched in the SINFONI data for compact (diameter ~1′′) gas emission associated with the SSCs, to identify which ones are still embedded in their parent clouds.

Fig. 1 K-band emission from the overlap region obtained with the CFHT (Herrera et al. 2011) overlaid with Brγ contours from the SINFONI observations. Small black squares highlight the position of the SSCs lying within the four SINFONI FOVs, which are marked with dashed boxes. An inserted image of the Antennae galaxies shows, in a black rectangle, the relative position of the overlap region in the merger.

Open with DEXTER

For each SSC, we obtained the K-band spectrum by circular aperture photometry. The aperture diameter is ~1.̋3 (140 pc), twice the SINFONI seeing size. We subtracted the nearby background measured within an annulus of 0.̋65 width, the seeing of our SINFONI data. Within the SINFONI spectra, we measured the molecular and ionized emission line fluxes, listed in Table 1. Upper limits for the fluxes were estimated as the standard deviation of aperture photometry of several positions in the integrated line image. In three of the selected SSCs – D, D2 and B1 – we detect Brγ. These clusters have ages ≲5 Myr and masses ≳10⁶M_⊙. The young and massive SSC C in SGMC 2, has an ionized nebula larger than the aperture size used to search for compact emission (Herrera et al. 2011). Other clusters may not be detected in the line emission due to the difficulty of isolating their emission from the background. This is a problem for the less massive clusters. Gilbert & Graham (2007) detected Brγ emission for all clusters using an aperture size of 2′′ (~ 220 pc), which includes not only compact H ii regions but also extended emission. In only one cluster, SSC B1, we do detect both H₂ and Brγ compact emission, suggesting that this cluster may still be embedded in its parent molecular cloud.

Figures 2 and 3 show the ratio between the line fluxes (both H₂ 1−0 S(1) and Brγ) and cluster masses versus cluster ages for six and five sources, respectively. The figures include all of the SSCs extracted from Gilbert & Graham (2007), Fig. 2 also includes the molecular compact source PCC 1 presented in Herrera et al. (2011). The H₂-to-stellar mass and Brγ-to-stellar mass ratio are normalized to the value for B1. Figure 2 shows that the H₂ emission per unit of stellar mass decreases sharply, by more than one order of magnitude for ages of approximately 3.5 Myr. Only PCC 1 and B1 show clear H₂ emission. Sources without an obvious H₂ detection are older than B1. This suggests that the disruption of the parent clouds by the newly formed stars occurs within a few Myr (<5 Myr). Figure 3 shows a large amount of scatter. The Brγ emission does not decrease as fast as the H₂ emission. The comparison between D and D1 in SGMC 1 gives us an idea of the time-scale for the disruption of H ii regions in the overlap region. Both clusters have similar masses (~10⁶M_⊙) but only D, which is 2 Myr younger, is associated with compact ionized gas emission. SSC D1, ~6 Myr old and with a mass of a few times 10⁵M_⊙, has already cleared out all or most of its parent cloud. Our findings agree with the upper limit for the parent cloud disruption timescale of <5 Myr estimated for SSCs in a sample of nearby galaxies and other regions of the Antennae pair (Bastian et al. 2014).

Fig. 2 H2 1−0 S(1) emission over the star-cluster mass, for SSCs with different ages and PCC 1 (Herrera et al. 2011), observed in the SINFONI fields. Since both the masses and H2 1−0 S(1) fluxes depend on the extinction in the same way, the H2-to-stellar mass ratio is independent of extinction. PCC 1 has a lower limit since it has not yet formed massive stars. The H2-to-stellar mass ratio is normalized to the value for B1.

Open with DEXTER

	Fig. 3 Brγ emission over the star cluster mass, for the SSCs in the SINFONI fields. As for H2 (see Fig. 2), the Brγ-to-stellar mass ratio is extinction independent. The ratio is normalized to the value for B1.
Open with DEXTER

In the rest of the paper, we focus on SSC B1 (also identified as “WS80” in other references), the only SSC across the overlap region associated with bright compact ionized and molecular gas emission. This cluster coincides with a bright emission peak in the mid-IR discovered by ISO, accounting for 15% of the mid-IR flux from the whole Antennae (Mirabel et al. 1998). SSC B1 is the brightest compact radio source at 4 and 6 cm in the overlap region (Neff & Ulvestad 2000). It is also one of the brightest sources in the K-band and Brγ emission, in the near-IR observations presented by Gilbert et al. (2000) and our SINFONI data. Gilbert et al. (2000) infers that the nebular part of this cluster is a compact, dense H ii region (n = 10⁴ cm^-3). Whitmore & Zhang (2002) find that SSC B1 is the intrinsically brightest cluster in the Antennae merger. The ALMA maps published by Herrera et al. (2012, see their Fig. 1) and by Whitmore et al. (2014, their Fig. 10) show that it coincides with a CO emission peak. SSC B1 is the more massive and youngest stellar cluster within our SINFONI fields (Table 1), which are the two fundamental properties to study embedded SSCs without confusion.

4. Physical properties of the SSC and GMC

SSC B1 in SGMC 4/5 is associated with compact ionized and molecular gas emission. This suggests that B1 may still be embedded in its parent molecular cloud. It is the target that we use to study how massive SSCs disperse their parent molecular clouds. In this section, we use the SINFONI and ALMA observations to characterize the physical properties of the cluster and its parent cloud.

4.1. SSC B1

We quantify the number of ionizing photons N_Lyc from the cluster. The radio flux of SSC B1 at 4 and 6 cm has a thermal spectral index (Neff & Ulvestad 2000). Both thermal radio emission and hydrogen recombination lines originate within ionized regions and are proportional to N_Lyc, but the radio emission, unlike the Brγ emission, is not affected by dust extinction. From the thermal radio emission, we estimate N_Lyc by using Eq. (1) in Roman-Lopes & Abraham (2004). The radio fluxes at 4 and 6 cm are listed in Tables 4 and 3 in Neff & Ulvestad (2000, brightest source in their region #2), respectively, which gives a mean value of N_Lyc = 2.2 ± 0.1 × 10⁵³ phot s^-1. Gilbert & Graham (2007) report a value ~1.4 lower, estimated from the extinction corrected Brγ flux. Uncertainties in the estimation of N_Lyc come from the unknown fraction of ionizing photons that do not escape the H ii region but are absorbed by dust. Such a fraction is difficult to quantify. Observations of Galactic H ii regions yield a photoionizing fraction between 0.3 and 1.0 (see Fig. 10 in Draine 2011).

The age and stellar mass of SSC B1 are hard to constrain and discrepancies are found when using different methods. From the K-band spectrum of SSC B1 (Fig. 5), we can set an upper limit to the age of ~6 Myr due to the absence of CO band-heads at 2.3 μm. This age limit agrees with the ages estimated by Whitmore & Zhang (2002) and Whitmore et al. (2010) of 2 Myr and 1 Myr, respectively, comparing UBVIHα photometry with population synthesis models. It also agrees with the cluster age derived by Gilbert & Graham (2007) of 3.5 Myr from the equivalent width (EW) of the Brγ line of 255 ± 8 Å, and our own estimation of the Brγ EW from our SINFONI observations of 295 ± 24 Å (3.5 ± 0.7 Myr). The error bar in our estimation of the age is a statistical value from the EW measurement. Systematical uncertainties from the unknown fraction of ionizing photons escaping the H ii region are larger. The stellar mass is estimated to be 6.6 × 10⁶ and 6.8 × 10⁶M_⊙ by Whitmore & Zhang (2002) and Whitmore et al. (2010), respectively.

We compute the cluster bolometric luminosity using the Starburst99 models, for a given age and mass of the cluster. Using the age and mass estimated by Whitmore et al. (2010), we find L_cl = 5.3 × 10⁹L_⊙. Note that the bolometric luminosity has a non-linear dependence with age, and the highest value is observed between 2 and 4 Myr. For 3.5 Myr, the age estimated by the Brγ EW, the luminosity is ~40% higher. The mass, age, N_Lyc, and L_cl values of the B1 cluster are listed in Table 2.

In this paper, we choose to use the mass estimated by Whitmore et al. (2010), which is obtained at high angular resolution (HST observations), and thus better resolving the stellar emission out of the nebular emission. Additionally, unlike our estimation of the stellar mass from the N_Lyc value, it does not have uncertainties on the photoionization fraction that escapes the H ii region.

4.2. Spatial extent of the ionized and molecular gas associated with SSC B1

We used the Brγ line to measure the extent of the ionized gas and the continuum emission at 345 GHz for that of the molecular gas. This emission is optically thin and, unlike the near-IR H₂ line emission, is a tracer of the column density. The size of the ionized emission associated with the cluster was measured by fitting Gaussian curves to the Brγ emission profiles. The observed size is 0.̋85 × 1.̋02 (FWHM), which, after correction by the seeing, gives a geometric mean of 0.̋67 (FWHM = 70 pc, listed in Table 2). A 2D Gauss fit of the K-band continuum emission gives a similar observed size of 0.̋79 × 1.̋02 (FWHM), yielding a seeing-corrected geometric mean of 0.̋62 (66 pc).

The size of the molecular cloud was measured by fitting a 2D Gaussian function to the continuum emission. We first subtracted the emission seen to be associated with a water maser detected by Brogan et al. (2010) at the position α: 12^h01^m549959, δ: –18°53′05.̋543, ~0.̋7 away from SSC B1 (see Fig. 6, peak emission in the continuum at 345 GHz). The absence of a Brγ, K-band continuum or H₂ counterpart suggests that the water maser is an extremely young star forming region. We measured the size (FWHM) of the emission associated with the water maser to be 0.̋76 (~80 pc). We created a 2D Gaussian array, centered at the position of the maser, with a peak emission of that seen in the continuum emission at the position of the maser. A multiplicative factor of 0.5 is found to minimise the dispersion of the residuals after subtracting the water maser contribution and the 2D Gaussian fit of the GMC emission. Figure 4 displays the result of the water maser subtraction, that is, the continuum emission associated with SSC B1. The black dashed contours correspond to the 2D Gaussian fit to this emission, black solid contours represent the Brγ emission, and the white cross shows the position of the water maser. The peak of the continuum emission coincides with the position of the SSC B1, offset (0,0). The observed source size (FWHM) is 0.̋96 × 1.̋22, which gives a geometrical mean, after beam correction, of 0.̋98 (FWHM = 104 pc). This is the spatial extent of the molecular cloud seen to be associated with SSC B1. It is larger than the Brγ extent, which suggests that the ionizing emission comes from within the cloud. Gaussian fits to the H₂ 1−0 S(1) emission profile yield a size of 2.̋0 × 2.̋1 (FWHM), which, after correction by the seeing, gives a geometric mean of 2′′ (FWHM = 210 pc, listed in Table 2). This value is larger than the GMC size. This indicates that part of the H₂ emission comes from the SGMC and is not directly associated with SSC B1.

Fig. 4 Continuum emission at 345 GHz. Emission from H2O maser was subtracted. The white cross marks the position of the H2O maser. Solid black lines are the 50%, 70%, and 90% of the peak of the Brγ line emission (see Fig. 6). Dashed black contours correspond to the 2D Gaussian fit to the continuum emission. Offset positions are relative to α: 12h01m5495, δ: –18°53′05.̋6 J2000.0. The offset position (0, 0) corresponds to the position of the SSC B1, which coincides with the continuum peak emission.

Open with DEXTER

4.3. Intensities of the ionized and molecular gas associated with SSC B1

In this section, we describe the measurements of the CO and near-IR line fluxes (H₂, Brγ, He i, and Brδ), listed in Table 3.

Fig. 5 Upper panel: K-band spectrum of SSC B1, obtained from aperture photometry based on the Brγ line. Ionized gas lines are the brightest lines across the spectrum. We mark the position of the H2 rovibrational lines with short red lines. Lower panels: zoom in the H2 1−0 S(1) and Brγ lines. Red lines represent the fitted Gaussian curves to each line, and residuals are plotted in the bottom part. Table 3 lists the parameters of the Gaussian fits.

Open with DEXTER

The K-band spectrum of SSC B1 is presented in Fig. 5. Several ionized and molecular lines are detected; we highlight the main ones. This spectrum is obtained by integrating the flux within an aperture of twice the FWHM observed in Brγ (1.̋9 diameter), and subtracting a background estimated within an annulus of 0.̋65 width (the size of the seeing disk of the SINFONI data). The same figure shows a zoom on the Brγ and H₂ 1−0 S(1) lines. Brγ is one order of magnitude brighter than the H₂ 1−0 S(1) line. Since the H₂ spatial extent is larger than that of Brγ, the H₂ fluxes measured from this K-band spectrum represent an approximate estimate of the H₂ emission associated with SSC B1. For a larger aperture of 4.̋2 diameter, twice the FWHM of the H₂ emission profile, with no background subtraction, since it is negligible, we measured a H₂ 1−0 S(1) flux four times larger. The aperture size of 4.̋2 is likely to be dominated by the SGMC (conservative upper limit). Measurements of the line intensities were taken by fitting Gaussians to each emission line and are listed in Table 3. For most of the lines, we estimated the Gaussian fit from a velocity range limited to the central channels. Figures 5b and c show the Gaussian fit to the H₂ 1−0 S(1) and Brγ lines, as well as the residuals of these fittings. Fluxes in Table 3 are not corrected for extinction. Errors in fluxes correspond to the errors obtained from the Gaussian fit, to which we include, as standard measurement errors, the quadratic sum of the standard deviation per channel σ_chan (close to the line emission) and the 10% uncertainty of the SINFONI absolute flux. Velocities are given with respect to the LSR with error-bars estimated from the fitting of OH sky lines (~1.5–3 km s^-1). Table 3 also lists the observed line widths (FWHM) derived from the the Gaussian fits, the spectral resolution for each line, and the intrinsic line width (observed value deconvolved by the spectral resolution).

We estimated the extinction in the K-band by comparing the N_Lyc values obtained with the thermal radio emission (Sect. 4.1) and that obtained with the Brγ emission. We computed N_Lyc from Brγ assuming case B recombination, an electron density n_e = 10⁴ cm^-3, and an electron temperature of T_e = 10⁴ K (see Table 6 of Hummer & Storey 1987). We estimated an extinction for Brγ of 0.8 mag (a factor 2), listed in Table 2. Gilbert et al. (2000) reported a slightly higher value of 1.2 mag. Their Brγ flux, measured from single-slit spectroscopy, is 1.5 times smaller than ours. Whitmore et al. (2010) reported E(B−V) to be 2.44, and assuming a Draine extinction law, we obtained A_K = 0.9, which agrees with our estimated value.

Fig. 6 Intensity maps for SGMC 4/5. From left to right, K-band continuum emission contours overlaid on the H2 1−0 S(1) emission, CO(3−2) emission, Brγ emission, and continuum emission at 345 GHz. Flux units of the color bars are × 10-17 erg s-1 cm-2 for H2 and Brγ, Jy beam-1 for CO(3−2), and mJy beam-1 for the continuum emission. The 1σ emission level of the K-band continuum emission is 2.05 × 10-20 erg s-1 cm-2. The peak emission of the K-band contours marks the position of the SSC B1, which corresponds to the offset position (0, 0). Offset positions are as in Fig. 4.

Open with DEXTER

Fig. 7 CO(3−2) channel maps of SGMC 4/5. The position of SSC B1 is marked with a cross and corresponds to the (0,0) offset position. Top panel: we present channel maps from 1440 km s-1 to 1550 km s-1 by 10 km s-1. Contours start at 10σ with a spacing of 10σ, with σ = 0.39 Jy beam-1 km s-1 measured in line-free channels. Bottom panel: CO channel maps from 1550 km s-1 to 1670 km s-1 by 10 km s-1. Contours start at 5σ with a spacing of 5σ. Offset positions are as in Fig. 4.

Open with DEXTER

	Fig. 8 CO(3−2) velocity components in SGMC 4/5. Left panels: CO emission integrated between 1370 km s-1 and 1540 km s-1 (top) and 1550 km s-1 and 1660 km s-1 (bottom). Contours represent the 10%, 30%, 50%, 70%, and 90% of the peak emission. Right panels: first moment of both velocity components. Contours are the same as in the left panel. Offset positions are as in Fig. 4.
Open with DEXTER

	Fig. 9 CO(3−2) and CO(2−1) spectra associated with SSC B1. For the CO J = 3−2 transition, we present the spectra measured with two apertures, 2.̋9 and 4.̋2, to characterize the two velocity components observed in the channel maps in Fig. 7.
Open with DEXTER

In Fig. 6, we have overlaid the K-band continuum emission on the H₂ 1−0 S(1), CO(3−2), Brγ line emission, and continuum emission at 345 GHz. The peak emission of the K-band continuum and Brγ emission coincide; they define the position of the SSC. The peak in the dust continuum emission at 345 GHz coincides with a water maser discovered by Brogan et al. (2010). Figure 4 displays the continuum emission with the water maser subtracted, its emission peak coincides with the position of the SSC. The molecular gas, that is, H₂ and CO emission, is more extended than the ionized gas. Neither H₂ or CO peaks exactly at the position of the SSC. Figure 7 presents the CO(3−2) channel maps from 1440 to 1670 km s^-1. From these channel maps, we can distinguish two spatially and spectrally separated components. Figure 8 presents the spatial extent (zero moment map) and first moment map of these components. The left top and bottom panels show the CO(3−2) emission integrated over the velocity ranges 1370–1540 km s^-1 and 1550–1660 km s^-1, respectively. There is a spatial offset between the two velocity components. The low-velocity emission is mainly associated with the SSC, while the higher velocity component seems to surround the cluster in a bubble-like shape structure, having the spatial extent of the SGMC. This bubble-like shape was already highlighted by Whitmore et al. (2014, see their Fig. 4). The first moment maps show that the low-velocity component has a velocity gradient while the high-velocity component is almost at the same velocity.

We measure the CO(3−2) line flux of the low-velocity component within an aperture diameter of 2.̋9, twice the observed FWHM size defined in the CO(3−2) emission integrated between 1370 km s^-1 and 1540 km s^-1 (see top-left figure in Fig. 8), and centered on the peak emission of this component. The resulting spectrum is plotted with a black line in Fig. 9. It has a non-Gaussian profile. Line parameters were estimated by fitting 2 Gaussian curves, but only those of the low-velocity component are listed in Table 3. Figure 8 shows that the high-velocity component has a larger spatial extent. We measure the CO(3−2) line flux of the high-velocity component within an aperture of 4.̋2 diameter, estimated from the size of the CO(3−2) emission integrated in the high-velocity channels. This size coincides with that measured in the H₂ 1−0 S(1). The result is plotted with a blue line in Fig. 9. The parameters listed in Table 3 for the high-velocity CO(3−2) emission component are derived from the fit of this spectrum with two Gaussian curves. To quantify the CO(2−1) emission, we used an aperture diameter of 4.̋7 (red line in Fig. 9), size computed from the total CO(2−1) integrated emission. We fit the emission with two Gaussian curves. Table 3 lists their line parameters. We do not subtract the water maser contribution to the CO emission, as we did for the continuum emission (Sect. 4.2), because the velocity structure is very complex, making it difficult to separate the contribution of the water maser from the total CO emission. However, this shortcoming does not affect the conclusions of this paper. Systematic uncertainties in the CO-to-H₂ conversion factor are much larger than the potential contribution of the maser to the CO emission.

Figure 9 displays the CO(3−2) spectra from both apertures mentioned above as well as the CO(2−1) spectra. Both CO transitions show the two velocity components observed in the channel maps. From the Gaussian fit, we find that the high-velocity component is almost twice wider than the low-velocity component. In Table 4, we list the CO(3−2)-to-CO(2−1) line ratios for both components. The ratio was measured within an aperture of 4.̋7 defined by the CO(2−1) size and after imaging the CO(3−2) observations using the CO(2−1) synthesized beam as restored beam during the imaging process (CLEAN algorithm in CASA). The intensity ratio between the J = 3–2 and 2–1 CO transitions depends on the gas temperature and density as discussed by Schulz et al. (2007) in the modeling of their single dish observations of the Antennae. We do not see a significant difference in the excitation of the gas between the low- and high-velocity components.

4.4. Mass of the GMC

We estimate the mass of the GMC, associated with SSC B1, from the continuum emission at 345 GHz. We approximate the spectral dependence on the dust emission with a gray body. We assume that the dust properties in the Antennae overlap region are the same as those in the Galaxy. Scoville et al. (2014) discuss the use of the Galactic values from Planck observations (Planck Collaboration XXV 2011) to estimate the millimetric masses in local starburst galaxies. They also show that dust temperatures between 15 and 30 K are typically found in starbursts (see also Dunne et al. 2011), even in extreme cases such as Arp 220, where the dust temperature does not exceed 45 K. In the Antennae, metallicities of approximately the solar one are observed towards several SSCs (Bastian et al. 2009). We assume a typical dust temperature of 25 K and the emissivity per hydrogen measured in the Taurus molecular cloud by Planck Collaboration XXV (2011), which is similar to that obtained for starburst galaxies (Scoville et al. 2014). Using Eq. (8) in Herrera et al. (2013), we derive the mass of the giant molecular cloud of M_GMC = 3.3 ± 0.1 × 10⁷M_⊙, listed in Table 2. The error is estimated as the 1σ standard deviation of the continuum image measured with the same aperture over nearby regions free of sources. Uncertainty on the dust temperature implies an uncertainty of a factor 2 in the mass (if T_d is 15 K, the mass will be twice the estimated value, while for T_d = 45 K, the mass will be half the estimated value). The gas column density N_H = 4.5 ± 0.3 × 10²³ H cm^-2 is estimated assuming spherical symmetry.

The virial mass, assuming a spherical cloud with a density profile ∝r^-1, can be written as M_vir = 190 Δv²R (MacLaren et al. 1988). We estimate the virial mass of the low-velocity component to be ~ 4 × 10⁷M_⊙. We also estimate the molecular gas mass from the observed CO(3−2) emission. We assume a CO-to-H₂ conversion factor of 0.6 × 10²⁰ H₂ cm^-2 (K km s^-1)^-1, estimated by Zhu et al. (2003) for the overlap region at scales of ~1.5 kpc and comparable to the value estimated by Kamenetzky et al. (2014) using CO modeling at 4.5 kpc scales (Herschel data) in the overlap region. This value is also comparable to the typical value found for starburst galaxies, ultraluminous infrared galaxies and towards the Galactic center (for a detailed discussion about the uncertainties on this factor see Bolatto et al. 2013). Other authors have estimated the conversion factor to be similar to the Galactic one (e.g., Wilson et al. 2003; Schulz et al. 2007). In this paper we choose to use a starburst-like value, which gives us a good agreement with our other measurements of the molecular gas (see below). We also assume a CO(3−2)-to-CO(1−0) ratio of ~0.6, as measured by Schulz et al. (2007) with single dish data in the overlap region and Ueda et al. (2012) with interferometric observations at 500 pc scales. Using Eq. (3) from Bolatto et al. (2013), we estimate M_mol to be 5.1 × 10⁷M_⊙. Using the Galactic CO-to-H₂ conversion factor of 2 × 10²⁰ H₂ cm^-2 (K km s^-1)^-1 (Solomon et al. 1987), the mass would be more than three times higher. The estimated M_mol value agrees with the value estimated from the continuum emission and the virial theorem, indicating that the CO-to-H₂ conversion factor estimated at 1.5 kpc scales by Zhu et al. (2003) is also adequate for SGMC 4/5. The measured flux from the high-velocity component yields a mass of 3.4 × 10⁷M_⊙. We list these values in Table 4.

5. Physical structure of the surrounding matter

In this section, we describe the physical structure of the molecular gas surrounding SSC B1 which we associate with the narrow velocity component of the CO spectrum in Fig. 9 and Table 3. We measure the molecular gas pressure and conclude that there is no trapping of IR light within the internal cavity (Sect. 5.1). In Sect. 5.2, we estimate the hot plasma pressure and find that the gas around SSC B1 is distributed in an inhomogeneous, already broken shell.

5.1. Gas pressure in molecular gas

Radiation pressure P_rad is the pressure exerted by the stellar radiation over the molecular gas that surrounds the central star cluster. We quantify the radiation pressure in the molecular gas surrounding SSC B1. We employ the expression from Murray et al. (2010), which includes the opacity of the shell, (1)where L_cl is the luminosity of the cluster, R_in is the H ii radius, and ⟨ τ_rad ⟩ accounts for the trapping of the IR photons in the cloud (⟨ τ_rad ⟩ ≥ 0). IR photons originate from UV and visible photons emitted by the cluster, which are absorbed by the dust grains in the cloud, and re-emitted in the IR. While the shell surrounding the cluster is closed and opaque, the IR photons are trapped within the shell, interacting more than once with the dust grains in the shell. The radiation pressure is estimated from the bolometric luminosity emitted by the SSC, L_cl = 5.3 × 10⁹L_⊙ and the H ii region radius R_in = 35 pc (Brγ size). We compute (2)There are large uncertainties for the age and mass of the cluster (Sect. 4.1) and the value given in Eq. (2) is a lower limit obtained for the lowest value of L_cl in Table 2. If we use the L_cl estimated for 3.5 Myr and 1.5 × 10⁷M_⊙, the age and mass computed from the SINFONI data, the radiation pressure is three times higher. In SSC B1, the moderate value of the effective extinction in the K-band (see Table 2) preclude large values of ⟨ τ_rad ⟩.

The observational results presented in Sect. 4 allow us to model the physical environment of SSC B1 (see Appendix A for more details). We model the H ii region surrounding the SSC as an ionized nebula with dust grains, with incident radiation field of that from the SSC. For different sizes of the internal cavity, we solve the radiative transfer equation within the H ii region and measure the outward radiation field and the pressure exerted at the surface of the PDR. We assume that the gas pressure is set by the radiation pressure and compare our grid of P_rad and radiation field values with the outputs of the PDR Meudon code (Le Bourlot et al. 2012), which predicts molecular line intensities from PDRs. We chose the standard isobaric model with fixed density. These models are detailed in Appendix A. We find that PDR models with radiation fields in a range of χ~10³–10⁴ times the mean radiation field in the solar neighborhood and a gas pressure of 3 × 10⁷−10⁸ K cm^-3 are needed to account for the observed H₂ 1−0/2−1 S(1) line emission (see Fig. A.1). The value of the radiation field is consistent with that estimated by Gilbert et al. (2000) from the comparison of PDR models and several near-IR H₂ line intensities (obtained with single-slit spectroscopy towards SSC B1). The estimated gas pressure of 3 × 10⁷−10⁸ K cm^-3 agrees with the measured molecular gas pressure (Eq. (2)), supporting low values of ⟨ τ_rad ⟩. There is no significant trapping of the IR photons within the cavity of the molecular gas surrounding SSC B1.

5.2. Gas pressure in the hot gas

SSC B1 was detected as a compact X-ray source by Chandra (source #88 in Table 5 of Zezas et al. 2006). Its X-ray luminosity, integrated between 0.1–10 keV, is L_X = 1.9 ± 0.2 × 10³⁸ erg s^-1. This X-ray luminosity can be used to estimate the electronic density n_e and thus the hot gas pressure. (3)where ϵ(T) = 3 × 10^-23 erg s^-1 cm³ is the emission coefficient (assuming a temperature of T = 10⁷ K, from Fig. B.1 in Guillard et al. 2009), is the emission measure, and is the area where the pressure is exerted. For the observed value of L_X, we find an electron density of n_e = 0.6 × ( [35 pc] /R_in)^3/2 cm^-3. The pressure of the hot plasma is P_hot,X = 1.9 n_ek_BT. For a temperature of T = 10⁷ K, the hot plasma pressure is (4)This value is a lower limit because we did not correct the X-ray luminosity for absorption by the shell. Figure 1 in Ryter (1996) shows that the absorption in the X-ray by the gas is comparable to the extinction in the K-band and it is thus a factor of approximately 2, increasing P_hot,X by a factor 1.4.

We find that the hot gas pressure is at least three times weaker than the pressure of the molecular gas. This difference indicates that the hot gas pressure is not dynamically important. The hot plasma produced by shocks driven by stellar winds (Castor et al. 1975) is not trapped within a closed shell. If the mechanical energy from stellar winds over 1 Myr (age of SSC B1) was confined within a closed shell of interstellar matter, then the gas pressure in the hot gas would be more than 50 times higher than that estimated in Eq. (4) (if we use the age and masses estimated from the Brγ EW and radio N_Lyc, then it would be three orders of magnitude higher). This result is consistent with our finding that ⟨ τ_rad ⟩ is negligible. It also implies that the molecular gas interfaces heated by the SSC must be distributed over a range of distances. The estimated molecular radius of ~50 pc (Table 2) is a mean value.

6. Radiation pressure feedback

In this section we discuss the role of radiation pressure as a stellar feedback process in SSC B1.

6.1. Acceleration of the gas by stellar feedback

In Sect. 5, we show that radiation pressure is the dominant pressure on the molecular gas closest to the cluster, in agreement with previous studies in massive clusters (e.g., Krumholz & Matzner 2009; Murray et al. 2010; Lopez et al. 2011). We compare this outward pressure with the gravitational attraction exerted by the central cluster and the mass of the cloud itself. For a spherical distribution of the matter, we estimate the inward force per unit surface as, (5)where r and M are the radius and mass of the GMC, and M_cl is the mass of the SSC taken from Table 2. We note that the factor 2 in the uncertainty on the GMC mass (see Sect. 4.4) introduces an uncertainty of a factor 3 in P_grav. The radiation pressure, as estimated in Eq. (2), is weaker than the gravitational force. It is not strong enough to push away the molecular gas associated with SSC B1.

In the following section, we speculate that most (or at least a significant part) of the pre-cluster gas has already been blown away and that the high-velocity component is tracing outflowing molecular gas. From the ALMA observations, Whitmore et al. (2014) also identify this component as a super-bubble. Earlier, this gas would have been accelerated by radiation pressure when it was much higher than today.

6.2. Outflowing gas

The momentum of the outflowing gas can be estimated from the molecular mass in the high-velocity component of CO(3−2) (taken from Table 4) and the expansion velocity of that gas v_exp. The latter value can be estimated using different assumptions that we discuss in Sect. 7.1. We find that the expansion velocity of the gas is v_exp ~ 80 km s^-1. It is larger than the escape velocity from the cluster (using the size measured in the K-band continuum emission of 66 pc), which is km s^-1. Uncertainties on the expansion velocity come from the unknown 3D geometry of the gas and from the mass of cluster.

To test the hypothesis that this gas was accelerated by radiation pressure over a time scale of t_feedback, we write, (6)where (7)with R_high ~ 100 pc, the mean size of the bubble-like shape structure observed in this component (Fig. 7), and τ_rad refers to the effective opacity of the cloud to radiation (Sect. 5.1). From Eq. (6), we can constrain the value of τ_rad as: (8)This value is higher than that observed today, but within plausible values earlier in the cluster evolution when the parent cloud was not yet disrupted. Indeed, in the models by Murray et al. (2010), τ_rad reaches values of several tens at the very beginning of the cluster evolution. Smaller values for τ_rad (~6) and v_exp (~60 km s^-1) are found if we use the age estimated from the Brγ EW (L_cl is higher, see Table 2).

We estimate a lower limit for the SFE within 100 pc from the cluster of 17%, considering that the high-velocity component may include not only natal gas but surrounding gas that did not participate in the cluster formation. We estimate the outflow rate of the parent molecular cloud to be Ṁ_outflow = M_high/t_feedback ≃ 30M_⊙/yr.

The broad velocity component is also observed in the H₂ 1−0 S(1) emission. The left panel in Fig. 10 presents the spectra for the v = 1–0 and 2−1 S(1) H₂ line emission integrated within an aperture of 4.̋2 (Sect. 4.2). The v = 2–1 line can be fitted with a single Gaussian curve (red-dashed line in Fig. 10) while v = 1–0 requires two Gaussian components (blue line in Fig. 10). We fit the v = 2–1 line and use the estimated center velocity as a fixed parameter to fit the low-velocity component of the v = 1−0 emission. The right panel of Fig. 10 shows the low- and high-velocity components, observed in v = 1–0, in red and black contours, respectively, overlaid on the CO(3−2) high-velocity component. The extent of the high-velocity components seen in CO(3−2) and H₂ 1−0 S(1) are similar. Comparing Fig. 10 with Fig. 8, we find that the extent of the low-velocity components seen in CO(3−2) and H₂ 1−0 S(1) are also similar.

Fig. 10 H2 2−1/1−0 S(1) ratio in SGMC 4/5. The left panel shows the spectra for the two near-IR H2 lines, H2 1−0 S(1) in solid-black and H2 2−1 S(1) in solid-red line scaled by a factor of 3.4. The dashed-red line is the Gaussian fit for the H2 2−1 S(1) line, while the blue-solid line is the fit for the H2 1−0 S(1) emission that corresponds to the addition of two Gaussian curves (the red and black dashed lines). The color image and white contours in the right panel show the integrated intensity of the high-velocity component observed in CO(3−2). Red and black contours are the H2 emission integrated in 1300–1512 km s-1 (low-velocity component) and 1580–1760 km s-1 (high-velocity component), respectively. Offset positions are as in Fig. 4.

Open with DEXTER

The R_S(1) = H₂ 2−1/1−0 S(1) ratio is commonly used to disentangle between H₂ gas heated by collisions or UV radiation. On the one hand, shock excitation is collisional and only low-ν H₂ levels are populated because the temperature of molecular gas is, at most, a few 1000 K; for higher temperatures the H₂ molecule is destroyed in collisions. On the other hand, H₂ excitation through UV pumping and fluorescence populates both high- and low-ν states (e.g., Herrera et al. 2011). We find R_S(1) = 0.32 ± 0.04 for the low-velocity component, and estimate an upper limit of 0.18 for the high-velocity component using the ratio between the 3σ emission for v = 2–1 (σ measured in the line-free channels) and the peak emission of the high-velocity component seen in v = 1–0. These values are listed in Table 4. The former ratio can be accounted for by UV heating of the gas in PDRs (see Sect. 5.1 and Appendix A). The latter ratio can be observed towards PDRs with high radiation fields (>10⁴) and gas pressure (≃10⁸ K cm^-3) (see Fig. A.1). However, since the high-velocity component is extended (a few 100 pc) and there is no evidence for extended massive stellar population, we favor shocks as the main gas heating mechanism. Values in the range 0.1–0.2 are reproduced in J- and C-shock models (Kristensen 2007) for gas with densities in the range 10⁴–10⁶ cm^-3 and shock velocities from 15 to 50 km s^-1. This supports the idea that the high-velocity component is tracing outflowing gas. This interpretation indicates that the two velocity components have different excitation mechanisms.

Evidence of high-velocity ionized gas can be observed in the Brγ line emission. Figure 11 presents the channel maps of the Brγ line emission overlaid on the CO(3−2) line emission. At the velocities of the low-velocity component (<1600 km s^-1), the Brγ emission is symmetrical. An excess of ionized gas emission, where the CO(3−2) has a low column density, is observed at velocities of 1600 ~ 1700 km s^-1, that we interpret as evidence of escaping ionized gas.

	Fig. 11 Contours are the channel maps of the ionized gas traced by Brγ overlaid on the CO(3−2) emission interpolated to the velocity of the SINFONI observations.
Open with DEXTER

7. Possible scenario

In our interpretation, it is too late to witness the first stages of the disruption of the parent molecular cloud by stellar feedback, when the parent cloud was still a bound cloud. Nevertheless, the high-velocity component may be tracing the outflowing gas from this cloud. In the following, we discuss this interpretation and its caveats.

7.1. Disruption of the parent cloud and outflowing gas

We propose that radiation pressure was the main mechanism disrupting the parent molecular cloud over a short timescale. At the beginning of the SSC evolution, radiation pressure was much higher than today due to the high density and opacity of the parent cloud. The IR photons from the dust heated by the UV photons from the cluster were highly trapped in the cloud implying a high τ_rad value and radiation pressure (see Eq. (1)). The momentum transfer rapidly accelerates the surrounding gas to velocities sufficiently high to escape from the cluster. Murray et al. (2010) modeled this disruption process for a source with similar characteristics (i.e., GMC mass and size, and star cluster mass), finding that the disruption occurs in less than 1 Myr.

Within this scenario, the high-velocity component is tracing outflowing gas from the parent cloud. The expansion time scale of 1.2 Myr, computed in Sect. 6.2, agrees with the model value from Murray et al. (2010). It is comparable to the age of the cluster (1–3.5 Myr), which makes this scenario of cloud disruption plausible. The bottom panel of Fig. 7 shows the channel maps of this component, which surrounds SSC B1 in a bubble-like shape structure with an inner radius of approximately 1′′ (~100 pc). Similar molecular gas morphology and kinematics are observed in the starburst galaxies M82 and NGC 253 (Weiß et al. 1999; Sakamoto et al. 2006), where the data have been interpreted as evidence of expanding molecular super-bubbles. Evidence of expanding bubbles has also been observed for ionized gas (e.g., in M83, Hollyhead et al. 2015).

In Fig. 12, we show the CO(3−2) emission averaged within concentric annuli centered on SSC B1 up to a distance of 1.̋8 from the cluster (azimuthally averaged CO(3−2) emission). For a closed, spherical and homogeneous expanding shell, the azimuthally averaged gas emission should show a semi-ellipse with constant flux. This is not the case for the super-bubble in SSC B1. The central velocity of the cluster is difficult to determine. We can approximate this velocity in two ways. First, we assume that the Brγ emission is symmetrical (no obscuration by dust) and therefore use the Brγ velocity as a proxy for the cluster velocity. In this case, we approximate the velocity expansion of the bubble as the difference between the Brγ velocity and the velocity of the high-velocity component, ~80 km s^-1. Figure 12 shows a cyan dashed-line semi-ellipse to illustrate this scenario. The size of the ellipse was chosen to fit the CO high-velocity emission. In this scenario, the CO data do not show any evidence of blue-shifted gas and we have to assume that the CO shell is not symmetric. Second, we estimate the velocity of the cluster fitting the CO azimuthally averaged emission. The red dashed-line in Fig. 12 shows the best fit for the CO emission (low- and high-velocity components), with a central velocity of ~1565 km s^-1, which would be tracing the expanding bubble. In this scenario, the central velocity of the ellipse is approximately 50 km s^-1 larger than that of the compact Brγ emission seen towards the SSC. This difference with the Brγ velocity may imply that the red-shifted Brγ gas is obscured and we only observe the gas that is the closest to us. An expanding ionized bubble with asymmetrical emission has been observed towards a young (6 Myr) and massive (10⁷M_⊙) SSC in a blue compact galaxy (Östlin et al. 2007).

Fig. 12 Azimuthally averaged CO(3−2) emission, around the position of SSC B1 up to a maximum radius of 1.̋8. The left panel shows the total emission and the right panel shows the emission clipped up to a maximum value of 20% of the peak emission. The semi-ellipse in dashed-red line highlights the expanding shell assuming a central velocity of ~1565 km s-1 marked with a red tick on the velocity axis, while for the dashed-cyan line, the central velocity is equal to that of the Brγ emission.

Open with DEXTER

7.2. The molecular gas surrounding SSC B1

If the parent cloud was disrupted early in the cluster evolution, we are curious as to why we see a CO emission peak coincident in velocity and position with the cluster. Within our scenario, we speculate that the low-velocity component traces gas that was near to the SSC when it formed but not part of its parent cloud. It may also include clumps that have migrated from the SGMC environment to the cluster neighborhood. This gas would have been too far from the SSC to be accelerated by radiation pressure. This possibility is exemplified by the 30 Doradus star-forming region. The gas surrounding the central massive cluster R136 is distributed in an irregular shell, where some molecular clouds that did not participate in the formation of the R136 cluster are observed very close to it (≲10 pc, Rubio et al. 2009).

We expect the low-velocity molecular gas to be dispersed by the action of stellar winds and future supernova explosions. Stellar winds from massive clusters can erode the surrounding gas and push away the clump material (Rogers & Pittard 2013).

In the scenario that the gas traced by the low-velocity component is indeed the parent molecular cloud, we do not really understand how the cluster will expel and disrupt its surrounding matter. This component is still very massive, comparable with the mass of the cluster, and, furthermore, it is self-bound. We need observations at higher resolution to resolve this component.

7.3. Star-formation efficiency

We are interested to know if SSC B1 is a bound cluster that will evolve as a globular cluster. As described in the introduction, a SFE lower limit of 30% is needed to create a bound stellar cluster. Within our scenario, for a region of ~100 pc surrounding SSC B1, we estimate a lower limit for the SFE of 17%. In the extreme, unlikely case that both the low- and high-CO velocity components were part of the parent molecular cloud, we estimate SFE> 10%. If we use the Galactic CO conversion factor, then the SFE decreases to ~5%. The estimated value is higher than what is measured in Galactic GMCs, but not high enough to ensure the formation of a bound star cluster. To reach the theoretical SFE limit, at least 30% of the high-velocity component should not be natal gas but, for example, surrounding gas that was swept-away by stellar feedback. This is possible for SSC B1, which is located in a complex interacting region rich in molecular gas. However, as discussed in the introduction, hydrodynamical simulations indicate that even if the global SFE is low, locally, the SFE can reach values large enough to form a bound cluster (Kruijssen et al. 2012).

Further observations are needed to quantify and test our proposed scenario. Line observations at high resolution of the stellar emission of the cluster would help settle the issues related to the cluster velocity discussed in Sect. 7.1. Higher-resolution observations will reveal the spatial distribution of the gas very close to the cluster, the low-velocity component. With the existing data, we cannot resolve structures smaller than 50 pc and, therefore, we cannot resolve and characterize the clumpy structure of this source. More sensitive observations will enable us to better trace the morphology of the high-velocity component. Also, by looking for high-velocity wings, we can elucidate whether or not the molecular gas is being carried out in a wind together with the hot gas. This can now be done with ALMA, which offers the required angular resolution and sensitivity to address these questions.

8. Conclusions

In this paper, we investigate the early evolution of SSC B1 located in SGMC 4/5 in the Antennae overlap region, and, specifically, the impact of the stellar feedback mechanisms on the surrounding matter.

We search for associations between SSCs and compact H ii and H₂ line emission in four FOVs observed with SINFONI/VLT across the Antennae overlap region. We focus on isolated SSCs previously studied by Gilbert & Graham (2007). Among our sample, there is only one SSC with both H ii and H₂ emission; SSC B1 in SGMC 4/5. This is the only SSC younger than 5 Myr. This is probably not a complete sample since it is difficult to isolate the molecular and ionized emission associated with SSCs. Nevertheless, our sample shows that the association between SSCs and their parent molecular clouds does not last much more than 5 Myr; similar to what is found in a sample of SSCs in nearby galaxies and other regions of the Antennae pair (Bastian et al. 2014).

CO emission in SSC B1 presents two velocity components, a compact, low-velocity one at the position of the SSC, and an extended, high-velocity one distributed in a bubble-like shape structure around the cluster. The virial and millimetric masses of the gas associated with SSC B1 low-velocity component are comparable. The comparison with its CO luminosity yields a CO-to-H₂ factor consistent with that found at larger scales (1.5 kpc) and in other starburst galaxies. The detection of molecular gas associated with SSC B1 suggests that the cluster may still be embedded in its parent molecular cloud. Although the actual velocity of the cluster is still unknown and thus it is not clear if the cluster is embedded in this cloud, this source is the only one in the overlap region that could tell us how the molecular gas is being pushed away from the cluster.

We analyzed the physical properties of SSC B1 and the molecular and ionized gas surrounding it. A model of the physical environment of the H ii region and the PDR surrounding SSC B1 shows that the observed H₂ 2−1/1−0 S(1) intensity line ratio is accounted for by a gas pressure of 3 × 10⁷−10⁸ K cm^-3, which agrees with the estimated radiation pressure from the cluster. This finding implies that there is no significant trapping of the IR photons within the cavity of the molecular gas surrounding the cluster. X-ray observations show that the plasma produced by shocks driven by stellar winds is not confined within a closed shell of cold gas. The plasma pressure is smaller than the radiation pressure. However, radiation pressure is weaker than the gravitational force from the cluster and the cloud itself. It is not strong enough to push away the observed molecular gas.

We propose that radiation pressure was early enhanced and was the main mechanism driving the disruption of the parent molecular cloud in ~1 Myr. The high-velocity component would be tracing outflowing molecular gas from this disruption. The association of this component with shock heated H₂ gas, and the presence of escaping high Brγ gas, are supporting evidence. The low-velocity component may be gas that was near to the SSC when it formed but not part of its parent cloud or clumps that migrated from the SGMC environment. This component would be disrupted by the action of stellar winds and supernova explosions. The existing data do not allow us to conclude whether or not SSC B1 is gravitationally bound and will evolve as a globular cluster. The lower limit for the SFE of 17%, estimated in a region of 100 pc size from the cluster, is slightly smaller than the theoretical limit for bound cluster formation.

Further observations are needed to quantify and test our proposed scenario. Additional observations of line emission from the cluster would allow us to directly measure the velocity of the cluster. Higher spatial and spectral-resolution observations are needed to probe the clumpy structure of the closest component to the cluster.

Acknowledgments

We would like to thank the anonymous referee for the detailed and constructive report that helped us to improve our manuscript. The Atacama Large Millimeter/submillimeter Array (ALMA), an international astronomy facility, is a partnership of Europe, North America and East Asia in cooperation with the Republic of Chile. This paper makes use of the following ALMA Science Verification data: ADS/JAO.ALMA#2011.0.00003.SV and Early Science data: ADS/JAO.ALMA#2011.0.00876.S.

References

Appendix A: Modeling the physical environment of SSC B1

In this section, we introduce a simple geometrical model that we use to quantify the action of the stellar feedback in the parent cloud of SSC B1. The geometrical model is combined with a PDR model to quantify the H₂ line emission. We introduce the model in two parts, we model the H ii region and then the PDR. This 1D model describes a segment of the shell. It does not assume that the cluster is fully embedded in a spherical shell. The matter around the SSC is likely to be clumpy, with some radiation escaping in some directions. This is in agreement with the fact that the mean extinction of the shell is moderate (a factor 2, Sect. 4.1). The H ii region is directly receiving the UV radiation from the SSC and the incident radiation on the PDR is the output radiation from the H ii gas layer. The internal radius of the H ii cavity R_in is a parameter of the model and the external radius is the GMC radius observed in the continuum emission, R_out = 52 pc. The stellar radiation was chosen to be that from the Starburst99 models, for an age of 3 Myr scaled by the cluster mass. This SED was used for the H ii region modeling as the ionizing radiation.

Appendix A.1: H II region

Fig. A.1 H2 2−1/1−0 S(1) intensity line ratio estimated from the PDR Meudon code for a grid of radiation fields and gas pressures. The measured value for the low-velocity component is 0.32, which is highlighted in the figure. The red area marks 1σ uncertainty associated with this value. Triangles represent the values given in Table A.1, estimated from the H ii/PDR model for radii from 30 to 100 pc.

Open with DEXTER

We model the H ii region as an ionized hydrogen nebula with dust grains. We compute the radiative transfer equation for this nebula in the plane-parallel approximation. In this nebula, the radiative transfer will be defined by the absorption produced by the ionization of the hydrogen atoms and light absorption and scattering by dust grains, (A.1)where s is the position in the nebula from R_in at which the absorption occurs, n_H and n_d, and σ_H and σ_d are the hydrogen and dust densities and cross sections, respectively. σ_H is zero for non-ionizing photons. We assume that the nebula extent is small compared to R_in. The ionization rate () equals the recombination rate (N_rec = n_en₊α₂), where is the mean photoionization cross section, S = S₀ e^{− τ} is the number of ionizing photons per unit time and unit area, n_e is the electron density, n₊ is the density of protons and α₂ is the recombination coefficient. Then, we can re-write Eq. (A.1) as: (A.2)where n is the gas density and x_d is the dust-to-gas mass ratio. We assume that all the protons and electrons come from the same atoms (n = n_e = n₊). The dust cross-sections are taken from the DUSTEM model (Desert et al. 1990; Compiegne et al. 2011). We include the extinction cross sections of the very small grains and big grains (there is no observational evidence for PAHs in H ii regions, i.e., Kassis et al. 2006).

We solve the radiative transfer equation (Eq. (A.2)) for different values of incident radiation on the H ii region, depending on the value of the unknown quantity R_in, which we vary from 30 to 100 by spacings of 10 pc. This range is distributed around the GMC radius (~50 pc), from the radius of the Brγ line emission to that of H₂. We compute the output flux in units of the Mathis radiation field (Mathis et al. 1983), integrating the emission between 91.2 to 200 nm. Table A.1 lists the radiation field for the different internal radii of the cavity. We also compute the radiation pressure exerted at the surface of the PDR from photons absorbed or scattered across the layer of the H ii gas, (A.3)where the luminosity L_HII is the luminosity incident onto the H ii region minus the output luminosity. P_rad(HII) is the pressure produced by the ionizing photons and the non-ionizing photons, which are absorbed by dust and gas in the H ii layer. This is the minimum value of the pressure at the surface of the PDR. Table A.1 lists the estimated values of P_rad.

Appendix A.2: PDR

To describe the PDR, we use the online tables of outputs from the PDR Meudon Code² (Le Bourlot et al. 2012). We chose to use results from isobaric rather than isochoric models to relate line intensities to the radiation pressure. Model results are available for a range of values of the incident radiation field and gas pressures. The model grid ranges are P/k_B = nT_K = 1 × 10⁵−1 × 10⁸ K cm^-3, for the gas pressure and χ = 1 to 1 × 10⁶ for the radiation field in Mathis units.

Figure A.1 shows the model results for the H₂2−1/1−0 S(1) line ratio. The red filled area indicates the 1σ uncertainty on the observed ratio for the low-velocity component (Table 4). The range of possible combinations of radiation fields and pressures given in Table A.1 are marked by triangles. This figure indicates that PDR models can account for the observed intensity ratios within the constrains on physical parameters set by the spatial extent of the emission lines. In the paper, we use a value R_in = 35 pc, which corresponds to the value observed in Brγ. From the model in Fig. A.1, the best R_in value is ~55 pc.

All Tables

All Figures

	Fig. 1 K-band emission from the overlap region obtained with the CFHT (Herrera et al. 2011) overlaid with Brγ contours from the SINFONI observations. Small black squares highlight the position of the SSCs lying within the four SINFONI FOVs, which are marked with dashed boxes. An inserted image of the Antennae galaxies shows, in a black rectangle, the relative position of the overlap region in the merger.
Open with DEXTER
In the text

	Fig. 2 H2 1−0 S(1) emission over the star-cluster mass, for SSCs with different ages and PCC 1 (Herrera et al. 2011), observed in the SINFONI fields. Since both the masses and H2 1−0 S(1) fluxes depend on the extinction in the same way, the H2-to-stellar mass ratio is independent of extinction. PCC 1 has a lower limit since it has not yet formed massive stars. The H2-to-stellar mass ratio is normalized to the value for B1.
Open with DEXTER
In the text

	Fig. 3 Brγ emission over the star cluster mass, for the SSCs in the SINFONI fields. As for H2 (see Fig. 2), the Brγ-to-stellar mass ratio is extinction independent. The ratio is normalized to the value for B1.
Open with DEXTER
In the text

	Fig. 4 Continuum emission at 345 GHz. Emission from H2O maser was subtracted. The white cross marks the position of the H2O maser. Solid black lines are the 50%, 70%, and 90% of the peak of the Brγ line emission (see Fig. 6). Dashed black contours correspond to the 2D Gaussian fit to the continuum emission. Offset positions are relative to α: 12h01m5495, δ: –18°53′05.̋6 J2000.0. The offset position (0, 0) corresponds to the position of the SSC B1, which coincides with the continuum peak emission.
Open with DEXTER
In the text

	Fig. 5 Upper panel: K-band spectrum of SSC B1, obtained from aperture photometry based on the Brγ line. Ionized gas lines are the brightest lines across the spectrum. We mark the position of the H2 rovibrational lines with short red lines. Lower panels: zoom in the H2 1−0 S(1) and Brγ lines. Red lines represent the fitted Gaussian curves to each line, and residuals are plotted in the bottom part. Table 3 lists the parameters of the Gaussian fits.
Open with DEXTER
In the text

	Fig. 6 Intensity maps for SGMC 4/5. From left to right, K-band continuum emission contours overlaid on the H2 1−0 S(1) emission, CO(3−2) emission, Brγ emission, and continuum emission at 345 GHz. Flux units of the color bars are × 10-17 erg s-1 cm-2 for H2 and Brγ, Jy beam-1 for CO(3−2), and mJy beam-1 for the continuum emission. The 1σ emission level of the K-band continuum emission is 2.05 × 10-20 erg s-1 cm-2. The peak emission of the K-band contours marks the position of the SSC B1, which corresponds to the offset position (0, 0). Offset positions are as in Fig. 4.
Open with DEXTER
In the text

	Fig. 7 CO(3−2) channel maps of SGMC 4/5. The position of SSC B1 is marked with a cross and corresponds to the (0,0) offset position. Top panel: we present channel maps from 1440 km s-1 to 1550 km s-1 by 10 km s-1. Contours start at 10σ with a spacing of 10σ, with σ = 0.39 Jy beam-1 km s-1 measured in line-free channels. Bottom panel: CO channel maps from 1550 km s-1 to 1670 km s-1 by 10 km s-1. Contours start at 5σ with a spacing of 5σ. Offset positions are as in Fig. 4.
Open with DEXTER
In the text

	Fig. 8 CO(3−2) velocity components in SGMC 4/5. Left panels: CO emission integrated between 1370 km s-1 and 1540 km s-1 (top) and 1550 km s-1 and 1660 km s-1 (bottom). Contours represent the 10%, 30%, 50%, 70%, and 90% of the peak emission. Right panels: first moment of both velocity components. Contours are the same as in the left panel. Offset positions are as in Fig. 4.
Open with DEXTER
In the text

	Fig. 9 CO(3−2) and CO(2−1) spectra associated with SSC B1. For the CO J = 3−2 transition, we present the spectra measured with two apertures, 2.̋9 and 4.̋2, to characterize the two velocity components observed in the channel maps in Fig. 7.
Open with DEXTER
In the text

	Fig. 10 H2 2−1/1−0 S(1) ratio in SGMC 4/5. The left panel shows the spectra for the two near-IR H2 lines, H2 1−0 S(1) in solid-black and H2 2−1 S(1) in solid-red line scaled by a factor of 3.4. The dashed-red line is the Gaussian fit for the H2 2−1 S(1) line, while the blue-solid line is the fit for the H2 1−0 S(1) emission that corresponds to the addition of two Gaussian curves (the red and black dashed lines). The color image and white contours in the right panel show the integrated intensity of the high-velocity component observed in CO(3−2). Red and black contours are the H2 emission integrated in 1300–1512 km s-1 (low-velocity component) and 1580–1760 km s-1 (high-velocity component), respectively. Offset positions are as in Fig. 4.
Open with DEXTER
In the text

	Fig. 11 Contours are the channel maps of the ionized gas traced by Brγ overlaid on the CO(3−2) emission interpolated to the velocity of the SINFONI observations.
Open with DEXTER
In the text

	Fig. 12 Azimuthally averaged CO(3−2) emission, around the position of SSC B1 up to a maximum radius of 1.̋8. The left panel shows the total emission and the right panel shows the emission clipped up to a maximum value of 20% of the peak emission. The semi-ellipse in dashed-red line highlights the expanding shell assuming a central velocity of ~1565 km s-1 marked with a red tick on the velocity axis, while for the dashed-cyan line, the central velocity is equal to that of the Brγ emission.
Open with DEXTER
In the text

	Fig. A.1 H2 2−1/1−0 S(1) intensity line ratio estimated from the PDR Meudon code for a grid of radiation fields and gas pressures. The measured value for the low-velocity component is 0.32, which is highlighted in the figure. The red area marks 1σ uncertainty associated with this value. Triangles represent the values given in Table A.1, estimated from the H ii/PDR model for radii from 30 to 100 pc.
Open with DEXTER
In the text