β-delayed one-neutron and two-neutron branching ratios (P1n and P2n) have been measured in the decay of A=84 to 87 Ga isotopes at the Radioactive-Isotope Beam Factory (RIBF) at the RIKEN Nishina Center using a high-efficiency array of He3 neutron counters (BRIKEN). Two-neutron emission was observed in the decay of Ga84,85,87 for the first time and the branching ratios were measured to be P2n=1.6(2)%,1.3(2)%, and 10.2(28)stat(5)sys%, respectively. One-neutron branching ratio of Ga87(P1n=81(9)stat(8)sys%) and half-life of 29(4) ms were measured for the first time. The branching ratios of Ga86 were also measured to be P1n=74(2)stat(8)sys% and 16.2(9)stat(6)sys% with better precision than a previous study. The observation that P1n>P2n for both Ga86,87 was unexpected and is interpreted as a signature of dominating one-neutron emission from the two-neutron unbound excited states in Ge86,87. In order to interpret the experimental results, shell-model and Hauser-Feshbach statistical model calculations of delayed particle and γ-ray emission probabilities were performed. This model framework reproduces the experimental results. The shell model alone predicts P2n significantly larger than P1n for the Ga87 decay, and it is necessary to invoke a statistical description to successfully explain the observation that P1n>P2n. Our new results demonstrate the relevance and importance of a statistical description of neutron emission for the prediction of the decay properties of multineutron emitters and that it must be included in the r-process modeling
