We present observations of CO(3-2) and $^{13}$CO(3-2) emission near the supernebula in the dwarf galaxy NGC 5253, which contains one of the best examples of a potential globular cluster in formation. The 0.3" resolution images reveal an unusual molecular cloud, "Cloud D1", coincident with the radio-infrared supernebula. The ~6-pc diameter cloud has a linewidth, $\Delta$ v = 21.7 km/s, that reflects only the gravitational potential of the star cluster residing within it. The corresponding virial mass is 2.5 x 10$^5$ M$_\odot$. The cluster appears to have a top-heavy initial mass function, with $M_{low}$~1-2 M$_\odot$. Cloud D1 is optically thin in CO(3-2) probably because the gas is hot. Molecular gas mass is very uncertain but constitutes < 35% of the dynamical mass within the cloud boundaries. In spite of the presence of an estimated ~1500-2000 O stars within the small cloud, the CO appears relatively undisturbed. We propose that Cloud D1 consists of molecular clumps or cores, possibly star-forming, orbiting with more evolved stars in the core of the giant cluster.

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.

Following the observational and theoretical evidence that points at core collapse supernovae as major producers of dust, here we calculate the hydrodynamics of the matter reinserted within young and massive super stellar clusters under the assumption of gas and dust radiative cooling. The large supernova rate expected in massive clusters allows for a continuous replenishment of dust immersed in the high temperature thermalized reinserted matter and warrants a stationary presence of dust within the cluster volume during the type II supernova era. We first show that such a balance determines the range of dust to gas mass ratio and this the dust cooling law. We then search for the critical line that separates stationary cluster winds from the bimodal cases in the cluster mechanical luminosity (or cluster mass) vs cluster size parameter space. In the latter, strong radiative cooling reduces considerably the cluster wind mechanical energy output and affects particularly the cluster central regions, leading to frequent thermal instabilities that diminish the pressure and inhibit the exit of the reinserted matter. Instead matter accumulates there and is expected to eventually lead to gravitational instabilities and to further stellar formation with the matter reinserted by former massive stars. The main outcome of the calculations is that the critical line is almost two orders of magnitude or more, depending on the assumed value of the adiabatic wind terminal speed, lower than when only gas radiative cooling is applied. And thus, many massive clusters are predicted to enter the bimodal regime.

By means of 3-D hydrodynamical simulations, here we evaluate the impact that supernova explosions occurring within wind-driven bubbles have on the survival or destruction of dust grains. We consider both, the dust generated within the ejecta and the dust initially present in the ambient gas and later locked-up in the surrounding wind-driven shell. The collision of the supernova blast wave with the wind-driven shell leads to a transmitted shock that moves into the shell and a reflected shock into the ejecta. The transmitted shock is capable of destroying large amounts of the dust locked in the shell, but only if the mass of the wind-driven shell is small, less than a few tens the ejected mass. Conversely, massive wind-driven shells, with several times the ejected mass, lead upon the interaction to strong radiative cooling, which inhibits the Sedov-Taylor phase and weakens the transmitted shock, making it unable to traverse the wind-driven shell. In such a case, the destruction/disruption of the ambient dust is largely inhibited. On the other hand, the SNRs grow rapidly in the very tenuous region excavated by the stellar winds, and thus a large fraction of the dust generated within the ejecta is not efficiently destroyed by the supernova reverse shock, nor by the reflected shock. Our calculations favor a scenario in which core-collapse supernovae within sufficiently massive wind-driven shells supply more dust to the ISM than what they are able to destroy.

Outflows are a pervasive feature of mechanical feedback from super star clusters (SSC) in starburst galaxies, playing a fundamental role in galaxy evolution. Observations are now starting to confirm that outflows can undergo catastrophic cooling, suppressing adiabatic superwinds. Here we present a suite of one-dimensional, hydrodynamic simulations that study the ionization structure of these outflows and the resulting line emission generated by the cooling gas. We use the non-equilibrium atomic chemistry package within MAIHEM, our modified version of FLASH, which evolves the ionization state of the gas and computes the total cooling rate on an ion-by-ion basis. We find that catastrophically cooling models produce strong nebular line emission compared to adiabatic outflows. We also show that such models exhibit non-equilibrium conditions, thereby generating more highly ionized states than equivalent equilibrium models. When including photoionization from the parent SSC, catastrophically cooling models show strong C IV {\lambda}1549 and O VI {\lambda}1037 emission. For density bounded photoionization, He II {\lambda}1640, {\lambda}4686, C III] {\lambda}1908, Si IV {\lambda}1206, and Si III {\lambda}1400 are also strongly enhanced. These lines are seen in extreme starbursts where catastrophic cooling is likely to occur, suggesting that they may serve as diagnostics of such conditions. The higher ionization generated by these flows may help to explain line emission that cannot be attributed to SSC photoionization alone.

Supernova explosions and their remnants (SNRs) drive important feedback mechanisms that impact considerably the galaxies that host them. Then, the knowledge of the SNRs evolution is of paramount importance in the understanding of the structure of the interstellar medium (ISM) and the formation and evolution of galaxies. Here we study the evolution of SNRs in homogeneous ambient media from the initial, ejecta-dominated phase, to the final, momentum-dominated stage. The numerical model is based on the Thin-Shell approximation and takes into account the configuration of the ejected gas and radiative cooling. It accurately reproduces well known analytic and numerical results and allows one to study the SNR evolution in ambient media with a wide range of densities $n_{0}$. It is shown that in the high density cases, strong radiative cooling alters noticeably the shock dynamics and inhibits the Sedov-Taylor stage, thus limiting significantly the feedback that SNRs provide to such environments. For $n_{0}>5 \times 10^{5}$ cm$^{-3}$, the reverse shock does not reach the center of the explosion due to the rapid fall of the thermal pressure in the shocked gas caused by strong radiative cooling.

By adopting empirical estimates of the Helium enhancement (Delta Y) between consecutive stellar generations for a sample of Galactic globular clusters (GGC), we uniquely constraint the star formation efficiency of each stellar generation in these stellar systems. In our approach, the star formation efficiency is the central factor that links stellar generations as it defines both their stellar mass and the remaining mass available for further star formation, fixing also the amount of matter required to contaminate the next stellar generation. In this way, the star formation efficiency is here shown to be fully defined by the He enhancement between successive stellar generations in a GC. Our approach has also an impact on the evolution of clusters and thus considers the possible loss of stars through evaporation, tidal interactions and stellar evolution. We focus on the present mass ratio between consecutive stellar generations and the present total mass of Galactic globular clusters. Such considerations suffice to determine the relative proportion of stars of consecutive generations that remain today in globular clusters. The latter is also shown to directly depend on the values of Delta Y and thus the He enhancement between consecutive stellar generations in GGC places major constraints on models of star formation and evolution of GC.

Following the current debate on the fate of SN-condensed dust grains, here we explore by means of three-dimensional hydrodynamical simulations the interaction of dusty supernova remnants (SNRs) with the shocked winds of neighboring massive stars within young massive stellar clusters (SSCs). As a comparison, we have also explored the evolution of supernova remnants in the diffuse ISM with constant density. Since the hydrodynamics of SNRs is intimately related to the properties of their immediate environment, the lifecycle of dust grains in SNRs within SSCs is radically different from that in the diffuse ISM. Moreover, off-centered SNRs evolving in the steep density gradient established due to a star cluster wind experience a blowout phase: shell fragmentation due to protruding Rayleigh-Taylor instabilities and the venting of SN ejecta. Our main finding is that clustered SN explosions will cause a net increase in the amount of dust in the surroundings of young massive stellar clusters. Our analysis considers the multiple dust processing resulting from the passage of the SN reverse shock, including its reflection at the SNR's center, the injection of shocked stellar winds within the respective remnant's volume and the effect of secondary forward shocks produced in sequential SN explosions. In the simulations, we have on-the-fly calculated the rates of thermal sputtering and dust-induced radiative cooling provided an initial distribution of grain sizes and dust content. Fast-moving elongated dusty SN ejecta resemble mushroom clouds violently ascending in a stratified atmosphere after volcanic super-eruptions, where the pyroclasts carried by the clouds are wind-driven and eventually accumulate into the vast surroundings.

