We determine the mass of the Milky Way bar and the torque it causes, using Gaia DR2, by applying the orbital arc method. Based on this, we have found that the gravitational acceleration is not directed towards the centre of our Galaxy but a few degrees away from it. We propose that the tangential acceleration component is caused by the bar of the Galaxy. Calculations based on our model suggest that the torque experienced by the region around the Sun is $\approx 2400\, km^2 s^{-2}$ per solar mass. The mass estimate for the bar is $\sim 1.6\pm0.3\times10^{10} M_\odot$. Using greatly improved data from Gaia DR2, we have computed the acceleration field to great accuracy by adapting the oPDF method (Han et al. 2016) locally and used the phase space coordinates of $\sim 4\times10^5$ stars within a distance of 0.5 kpc from the Sun. In the orbital arc method, the first step is to guess an acceleration field and then reconstruct the stellar orbits using this acceleration for all the stars within a specified region. Next, the stars are redistributed along orbits to check if the overall phase space distribution has changed. We repeat this process until we find an acceleration field that results in a new phase space distribution that is the same as the one that we started with; we have then recovered the true underlying acceleration.

We study the effect of dynamical friction on globular clusters and on the stars evaporated from the globular clusters (stellar streams) moving in a galactic halo. Due to dynamical friction, the position of a globular cluster (GC) as a stream progenitor starts to shift with respect to its original position in the reference frame of initial GC orbit. Therefore the stars that have evaporated at different times have different mean position with respect to the GC position. This shifting results in a certain asymmetry in stellar density distribution between the leading and trailing arms of the stream. The degree of the asymmetry depends on the characteristics of the environment in which the GC and the stream stars move. As GCs are located mainly in outer parts of a galaxy, this makes dynamical friction a unique probe to constrain the underlying dark matter spatial density and velocity distributions. For a GC NGC 3201 we compared our theoretical shift estimates with available observations. Due to large uncertainties in current observation data, we can only conclude that the derived estimates have the same order of magnitude.

Galaxy filaments are the dominant feature in the overall structure of the cosmic web. The study of the filamentary web is an important aspect in understanding galaxy evolution and the evolution of matter in the Universe. A map of the filamentary structure is an adequate probe of the web. We propose that photometric redshift galaxies are significantly positively associated with the filamentary structure detected from the spatial distribution of spectroscopic redshift galaxies. The catalogues of spectroscopic and photometric galaxies are seen as point-process realisations in a sphere, and the catalogue of filamentary spines is proposed to be a realisation of a random set in a sphere. The positive association between these sets was studied using a bivariate $J-$function, which is a summary statistics studying clustering. A quotient $D$ was built to estimate the distance distribution of the filamentary spine to galaxies in comparison to the distance distribution of the filamentary spine to random points in $3-$dimensional Euclidean space. This measure gives a physical distance scale to the distances between filamentary spines and the studied sets of galaxies. The bivariate $J-$function shows a statistically significant clustering effect in between filamentary spines and photometric redshift galaxies. The quotient $D$ confirms the previous result that smaller distances exist with higher probability between the photometric galaxies and filaments. The trend of smaller distances between the objects grows stronger at higher redshift. Additionally, the quotient $D$ for photometric galaxies gives a rough estimate for the filamentary spine width of about $1$~Mpc. Photometric redshift galaxies are positively associated with filamentary spines detected from the spatial distribution of spectroscopic galaxies.

We present a method to calculate gravitational potential gradients within regions containing few tens of thousands stars with known phase space coordinates. The central idea of the method is to calculate orbital arcs for each star within a given region for a certain parametrised potential (gravitational acceleration) and to assume that position of each star on its orbital arc is a random variable with a uniform probability density in time. Thereafter, by combining individual probability densities of stars it is possible to calculate the overall probability density distribution and likelihood for a given region as a function of gravitational acceleration parameters. The likelihood has a maximum if the calculated probability distribution and the observed distribution of stars in phase space are consistent. This allows us to constrain gravitational accelerations and potential gradient values. The method assumes that phases of stars are mixed within the regions where stellar orbits are calculated. We tested the method for 12 small rectangular regions within simulated disc galaxy from Gaia Wiki. Tests show that even with a rather simple acceleration form the calculated accelerations in galactic plane coincide with their true values from simulation about 5 per cent, misalignment between the calculated and true acceleration vector directions is less than 1 degree (median values). The model can be used with the Milky Way Gaia complete solution data.

