Stellar shells around galaxies could provide precious insights into their assembly history. However, their formation mechanism remains poorly empirically constrained. We present MUSE@VLT data of the most prominent outer shell of NGC474, to constrain its formation history. The stellar shell spectrum is clearly detected, with a signal-to-noise ratio of around 65 pix$^{-1}$. We use a full spectral fitting method to determine the line-of-sight velocity and the age and metallicity of the shell and associated point-like sources within the MUSE field of view. We detect six globular cluster (GC) candidates and eight planetary nebula (PN) candidates which are all kinematically associated to the stellar shell. We show that the shell has an intermediate metallicity, [M/H] = $-0.83^{+0.12}_{-0.12}$ and a possible alpha-enrichement, [alpha/Fe] ~ 0.3. This metallicity and the number of PNe are consistent with a progenitor galaxy more massive than 10^9 Msun. The photometry of the shell suggest that it is possibly composed of stars from the low metallicity outskirts of an intermediate mass (~1:10) companion. We show that at least two globular cluster candidates are quite young, with ages below 1.5 Gyr, which gives a lower-limit for the merger age. We note that spectroscopic data on the center of NGC~474 shows the presence of a young (~1~Gyr) stellar population. We conclude that the shell formation event is likely to be an intermediate-mass merger which happened around 1.5~Gyr ago and which has collaterally triggered the formation of massive, up to ~ 6 x 10^4 Msun, star cluster, and a nuclear starburst in NGC474.

NGC 4993 is the shell galaxy host of the GRB170817A short gamma-ray burst and the GW170817 gravitational-wave event produced during a binary-neutron-star coalescence. The galaxy shows signs, including the stellar shells, that it has recently accreted a smaller, late-type galaxy. The accreted galaxy might be the original host of the binary neutron star. We measured the positions of the stellar shells of NGC 4993 in an HST/ACS archival image and use the shell positions to constrain the time of the galactic merger. According to the analytical model of the evolution of the shell structure in the expected gravitational potential of NGC 4993, the galactic merger happened at least 200 Myr ago, with a probable time roughly around 400 Myr and the estimates higher than 600 Myr being improbable. This constitutes the lower limit on the age of the binary neutron star, because the host galaxy was probably quenched even before the galactic merger, and the merger has likely shut down the star formation in the accreted galaxy. We roughly estimate the probability that the binary neutron star originates in the accreted galaxy to be around 30%.

We present the photometric properties of 2210 newly identified dwarf galaxy candidates in the MATLAS fields. The Mass Assembly of early Type gaLAxies with their fine Structures (MATLAS) deep imaging survey mapped $\sim$142 deg$^2$ of the sky around nearby isolated early type galaxies using MegaCam on the Canada-France-Hawaii Telescope, reaching surface brightnesses of $\sim$ 28.5 - 29 in the g-band. The dwarf candidates were identified through a direct visual inspection of the images and by visually cleaning a sample selected using a partially automated approach, and were morphologically classified at the time of identification. Approximately 75% of our candidates are dEs, indicating that a large number of early type dwarfs also populate low density environments, and 23.2% are nucleated. Distances were determined for 13.5% of our sample using pre-existing $z_{spec}$ measurements and HI detections. We confirm the dwarf nature for 99% of this sub-sample based on a magnitude cut $M_g$ = -18. Additionally, most of these ($\sim$90%) have relative velocities suggesting that they form a satellite population around nearby massive galaxies rather than an isolated field sample. Assuming that the candidates over the whole survey are satellites of the nearby galaxies, we demonstrate that the MATLAS dwarfs follow the same scaling relations as dwarfs in the Local Group as well as the Virgo and Fornax clusters. We also find that the nucleated fraction increases with $M_g$, and find evidence of a morphology-density relation for dwarfs around isolated massive galaxies.

Stellar streams are regarded as crucial objects to test galaxy formation models, with their morphology tracing the underlying potentials and their occurrence tracking the assembly history of the galaxies. The existence of one of the most iconic stellar streams, the double loop around NGC5907, has recently been questioned by new observations with the Dragonfly telescope. This new work only finds parts of the stream, even though they reach a 1 sigma surface brightness limit of 30.3 mag per sq. arcsec in the g-band. Using 7.2 hours of Luminance L-band imaging with the Milankovi\'c 1.4 meter telescope, we have re-observed the putative double loop part to confirm or reject this assessment. We do not find signs of the double loop, but see only a single, knee-shaped stellar stream. Comparing our to the data by the Dragonfly team, we find the same features. Our observations reach a 1 sigma surface brightness limit of 29.7 mag per sq. arcsec in the g-band. These findings emphasize the need for independent confirmation of detections of very low-surface brightness features.