Gas expulsion or gas retention is a central issue in most of the models for multiple stellar populations and light element anti-correlations in globular clusters. The success of the residual matter expulsion or its retention within young stellar clusters has also a fundamental importance in order to understand how star formation proceeds in present-day and ancient star-forming galaxies and if proto-globular clusters with multiple stellar populations are formed in the present epoch. It is usually suggested that either the residual gas is rapidly ejected from star-forming clouds by stellar winds and supernova explosions, or that the enrichment of the residual gas and the formation of the second stellar generation occur so rapidly, that the negative stellar feedback is not significant. Here we continue our study of the early development of star clusters in the extreme environments and discuss the restrictions that strong radiative cooling and stellar mass segregation provide on the gas expulsion from dense star-forming clouds. A large range of physical initial conditions in star-forming clouds which include the star-forming cloud mass, compactness, gas metallicity, star formation efficiency and effects of massive stars segregation are discussed. It is shown that in sufficiently massive and compact clusters hot shocked winds around individual massive stars may cool before merging with their neighbors. This dramatically reduces the negative stellar feedback, prevents the development of the global star cluster wind and expulsion of the residual and the processed matter into the ambient interstellar medium. The critical lines which separate the gas expulsion and the gas retention regimes are obtained.

The High Altitude Water Cherenkov observatory (HAWC) is an air shower array devised for TeV gamma-ray astronomy. HAWC is located at an altitude of 4100 m a.s.l. in Sierra Negra, Mexico. HAWC consists of 300 Water Cherenkov Detectors, each instrumented with 4 photomultiplier tubes (PMTs). HAWC re-uses the Front-End Boards from the Milagro experiment to receive the PMT signals. These boards are used in combination with Time to Digital Converters (TDCs) to record the time and the amount of light in each PMT hit (light flash). A set of VME TDC modules (128 channels each) is operated in a continuous (dead time free) mode. The TDCs are read out via the VME bus by Single-Board Computers (SBCs), which in turn are connected to a gigabit Ethernet network. The complete system produces ~ 500 MB/s of raw data. A high-throughput data processing system has been designed and built to enable real-time data analysis. The system relies on off-the-shelf hardware components, an open-source software technology for data transfers (ZeroMQ) and a custom software framework for data analysis (AERIE). Multiple trigger and reconstruction algorithms can be combined and run on blocks of data in a parallel fashion, producing a set of output data streams which can be analyzed in real time with minimal latency (< 5 s). This paper provides an overview of the hardware set-up and an in-depth description of the software design, covering both the TDC data acquisition system and the real-time data processing system. The performance of these systems is also discussed.

We report ALMA observations of $^{12}$CO(3-2) and $^{13}$CO(3-2) in the gas-poor dwarf galaxy NGC 5253. These 0.3"(5.5 pc) resolution images reveal small, dense molecular gas clouds that are located in kinematically distinct, extended filaments. Some of the filaments appear to be falling into the galaxy and may be fueling its current star formation. The most intense CO(3-2) emission comes from the central $\sim$100 pc region centered on the luminous radio-infrared HII region known as the supernebula. The CO(3-2) clumps within the starburst region are anti-correlated with H$\alpha$ on $\sim$5 pc scales, but are well-correlated with radio free-free emission. Cloud D1, which enshrouds the supernebula, has a high $^{12}$CO/$^{13}$CO ratio, as does another cloud within the central 100 pc starburst region, possibly because the clouds are hot. CO(3-2) emission alone does not allow determination of cloud masses as molecular gas temperature and column density are degenerate at the observed brightness, unless combined with other lines such as $^{13}$CO.

We report the detection of CO(J=2-1) coincident with the super star cluster (SSC) Mrk 71-A in the nearby Green Pea analog galaxy, NGC 2366. Our NOEMA observations reveal a compact, ~7 pc, molecular cloud whose mass (10^5 M_sun) is similar to that of the SSC, consistent with a high star-formation efficiency, on the order of 0.5. There are two, spatially distinct components separated by 11 km/s. If expanding, these could be due to momentum-driven, stellar wind feedback. Alternatively, we may be seeing the remnant infalling, colliding clouds responsible for triggering the SSC formation. The kinematics are also consistent with a virialized system. These extreme, high-density, star-forming conditions inhibit energy-driven feedback; the co-spatial existence of a massive, molecular cloud with the SSC supports this scenario, and we quantitatively confirm that any wind-driven feedback in Mrk 71-A is momentum-driven, rather than energy-driven. Since Mrk 71-A is a candidate Lyman continuum emitter, this implies that energy-driven superwinds may not be a necessary condition for the escape of ionizing radiation. In addition, the detection of the nebular continuum emission yields an accurate astrometric position for the Mrk 71-A. We also detect four other massive, molecular clouds in this giant star-forming complex.

The ability of young stellar clusters to expel or retain the gas left over after a first episode of star formation is a central issue in all models aiming to explain multiple stellar populations and the peculiar light element abundance patterns in globular clusters. Recent attempts to detect the gas left over from star formation in present day clusters with masses similar to those of globular clusters did not reveal a significant amount of gas in the majority of them, which strongly restricts the scenarios of multiple stellar population formation. Here the conditions required to retain the gas left over from star formation within the natal star forming cloud are revised. It is shown that the usually accepted concept regarding the thermalization of the star cluster kinetic energy due to nearby stellar winds and SNe ejecta collisions must be taken with care in the case of very compact and dense star forming clouds where three star formation regimes are possible if one considers different star formation efficiencies and mass concentrations. The three possible regimes are well separated in the half-mass radius and in the natal gas central density vs pre-stellar cloud mass parameter space. The two gas free clusters in the Antennae galaxies and the gas rich cluster with a similar mass and age in the galaxy NGC 5253 appear in different zones in these diagrams. The critical lines obtained for clusters with a solar and a primordial gas metallicity are compared.

By adopting the empirical constraints related to the estimates of Helium enhancement ($\Delta Y$), present mass ratio between first and second stellar generations ($M_{1G}/M_{2G}$) and the actual mass of Galactic globular clusters ($M_{GC}$), we envisage a possible scenario for the formation of these stellar systems. Our approach allows for the possible loss of stars through evaporation or tidal interactions and different star formation efficiencies. In our approach the star formation efficiency of the first generation ($\epsilon_{1G}$) is the central factor that links the stellar generations as it not only defines both the mass in stars of the first generation and the remaining mass available for further star formation, but it also fixes the amount of matter required to contaminate the second stellar generation. In this way, $\epsilon_{1G}$ is fully defined by the He enhancement between successive generations in a GC. We also show that globular clusters fit well within a $\Delta Y$ {\it vs} $M_{1G}/M_{2G}$ diagram which indicates three different evolutionary paths. The central one is for clusters that have not loss stars, through tidal interactions, from either of their stellar generations, and thus their present $M_{GC}$ value is identical to the amount of low mass stars ($M_* \le$ 1 M$_\odot$) that resulted from both stellar generations. Other possible evolutions imply either the loss of first generation stars or the combination of a low star formation efficiency in the second stellar generation and/or a loss of stars from the second generation. From these considerations we derive a lower limit to the mass ($M_{tot}$) of the individual primordial clouds that gave origin to globular clusters.

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.

