We perform a search for the rare decay B0s→μ+μ− using data collected by the D0 experiment at the Fermilab Tevatron Collider. This result is based on the full D0 Run II data set corresponding to 10.4 fb−1 of p¯p collisions at √s=1.96 TeV. We use a multivariate analysis to increase the sensitivity of the search. In the absence of an observed number of events above the expected background, we set an upper limit on the decay branching fraction of B(B0s→μ+μ−)<15×10−9 (12×10−9) at the 95% (90%) C.L.5 MoreReceived 23 January 2013DOI:https://doi.org/10.1103/PhysRevD.87.072006© 2013 American Physical Society