Possible bound nuclei beyond the two-neutron drip line in the <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mn>50</mml:mn><mml:mo>≤</mml:mo><mml:mi>Z</mml:mi><mml:mo>≤</mml:mo><mml:mn>70</mml:mn></mml:mrow></mml:math> region

Possible bound nuclei beyond the two-neutron drip line in the $50\ensuremath{\le}Z\ensuremath{\le}70$ region are investigated by using the deformed relativistic Hartree-Bogoliubov theory in continuum with density functional PC-PK1. Bound nuclei beyond the drip lines of $_{56}\mathrm{Ba}$, $_{58}\mathrm{Ce}$, $_{62}\mathrm{Sm}$, $_{64}\mathrm{Gd}$, and $_{66}\mathrm{Dy}$ are predicted, forming peninsulas of stability in the nuclear landscape. Near these peninsulas, several multineutron emitters are predicted. The underlying mechanism of the peninsulas of stability is investigated by studying the total energy, Fermi surface, quadrupole deformation, and the single-neutron spectrum in the canonical basis. It is found that the deformation effect is crucial for forming the peninsulas of stability, and pairing correlations are also essential in specific cases. The dependence on the deformation evolution is also discussed. The decay rates of multineutron radioactivity in Ba and Sm isotopic chains are estimated by using the direct decay model.