Our arguments deal with the early evolution of Galactic globular clusters and show why only a few of the supernovae products were retained within globular clusters and only in the most massive cases ($M \ge 10^6$ Msol), while less massive clusters were not contaminated at all by supernovae. Here we show that supernova blast waves evolving in a steep density gradient undergo blowout and end up discharging their energy and metals into the medium surrounding the clusters. This inhibits the dispersal and the contamination of the gas left over from a first stellar generation. Only the ejecta from well centered supernovae, that evolve into a high density medium available for a second stellar generation in the most massive clusters would be retained. These are likely to mix their products with the remaining gas, leading in these cases eventually to an Fe contaminated second stellar generation.

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.

The High Altitude Water Cherenkov (HAWC) observatory is a wide field-of-view detector sensitive to gamma rays of 100 GeV to a few hundred TeV. Located in central Mexico at 19 degrees North latitude and 4100 m above sea level, HAWC will observe gamma rays and cosmic rays with an array of water Cherenkov detectors. The full HAWC array is scheduled to be operational in Spring 2015. In this paper, we study the HAWC sensitivity to the gamma-ray signatures of high-mass (multi- TeV) dark matter annihilation. The HAWC observatory will be sensitive to diverse searches for dark matter annihilation, including annihilation from extended dark matter sources, the diffuse gamma-ray emission from dark matter annihilation, and gamma-ray emission from non-luminous dark matter subhalos. Here we consider the HAWC sensitivity to a subset of these sources, including dwarf galaxies, the M31 galaxy, the Virgo cluster, and the Galactic center. We simulate the HAWC response to gamma rays from these sources in several well-motivated dark matter annihilation channels. If no gamma-ray excess is observed, we show the limits HAWC can place on the dark matter cross-section from these sources. In particular, in the case of dark matter annihilation into gauge bosons, HAWC will be able to detect a narrow range of dark matter masses to cross-sections below thermal. HAWC should also be sensitive to non-thermal cross-sections for masses up to nearly 1000 TeV. The constraints placed by HAWC on the dark matter cross-section from known sources should be competitive with current limits in the mass range where HAWC has similar sensitivity. HAWC can additionally explore higher dark matter masses than are currently constrained.

Primordial Black Holes (PBHs) are gravitationally collapsed objects that may have been created by density fluctuations in the early universe and could have arbitrarily small masses down to the Planck scale. Hawking showed that due to quantum effects, a black hole has a temperature inversely proportional to its mass and will emit all species of fundamental particles thermally. PBHs with initial masses of ~5.0 x 10^14 g should be expiring in the present epoch with bursts of high-energy particles, including gamma radiation in the GeV - TeV energy range. The Milagro high energy observatory, which operated from 2000 to 2008, is sensitive to the high end of the PBH evaporation gamma-ray spectrum. Due to its large field-of-view, more than 90% duty cycle and sensitivity up to 100 TeV gamma rays, the Milagro observatory is well suited to perform a search for PBH bursts. Based on a search on the Milagro data, we report new PBH burst rate density upper limits over a range of PBH observation times. In addition, we report the sensitivity of the Milagro successor, the High Altitude Water Cherenkov (HAWC) observatory, to PBH evaporation events.

VAMOS was a prototype detector built in 2011 at an altitude of 4100m a.s.l. in the state of Puebla, Mexico. The aim of VAMOS was to finalize the design, construction techniques and data acquisition system of the HAWC observatory. HAWC is an air-shower array currently under construction at the same site of VAMOS with the purpose to study the TeV sky. The VAMOS setup included six water Cherenkov detectors and two different data acquisition systems. It was in operation between October 2011 and May 2012 with an average live time of 30%. Besides the scientific verification purposes, the eight months of data were used to obtain the results presented in this paper: the detector response to the Forbush decrease of March 2012, and the analysis of possible emission, at energies above 30 GeV, for long gamma-ray bursts GRB111016B and GRB120328B.

We present HST/STIS optical and Gemini/NIFS near-IR IFU spectroscopy, and archival HST imaging of the triplet of super star clusters (A1, A2 and A3) in the core of the M82 starburst. Using model fits to the STIS spectra, and the weakness of red supergiant CO absorption features (appearing at ~6 Myr) in the NIFS H-band spectra, the ages of A2 and A3 are $4.5\pm1.0$~Myr. A1 has strong CO bands, consistent with our previously determined age of $6.4\pm0.5$~Myr. The photometric masses of the three clusters are 4--$7\times10^5$~\Msol, and their sizes are $R_{\rm eff}=159$, 104, 59~mas ($\sim$2.8, 1.8, 1.0~pc) for A1,2 and 3. The STIS spectra yielded radial velocities of $320\pm2$, $330\pm6$, and $336\pm5$~\kms\ for A1,2, and 3, placing them at the eastern end of the $x_2$ orbits of M82's bar. Clusters A2 and A3 are in high density (800--1000~\cmt) environments, and like A1, are surrounded by compact H\two\ regions. We suggest the winds from A2 and A3 have stalled, as in A1, due to the high ISM ambient pressure. We propose that the 3 clusters were formed \textit{in-situ} on the outer $x_2$ orbits in regions of dense molecular gas subsequently ionized by the rapidly evolving starburst. The similar radial velocities of the 3 clusters and their small projected separation of $\sim 25$~pc suggest that they may merge in the near future unless this is prevented by velocity shearing.

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.

We describe measurements of GeV and TeV cosmic rays with the High-Altitude Water Cherenkov Gamma-Ray Observatory, or HAWC. The measurements include the observation of the shadow of the moon; the observation of small-scale and large-scale angular clustering of the TeV cosmic rays; the prospects for measurement of transient solar events with HAWC; and the observation of Forbush decreases with the HAWC engineering array and HAWC-30.

The High-Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory is designed to record air showers produced by cosmic rays and gamma rays between 100 GeV and 100 TeV. Because of its large field of view and high livetime, HAWC is well-suited to measure gamma rays from extended sources, diffuse emission, and transient sources. We describe the sensitivity of HAWC to emission from the extended Cygnus region as well as other types of galactic diffuse emission; searches for flares from gamma-ray bursts and active galactic nuclei; and the first measurement of the Crab Nebula with HAWC-30.

The High-Altitude Water Cherenkov Gamma Ray Observatory (HAWC) is under construction 4100 meters above sea level at Sierra Negra, Mexico. We describe the design and cabling of the detector, the characterization of the photomultipliers, and the timing calibration system. We also outline a next-generation detector based on the water Cherenkov technique.

The High-Altitude Water Cherenkov Gamma Ray Observatory (HAWC) is designed to perform a synoptic survey of the TeV sky. The high energy coverage of the experiment will enable studies of fundamental physics beyond the Standard Model, and the large field of view of the detector will enable detailed studies of cosmologically significant backgrounds and magnetic fields. We describe the sensitivity of the full HAWC array to these phenomena in five contributions shown at the 33rd International Cosmic Ray Conference in Rio de Janeiro, Brazil (July 2013).

