Searches for B0 decays to combinations of charmless isoscalar mesons

We search for B meson decays into two-body combinations of eta, eta', omega, and phi mesons from 89 million B B-bar pairs collected with the BaBar detector at the PEP-II asymmetric-energy e+e- collider at SLAC. We find the branching fraction BF(B0 -> eta omega) = (4.0^{+1.3}_{-1.2} +- 0.4) x 10^-6 with a significance of 4.3 sigma. For all the other decay modes we set the following 90% confidence level upper limits on the branching fractions, in units of 10^-6 : BF(B0 -> eta eta)<2.8, BF(B0 -> eta eta')<4.6, BF(B0 -> eta' eta')<10, BF(B0 -> eta'omega)<2.8, BF(B0 -> eta phi)<1.0, BF(B0 -> eta' phi)<4.5, BF(B0 -> phi phi)<1.5.