  We explain the 2.3 sigma deviation in the recent measurements of the neutral
B mesons decay into muon pairs from the standard model prediction in the
framework of supersymmetric grand unified models using anti-symmetric coupling
as a new source of flavor violation. We show a correlation between the B_d ->
mu^+mu^- decay and the CP phase in the B_d -> J/psi K decay and that their
deviations from the standard model predictions can be explained after
satisfying constraints arising from various hadronic and leptonic rare decay
processes, B-bar{B}, K-bar{K} oscillations data and electric dipole moments of
electron and neutron. The allowed parameter space is typically represented by
pseudoscalar Higgs mass m_A < 1 TeV and tan beta_H (= {v_u}/{v_d}) < 20 for
squark and gluino masses around 2 TeV.