The High Altitude Water Cherenkov (HAWC) observatory is an array of large water Cherenkov detectors sensitive to gamma rays and hadronic cosmic rays in the energy band between 100 GeV and 100 TeV. The observatory will be used to measure high-energy protons and cosmic rays via detection of the energetic secondary particles reaching the ground when one of these particles interacts in the atmosphere above the detector. HAWC is under construction at a site 4100 meters above sea level on the northern slope of the volcano Sierra Negra, which is located in central Mexico at 19 degrees N latitude. It is scheduled for completion in 2014. In this paper we estimate the sensitivity of the HAWC instrument to point-like and extended sources of gamma rays. The source fluxes are modeled using both unbroken power laws and power laws with exponential cutoffs. HAWC, in one year, is sensitive to point sources with integral power-law spectra as low as 5x10^-13 cm^-2 sec^-1 above 2 TeV (approximately 50 mCrab) over 5 sr of the sky. This is a conservative estimate based on simple event parameters and is expected to improve as the data analysis techniques are refined. We discuss known TeV sources and the scientific contributions that HAWC can make to our understanding of particle acceleration in these sources.

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.

This paper presents a hydrodynamical model describing the evolution of the gas reinserted by stars within a rotating young nuclear star cluster (NSC). We explicitly consider the impact of the stellar component to the flow by means of a uniform insertion of mass and energy within the stellar cluster. The model includes the gravity force of the stellar component and a central supermassive black hole (SMBH), and accounts for the heating from the central source of radiation and the radiative cooling of the thermalized gas. By using a set of parameters typical for NSCs and SMBHs in Seyfert galaxies our simulations show that a filamentary/clumpy structure is formed in the inner part of the cluster. This "torus" is Compton thick and covers a large fraction of the sky (as seen from the SMBH). In the outer parts of the cluster a powerful wind is produced, that inhibits the infall of matter from larger scales and thus the NSC-SMBH interplay occurs in isolation.

We report the discovery of extended X-ray emission within the young star cluster NGC 602 in the Wing of the Small Magellanic Cloud (SMC) based on observations obtained with the Chandra X-ray Observatory. X-ray emission is detected from the cluster core area with the highest stellar density and from a dusty ridge surrounding the HII region. We use a census of massive stars in the cluster to demonstrate that a cluster wind or wind-blown bubble is unlikely to provide a significant contribution to the X-ray emission detected from the central area of the cluster. We therefore suggest that X-ray emission at the cluster core originates from an ensemble of low- and solar-mass pre-main-sequence (PMS) stars, each of which would be too weak in X-rays to be detected individually. We attribute the X-ray emission from the dusty ridge to the embedded tight cluster of the new-born stars known in this area from infrared studies. Assuming that the levels of X-ray activity in young stars in the low-metallicity environment of NGC 602a are comparable to their Galactic counterparts, then the detected spatial distribution, spectral properties, and level of X-ray emission are largely consistent with those expected from low- and solar-mass PMS stars and young stellar objects (YSOs). This is the first discovery of X-ray emission attributable to PMS stars and YSOs in the SMC, which suggests that the accretion and dynamo processes in young, low-mass objects in the SMC resemble those in the Galaxy.

The impact of radiation pressure on the dynamics of the gas in the vicinity of young stellar clusters is thoroughly discussed. The radiation over the thermal/ram pressure ratio time evolution is calculated explicitely and the crucial role of the cluster mechanical power and of the strong time evolution of the ionizing photon flux and of the bolometric luminosity of the exciting cluster is stressed. It is shown that radiation has only a narrow window of opportunity to dominate the wind-driven shell dynamics. This may occur only at early stages of the bubble evolution and if the shell expands into a dusty and/or a very dense proto-cluster medium. The impact of radiation pressure on the wind-driven shell becomes always negligible after about 3 Myr. Finally, the wind-driven model results allow one to compare the model predictions with the distribution of thermal pressure derived from X-ray observations. The shape of the thermal pressure profile allows then to distinguish between the energy and the momentum dominated regimes of expansion and thus conclude whether radiative losses of energy or the leakage of hot gas from the bubble interior have been significant during the bubble evolution.

We present the sensitivity of HAWC to Gamma Ray Bursts (GRBs). HAWC is a very high-energy gamma-ray observatory currently under construction in Mexico at an altitude of 4100 m. It will observe atmospheric air showers via the water Cherenkov method. HAWC will consist of 300 large water tanks instrumented with 4 photomultipliers each. HAWC has two data acquisition (DAQ) systems. The main DAQ system reads out coincident signals in the tanks and reconstructs the direction and energy of individual atmospheric showers. The scaler DAQ counts the hits in each photomultiplier tube (PMT) in the detector and searches for a statistical excess over the noise of all PMTs. We show that HAWC has a realistic opportunity to observe the high-energy power law components of GRBs that extend at least up to 30 GeV, as it has been observed by Fermi LAT. The two DAQ systems have an energy threshold that is low enough to observe events similar to GRB 090510 and GRB 090902b with the characteristics observed by Fermi LAT. HAWC will provide information about the high-energy spectra of GRBs which in turn could help to understanding about e-pair attenuation in GRB jets, extragalactic background light absorption, as well as establishing the highest energy to which GRBs accelerate particles.

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.

We study the evolution of Super Star Cluster (SSC) winds driven by stellar winds and supernova (SN) explosions. Time-dependent rates at which mass and energy are deposited into the cluster volume, as well as the time-dependent chemical composition of the re-inserted gas, are obtained from the population synthesis code Starburst99. These results are used as input for a semi-analytic code which determines the hydrodynamic properties of the cluster wind as a function of cluster age. Two types of winds are detected in the calculations. For the quasi-adiabatic solution, all of the inserted gas leaves the cluster in the form of a stationary wind. For the bimodal solution, some of the inserted gas becomes thermally unstable and forms dense warm clumps which accumulate inside the cluster. We calculate the evolution of the wind velocity and energy flux and integrate the amount of accumulated mass for clusters of different mass, radius and initial metallicity. We consider also conditions with low heating efficiency of the re-inserted gas or mass loading of the hot thermalized plasma with the gas left over from star formation. We find that the bimodal regime and the related mass accumulation occur if at least one of the two conditions above is fulfilled.

We present 1D numerical simulations, which consider the effects of radiative cooling and gravity on the hydrodynamics of the matter reinserted by stellar winds and supernovae within young nuclear starbursts with a central supermassive black hole (SMBH). The simulations confirm our previous semi-analytic results for low energetic starbursts, evolving in a quasi-adiabatic regime, and extend them to more powerful starbursts evolving in the catastrophic cooling regime. The simulations show a bimodal hydrodynamic solution in all cases. They present a quasi-stationary accretion flow onto the black hole, defined by the matter reinserted by massive stars within the stagnation volume and a stationary starburst wind, driven by the high thermal pressure acquired in the region between the stagnation and the starburst radii. In the catastrophic cooling regime, the stagnation radius rapidly approaches the surface of the starburst region, as one considers more massive starbursts. This leads to larger accretion rates onto the SMBH and concurrently to powerful winds able to inhibit interstellar matter from approaching the nuclear starburst. Our self-consistent model thus establishes a direct physical link between the SMBH accretion rate and the nuclear star formation activity of the host galaxy and provides a good upper limit to the accretion rate onto the central black hole.

We present a detailed study of the hydrodynamics of the matter reinserted by massive stars via stellar winds and supernovae explosions in young assembling galaxies. We show that the interplay between the thermalization of the kinetic energy provided by massive stars, radiative cooling of the thermalized plasma and the gravitational pull of the host galaxy, lead to three different hydrodynamic regimes. These are: a) The quasi-adiabatic supergalactic winds. b) The bimodal flows, with mass accumulation in the central zones and gas expulsion from the outer zones of the assembling galaxy. c) The gravitationally bound regime, for which all of the gas returned by massive stars remains bound to the host galaxy and is likely to be reprocessed into futher generations of stars. Which of the three possible solutions takes place, depends on the mass of the star forming region its mechanical luminosity (or star formation rate) and its size. The model predicts that massive assembling galaxies with large star formation rates similar to those detected in SCUBA sources ($\sim 1000$ M$_\odot$ yr$^{-1}$) are likely to evolve in a positive star-formation feedback condition, either in the bimodal, or in the gravitationally bound regime. This implies that star formation in these sources may have little impact on the intergalactic medium and result instead into a fast interstellar matter enrichment, as observed in high redshift quasars.

