We study the charmless decays B to Lambda Lambdabar h, where h stands for pi^+, K^+,K^0,K^{*+}, or K^{*0}, using a 605 fb^{-1} data sample collected at the Upsilon(4S)resonance with the Belle detector at the KEKB asymmetric energy e^+ e^- collider. W
e observe B0 to Lambda Lambdabar K0 and B0 to Lambda Lambdabar K*0 with branching fractions of (4.76^{+0.84}_{-0.68} (stat.) pm 0.61 (syst.)) times 10^{-6} and (2.46^{+0.87}_{-0.72} pm 0.34) times 10^{-6}, respectively. The significances of these signals in the threshold-mass enhanced mass region are 12.4 sigma and 9.3 sigma, respectively. We also update the branching fraction BF(B+ to Lambda Lambdabar K+) = (3.38^{+0.41}_{-0.36} pm 0.41)times 10^{-6} with better accuracy,and report the following measurement or 90% confidence level upper limit in the threshold-mass-enhanced region: BF(B+ to Lambda Lambdabar K*+) = (2.19^{+1.13}_{-0.89} pm 0.33) times 10^{-6} with 3.7 sigma significance; BF(B+ to Lambda Lambdabar pi+) < 0.94 times 10^{-6}. A related search for B0 to Lambda Lambdabar D0bar yields a branching fraction BF(B0 to Lambda Lambdabar D0bar) = (1.05_{-0.44}^{+0.57} pm 0.13)times 10^{-5}. This may be compared with the large, ~10^{-4}, branching fraction observed for B0 to p pbar D0bar. The m_{Lambda Lambdabar} enhancements near threshold and related angular distributions for the observed modes are also reported. PACS:13.25.Hw, 14.40.Nd
