  We study radial oscillations of non-rotating neutron stars (NSs) in
four-dimensional General Relativity. The interior of the NS was modelled within
a recently proposed multicomponent realistic equation of state (EoS) with the
induced surface tension (IST). In particular, we considered the IST EoS with
two sets of model parameters, that both reproduce all the known properties of
normal nuclear matter, give a high quality description of the proton flow
constraint, hadron multiplicities created in nuclear-nuclear collisions,
consistent with astrophysical observations and the observational data from the
NS-NS merger.
  We computed the 12 lowest radial oscillation modes, their frequencies and
corresponding eigenfunctions, as well as the large frequency separation for six
selected fiducial NSs (with different radii and masses of 1.2, 1.5 and 1.9
solar masses) of the two distinct model sets. The calculated frequencies show
their continuous growth with an increase of the NS central baryon density.
Moreover, we found correlations between the behaviour of first eigenfunction
calculated for the fundamental mode, the adiabatic index and the speed of sound
profile, which could be used to probe the internal structure of NSs with the
asteroseismology data.