The origin of supersonic infrared and radio recombination nebular lines often detected in young and massive superstar clusters are discussed. We suggest that these arise from a collection of repressurizing shocks (RSs), acting effectively to re-establish pressure balance within the cluster volume and from the cluster wind which leads to an even broader although much weaker component. The supersonic lines are here shown to occur in clusters that undergo a bimodal hydrodynamic solution (Tenorio-Tagle et al. 2007), that is within clusters that are above the threshold line in the mechanical luminosity or cluster mass vs the size of the cluster (Silich et al. 2004). The plethora of repressurizing shocks is due to frequent and recurrent thermal instabilities that take place within the matter reinserted by stellar winds and supernovae. We show that the maximum speed of the RSs and of the cluster wind, are both functions of the temperature reached at the stagnation radius. This temperature depends only on the cluster heating efficiency ($\eta$). Based on our two dimensional simulations (Wunsch et al. 2008) we calculate the line profiles that result from several models and confirm our analytical predictions. From a comparison between the predicted and observed values of the half-width zero intensity of the two line components we conclude that the thermalization efficiency in SSC's above the threshold line must be lower than 20%.

Here we discuss the mechanical feedback that massive stellar clusters provide to the interstellar medium of their host galaxy. We apply an analytic theory developed in a previous study for M82-A1 to a sample of 10 clusters located in the central zone of the starburst galaxy M82, all surrounded by compact and dense HII regions. We claim that the only way that such HII regions can survive around the selected clusters, is if they are embedded into a high pressure ISM and if the majority of their mechanical energy is lost within the star cluster volume via strong radiative cooling. The latter implies that these clusters have a low heating efficiency, $\eta$, and evolve in the bimodal hydrodynamic regime. In this regime the shock-heated plasma in the central zones of a cluster becomes thermally unstable, loses its pressure and is accumulated there, whereas the matter injected by supernovae and stellar winds outside of this volume forms a high velocity outflow - the star cluster wind. We calculated the heating efficiency for each of the selected clusters and found that in all cases it does not exceed 10% . Such low heating efficiency values imply a low mechanical energy output and the impact that the selected clusters provide to the ISM of M82 is thus much smaller than what one would expect using stellar cluster synthetic models.

In this work, we present a comprehensive X-ray picture of the interaction between a super star cluster and the ISM. In order to do that, we compare and combine the X-ray emission from the superwind driven by the cluster with the emission from the wind-blown bubble. Detailed analytical models for the hydrodynamics and X-ray luminosity of fast polytropic superwinds are presented. The superwind X-ray luminosity models are an extension of the results obtained in Paper I of this series. Here, the superwind polytropic character allows to parameterize a wide variety of effects, for instance, radiative cooling. Additionally, X-ray properties that are valid for all bubble models taking thermal evaporation into account are derived. The final X-ray picture is obtained by calculating analytically the expected surface brightness and weighted temperature of each component. All of our X-ray models have an explicit dependence on metallicity and admit general emissivities as functions of the hydrodynamical variables. We consider a realistic X-ray emissivity that separates the contributions from hydrogen and metals. The paper ends with a comparison of the models with observational data.

The mounting evidence for violent nuclear star formation in Seyfert galaxies has led us to consider the hydrodynamics of the matter reinserted by massive stars through strong stellar winds and supernovae, under the presence of a central massive BH. We show that in all cases there is a bimodal solution strongly weighted by the location of the stagnation radius (Rst), which splits the star cluster into two different zones. Matter reinserted within the stagnation volume is to be accreted by the BH while its outer counterpart would composed a star cluster wind. The mechanical power of the latter, ensures that there is no accretion of the ISM into the BH and thus the BH accretion and its luminosity is regulated by the star formation feedback. The location of the stagnation radius is a function of three parameters: the BH mass, the mechanical power (or mass) of the star formation event and the size of the star forming region. Here we present our self-consistent, stationary solution, discuss the accretion rates and BH luminosities and show that our model predicts the intrinsic link between the BH activity and the starburst parameters.

The mass reinserted by young stars of an emerging massive compact cluster shows a bimodal hydrodynamic behaviour. In the inner part of the cluster, it is thermally unstable, while in its outer parts it forms an out-blowing wind. The chemical homogeneity/inhomogeneity of low/high mass clusters demonstrates the relevance of this solution to the presence of single/multiple stellar populations. We show the consequences that the thermal instability of the reinserted mass has to the galactic super-winds.

We present a self-consistent, bimodal stationary solution for spherically symmetric flows driven by young massive stellar clusters with a central supermassive black hole. We demonstrate that the hydrodynamic regime of the flow depends on the location of the cluster in the 3D (star cluster mechanical luminosity - BH mass - star cluster radius) parameter space. We show that a threshold mechanical luminosity (L_crit) separates clusters which evolve in the BH dominated regime frome those whose internal structure is strongly affected by the radiative cooling. In the first case(below the threshold energy) gravity of the BH separates the flow into two distinct zones: the inner accretion zone and the outer zone where the star cluster wind is formed. In the second case (above the critical luminosity), catastrophic cooling sets in inside the cluster. In this case the injected plasma becomes thermally unstable that inhibits a complete stationary solution. We compared the calculated accretion rates and the BH luminosities with those predicted by the classic Bondi accretion theory and found that Bondi's theory is in good agreement with our results in the case of low mass clusters. However, it substantially underestimates the accretion rates and BH luminosities if the star cluster mechanical luminosity, L_sc, approaches the threshold value (L_sc > 0.1 L_crit).

Using the hydrodynamic code ZEUS, we perform 2D simulations to determine the fate of the gas ejected by massive stars within super star clusters. It turns out that the outcome depends mainly on the mass and radius of the cluster. In the case of less massive clusters, a hot high velocity ($\sim 1000$ km s$^{-1}$) stationary wind develops and the metals injected by supernovae are dispersed to large distances from the cluster. On the other hand, the density of the thermalized ejecta within massive and compact clusters is sufficiently large as to immediately provoke the onset of thermal instabilities. These deplete, particularly in the central densest regions, the pressure and the pressure gradient required to establish a stationary wind, and instead the thermally unstable parcels of gas are rapidly compressed, by a plethora of re-pressurizing shocks, into compact high density condensations. Most of these are unable to leave the cluster volume and thus accumulate to eventually feed further generations of star formation. The simulations cover an important fraction of the parameter-space, which allows us to estimate the fraction of the reinserted gas which accumulates within the cluster and the fraction that leaves the cluster as a function of the cluster mechanical luminosity, the cluster size and heating efficiency.

Here we address the question whether the ionized shells associated with giant HII regions can be progenitors of the larger HI shell-like objects found in the Milky Way and other spiral and dwarf irregular galaxies. We use for our analysis a sample of 12 HII shells presented recently by Rela\~no et al. (2005, 2007). We calculate the evolutionary tracks that these shells would have if their expansion is driven by multiple supernovae explosions from the parental stellar clusters. We find, contrary to Rela\~no et al. (2007), that the evolutionary tracks of their sample HII shells are inconsistent with the observed parameters of the largest and most massive neutral hydrogen supershells. We conclude that HII shells found inside giant HII regions may represent the progenitors of small or intermediate HI shells, however they cannot evolve into the largest HI objects unless, aside from the multiple supernovae explosions, an additional energy source contributes to their expansion.