Stellar shells are low surface brightness features in the form of open, concentric arcs, formed in close-to-radial collisions of galaxies. They occur in at least 10% of early-type galaxies and a small portion of spirals and their unique kinematics carry valuable information about the host galaxies. We discuss a method using measurements of the number and distribution of shells to estimate the mass distribution of the galaxies and the time since the merger. The method is applied on the shells of NGC 4993 - a galaxy hosting the electromagnetic counterpart of the gravitational wave event GW170817, to estimate the probable time since the galactic merger. We used analytical calculations and particle simulations to show that, in special cases, when kinematic data are available, further constraints on mass distribution and merger time can be derived. Applying the methods to the rapidly growing sample of known shell galaxies will constrain the dark-matter content in the galaxies and reveal detailed information on the recent merger history of the Universe.

The missing mass problem has not been solved decisively yet. Observations show that if gravity is to be modified, then the MOND theory is its excellent approximation on galactic scales. MOND suggests an adjustments of the laws of physics in the limit of low accelerations. Comparative simulations of interacting galaxies in MOND and Newtonian gravity with dark matter revealed two principal differences: 1) galaxies can have close flybys without ending in mergers in MOND because of weaker dynamical friction, and 2) tidal dwarf galaxies form very easily in MOND. When this is combined with the fact that many interacting galaxies are observed at high redshift, we obtain a new perspective on tidal features: they are often formed by non-merging encounters and tidal disruptions of tidal dwarf galaxies. Here we present the results from our self-consistent MOND $N$-body simulation of a close flyby of two galaxies similar to the Milky Way. It turns out that most types of the structures that are traditionally assigned to galaxy mergers can be formed by non-merging encounters, including tidal arms, bridges, streams, shells, disk warps, thick disks, and most probably also disks of satellites. The success of MOND in explaining the dynamics of galaxies hints us that this way of formation of tidal structures should be considered seriously.

We report that the density profiles of globular cluster (GC) systems in a sample of 17 early-type galaxies (ETGs) show breaks at the radii where the gravitational acceleration exerted by the stars equals the galactic acceleration scale $a_0$ known from the radial acceleration relation or MOND. The match with the other characteristic radii in the galaxy is not that close. We propose possible explanations in the frameworks of the $\Lambda$CDM model and MOND. We find tentative evidence that in the $\Lambda$CDM context, GCs reveal not only the masses of the dark halos through the richness of the GC systems but also the concentrations through the break radii of the GC systems.

Using Jeans modelling, we calculate the velocity dispersion profile of the ultra-diffuse galaxy (UDG) Dragonfly 44 in MOND. For the nominal mass-to-light ratio from the literature and an isotropic profile, the agreement with the data is excellent near the center of the galaxy. However, in modified gravity, close to the cluster core, the gravitational environment should bring the galaxy back towards the Newtonian behavior. This success of the isolated MOND prediction for the central velocity dispersion could then mean that the galaxy is at large distances ($\gg5$ Mpc) from the cluster core, as hinted by the fact that close-by UDGs share similar velocities with a dispersion well below that of the cluster itself. There is however a two-sigma tension in the outer part of the UDG due to an increase of the observed dispersion profile with respect to the flat MOND prediction. This deviation could simply be a measurement error. Other possibilities could be, for a UDG far from the cluster, a higher-than-nominal baryonic mass with a tangentially anisotropic dispersion profile, or even a sign of a dark baryonic halo. If the UDG is closer to the cluster core, the deviation could be a sign that it is in the process of disruption.

The possible existence of two dark matter free galaxies (NGC1052-DF2 and NGC1052-DF4) in the field of the early-type galaxy NGC1052 presents a challenge to theories of dwarf galaxy formation according to the current cosmological paradigm. We carried out a search for signatures of past interactions connected to the putative hosts of NGC1052-DF2 and NGC1052-DF4 using a very deep L-band image obtained with the 0.7m Jeanne Rich telescope that reach a surface brightness limit of 28.5 mag arcsec-2 in the r band. We found several low-surface brightness features, possibly consistent with an ongoing merger history in this group. We find a tidal interaction between NGC1052 and NGC1047, confirming a physical association. Furthermore we find a stellar loop around NGC1052 in the direction of NGC1042 and a stellar stream pointing in the direction of NGC1052-DF2 - but both are not directly connected. We find no evidence for a recent tidal interaction for NGC1052-DF2 and NGC1052-DF4. No low surface brightness features have been uncovered around the spiral galaxy NGC1042, which leaves the association (physical or projected) between NGC1052 and NGC1042 ambiguous, though they have similar radial velocities. Their association will only be established once accurate distances to both objects have been measured.

