Neutron drops confined in an external field are studied in the framework of relativistic Brueckner-Hartree-Fock theory using the bare nucleon-nucleon interaction. The ground state energies and radii of neutron drops with even numbers from N = 4 to N = 50 are calculated and compared with results obtained from other nonrelativistic ab initio calculations and from relativistic density functional theory. Special attention has been paid to the magic numbers and to the sub-shell closures. The single-particle energies are investigated and the monopole effect of the tensor force on the evolutions of the spin-orbit and the pseudospin-orbit splittings is discussed. The results provide interesting insight of neutron rich systems and can form an important guide for future density functionals.Recently, the self-consistent relativistic Brueckner-Hartree-Fock (RBHF) theory for finite nuclei has been established, and the results are in much better agreement with experimental data than the nonrelativistic calculations with the 2N interaction only [22,23]. Indeed, it is known since more than 30 years that relativistic Brueckner-Hartree-Fock theory gives a much better description of the nuclear matter saturation properties than nonrelativistic BHF theories [38][39][40]. In nonrelativistic many-body investigations on the influence of various types of 3N-interactions, it was found that a relativistic effect, the so-called Z-diagram, plays a major role [41].Having these progresses in mind, it is important to study the neutron drops in more detail in the framework of RBHF theory and compare the results with other nonrelativistic ab initio calculations using various 2N or 2N + 3N interactions, as well as calculations using various density functionals. This can also provide valuable insight to improve current relativistic density functionals. In Ref. [28], a systematic and specific pattern due to the tensor forces in the evolution of spin-orbit splittings based on RBHF theory is reported.In this work, we investigate neutron drops confined in an external harmonic oscillator potential using relativistic Brueckner-Hartree-Fock theory, and present the numerical details and calculated results in detail. In Sec. II, we give a brief outline of the RBHF framework for neutron drops. The numerical details are discussed in Sec. III. Results and discussion for neutron drops with an even number of neutrons from N = 4 to 50 will be presented in Sec. IV. Finally, a summary and perspectives for future investigations will be given in Sec.V.