The observed parameters of the young superstar cluster M82-A1 and its associated compact HII region are here shown to indicate a low heating efficiency or immediate loss, through radiative cooling, of a large fraction of the energy inserted by stellar winds and supernovae during the early evolution of the cluster. This implies a bimodal hydrodynamic solution which leads to a reduced mass deposition rate into the ISM, with a much reduced outflow velocity. Furthermore, to match the observed parameters of the HII region associated to M82-A1, the resultant star cluster wind is here shown to ought to be confined by a high pressure interstellar medium. The cluster wind parameters, as well as the location of the reverse shock, its cooling length and the radius of the standing outer HII region are derived analytically. All of these properties are then confirmed with a semi-analytical integration of the flow equations, which provides us also with the run of the hydrodynamic variables as a function of radius. The impact of the results is discussed and extended to other massive and young superstar clusters surrounded by a compact HII region.

We look for a simple analytic model to distinguish between stellar clusters undergoing a bimodal hydrodynamic solution from those able to drive only a stationary wind. Clusters in the bimodal regime undergo strong radiative cooling within their densest inner regions, which results in the accumulation of the matter injected by supernovae and stellar winds and eventually in the formation of further stellar generations, while their outer regions sustain a stationary wind. The analytic formulae are derived from the basic hydrodynamic equations. Our main assumption, that the density at the star cluster surface scales almost linearly with that at the stagnation radius, is based on results from semi-analytic and full numerical calculations. The analytic formulation allows for the determination of the threshold mechanical luminosity that separates clusters evolving in either of the two solutions. It is possible to fix the stagnation radius by simple analytic expressions and thus to determine the fractions of the deposited matter that clusters evolving in the bimodal regime blow out as a wind or recycle into further stellar generations.

We present semi-analytical and numerical models, accounting for the impact of radiative cooling on the hydrodynamics of the matter reinserted as strong stellar winds and supernovae within the volume occupied by young, massive and compact superstellar clusters. First of all we corroborate the location of the threshold line in the mechanical energy input rate vs the cluster size plane, found by Silich et al. (2004). Such a line separates clusters able to drive a quasi-adiabatic or a strongly radiative wind from clusters in which catastrophic cooling occurs within the star cluster volume. Then we show that the latter, clusters above the threshold line, undergo a bimodal behavior in which the central densest zones cool rapidly and accumulate the injected matter to eventually feed further generations of star formation, while the outer zones are still able to drive a stationary wind. The results are presented into a series of universal dimensionless diagrams from which one can infer: the size of the two zones, the fraction of the deposited mass that goes into each of them and the luminosity of the resultant winds, for clusters of all sizes and energy input rates, regardless the assumed adiabatic terminal speed V_A.

Different hydrodynamic regimes for the gaseous outflows generated by multiple supernovae explosions and stellar winds occurring within compact and massive star clusters are discussed. It is shown that there exists the threshold energy that separates clusters whose outflows evolve in the quasi-adiabatic or radiative regime from those within which catastrophic cooling and a positive feedback star-forming mode sets in. The role of the surrounding ISM and the observational appearance of the star cluster winds evolving in different hydrodynamic regimes are also discussed.

We review the structural properties of giant extragalactic HII regions and HII galaxies based on two dimensional hydrodynamic calculations, and propose an evolutionary sequence that accounts for their observed detailed structure. The model assumes a massive and young stellar cluster surrounded by a large collection of clouds. These are thus exposed to the most important star-formation feedback mechanisms: photoionization and the cluster wind. The models show how the two feedback mechanisms compete with each other in the disruption of clouds and lead to two different hydrodynamic solutions: The storage of clouds into a long lasting ragged shell that inhibits the expansion of the thermalized wind, and the steady filtering of the shocked wind gas through channels carved within the cloud stratum that results into the creation of large-scale superbubbles. Both solutions are here claimed to be concurrently at work in giant HII regions and HII galaxies, causing their detailed inner structure. A full description of the calculations can be found in The Astrophysical Journal 643, 186. Animated version of the models presented can be found at http://www.iaa.csic.es/\~{}eperez/ssc/ssc.html.

We review the structural properties of giant extragalactic HII regions and HII galaxies based on 2D hydrodynamic calculations, and propose an evolutionary sequence that accounts for their observed detailed structure. The model assumes a massive and young stellar cluster surrounded by a large collection of clouds. These are thus exposed to the most important star-formation feedback mechanisms: photoionization and the cluster wind. The models show how the two feedback mechanisms compete in the disruption of clouds and lead to two different hydrodynamic solutions: The storage of clouds into a long lasting ragged shell that inhibits the expansion of the thermalized wind, and the steady filtering of the shocked wind gas through channels carved within the cloud stratum. Both solutions are claimed to be concurrently at work in giant HII regions and HII galaxies, causing their detailed inner structure. This includes multiple large-scale shells, filled with an X-ray emitting gas, that evolve to finally merge with each other, giving the appearance of shells within shells. The models also show how the inner filamentary structure of the giant superbubbles is largely enhanced with matter ablated from clouds and how cloud ablation proceeds within the original cloud stratum. The calculations point at the initial contrast density between the cloud and the intercloud media as the factor that defines which of the two feedback mechanisms becomes dominant throughout the evolution. Animated version of the models can be found at http://www.iaa.csic.es/\~{}eperez/ssc/ssc.html.

We study the neutral hydrogen distribution and kinematics in the Local Group dwarf irregular galaxy IC 1613 and compare them with the ionized gas distribution and stellar content of the galaxy. We discuss several mechanisms which may be responsible for the origin of the observed complicated HI morphology and compare parameters of the most prominent kpc-scale HI structure with the multiple SNe scenario. The observations were performed with the Vary Large Array of NRAO with a linear resolution ~23 pc at the adopted distance of 725 kpc and the spectral channel width of 2.57 km/s. The numerical calculations have been provided with our 2.5D Lagrangian scheme based on the thin layer approximation. We found that the ISM of the galaxy is highly inhomogeneous and identified a number of intermediate-scale (200 pc - 300 pc in diameter) HI arcs and shells having expansion velocities of 10 to 20 km/s. Besides these shells, several giant holes and arc-shaped structures have been revealed, whose radii exceed several hundred parsecs. We found that parameters of the most prominent (M_HI = 2.8 x 10^7 Msol) kpc-scale structure and the level of the detected star formation activity are inconsistent with the multiple SNe hypothesis.

The X-ray emission properties from the hot thermalized plasma that results from the collisions of individual stellar winds and supernovae ejecta within rich and compact star clusters are discussed. We propose a simple analytical way of estimating the X-ray emission generated by super star clusters and derive an expression that indicates how this X-ray emission depends on the main cluster parameters. Our model predicts that the X-ray luminosity from the star cluster region is highly dependent on the star cluster wind terminal speed, a quantity related to the temperature of the thermalized ejecta.We have also compared the X-ray luminosity from the SSC plasma with the luminosity of the interstellar bubbles generated from the mechanical interaction of the high velocity star cluster winds with the ISM.We found that the hard (2.0 keV - 8.0 keV) X-ray emission is usually dominated by the hotter SSC plasma whereas the soft (0.3 keV - 2.0 keV) component is dominated by the bubble plasma. This implies that compact and massive star clusters should be detected as point-like hard X-ray sources embedded into extended regions of soft diffuse X-ray emission. We also compared our results with predictions from the population synthesis models that take into consideration binary systems and found that in the case of young,massive and compact super star clusters the X-ray emission from the thermalized star cluster plasma may be comparable or even larger than that expected from the HMXB population.

