Not yet reviewed by Pith; the record is open.

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.