Galaxy groups and clusters are formidable cosmological probes. They permit the studying of the environmental effects on galaxy formation. A reliable detection of galaxy groups is an open problem and is important for ongoing and future cosmological surveys. We propose a probabilistic galaxy group detection algorithm based on marked point processes with interactions. The pattern of galaxy groups in a catalogue is seen as a random set of interacting objects. The positions and the interactions of these objects are governed by a probability density. The estimator of the unknown cluster pattern is given by the configuration of objects maximising the proposed probability density. Adopting the Bayesian framework, the proposed probability density is maximised using a simulated annealing algorithm. The method provides "for free" additional information such as the probabilities that a point or two points in the observation domain belong to the cluster pattern, respectively. These supplementary tools allow the construction of tests and techniques to validate and to refine the detection result. To test the feasibility of the proposed methodology, we applied it to the well-studied 2MRS data set. Compared to previously published Friends-of-Friends (FoF) group finders, the proposed Bayesian group finder gives overall similar results. The proposed Bayesian group finder is tested on a galaxy redshift survey, but more detailed analyses are needed to understand the actual capabilities of the algorithm regarding upcoming cosmological surveys. The presented mathematical framework permits adapting it easily for other data sets (in astronomy and in other fields of sciences). In cosmology, one promising application is the detection of galaxy groups in photometric galaxy redshift surveys, while taking into account the full photometric redshift posteriors. (abridged)

We study galaxy populations and search for possible merging substructures in the rich galaxy cluster A2142. Normal mixture modelling revealed in A2142 several infalling galaxy groups and subclusters. The projected phase space diagram was used to analyse the dynamics of the cluster and study the distribution of various galaxy populations in the cluster and subclusters. The cluster, supercluster, BCGs, and one infalling subcluster are aligned. Their orientation is correlated with the alignment of the radio and X-ray haloes of the cluster. Galaxies in the centre of the main cluster at the clustercentric distances $0.5~h^{-1}Mpc$ have older stellar populations (with the median age of $10 - 11$~Gyrs) than galaxies at larger clustercentric distances. Star-forming and recently quenched galaxies are located mostly in the infall region at the clustercentric distances $D_{\mathrm{c}} \approx 1.8~h^{-1}Mpc$, where the median age of stellar populations of galaxies is about $2$~Gyrs. Galaxies in A2142 have higher stellar masses, lower star formation rates, and redder colours than galaxies in other rich groups. The total mass in infalling groups and subclusters is $M \approx 6\times10^{14}h^{-1}M_\odot$, approximately half of the mass of the cluster, sufficient for the mass growth of the cluster from redshift $z = 0.5$ (half-mass epoch) to the present. The cluster A2142 may have formed as a result of past and present mergers and infallen groups, predominantly along the supercluster axis. Mergers cause complex radio and X-ray structure of the cluster and affect the properties of galaxies in the cluster, especially in the infall region. Explaining the differences between galaxy populations, mass, and richness of A2142, and other groups and clusters may lead to better insight about the formation and evolution of rich galaxy clusters.