The force of gravity acting within the volume occupied by young, compact and massive superstar clusters, is here shown to drive in situ all the matter deposited by winds and supernovae into several generations of star formation. These events are promoted by radiative cooling which drains the thermal energy of the ejected gas causing its accumulation to then rapidly exceed the gravitational instability criterion. A detailed account of the integrated ionizing radiation and mechanical luminosity as a function of time is here shown to lead to a new stationary solution. In this, the mass deposition rate $\dot M$, instead of causing a wind as in the adiabatic solution, turns into a positive feedback star-forming mode equal to the star formation rate. Some of the implications of this extreme positive feedback mode are discussed.

We review the properties of young superstellar clusters and the impact that their evolution has in their host galaxies. In particular we look at the two different star-forming feedback modes: positive and negative feedback. The development of strong isotropic winds emanating from massive clusters, capable of disrupting the remains of the parental cloud as well as causing the large-scale restructuring of the surrounding ISM, has usually been taken as a negative feedback agent. Here we show the impact that radiative cooling has on the resultant outflows and then, as an extreme example, we infer from the observations of M82 the detailed inner structure of supergalactic winds and define through numerical simulations the ingredients required to match such structures.We also show how when radiative cooling becomes significant within the star cluster volume itself (~ 30% of the deposited energy), the force of gravity takes over and drives in situ all the matter deposited by winds and supernovae into several generations of star formation. A situation in which the mass deposition rate from the star cluster instead of causing a wind as in the adiabatic solution, turns into a positive feedback star-forming mode equal to the star formation rate.

The properties of star cluster winds in the supercritical, catastrophic cooling regime are discussed. We demonstrate that strong radiative cooling may inhibit superwinds and, after a rapid phase of accumulation of the ejected material within the star-forming volume, a new stationary isothermal regime, supported by the ionizing radiation from the central cluster, is established. The expected appearance of this core/halo supernebula in the visible line regime and possible late evolutionary tracks for super-star cluster winds, in the absence of ionizing radiation, are thoroughly discussed.

Here we present a self-consistent stationary solution for spherically symmetric winds driven by massive star clusters under the impact of radiative cooling. We demonstrate that cooling may modify drastically the distribution of temperature if the rate of injected energy approaches a critical value. We also prove that the stationary wind solution does not exist whenever the energy radiated away at the star cluster center exceeds ~ 30% of the energy deposition rate. Finally we thoroughly discuss the expected appearance of super-star cluster winds in the X-ray and visible line regimes. The three solutions here found: the quasi-adiabatic, the strongly radiative wind and the inhibited stationary solution, are then compared to the winds from Arches cluster, NGC 4303 central cluster and to the supernebula in NGC 5253.

Using a novel approach we have reanalized the question of whether the extreme star forming galaxies known as HII galaxies are truly young or rejuvenated old systems. We first present a method of inversion that applies to any monotonic function of time describing the evolution of independent events. We show that, apart from a normalization constant, the ``true'' time dependence can be recovered from the inversion of its probability density function. We applied the inversion method to the observed equivalent width of Hbeta (EW(Hbeta)) distribution for objects in the Terlevich and collaborators Spectrophotometric Catalogue of HII galaxies and found that their global history of star formation behaves much closer to the expectations of a continuos star formation model than to an instantaneous one. On the other hand, when the inversion method is applied to samples within a restricted metallicity range we find that their history of star formation behaves much closer to what the instantaneous model predicts. Our main conclusion is that, globally, the evolution of HII galaxies seems consistent with a succession of short starbursts separated by quiescent periods and that, while the emission lines trace the properties of the present burst, the underlying stellar continuum traces the whole star formation history of the galaxy. Thus, observables like the EW(Hbeta) that combine an emission line flux, i.e. a parameter pertaining to the present burst, with the continuum flux, i.e. a parameter that traces the whole history of star formation, should not be used alone to characterize the present burst.

We report the discovery of a double ring of emission-line regions around the nucleus of the Blue Compact Dwarf (BCD) galaxy Mrk409 as seen by deep, ground-based Halpha images. Echelle spectroscopy obtained at Magellan-I with MIKE shows the presence of ionized gas flowing out of the galaxy from a very massive (>7E6 Msun) nuclear starburst with projected expansion velocities of \~50 km/s. Different scenarios for the formation of these rings are discussed. While the innermost, nuclear ring is most probably formed by the interaction of a starburst-driven shock with the surrounding interstellar medium, the origin of the outer ring is less clear.

Here we model a star forming factory in which the continuous creation of stars results in a highly concentrated, massive (globular cluster-like) stellar system. We show that under very general conditions a large-scale gravitational instability in the ISM, which triggers the collapse of a massive cloud, leads with the aid of a spontaneous first generation of massive stars, to a standing, small-radius, cold and dense shell. Eventually, as more of the collapsing matter is processed and incorporated, the shell becomes gravitationally unstable and begins to fragment, allowing the formation of new stars, while keeping its location. This is due to a detailed balance established between the ram pressure from the collapsing cloud which, together with the gravitational force exerted on the shell by the forming cluster, acts against the mechanical energy deposited by the collection of new stars. We analyze the mass spectrum of fragments that result from the continuous fragmentation of the standing shell and show that its shape is well approximated at the high mass end by a power law with slope -2.25, very close to the value that fits the universal IMF. The self-contamination resultant from the continuous generation of stars is shown to lead to a large metal spread in massive ($\sim$ 10$^6$ M$_\odot$) clusters, while clusters with a mass similar to 10$^5$ M$_\odot$ or smaller, simply reflect the initial metalicity of the collapsing cloud. This is in good agreement with the data available for globular clusters in the Galaxy.

We present two dimensional hydrodynamic calculations of free expanding supergalactic winds, taking into consideration strong radiative cooling. Our main premise is that supergalactic winds are powered by collections of superstar clusters. Every individual superstar cluster is a source of a high metallicity radiative supersonic outflow (paper I, 2003, ApJ, 590, 791). The interaction of winds from neighboring knots of star formation is shown to lead to a collection of stationary oblique shocks and crossing shocks, able to structure the general outflow into a network of dense and cold, kpc long filaments that originate near the base of the outflow. The shocks also lead to extended regions of diffuse soft X-ray emission and furthermore, to channel the outflow with a high degree of collimation into the intergalactic medium.

Here we present semi-analytical solutions of the hydrodynamics of compact supernova remnants (cSNRs) and their chemical evolution. Two models of $15 M_{\odot}$ and $25 M_{\odot}$ supernova explosions evolving in a high density medium ($n_0 \sim 10^7 cm^{-3}$) are throughly discussed.

The steady state numerical solution and two dimensional hydrodynamic calculations of supergalactic winds generated by super star clusters are presented. Two possible regimes: quasi-adiabatic and strongly radiative are thoroughly discussed. It is shown that superwinds driven by compact and powerful starbursts undergo catastrophic cooling close to their sources, and establish a temperature distribution radically different from those obtained in the adiabatic solution. At the same time, both velocity and density radial distributions remain almost unchanged given the speed of the rapidly diverging flow. The rapid fall in temperature as a function of radius reduces the size of the zone radiating in X-rays and decreases the superwind X-ray luminosity. At the same time, cooling brings the 10$^4$ K boundary closer to the wind center, and promotes the establishment of an extended ionized fast moving envelope. This should show up as a weak and broad ($\sim 1000$ km s$^{-1}$) line emission component at the base of the much narrower line caused by the central HII region.

