We present new high spectral resolution observations of 15 high-z ($1.3 \leq$ z $\leq 2.5$) HII Galaxies (HIIG) obtained with MOSFIRE at the Keck Observatory. These data, combined with already published data for another 31 high-z and 107 z $\leq 0.15$ HIIG, are used to obtain new independent cosmological results using the distance estimator based on the established correlation between the Balmer emission line velocity dispersion and luminosity for HIIG. Our results are in excellent agreement with the latest cosmological concordance model ($\Lambda$CDM) published results. From our analysis, we find a value for the mass density parameter of $\Omega_m=0.290^{+0.056}_{-0.069}$ (stat). For a flat Universe we constrain the plane $\lbrace\Omega_m;w_0\rbrace = \lbrace 0.280^{+0.130}_{-0.100} ; -1.12^{+0.58}_{-0.32}\rbrace $ (stat). The joint likelihood analysis of HIIG with other complementary cosmic probes (Cosmic Microwave Background and Baryon Acoustic Oscillations) provides tighter constraints for the parameter space of the Equation of State of Dark Energy that are also in excellent agreement with those of similar analyses using Type Ia Supernovae instead as the geometrical probe.

Massive starforming regions like Giant HII Regions (GHIIR) and HII Galaxies (HIIG) are emission line systems ionized by compact young massive star clusters (YMC) with masses ranging from $10^4$M$_\odot$ to $10^8$M$_\odot$. We model the photometric and dynamical evolution over a Hubble time of the massive gravitationally bound systems that populate the tight relation between absolute blue magnitude and velocity dispersion ($M_{B}-\sigma$) of GHIIR and HIIG and compare the resulting relation with that one of old stellar systems: globular clusters, elliptical galaxies, bulges of spirals. After 12~Gyr of evolution their position on the $\sigma$ vs. M$_B$ plane coincides -- depending on the initial mass -- either with the globular clusters for systems with initial mass $M < 10^6$M$_\odot$ or with a continuation of the ellipticals, bulges of spirals and ultracompact dwarfs for YMC with $M >10^6$M$_\odot$. The slope change in the $M_{B}-\sigma$ and $M_B$-size relations at cluster masses around $10^6$M$_\odot$ is due to the larger impact of the dynamical evolution on the lower mass clusters. We interpret our result as an indication that the YMC that ionize GHIIR and HIIG can evolve to form globular clusters and ultra compact dwarf ellipticals in about 12 Gyr so that present day globular clusters and ultra compact dwarf ellipticals may have formed in conditions similar to those observed in today GHIIR and HIIG.

We use the HII galaxies $L - \sigma$ relation and the resulting Hubble expansion cosmological probe of a sample of just 25 high-$z$ (up to $z \sim 2.33$) HII galaxies, in a joint likelihood analysis with other well tested cosmological probes (CMB, BAOs) in an attempt to constrain the dark energy equation of state (EoS). The constraints, although still weak, are in excellent agreement with those of a similar joint analysis using the well established SNIa Hubble expansion probe. Interestingly, even with the current small number of available high redshift HII galaxies, the HII/BAO/CMB joint analysis gives a 13% improvement of the quintessence dark energy cosmological constraints compared to the BAO/CMB joint analysis. We have further performed extensive Monte Carlo simulations, with a realistic redshift sampling, to explore the extent to which the use of the $L - \sigma$ relation, observed in HII galaxies, can constrain effectively the parameter space of the dark energy EoS. The simulations predict substantial improvement in the constraints when increasing the sample of high-$z$ HII galaxies to 500, a goal that can be achieved in reasonable observing times with existing large telescopes and state-of-the-art instrumentation.