New and more reliable distances and proper motions of a large number of stars in the Tycho-Gaia Astrometric Solution (TGAS) catalogue allow to calculate the local matter density distribution more precisely than earlier. We devised a method to calculate the stationary gravitational potential distribution perpendicular to the Galactic plane by comparing the vertical probability density distribution of a sample of observed stars with the theoretical probability density distribution computed from their vertical coordinates and velocities. We applied the model to idealised test stars and to the real observational samples. Tests with two mock datasets proved that the method is viable and provides reasonable results. Applying the method to TGAS data we derived that the total matter density in the Solar neighbourhood is $0.09\pm 0.02 \text{M}_\odot\text{pc}^{-3}$ being consistent with the results from literature. The matter surface density within $|z|\le 0.75 \text{kpc}$ is $42\pm 4 \text{M}_\odot\text{pc}^{-2}$. This is slightly less than the results derived by other authors but within errors is consistent with previous estimates. Our results show no firm evidence for significant amount of dark matter in the Solar neighbourhood. However, we caution that our calculations at $|z| \leq 0.75$ kpc rely on an extrapolation from the velocity distribution function calculated at $|z| \leq 25$ pc. This extrapolation can be very sensitive to our assumption that the stellar motions are perfectly decoupled in R and z, and to our assumption of equilibrium. Indeed, we find that $\rho (z)$ within $|z|\le 0.75$ kpc is asymmetric with respect to the Galactic plane at distances $|z| = 0.1-0.4$ kpc indicating that the density distribution may be influenced by density perturbations.

Galaxy groups and clusters are the main tools used to test cosmological models and to study the environmental effect of galaxy formation. This work provides a catalogue of galaxy groups and clusters, as well as potentially merging systems based on the SDSS main galaxy survey. We identified galaxy groups and clusters using the modified friends-of-friends (FoF) group finder designed specifically for flux-limited galaxy surveys. The FoF group membership is refined by multimodality analysis to find subgroups and by using the group virial radius and escape velocity to expose unbound galaxies. We look for merging systems by comparing distances between group centres with group radii. The analysis results in a catalogue of 88662 galaxy groups with at least two members. Among them are 6873 systems with at least six members which we consider to be more reliable groups. We find 498 group mergers with up to six groups. We performed a brief comparison with some known clusters in the nearby Universe, including the Coma cluster and Abell 1750. The Coma cluster in our catalogue is a merging system with six distinguishable subcomponents. In the case of Abell 1750 we find a clear sign of filamentary infall toward this cluster. Our analysis of mass-to-light ratio (M/L) of galaxy groups reveals that M/L slightly increases with group richness.

Aims: Density waves are often considered as the triggering mechanism of star formation in spiral galaxies. Our aim is to study relations between different star formation tracers (stellar UV and near-IR radiation and emission from HI, CO and cold dust) in the spiral arms of M31, to calculate stability conditions in the galaxy disc and to draw conclusions about possible star formation triggering mechanisms. Methods: We select fourteen spiral arm segments from the de-projected data maps and compare emission distributions along the cross sections of the segments in different datasets to each other, in order to detect spatial offsets between young stellar populations and the star forming medium. By using the disc stability condition as a function of perturbation wavelength and distance from the galaxy centre we calculate the effective disc stability parameters and the least stable wavelengths at different distances. For this we utilise a mass distribution model of M31 with four disc components (old and young stellar discs, cold and warm gaseous discs) embedded within the external potential of the bulge, the stellar halo and the dark matter halo. Each component is considered to have a realistic finite thickness. Results: No systematic offsets between the observed UV and CO/far-IR emission across the spiral segments are detected. The calculated effective stability parameter has a minimal value Q_{eff} ~ 1.8 at galactocentric distances 12 - 13 kpc. The least stable wavelengths are rather long, with the minimal values starting from ~ 3 kpc at distances R > 11 kpc. Conclusions: The classical density wave theory is not a realistic explanation for the spiral structure of M31. Instead, external causes should be considered, e.g. interactions with massive gas clouds or dwarf companions of M31.

We probe the feasibility of describing the structure of a multi-component axisymmetric galaxy with a dynamical model based on the Jeans equations while taking into account a third integral of motion. We demonstrate that using the third integral in the form derived by G. Kuzmin, it is possible to calculate the stellar kinematics of a galaxy from the Jeans equations by integrating the equations along certain characteristic curves. In cases where the third integral of motion does not describe the system exactly, the derived kinematics would describe the galaxy only approximately. We apply our method to the Andromeda galaxy, for which the mass distribution is relatively firmly known. We are able to reproduce the observed stellar kinematics of the galaxy rather well. The calculated model suggests that the velocity dispersion ratios ${\sigma}_z^2/{\sigma}_R^2$ of M31 decrease with increasing R. Moving away from the galactic plane, ${\sigma}_z^2/{\sigma}_R^2$ remains the same. The velocity dispersions ${\sigma}_{\theta}^2$ and ${\sigma}_R^2$ are roughly equal in the galactic plane.