We stress some of the major differences between supergalactic winds and giant superbubbles evolving into the giant low density haloes of galaxies. We review major properties of bursts of star formation and the critical energy required for mass ejection both in the case of an ISM strongly flattened by rotation and that imposed by a more extended ISM distribution arising from a smaller rotation. The limits are thus establish for galaxies with an ISM mass in the range 10$^6$ M$_\odot$ to more than 10$^9$ M$_\odot$, and are compared with a sample of local galaxies. Some of these galaxies seem to be above the critical limit despite the fact that their structure is clearly that of a superbubble. True supergalactic winds, as evidence by M82, are shown to exceed the critical limit by more than an order of magnitude and thus the limit derived by Silich & Tenorio-Tagle (2001) for mass ejection should be regarded as a lower limit.

The hydrodynamics and intrinsic properties of galactic-scale gaseous outflows generated by violent starbursts are thoroughly discussed, taking into account the hot gas chemical evolution and radiative cooling. It is shown that the metallicity of superbubbles vary with time and can easily exceed the solar value. Galactic superwinds driven by compact and powerful starbursts undergo catastrophic cooling and establish a temperature distribution radically different to that predicted by adiabatic solution. The observational properties of supergalactic winds in X-rays and visible line regimes derived from our calculations are compared with the standard adiabatic model predictions.

We have mapped the 12CO J=1-0 and J=2-1 line emission in Mrk86, one of the most metal-deficient Blue Compact Dwarf galaxies so far detected in 12CO. The 12CO emission is distributed in a horseshoe-like structure that follows the locus of the most recent star formation regions. The minimum in molecular-line emission corresponds to the position of an older, massive nuclear starburst. The H2 mass of the galaxy (in the range 0.4-5x10^7 Msun) and its morphology have been compared with the predictions of hydrodynamic simulations of the evolution of the interstellar medium surrounding a nuclear starburst. These simulations suggest that the physical conditions in the gas swept out by the starburst could have led to the formation of the ring of molecular gas reported here. This result provides an attractive scenario for explaining the propagation (in a galactic scale) of the star formation in dwarf galaxies.

The recent history of star formation in the BCD galaxy VII Zw403, is analysed with accounts for the dynamics of the remnant generated either by an instantaneous burst or by a continuous star formation event. We have re-observed VII Zw403 with a better sensitivity corresponding to the threshold H$_{\alpha}$ flux $8.15 \times 10^{-17}$ erg cm$^{-2}$ s$^{-1}$. The total H$_{\alpha}$ luminosity derived from our data is much larger than reported before, and presents a variety of ionized filaments and incomplete shells superimposed on the diffuse H$_{\alpha}$ emission. This result has a profound impact on the predicted properties of the starburst blown superbubble. Numerical calculations based on the previous H$_{\alpha}$ data, predict two different scenarios of star formation able to match simultaneously all observed parameters. These are an instantaneous burst of star formation with a total mass of $5 \times 10^5$ M$_{\odot}$ and a star-forming event with a constant SFR = $4 \times 10^{-3}$ M$_{\odot}$ yr$^{-1}$, which lasts for 35 Myrs. The numerical calculations based on the energy input rate derived from our observations predict a short episode of star formation lasting less than 10 Myrs with a total star cluster mass $\sim (1 - 3) \times 10^6$ M$_{\odot}$.

We construct three extreme different scenarios of the star formation histories applicable to a sample of dwarf galaxies, based either on their present metallicity or their luminosity. The three possible scenarios imply different mechanical energy input rates and these we compare with the theoretical lower limits established for the ejection of processed matter out of dwarf galaxies. The comparison strongly points at the existence of extended gaseous haloes in these galaxies, acting as the barrier that allows galaxies to retain their metals and enhance their abundance. At the same time our findings strongly point at a continuous star-forming process, rather than to coeval bursts, as the main contributors to the overall metallicity in our galaxy sample.

The recent history of star formation in VIIZw 403 is analyzed taking into account the dynamics of the starburst blown superbubble and the restrictions that follow from X-ray observations, and from our H$\alpha$ data. Our results show that the starburst energizing VIIZw 403 is NOT coeval, but rather the star forming phase has lasted for more than 30Myr and at such a low pace that most probably the newly processed metals will not be ejected into the intergalactic medium.

We evaluate the minimum energy input rate that starbursts require for expelling their newly processed matter from their host galaxies. Special attention is given to the pressure caused by the environment in which a galaxy is situated, as well as to the intrinsic rotation of the gaseous component. We account for these factors and for a massive dark matter distribution, and develop a self-consistent solution for the interstellar matter gas distribution. Our results are in excellent agreement with the results of Mac Low & Ferrara (1999) for galaxies with a flattened disk-like ISM density distribution and a low intergalactic gas pressure ($P_{IGM}/k$ $\leq $ 1 cm$^{-3}$ K). However, our solution also requires a much larger energy input rate threshold when one takes into consideration both a larger intergalactic pressure and the possible existence of a low-density, non-rotating, extended gaseous halo component.

We calculate the time-dependent metal production expected from starbursts and use them as boundary conditions in our 2D simulations of evolving superbubbles. We assume that the produced metals (oxygen and iron) thoroughly mix with the ejected stellar envelopes, and/or with the matter thermally evaporated from the superbubble cold outer shell. The metal production process determines the time-dependent metallicity in hot superbubble interiors, and leads to values of Z greater or equal than solar, when oxygen is used as tracer, and under-solar when the metallicity is measured with respect to iron. In either case, the enhanced metallicity boosts the X-ray emissivity of superbubbles, bringing theory and observations closer together.

The detection of Ly alpha emission in star-forming galaxies in different shapes and intensities (always smaller than predicted for case B recombination) has puzzled the astronomical community for more than a decade. Here we use two dimensional calculations to follow the evolution of superbubbles and of the H II regions generated by the output of UV photons from massive stars. We show the impact caused by massive star formation in the ISM of different galaxies and we look at the conditions required to detect Ly alpha emission from a nuclear H II region, and the variety of profiles that may be expected as a function of time.

Here we present a possible solution to the apparent discrepancy between the observed properties of LMC bubbles and the standard, constant density bubble model. A two-dimensional model of a wind-driven bubble expanding from a flattened giant molecular cloud is examined. We conclude that the expansion velocities derived from spherically symmetric models are not always applicable to elongated young bubbles seen almost face-on due to the LMC orientation. In addition, an observational test to differentiate between spherical and elongated bubbles seen face-on is discussed.

Two dimensional calculations of the evolution of remnants generated by the strong mechanical energy deposited by stellar clusters in dwarf galaxies (M \si $10^9 - 10^{10}$ \msun) are presented. The evolution is followed for times longer than both the blowout time and the presumed span of energy injection generated by a coeval massive stellar cluster. The remnants are shown to end up wrapping around the central region of the host galaxy, while growing to kpc-scale dimensions. Properties of the remnants such as luminosity, size, swept up mass, and expansion speed are given as a function of time for all calculated cases. The final fate of the swept-up galactic gas and of the matter processed by the central starburst is shown to be highly-dependent on the properties of the low density galactic halo. Superbubbles powered by star clusters, with properties similar to those inferred from the observations, slow down in the presence of an extended halo to expansion speeds smaller than the host galaxy escape velocity. Values of the critical luminosity required for the superbubbles to reach the edge of the galaxies with a speed comparable to the escape speed are derived analytically and numerically. The critical luminosities are larger than those in the detected sources and thus, the superbubbles in amorphous dwarf galaxies must have already undergone blowout and are presently evolving into an extended low density halo. This will inhibit the loss of the swept-up and processed matter from the galaxy.

Equations are derived which describe a propagation of strong shocks in the interstellar matter, without any demands for a symmetry, in a thin layer approximation (2.5 dimensions). Using these equations permits to calculate a propagation of shock waves from nonsymmetric supernovae explosions in the medium with arbitrary density distribuion, formation of superbubbles in galaxies.