Context. The elliptical galaxy NGC 3923 is known to be surrounded by a number of stellar shells, probable remnants of an accreted galaxy. Despite its uniqueness, the deepest images of its outskirts come from the 1980s. On the basis of the modified Newtonian dynamics (MOND), it has recently been predicted that a new shell lies in this region. Aims. We obtain the deepest image ever of the galaxy, map the tidal features in it, and search for the predicted shell. Methods. The image of the galaxy was taken by the MegaCam camera at the Canada-France-Hawaii Telescope in the g' band. It reached the surface-brightness limit of 29 mag/arcsec2. In addition, we reanalyzed an archival HST image of the galaxy. Results. We detected up to 42 shells in NGC 3923. This is by far the highest number among all shell galaxies. We present the description of the shells and other tidal features in the galaxy. A probable progenitor of some of these features was discovered. The shell system likely originates from two or more progenitors. The predicted shell was not detected, but the new image revealed that the prediction was based on incorrect assumptions and poor data.

Context. Many ellipticals are surrounded by round stellar shells probably stemming from minor mergers. A new method for constraining gravitational potential in elliptical galaxies has recently been suggested. It uses the spectral line profiles of these shells to measure the circular velocity at the edge of the shell and the expansion velocity of the shell itself. MOND is an alternative to the dark matter framework aiming to solve the missing mass problem. Aims. We study how the circular and expansion velocities behave in MOND for large shells. Methods. The asymptotic behavior for infinitely large shells is derived analytically. The applicability of the asymptotic results for finitely sized shells is studied numerically on a grid of galaxies modeled with S\'ersic spheres. Results. Circular velocity settles asymptotically at a value determined by the baryonic mass of the galaxy forming the baryonic Tully-Fisher relation known for disk galaxies. Shell expansion velocity also becomes asymptotically constant. The expansion velocities of large shells form a multibranched analogy to the baryonic Tully-Fisher relation, together with the galactic baryonic masses. For many - but not all - shell galaxies, the asymptotic values of these two types of velocities are reached under the effective radius. If MOND is assumed to work in ellipticals, then the shell spectra allow many details of the history to be revealed about the formation of the shell system, including its age. The results pertaining to circular velocities apply to all elliptical galaxies, not only those with shells.

Tests of MOND in ellipticals are relatively rare because these galaxies often lack kinematic tracers in the regions where the MOND effects are significant. Stellar shells observed in many elliptical galaxies offer a promising way to constrain their gravitational field. Shells appear as glowing arcs around their host galaxy. They are observed up to ~100 kpc. The stars in axially symmetric shell systems move in nearly radial orbits. The radial distributions of shell locations and the spectra of stars in shells can be used to constrain the gravitational potential of their host galaxy. The symmetrical shell systems, being especially suitable for these studies, occur in approximately 3% of all early-type galaxies. Hence the shells substantially increase the number of ellipticals in which MOND can be tested up to large radii. In this paper, we review our work on shell galaxies in MOND. We summarize the paper B\'{i}lek et al. (2013), where we demonstrated the consistency of shell radii in an elliptical NGC 3923 with MOND, and the work B\'{i}lek et al. (2014), in which we predicted a giant (~200 kpc), as yet undiscovered shell of NGC 3923. We explain the shell identification method, which was used in these two papers. We further describe the expected shape of line profiles in shell spectra in MOND which is very special due to the direct relation of the gravitational field and baryonic matter distribution (B\'{i}lek et al., 2014, in preparation).

Context. Stellar shells, which form axially symmetric systems of arcs in some elliptical galaxies, are most likely remnants of radial minor mergers. They are observed up a radius of $\sim$100 kpc. The stars in them oscillate in radial orbits. The radius of a shell depends on the free-fall time at the position of the shell and on the time since the merger. We previously verified the consistency of shell radii in the elliptical galaxy NGC 3923 with its most probable MOND potential. Our results implied that an as~yet undiscovered shell exists at the outskirts of the galaxy. Aims. We here extend our study by assuming more general models for the gravitational potential to verify the prediction of the new shell and to estimate its position. Methods. We tested the consistency of the shell radial distribution observed in NGC 3923 with a wide variety of MOND potentials of the galaxy. The potentials differed in the mass-to-light ratio and in distance to the galaxy. We considered different MOND interpolation functions, values of the acceleration constant $a_0$, and density profiles of the galaxy. We verified the functionality of our code on a Newtonian self-consistent simulation of the formation of a shell galaxy. Results. Our method reliably predicts that exactly one new outermost shell exists at a galactocentric radius of about 1900$^{\prime\prime}$ ($\sim$210 kpc) on the southwestern side of the galaxy. Its estimated surface brightness is about 28 mag arcsec$^{-2}$ in $B$ - a value accessible by current instruments. This prediction enables a rare test of MOND in an elliptical down to an acceleration of $a_0/10$. The predictive power of our method is verified by reconstructing the position of the largest known shell from the distribution of the remaining shells.

As the missing complement to existing studies of shell galaxies, we carried out a set of self-consistent N-body simulations of a minor merger forming a stellar shell system within a giant elliptical galaxy. We discuss the effect of a phenomenon possibly associated with the galaxy merger simulations --- a presence of multiple generations of shells.