We study the origin of bulge stars and their angular momentum (AM) evolution
in 10 spiral galaxies with baryonic masses above 1010M⊙ in the
NIHAO galaxy formation simulations. The simulated galaxies are in good
agreement with observations of the relation between specific AM and mass of the
baryonic component and the stellar bulge-to-total ratio (B/T). We divide the
star particles at z=0 into disc and bulge components using a hybrid
photometric/kinematic decomposition method that identifies all central mass
above an exponential disc profile as the `bulge'. By tracking the bulge star
particles back in time, we find that on average 95\% of the bulge stars formed
{\it in situ}, 3\% formed {\it ex situ} in satellites of the same halo, and
only 2\% formed {\it ex situ} in external galaxies. The evolution of the AM
distribution of the bulge stars paints an interesting picture: the higher the
final B/T ratio, the more the specific AM remains preserved during the bulge
formation. In all cases, bulge stars migrate significantly towards the central
region, reducing their average galactocentric radius by roughly a factor 2,
independently of the final B/T value. However, in the higher B/T
(≳0.2) objects, the velocity of the bulge stars increases and the AM of
the bulge is almost conserved, whereas at lower B/T values, the velocity of
the bulge stars decreases and the AM of bulge reduces. The correlation between
the evolution of the AM and B/T suggests that bulge and disc formation are
closely linked and cannot be treated as independent processes.Comment: 17 pages, 16 Figures, 1 table; accepted for publication in MNRA