Measurements of various lepton flavor universal observables in $b \to\, s\,l^+ l^-$ transition decays continue to disagree

with the standard model expectations. The recent update of $R_K$ measurement from LHCb still indicates $3.1\sigma$

deviation from the standard model.

Similarly, the measurements of other observables such as $R_{K^*}$, $P_{5}^{\prime}$ and

$\mathcal{B}(B_s\, \to\, \phi \,\mu^+ \,\mu^-)$ continue to show disagreement with standard model predictions.

It is well known that there exists a very close relation between $b \to\, s\,l^+ l^-$ and

$b \to\, s\,\nu \bar{\nu}$ decays not only in standard model but also in beyond the standard model physics.

In beyond the standard model physics these decay processes are related via $SU(2)_L$ gauge symmetry which relates

neutrinos to the left handed charged leptons.

Moreover, the $B$ decays with $\nu \bar{\nu}$ final state

are theoretically cleaner channels than the corresponding $b \to\, s\,l^+ l^-$ neutral transitions

as they do not suffer from hadronic uncertainties beyond the form factors such as the non-factorizable

corrections and photonic penguin contributions.

Hence, we explore $B_s\, \to\, (\phi\,,\eta\,,\eta^\prime) \,\nu \,\bar{\nu}$ decays mediated via

$b\, \to\, s\, \nu\,\bar{\nu}$ transitions using the standard model effective field theory formalism, $Z'$ and scalar and vector leptoquarks new physics models.

We give predictions of branching ratio in standard model and in the presence of various new physics couplings.