We present a dense model grid with tailored input chemical composition appropriate for the Large Magellanic Cloud. We use a one-dimensional hydrodynamic stellar evolution code, which accounts for rotation, transport of angular momentum by magnetic fields, and stellar wind mass loss to compute our detailed models. We calculate stellar evolution models with initial masses of 70-500 Msun and with initial surface rotational velocities of 0-550 km/s, covering the core-hydrogen burning phase of evolution. We find our rapid rotators to be strongly influenced by rotationally induced mixing of helium, with quasi-chemically homogeneous evolution occurring for the fastest rotating models. Above 160 Msun, homogeneous evolution is also established through mass loss, producing pure helium stars at core hydrogen exhaustion independent of the initial rotation rate. Surface nitrogen enrichment is also found for slower rotators, even for stars that lose only a small fraction of their initial mass. For models above 150 MZAMS, and for models in the whole considered mass range later on, we find a considerable envelope inflation due to the proximity of these models to their Eddington limit. This leads to a maximum zero-age main sequence surface temperature of 56000 K, at 180 Msun, and to an evolution of stars in the mass range 50-100 Msun to the regime of luminous blue variables in the HR diagram with high internal Eddington factors. Inflation also leads to decreasing surface temperatures during the chemically homogeneous evolution of stars above 180 Msun. The cool surface temperatures due to the envelope inflation in our models lead to an enhanced mass loss, which prevents stars at LMC metallicity from evolving into pair-instability supernovae. The corresponding spin-down will also prevent very massive LMC stars to produce long-duration gamma-ray bursts, which might, however, originate from lower masses.

Rotational mixing in massive main sequence stars is predicted to monotonically increase their surface nitrogen abundance with time. We use this effect to design a method for constraining the age and the inclination angle of massive main sequence stars, given their observed luminosity, effective temperature, projected rotational velocity and surface nitrogen abundance. This method relies on stellar evolution models for different metallicities, masses and rotation rates. We use the population synthesis code STARMAKER to show the range of applicability of our method. We apply this method to 79 early B-type main sequence stars near the LMC clusters NGC 2004 and N 11 and the SMC clusters NGC 330 and NGC 346. From all stars within the sample, 17 were found to be suitable for an age analysis. For ten of them, which are rapidly rotating stars without a strong nitrogen enhancement, it has been previously concluded that they did not evolve as rotationally mixed single stars. This is confirmed by our analysis, which flags the age of these objects as highly discrepant with their isochrone ages. For the other seven stars, their nitrogen and isochrone ages are found to agree within error bars, what validates our method. Constraints on the inclination angle have been derived for the other 62 stars,with the implication that the nitrogen abundances of the SMC stars, for which mostly only upper limits are known, fall on average significantly below those limits. Nitrogen chronology is found to be a new useful tool for testing stellar evolution and to constrain fundamental properties of massive main sequence stars. A web version of this tool is provided.