The cosmic web is a highly complex geometrical pattern, with galaxy clusters at the intersection of filaments and filaments at the intersection of walls. Identifying and describing the filamentary network is not a trivial task due to the overwhelming complexity of the structure, its connectivity and the intrinsic hierarchical nature. To detect and quantify galactic filaments we use the Bisous model, which is a marked point process built to model multi-dimensional patterns. The Bisous filament finder works directly with the galaxy distribution data and the model intrinsically takes into account the connectivity of the filamentary network. The Bisous model generates the visit map (the probability to find a filament at a given point) together with the filament orientation field. Using these two fields, we can extract filament spines from the data. Together with this paper we publish the computer code for the Bisous model that is made available in GitHub. The Bisous filament finder has been successfully used in several cosmological applications and further development of the model will allow to detect the filamentary network also in photometric redshift surveys, using the full redshift posterior. We also want to encourage the astro-statistical community to use the model and to connect it with all other existing methods for filamentary pattern detection and characterisation.

Context. Groups form the most abundant class of galaxy systems. They act as the principal drivers of galaxy evolution and can be used as tracers of the large-scale structure and the underlying cosmology. However, the detection of galaxy groups from galaxy redshift survey data is hampered by several observational limitations. Aims. We improve the widely used friends-of-friends (FoF) group finding algorithm with membership refinement procedures and apply the method to a combined dataset of galaxies in the local Universe. A major aim of the refinement is to detect subgroups within the FoF groups, enabling a more reliable suppression of the fingers-of-God effect. Methods. The FoF algorithm is often suspected of leaving subsystems of groups and clusters undetected. We used a galaxy sample built of the 2MRS, CF2, and 2M++ survey data comprising nearly 80000 galaxies within the local volume of 430 Mpc radius to detect FoF groups. We conducted a multimodality check on the detected groups in search for subgroups. We furthermore refined group membership using the group virial radius and escape velocity to expose unbound galaxies. We used the virial theorem to estimate group masses. Results. The analysis results in a catalogue of 6282 galaxy groups in the 2MRS sample with two or more members, together with their mass estimates. About half of the initial FoF groups with ten or more members were split into smaller systems with the multimodality check. An interesting comparison to our detected groups is provided by another group catalogue that is based on similar data but a completely different methodology. Two thirds of the groups are identical or very similar. Differences mostly concern the smallest and largest of these other groups, the former sometimes missing and the latter being divided into subsystems in our catalogue.

Because of the 3D nature of galaxies, an algorithm for constructing spatial density distribution models of galaxies on the basis of galaxy images has many advantages over surface density distribution approximations. We present a method for deriving spatial structure and overall parameters of galaxies from images and estimate its accuracy and derived parameter degeneracies on a sample of idealised model galaxies. The test galaxies consist of a disc-like component and a spheroidal component with varying proportions and properties. Both components are assumed to be axially symmetric and coplanar. We simulate these test galaxies as if observed in the SDSS project through ugriz filters, thus gaining a set of realistically imperfect images of galaxies with known intrinsic properties. These artificial SDSS galaxies were thereafter remodelled by approximating the surface brightness distribution with a 2D projection of a bulge+disc spatial distribution model and the restored parameters were compared to the initial ones. Down to the r-band limiting magnitude 18, errors of the restored integral luminosities and colour indices remain within 0.05 mag and errors of the luminosities of individual components within 0.2 mag. Accuracy of the restored bulge-to-disc ratios (B/D) is within 40% in most cases, and becomes worse for galaxies with low B/D, but the general balance between bulges and discs is not shifted systematically. Assuming that the intrinsic disc axial ratio is < 0.3, the inclination angles can be estimated with errors < 5deg for most of the galaxies with B/D < 2 and with errors < 15deg up to B/D = 6. Errors of the recovered sizes of the galactic components are below 10% in most cases. In general, models of disc components are more accurate than models of spheroidal components for geometrical reasons.

