We discuss a theoretical model for the early evolution of massive star clusters and confront it with the ALMA, radio and infrared observations of the young stellar cluster highly obscured by the molecular cloud D1 in the nearby dwarf spheroidal galaxy NGC 5253. We show that a large turbulent pressure in the central zones of D1 cluster may cause individual wind-blown bubbles to reach pressure confinement before encountering their neighbors. In this case stellar winds are added to the hot shocked wind pockets of gas around individual massive stars that leads them to meet and produce a cluster wind in time-scales less than $10^5$ yrs. In order to inhibit the possibility of cloud dispersal, or the early negative star formation feedback, one should account for mass loading that may come, for example, from pre-main sequence (PMS) low-mass stars through photo-evaporation of their proto-stellar disks. Mass loading at a rate in excess of 8$\times 10^{-9}$ M$_{\odot}$ yr$^{-1}$ per each PMS star is required to extend the hidden star cluster phase in this particular cluster. In this regime, the parental cloud remains relatively unperturbed, while pockets of molecular, photoionized and hot gas coexist within the star forming region. Nevertheless, the most likely scenario for cloud D1 and its embedded cluster is that the hot shocked winds around individual massive stars should merge at an age of a few millions of years when the PMS star proto-stellar disks vanish and mass loading ceases that allows a cluster to form a global wind.

We present a physically-motivated model involving the different processes affecting supernova dust grains as they are incorporated into the thermalized medium within young massive star clusters. The model is used to explain the near- to mid-infrared (NIR-MIR) excess found in such clusters and usually modeled as a blackbody with temperature $\sim (400-1000)$ K. In our approach, dust grains are efficiently produced in the clumpy ejecta of core-collapse supernovae, fragmented into small pieces ($\lesssim 0.05$ $\mu$m) as they are incorporated into the hot and dense ISM, heated via frequent collisions with electrons and the absorption of energetic photons. Grains with small sizes can more easily acquire the high temperatures ($\sim 1000$ K) required to produce a NIR-MIR excess with respect to the emission of foreground PAHs and starlight. However, the extreme conditions inside young massive clusters make difficult for these small grains to have a persistent manifestation at NIR-MIR wavelengths as they are destroyed by efficient thermal sputtering. Nevertheless, the chances for a persistent manifestation are increased by taking into account that small grains become increasingly transparent to their impinging ions as their size decreases. For an individual SN event, we find that the NIR-MIR excess last longer if the time required to incorporate all the grains into the thermalized medium is also longer, and in some cases, comparable to the characteristic interval between supernova explosions. Our models, can successfully explain the near-infrared excesses found in the star clusters observed in M33 (Rela\~no et al. 2016) assuming a low heating efficiency and mass-loading. In this scenario, the presence of the NIR-MIR excess is an indication of efficient dust production in SNe and its subsequent destruction

The growing evidence pointing at core-collapse supernovae as large dust producers makes young massive stellar clusters ideal laboratories to study the evolution of dust immersed into a hot plasma. Here we address the stochastic injection of dust by supernovae and follow its evolution due to thermal sputtering within the hot and dense plasma generated by young stellar clusters. Under these considerations, dust grains are heated by means of random collisions with gas particles which results on the appearance of infrared spectral signatures. We present time-dependent infrared spectral energy distributions which are to be expected from young stellar clusters. Our results are based on hydrodynamic calculations that account for the stochastic injection of dust by supernovae. These also consider gas and dust radiative cooling, stochastic dust temperature fluctuations, the exit of dust grains out of the cluster volume due to the cluster wind and a time-dependent grain size distribution.

From the structure of PHL 293B and the physical properties of its ionizing cluster and based on results of hydrodynamic models, we point at the various events required to explain in detail the emission and absorption components seen in its optical spectrum. We ascribe the narrow and well centered emission lines, showing the low metallicity of the galaxy, to an HII region that spans through the main body of the galaxy. The broad emission line components are due to two off-centered supernova remnants evolving within the ionizing cluster volume and the absorption line profiles are due to a stationary cluster wind able to recombine at a close distance from the cluster surface, as originally suggested by Silich et al. (2004). Our numerical models and analytical estimates confirm the ionized and neutral column density values and the inferred X-ray emission derived from the observations.

Massive young stellar clusters are strong sources of radiation and mechanical energy. Their powerful winds and radiation pressure sweep-up interstellar gas into thin expanding shells which trap the ionizing radiation produced by the central clusters affecting the dynamics and the distribution of their ionized gas. Here we continue our comparison of the star cluster winds and radiation pressure effects on the dynamics of shells around young massive clusters. We calculate the impact that radiation pressure has on the distribution of matter and thermal pressure within such shells as well as on the density weighted ionization parameter $U_w$ and put our results on the diagnostic diagram which allows one to discriminate between the wind-dominated and radiation-dominated regimes. We found that model predicted values of the ionization parameter agree well with typical values found in local starburst galaxies. Radiation pressure may affect the inner structure and the dynamics of wind-driven shells only at the earliest stages of evolution or if a major fraction of the star cluster mechanical luminosity is dissipated or radiated away within the star cluster volume and thus the star cluster mechanical energy output is significantly smaller than star cluster synthetic models predict. However, even in these cases radiation dominates over the wind dynamical pressure only if the exciting cluster is embedded into a high density ambient medium.

Hydrodynamic models for spherically-symmetric winds driven by young stellar clusters with a generalized Schuster stellar density profile are explored. For this we use both semi-analytic models and 1D numerical simulations. We determine the properties of quasi-adiabatic and radiative stationary winds and define the radius at which the flow turns from subsonic into supersonic for all stellar density distributions. Strongly radiative winds diminish significantly their terminal speed and thus their mechanical luminosity is strongly reduced. This also reduces their potential negative feedback into their host galaxy ISM. The critical luminosity above which radiative cooling becomes dominant within the clusters, leading to thermal instabilities which make the winds non-stationary, is determined, and its dependence on the star cluster density profile, core radius and half mass radius is discussed.

A hydrodynamic model for steady state, spherically-symmetric winds driven by young stellar clusters with an exponential stellar density distribution is presented. Unlike in most previous calculations, the position of the singular point R_sp, which separates the inner subsonic zone from the outer supersonic flow, is not associated with the star cluster edge, but calculated self-consistently. When the radiative losses of energy are negligible, the transition from the subsonic to the supersonic flow occurs always at R_sp ~ 4 R_c, where R_c is the characteristic scale for the stellar density distribution, irrespective of other star cluster parameters. This is not the case in the catastrophic cooling regime, when the temperature drops abruptly at a short distance from the star cluster center and the transition from the subsonic to the supersonic regime occurs at a much smaller distance from the star cluster center. The impact from the major star cluster parameters to the wind inner structure is thoroughly discussed. Particular attention is paid to the effects which radiative cooling provides to the flow. The results of the calculations for a set of input parameters, which lead to different hydrodynamic regimes, are presented and compared to the results from non-radiative 1D numerical simulations and to those from calculations with a homogeneous stellar mass distribution.