Multiscalar B-L extension based on S4 flavor symmetry for neutrino masses and mixing

A multiscalar and nonrenormalizable B - L extension of the standard model (SM) with S-4 symmetry which successfully explains the recently observed neutrino oscillation data is proposed. The tiny neutrino masses and their hierarchies are generated via the type-I seesaw mechanism. The model reproduces the recent experiments of neutrino mixing angles and Dirac CP violating phase in which the atmospheric angle (theta(23)) and the reactor angle (theta(13)) get the best-fit values while the solar angle (theta(12)) and Dirac CP violating phase (delta) are in 3 sigma range of the best-fit value for the normal hierarchy (NH). For the inverted hierarchy (IH), theta(13) gets the best-fit value and theta(23) together with delta are in the 1 sigma range, while theta(12) is in 3 sigma range of the best-fit value. The effective neutrino masses are predicted to be < m(ee)> = 6.81 meV for the NH and < m(ee)> = 48.48 meV for the IH, in good agreement with the most recent experimental data.