The first stars are assumed to be predominantly massive. Although, due to the low initial abundances of heavy elements the line-driven stellar winds are supposed to be inefficient in the first stars, these stars may loose a significant amount of their initial mass by other mechanisms. In this work, we study the evolution with a prescribed mass loss rate of very massive, galactic and pregalactic, Population III stars, with initial metallicities $Z=10^{-6}$ and $Z=10^{-9}$, respectively, and initial masses 100, 120, 150, 200, and 250$\,M_{\odot}$ during the hydrogen and helium burning phases. The evolution of these stars depends on their initial mass, metallicity and the mass loss rate. Low metallicity stars are hotter, compact and luminous, and they are shifted to the blue upper part in the Hertzprung-Russell diagram. With mass loss these stars provide an efficient mixing of nucleosynthetic products, and depending on the He-core mass their final fate could be either pair-instability supernovae or energetic hypernovae. These stars contributed to the reionization of the universe and its enrichment with heavy elements, which influences the subsequent star formation properties.

We examine the problem of the collapse and fragmentation of molecular clouds with a Gaussian density distribution with high resolution, double precision numerical simulations using the GADGET-2 code. To describe the thermodynamic properties of the cloud during the collapse -to mimic the rise of temperature predicted by radiative transfer- we use a barotropic equation of state that introduces a critical density to separate the isothermal and adiabatic regimes. We discuss the effects of this critical density in the formation of multiple systems. We confirm the tendency found for Plummer and Gaussian models that if the collapse changes from isothermal to adiabatic at earlier times that occurs for the models with a lower critical density, the collapse is slowed down, and this enhances the fragments' change to survive. However, this effect happens up to a threshold density below which single systems tend to form. On the other hand, by setting a bigger initial perturbation amplitude, the collapse is faster and in some cases a final single object is formed.

The first generation of stars was formed from primordial gas. Numerical simulations suggest that the first stars were predominantly very massive, with typical masses M > 100 Mo. These stars were responsible for the reionization of the universe, the initial enrichment of the intergalactic medium with heavy elements, and other cosmological consequences. In this work, we study the structure of Zero Age Main Sequence stars for a wide mass and metallicity range and the evolution of 100, 150, 200, 250 and 300 Mo galactic and pregalactic Pop III very massive stars without mass loss, with metallicity Z=10E-6 and 10E-9, respectively. Using a stellar evolution code, a system of 10 equations together with boundary conditions are solved simultaneously. For the change of chemical composition, which determines the evolution of a star, a diffusion treatment for convection and semiconvection is used. A set of 30 nuclear reactions are solved simultaneously with the stellar structure and evolution equations. Several results on the main sequence, and during the hydrogen and helium burning phases, are described. Low metallicity massive stars are hotter and more compact and luminous than their metal enriched counterparts. Due to their high temperatures, pregalactic stars activate sooner the triple alpha reaction self-producing their own heavy elements. Both galactic and pregalactic stars are radiation pressure dominated and evolve below the Eddington luminosity limit with short lifetimes. The physical characteristics of the first stars have an important influence in predictions of the ionizing photon yields from the first luminous objects; also they develop large convective cores with important helium core masses which are important for explosion calculations.