Applying the semi-realistic N N interactions that include realistic tensor force to the Hartree-Fock calculations, we investigate tensor-force effects on the single-particle levels in the Ca isotopes. The semi-realistic interaction successfully describes the experimental difference between ε(p1s 1/2 ) and ε(p0d 3/2 ) (denoted by ∆ε13) both at 40 Ca and 48 Ca, confirming importance of the tensor force. The tensor force plays a role in the N -dependence of ∆ε13 also in neutron-rich Ca nuclei. While the p1s 1/2 -p0d 3/2 inversion is predicted in heavier Ca nuclei as in 48 Ca, it takes place only in N ≥ 46, delayed by the tensor force. We further investigate possibility of proton bubble structure in Ar, which is suggested by the p1s 1/2 -p0d 3/2 inversion in 48 Ca and more neutron-rich Ca nuclei, by the spherical Hartree-Fock-Bogolyubov calculations. Even with the inversion at 48 Ca the pair correlation prohibits prominent bubble distribution in 46 Ar. Bubble in Ar is unlikely also near neutron drip line because either of unboundness or of deformation. However, 34 Si remains a candidate for proton bubble structure, owing to large shell gap between p1s 1/2 and p0d 5/2 .PACS numbers: 21.10. Pc, 21.10.Ft, 21.60.Jz, 27.40.+z, 27.30.+t Introduction. Owing to the progress in experiments on unstable nuclei, it has been recognized [1] that the nuclear shell structure depends on Z or N as often called shell evolution. Moreover, it is now known [2-4] that the tensor force, which should be contained in the nucleonnucleon (N N ) interaction, plays a crucial role in the shell evolution. While the tensor force had been ignored in the conventional mean-field (MF) or the energy-densityfunctional (EDF) approaches, there have been several attempts to incorporate the tensor force into those approaches; e.g. the calculations with the Skyrme [5][6][7][8] or the Gogny [9, 10] interactions. However, without wellestablished strengths and/or radial-dependence of the tensor force, it is not straightforward to pin down tensorforce effects quantitatively from those interactions. One of the authors (H.N.) has developed the so-called semirealistic N N interactions [11], which is applicable to the self-consistent MF calculations. The recent parametersets, M3Y-Pn with n ≥ 5 [3,4,12], include the tensor force originating from the G-matrix at the nuclear surface. Because of the realistic nature of the tensor force in them, these semi-realistic interactions are suitable to investigate the tensor-force effects on the shell evolution.