The radial wave functions of the Bear–Hodgson potential have been used to study the ground state features such as the proton, neutron and matter densities and the as- sociated rms radii of two neutrons halo 6He, 11Li, 14Be and 17B nuclei. These halo nuclei are treated as a three-body system composed of core and outer two-neutron (Core + n + n). The radial wave functions of the Bear–Hodgson potential are used to describe the core and halo density distributions. The interaction of core-neutron takes the Bear–Hodgson potential form. The outer two neutrons of 6He and 11Li interact by the realistic interaction REWIL whereas those of 14Be and 17B interact by the realistic interaction of HASP. The obtained results show that this model succeeds in reproducing the neutron halo in these nuclei. From the calculated densities, it is found that 6He, 11Li, 14Be and 17B have a long tail in neutron and matter densities which is consistent with the experimental data. Elastic charge form factors for these halo nuclei are analyzed via the plane wave Born approximation.