Stochastic Variational Search for $^{4}_{\Lambda\Lambda}$H

A four-body calculation of the $pn\Lambda\Lambda$ bound state, $^{\ 4}_{\Lambda\Lambda}$H, is performed using the stochastic variational method and phenomenological potentials. The $NN$, $\Lambda N$, and $\Lambda\Lambda$ potentials are taken from a recent Letter by Filikhin and Gal, PRL{\bf 89}, 172502 (2002). Although their Faddeev-Yakubovsky calculation found no bound-state solution over a wide range of $\Lambda\Lambda$ interaction strengths, the present variational calculation gives a bound-state energy, which is clearly lower than the $_\Lambda^3{H}+\Lambda$ threshold, even for a weak $\Lambda\Lambda$ interaction strength deduced from a recent experimental $B_{\Lambda\Lambda}(^{6}_{\Lambda\Lambda}{He})$ value. The binding energies obtained are close to, and slightly larger than, the values obtained from the three-body $d\Lambda\Lambda$ model in the Letter.