(abridged) Context. Gravitational fields at the outskirts of early-type galaxies (ETGs) are difficult to constrain observationally. It thus remains poorly explored how well the LCDM and MOND hypotheses agree with ETGs. Aims. This led us to gather a large sample of ETGs and examine homogeneously which dark matter halos they occupy, whether the halos follow the theoretically predicted stellar-to-halo mass relation (SHMR) and the halo mass-concentration relation (HMCR), whether ETGs obey MOND and the radial acceleration relation (RAR) observed for late-type galaxies (LTGs), and finally whether LCDM or MOND perform better in ETGs. Methods. We employed Jeans analysis of radial velocities of globular clusters (GCs). We analysed nearly all ETGs having more than about 100 archival GC radial velocity measurements. The GC systems of our 17 ETGs extend mostly over ten effective radii. A LCDM simulation of GC formation helped us to interpret the results. Results. Successful LCDM fits are found for all galaxies, but compared to the theoretical HMCR and SHMR, the best-fit halos usually have concentrations that are too low and stellar masses that are too high for their masses. This might be because of tidal stripping of the halos or because ETGs and LTGs occupy different halos. Most galaxies can be fitted by the MOND models successfully as well, but for some of the galaxies, especially those in centers of galaxy clusters, the observed GC velocity dispersions are too high. This might be a manifestation of the additional dark matter that MOND requires in galaxy clusters. Additionally, we find many signs that the GC systems were perturbed by galaxy interactions. Formal statistical criteria prefer the best-fit LCDM models over the MOND models, but this might be due to the higher flexibility of the LCDM models. The MOND approach can predict the GC velocity dispersion profiles better.

(abridged) The Milky Way (MW) and Andromeda (M31) galaxies possess rotating planes of satellites. Their formation has not been explained satisfactorily yet. It was suggested that the MW and M31 satellites are ancient tidal dwarf galaxies, which could explain their configuration. This suggestion gained support by an analytic backward calculation of the relative MW-M31 orbit in the MOND modified dynamics paradigm by Zhao et al. (2013) implying their close flyby 7-11 Gyr ago. Here we explore the Local Group history in MOND in more detail using a simplified first-ever self-consistent simulation. We note the features induced by the encounter in the simulation and identify their possible real counterparts. The simulation was set to approximately reproduce the observed MW and M31 masses, effective radii, separation, relative velocity and disk inclinations. We used the publicly available adaptive-mesh-refinement code Phantom of RAMSES. In the simulation, matter was transferred from the MW to M31 along a tidal tail. The encounter induced formation of several structures resembling the peculiarities of the Local Group. Most notably: 1) A rotating planar structure formed around M31 from the transferred material. It had a size similar to the observed satellite plane and was oriented edge-on to the simulated MW, just as the real one. 2) The same structure also resembled the tidal features observed around M31 by its size and morphology. 3) A warp in the MW developed with an amplitude and orientation similar to that observed. 4) A cloud of particles formed around the simulated MW, with the extent of the actual MW satellite system. The encounter did not end by merging in a Hubble time. The simulation thus demonstrated that MOND can possibly explain many peculiarities of the Local Group and, moreover, that non-merging galaxy encounters in MOND can produce tidal features in galaxies.

Investigating the gravitational field in the early-type galaxies (ETGs, i.e. ellipticals and lenticulars) up to large radii is observationally difficult. It is questionable how the radial acceleration (RAR) in the ETGs looks like, i.e. the relation between the dynamically inferred gravitational acceleration and the acceleration expected from the distribution of the visible matter. This relation is nearly universal for the spiral galaxies, in agreement with the MOND modified dynamics paradigm. In this contribution, we investigate a sample of 15 ETGs. We extract their full kinematic profiles out to several effective radii from their globular cluster systems and estimate their gravitational field using the Jeans equation. We parametrize the gravitational field by that produced by the stars and a Navarro-Frenk-White DM halo. We find that only 4-5 of our ETGs follow the RAR for the spiral galaxies. All these galaxies are fast rotators, have disky isophotes, appear mostly very elongated and the have bluest colors in our sample. This suggests that they might be spiral galaxies which lost their gas. Our galaxies deviating from the RAR for the spirals either disprove MOND, contain unobserved matter, or indicate a flaw in the method.

We present the first release of the Czech Variable star catalogue that currently contains 1228 stars whose variability was discovered by 60 Czech observers. The catalogue contains confirmed variable stars of various types, but also candidates. We give precise coordinates, cross identification with other catalogues, information about constellation, variability type, brightness, light elements, name of the discoverer and year of discovery. In eighty-eight percent of stars the variability type is estimated, for more than 60 % of the stars the light ephemerides are given.

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.

(doctoral thesis of Michal Bilek, finished on June 19, 2015) MOND is an observational rule for predicting the acceleration of stars and galaxies from the distribution of the visible matter. It possibly stems from a new law of physics. I list the theoretical aspects of MOND, its achievements and problems. MOND has been tested mainly in disc galaxies so far. Its tests in elliptical galaxies are rare because the MOND effects are small for them in the parts observable by the conventional methods. In the thesis, I explain the methods and ideas I developed for testing MOND in the ellipticals using stellar shells. Moreover, the shells enable us to test MOND for stars in radial orbits for the first time. The shells are results of galactic interactions. I discuss the shell formation mechanisms and summarize the findings from shell observations and simulations. The thesis contains as yet unpublished results mainly in: 1) the introduction of Sect. 3 (the expected differences in the shell morphology in the Newtonian dynamics with dark matter and in MOND), and 2) Sect. 2.4.1 (formulas for modeling the evolution of shell radii).

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.