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 .

Keywords : neutron stars , equation of state , oscillations , asteroseismology