The accretion of satellites onto central galaxies along vast cosmic filaments is an apparent outcome of the anisotropic collapse of structure in our Universe. Numerical work (based on gravitational dynamics of N-body simulations) indicates that satellites are beamed towards hosts along preferred directions imprinted by the velocity shear field. Here we use the Sloan Digital Sky Survey to observationally test this claim. We construct 3D filaments and sheets and examine the relative position of satellite galaxies. A statistically significant alignment between satellite galaxy position and filament axis in observations is confirmed. We find a qualitatively compatible alignments by examining satellites and filaments similarly identified in the Millennium simulation, semi-analytical galaxy catalogue. We also examine the dependence of the alignment strength on galaxy properties such as colour, magnitude and (relative) satellite magnitude, finding that the alignment is strongest for the reddest and brightest central and satellite galaxies. Our results confirm the theoretical picture and the role of the cosmic web in satellite accretion. Furthermore our results suggest that filaments identified on larger scales can be reflected in the positions of satellite galaxies that are quite close to their hosts.

Context. Galaxies in the Universe form chains (filaments) that connect groups and clusters of galaxies. The filamentary network includes nearly half of the galaxies and is visually the most striking feature in cosmological maps. Aims. We study the distribution of galaxies along the filamentary network, trying to find specific patterns and regularities. Methods. Galaxy filaments are defined by the Bisous model, a marked point process with interactions. We use the two-point correlation function and the Rayleigh Z-squared statistic to study how galaxies and galaxy groups are distributed along the filaments. Results. We show that galaxies and groups are not uniformly distributed along filaments, but tend to form a regular pattern. The characteristic length of the pattern is around 7 Mpc/h. A slightly smaller characteristic length 4 Mpc/h can also be found, using the Z-squared statistic. Conclusions. We find that galaxy filaments in the Universe are like pearl necklaces, where the pearls are galaxy groups distributed more or less regularly along the filaments. We propose that this well defined characteristic scale could be used to test various cosmological models and to probe environmental effects on the formation and evolution of galaxies.

We provide flux-limited and volume-limited galaxy group and cluster catalogues, based on the spectroscopic sample of the SDSS data release 10 galaxies. We used a modified friends-of-friends (FoF) method with a variable linking length in the transverse and radial directions to identify as many realistic groups as possible. The flux-limited catalogue incorporates galaxies down to m_r = 17.77 mag. It includes 588193 galaxies and 82458 groups. The volume-limited catalogues are complete for absolute magnitudes down to M_r = -18.0, -18.5, -19.0, -19.5, -20.0, -20.5, and -21.0; the completeness is achieved within different spatial volumes, respectively. Our analysis shows that flux-limited and volume-limited group samples are well compatible to each other, especially for the larger groups/clusters. Dynamical mass estimates, based on radial velocity dispersions and group extent in the sky, are added to the extracted groups. The catalogues can be accessed via http://cosmodb.to.ee and the Strasbourg Astronomical Data Center (CDS).

Evolution of galaxies is one of the most actual topics in astrophysics. Among the most important factors determining the evolution are two galactic components which are difficult or even impossible to detect optically: the gaseous disks and the dark matter halo. We use deep Hubble Space Telescope images to construct a two-component (bulge + disk) model for stellar matter distribution of galaxies. Properties of the galactic components are derived using a three-dimensional galaxy modeling software, which also estimates disk thickness and inclination angle. We add a gas disk and a dark matter halo and use hydrodynamical equations to calculate gas rotation and dispersion profiles in the resultant gravitational potential. We compare the kinematic profiles with the Team Keck Redshift Survey observations. In this pilot study, two galaxies are analyzed deriving parameters for their stellar components; both galaxies are found to be disk-dominated. Using the kinematical model, the gas mass and stellar mass ratio in the disk are estimated.