We study the properties, connectivity, and galaxy content of groups and filaments in the low-density region (cocoon) around A2142 supercluster (SClA2142). We traced the SClA2142 cocoon boundaries by the lowest luminosity-density regions that separate SClA2142 from other superclusters. We determined galaxy filaments and groups in the cocoon and analysed the connectivity of groups, the high density core (HDC) of the supercluster, and the whole of the supercluster. We compared the distribution and properties of galaxies with different star-formation properties in the supercluster and in the cocoon. SClA2142 and the long filament that is connected to it forms the longest straight structure in the Universe detected so far, with a length of $75$ $h^{-1}$ Mpc. The connectivity of the supercluster is C = 6 - 7; poor groups have C = 1 - 2. Long filaments around the supercluster's main body are detached from it at the turnaround region. Galaxies with very old stellar populations lie in systems across a wide range of richness from the richest cluster to poorest groups and single galaxies. They lie even at local densities as low as $D1 < 1$ in the cocoon and up to $D1 > 800$ in the supercluster. Recently quenched galaxies lie in the cocoon mainly in one region and their properties are different in the cocoon and in the supercluster. The star-formation properties of single galaxies are similar across all environments. The collapsing main body of SClA2142 with the detached long filaments near it are evidence of an important epoch in the supercluster evolution. Further studies are needed to understand the reasons of similarity of galaxies with very old stellar populations in extremely different environments. The presence of long, straight structures in the cosmic web may serve as a test for cosmological models.

The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

Large multi-object spectroscopic surveys require automated algorithms to optimise their observing strategy. One of the most ambitious upcoming spectroscopic surveys is the 4MOST survey. The 4MOST survey facility is a fibre-fed spectroscopic instrument on the VISTA telescope with a large enough field of view to survey a large fraction of the southern sky within a few years. Several Galactic and extragalactic surveys will be carried out simultaneously, so the combined target density will strongly vary. In this paper, we describe a new tiling algorithm that can naturally deal with the large target density variations on the sky and which automatically handles the different exposure times of targets. The tiling pattern is modelled as a marked point process, which is characterised by a probability density that integrates the requirements imposed by the 4MOST survey. The optimal tilling pattern with respect to the defined model is estimated by the tiles configuration that maximises the proposed probability density. In order to achieve this maximisation a simulated annealing algorithm is implemented. The algorithm automatically finds an optimal tiling pattern and assigns a tentative sky brightness condition and exposure time for each tile, while minimising the total execution time that is needed to observe the list of targets in the combined input catalogue of all surveys. Hence, the algorithm maximises the long-term observing efficiency and provides an optimal tiling solution for the survey. While designed for the 4MOST survey, the algorithm is flexible and can with simple modifications be applied to any other multi-object spectroscopic survey.

Context. The cosmic web, a complex network of galaxy groups and clusters connected by filaments, is a dynamical environment in which galaxies form and evolve. However, the impact of cosmic filaments on the properties of galaxies is difficult to study because of the much more influential local (galaxy-group scale) environment. Aims. The aim of this paper is to investigate the dependence of intrinsic galaxy properties on distance to the nearest cosmic web filament, using a sample of galaxies for which the local environment is easily assessable.} Methods. Our study is based on a volume-limited galaxy sample with $M_\mathrm{r}$ $\leq -19$ mag, drawn from the SDSS DR12. We chose brightest group galaxies (BGGs) in groups with two to six members as our probes of the impact of filamentary environment because their local environment can be determined more accurately. We use the Bisous marked point process method to detect cosmic-web filaments with radii of $0.5-1.0$ Mpc and measure the perpendicular filament spine distance ($D_{\mathrm{fil}}$) for the BGGs. We limit our study to $D_{\mathrm{fil}}$ values up to 4 Mpc. We use the luminosity density field as a tracer of the local environment. To achieve uniformity of the sample and to reduce potential biases we only consider filaments longer than 5 Mpc. Our final sample contains 1427 BGGs. Results. We note slight deviations between the galaxy populations inside and outside the filament radius in terms of stellar mass, colour, the 4000AA break, specific star formation rates, and morphologies. However, all these differences remain below 95% confidence and are negligible compared to the effects arising from local environment density. Conclusions. Within a 4 Mpc radius of the filament axes, the effect of filaments on BGGs is marginal. The local environment is the main factor in determining BGG properties.

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.

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

We explore the high spectral resolution X-ray data towards the quasar 3C273 to search for signals of hot ($\sim10^{6-7}$ K) X-ray-absorbing gas co-located with two established intergalactic FUV OVI absorbers. We analyze the soft X-ray band grating data of all XMM-Newton and Chandra instruments to search for the hot phase absorption lines at the FUV predicted redshifts. The viability of potential line detections is examined by adopting the constraints of a physically justified absorption model. The WHIM hypothesis is investigated with a complementary 3D galaxy distribution analysis, and by comparison of the measurement results to the WHIM properties in the EAGLE cosmological, hydrodynamical simulation. At FUV redshift z=0.09017, we measured signals of two hot ion species, OVIII and NeIX, with a $3.9\sigma$ combined significance level. Considering the line features in all instruments collectively and assuming collisional equilibrium for absorbing gas, we were able to constrain the temperature ($kT=0.26\pm0.03$ keV) and the column density ($N_H\times{Z_\odot/Z}=1.3_{-0.5}^{+0.6}\times10^{19}$ cm$^{-2}$) of the absorber. Thermal analysis indicates that FUV and X-ray absorption relate to different phases, with estimated temperatures $T_{FUV}\approx3\times10^5$ and $T_{X-ray}\approx3\times10^6$ K, which match the EAGLE predictions for WHIM at the FUV/X-ray measured $N_{ion}$-ranges. We detected a large scale galactic filament crossing the sightline at the redshift of the absorption, linking the absorption to this structure. This study provides insights into co-existing warm and hot gas within a WHIM filament and estimates the ratio of the hot and warm phases. Because the hot phase is thermally distinct from the OVI gas, the estimated baryon content of the absorber is increased, conveying the promise of X-ray follow-up studies of FUV detected WHIM in refining the picture of the missing baryons.

The intergalactic medium (IGM) plays an important role in the formation and evolution of galaxies. Recent developments in upcoming radio telescopes are starting to open up the possibility of making a first direct detection of the 21 cm signal of neutral hydrogen (HI) from the warm gas of the IGM in large-scale filaments. The cosmological hydrodynamical EAGLE simulation is used to estimate the typical IGM filament signal. Assuming the same average signal for all filaments, a prediction is made for the detectability of such a signal with the upcoming mid-frequency array of the Square Kilometer Array (SKA1-mid) or the future upgrade to SKA2. The signal-to-noise (S/N) then only depends on the size and orientation of each filament. With filament spines inferred from existing galaxy surveys as a proxy for typical real filaments, we find hundreds of filaments in the region of the sky accessible to the SKA that can be detected. Once the various phases of the SKA telescope become operational, their own surveys will be able to find the galaxies required to infer the position of even more filaments within the survey area. We find that in 120 h, SKA1-mid/SKA2 will detect HI emission from the strongest filaments in the field with a S/N of the order of 10 to $\sim$150 for the most pessimistic model considered here. Some of the brighter filaments can be detected with an integration time of a few minutes with SKA1-mid and a few seconds with SKA2. Therefore, SKA2 will be capable of not only detecting but also mapping a large part of the IGM in these filaments.

In a framework where galaxies mostly migrate on the colour-magnitude diagram from star-forming to quiescent, the green valley is considered a transitional galaxy stage. The details of the processes that drive galaxies from star-forming to passive systems still remain unknown. We developed a method that estimates empirically the star formation quenching times-scales of green valley galaxies, assuming an exponential decay model of the SFH and through a combination of narrow and broad bands from J-PLUS and GALEX. We correlate these quenching time-scales with the presence of bars. We find that the J-PLUS colours F0395-g and F0415-g are sensitive to different SFH, showing, a clear correlation with the Dn(4000) and H-delta,A spectral indices. We find that quenching time-scales obtained with our new approach are in agreement with those determined using spectral indices. We also find that galaxies with high bar probability tend to quench their star formation slowly. We conclude that: 1) J-PLUS filters can be used to measure quenching timescales in nearby green valley galaxies; and 2) the resulting star formation quenching time-scales are longer for barred green valley galaxies. Considering that the presence of a bar indicates that more violent processes (e.g., major mergers) are absent in host galaxies, we conclude that the presence of a bar can be used as a morphological signature for slow star formation quenching.

It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the $\Lambda$CDM model. One type of anisotropy recently detected in the SDSS (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called "lopsided". Namely, in pairs of galaxies (like Andromeda and the Milky Way) the satellites are more likely to inhabit the region in between the pair, rather than on opposing sides. Although recent studies found a similar set up when comparing pairs of galaxies in $\Lambda$CDM simulations indicating that such a set up is not inconsistent with $\Lambda$CDM, the origin has yet to be explained. Here we examine the origin of such lopsided setups by first identifying such distributions in pairs of galaxies in numerical cosmological simulations, and then tracking back the orbital trajectories of satellites (which at $z=0$ display the effect). We report two main results: first, the lopsided distribution was stronger in the past and weakens towards $z=0$. Second, the weakening of the signal is due to the interaction of satellite galaxies with the pair. Finally, we show that the $z=0$ signal is driven primarily by satellites that are on first approach, who have yet to experience a "flyby". This suggests that the signal seen in the observations is also dominated by dynamically young accretion events.

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.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

The 4MOST Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects over the redshift range $z = 0.15$ to $3.5$, will allow definitive tests of gravitational physics. Many key science questions will be addressed by combining CRS spectra of these targets with data from current or future facilities such as the Large Synoptic Survey Telescope, the Square Kilometre Array and the Euclid mission.

Groups and clusters of galaxies are a current focus of astronomical research owing to their role in determining the environmental effects on galaxies and the constraints they provide to cosmology. The eROSITA X-ray telescope on board the Spectrum Roentgen Gamma observatory will be launched in 2019 and will have completed eight scans of the full sky when 4MOST starts operating. The experiment will detect groups and clusters of galaxies through X-ray emission from the hot intergalactic medium. The purpose of the 4MOST eROSITA Galaxy Cluster Redshift Survey is to provide spectroscopic redshifts of the optical counterparts to the X-ray emission from 40,000 groups and clusters of galaxies so as to perform dynamical estimates of the total mass and to measure the properties of the member galaxies. The survey aims to obtain precise redshift measurements of the photometrically identified brightest cluster galaxies at redshift $z > 0.7$. At lower redshifts ($z < 0.7$), the programme aims to sample over 15 member galaxies per cluster and enable dynamical mass measurements to calibrate the clusters for cosmological experiments. At $z < 0.2$, eROSITA will also detect X-ray emission from galaxy groups and filaments. 4MOST spectroscopic data from the survey will be used for optical identification of galaxy groups down to eROSITA's mass detection limits of $10^{13} M_\odot$, as well as the detection of the largest filaments for pioneering studies of their X-ray emission.

The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surveys that have vastly different target brightnesses and densities, sample completeness levels, and signal-to-noise requirements. We introduce here a number of concepts that we will use to ensure all surveys are optimised. Astronomers who are planning to submit a Participating Survey proposal are strongly encouraged to read this article and any relevant 4MOST Survey articles in this issue of The Messenger such that they can optimally complement and benefit from the planned surveys of the 4MOST Consortium.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

It has been shown, both in simulations and observationally, that the tidal field of a large galaxy can torque its satellites such that the major axis of satellite galaxies points towards their hosts. This so-called `shape alignment' has been observed in isolated Milky Way-like galaxies but not in `Local Group'-like pairs. In this study, we investigate the shape alignment of satellite galaxies in galaxy pairs similar to the Local Group identified in the Sloan Digital Sky Survey Data Release 13 (SDSS DR13). By stacking tens of thousands of satellite galaxies around primary galaxy pairs, we find two statistically strong alignment signals. (1) The major axes of satellite galaxies located in the (projected) area between two primaries (the {\it facing} region) tend to be perpendicular to the line connecting the satellite to its host (tangential alignment), while (2) the major axes of satellite galaxies located in regions away from the other host (the {\it away} region) tend to be aligned with the line connecting the satellite to its host (radial alignment). These alignments are confirmed at $\sim5\sigma$ levels. The alignment signal increases with increasing primary brightness, decreasing pair separation, and decreasing satellite distance. The alignment signal is also found to be stronger in filamentary environments. These findings will shed light on understanding the mechanisms of how satellite galaxies are affected by the tidal field in galaxy pairs and will be useful for investigating galaxy intrinsic alignment in the analyses of weak gravitational lensing.

We investigate the alignment of galaxies and haloes relative to cosmic web filaments using the EAGLE hydrodynamical simulation. We identify filaments by applying the NEXUS+ method to the mass distribution and the Bisous formalism to the galaxy distribution. Both web finders return similar filamentary structures that are well aligned and that contain comparable galaxy populations. EAGLE haloes have an identical spin alignment with filaments as their counterparts in dark matter only simulations: a complex mass dependent trend with low mass haloes spinning preferentially parallel to and high mass haloes spinning preferentially perpendicular to filaments. In contrast, galaxy spins do not show such a spin transition and have a propensity for perpendicular alignments at all masses, with the degree of alignment being largest for massive galaxies. This result is valid for both NEXUS+ and Bisous filaments. When splitting by morphology, we find that elliptical galaxies show a stronger orthogonal spin--filament alignment than spiral galaxies of similar mass. The same is true of their haloes, with the host haloes of elliptical galaxies having a larger degree of orthogonal alignment than the host haloes of spirals. Due to the misalignment between galaxy shape and spin, galaxy minor axes are oriented differently with filaments than galaxy spins. We find that the galaxies whose minor axis is perpendicular to a filament are much better aligned with their host haloes. This suggests that many of the same physical processes determine both the galaxy--filament and the galaxy--halo alignments.

The cosmological missing baryons at z<1 most likely hide in the hot (T$\gtrsim10^{5.5}$ K) phase of the Warm Hot Intergalactic Medium (WHIM). While the hot WHIM is hard to detect due to its high ionisation level, the warm (T$\lesssim10^{5.5}$ K) phase of the WHIM has been very robustly detected in the FUV band. We adopted the assumption that the hot and warm WHIM phases are co-located and thus used the FUV-detected warm WHIM as a tracer for the cosmologically interesting hot WHIM. We utilised the assumption by performing an X-ray follow-up in the sight line of a blazar PKS 2155-304 at the redshifts where previous FUV measurements of OVI, SiIV and BLA absorption have indicated the existence of the warm WHIM. We looked for the OVII He$\alpha$ and OVIII Ly$\alpha$ absorption lines, the most likely hot WHIM tracers. Despite of the very large exposure time ($\approx$ 1 Ms), the XMM-Newton/RGS1 data yielded no significant detection which corresponds to upper limits of $\log{N({\rm OVII})({\rm cm}^{-2}))} \le 14.5-15.2$ and $\log{N({\rm OVIII})({\rm cm}^{-2}))} \le 14.9-15.2$. An analysis of LETG/HRC data yielded consistent results. However, the LETG/ACIS data yielded a detection of an absorption line - like feature at $\lambda \approx$ 20 \AA\ at simple one parameter uncertainty - based confidence level of 3.7 $\sigma$, consistently with several earlier LETG/ACIS reports. Given the high statistical quality of the RGS1 data, the possibility of RGS1 accidentally missing the true line at $\lambda \sim$ 20 \AA\ is very low, 0.006%. Neglecting this, the LETG/ACIS detection can be interpreted as Ly$\alpha$ transition of OVIII at one of the redshifts (z$\approx$ 0.054) of FUV-detected warm WHIM. Given the very convincing X-ray spectral evidence for and against the existence of the $\lambda \sim$ 20 \AA\ feature, we cannot conclude whether or not it is a true astrophysical absorption line.

We study the dynamical state and properties of galaxies and groups in the supercluster SClA2142 that has a collapsing core, to understand its possible formation and evolution. We find the substructure of galaxy groups using normal mixture modelling. We have used the projected phase space (PPS) diagram, spherical collapse model, clustercentric distances, and magnitude gap between the brightest galaxies in groups to study the dynamical state of groups and to analyse group and galaxy properties. We compared the alignments of groups and their brightest galaxies with the supercluster axis. The supercluster core has a radius of about $8 h^{-1}$Mpc and total mass $M_{\mathrm{tot}} \approx 2.3\times~10^{15}h^{-1}M_\odot$ and is collapsing. Galaxies in groups on the supercluster axis have older stellar populations than off-axis groups, with median stellar ages $4 - 6$ and $< 4$Gyr, correspondingly. The cluster A2142 and the group Gr8 both host galaxies with the oldest stellar populations among groups in SClA2142 having the median stellar age $t > 8$Gyr. Recently quenched galaxies and active galactic nuclei (AGNs) are mostly located at virial radii or in merging regions of groups, and at clustercentric distances $D_c \approx 6 h^{-1}$ Mpc. The most elongated groups lie along the supercluster axis and are aligned with it. Magnitude gaps between the brightest galaxies of groups are less than one magnitude, suggesting that groups in SClA2142 are dynamically young. The collapsing core of the supercluster, infall of galaxies and groups, and possible merging groups, which affect galaxy properties and may trigger the activity of AGNs, show how the whole supercluster is evolving.

We investigate the alignment of haloes with the filaments of the cosmic web using an unprecedently large sample of dark matter haloes taken from the P-Millennium $\Lambda$CDM cosmological N-body simulation. We use the state-of-the-art NEXUS morphological formalism which, due to its multiscale nature, simultaneously identifies structures at all scales. We find strong and highly significant alignments, with both the major axis of haloes and their peculiar velocity tending to orient along the filament. However, the spin - filament alignment displays a more complex trend changing from preferentially parallel at low masses to preferentially perpendicular at high masses. This "spin flip" occurs at an average mass of $5\times10^{11}~h^{-1}M_\odot$. This mass increases with increasing filament diameter, varying by more than an order of magnitude between the thinnest and thickest filament samples. We also find that the inner parts of haloes have a spin flip mass that is several times smaller than that of the halo as a whole. These results confirm that recent accretion is responsible for the complex behaviour of the halo spin - filament alignment. Low-mass haloes mainly accrete mass along directions perpendicular to their host filament and thus their spins tend to be oriented along the filaments. In contrast, high-mass haloes mainly accrete along their host filaments and have their spins preferentially perpendicular to them. Furthermore, haloes located in thinner filaments are more likely to accrete along their host filaments than haloes of the same mass located in thicker filaments.

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.

The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts.

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)

The alignment between satellites and central galaxies has been studied in detail both in observational and theoretical works. The widely accepted fact is that the satellites preferentially reside along the major axis of their central galaxy. However, the origin and large-scale environment dependence of this alignment are still unknown. In an attempt to figure out those, we use data constructed from SDSS DR7 to investigate the large-scale environmental dependence of this alignment with emphasis on examining the alignments' dependence on the colour of the central galaxy. We find a very strong large-scale environmental dependence of the satellite-central alignment in groups with blue centrals. Satellites of blue centrals in knots are preferentially located perpendicular to the major axis of the centrals, and the alignment angle decreases with environment namely when going from knots to voids. The alignment angle strongly depend on the ${}^{0.1}(g-r)$ colour of centrals. We suggest that the satellite-central alignment is the result of a competition between satellite accretion within large scale-structure and galaxy evolution inside host haloes. For groups containing red central galaxies, the satellite-central alignment is mainly determined by the evolution effect, while for blue central dominated groups, the effect of large-scale structure plays a more important role, especially in knots. Our results provide an explanation for how the satellite-central alignment forms within different large-scale environments. The perpendicular case in groups and knots with blue centrals may also provide insight into understanding similar polar arrangements such the formation of the Milky Way and Centaurus A's satellite system.

To study the full formation and evolution history of galaxy clusters and their population, high resolution simulations of the latter are flourishing. However comparing observed clusters to the simulated ones on a one-to-one basis to refine the models and theories down to the details is non trivial. The large variety of clusters limits the comparisons between observed and numerical clusters. Simulations resembling the local Universe down to the cluster scales permit pushing the limit. Simulated and observed clusters can be matched on a one-to-one basis for direct comparisons provided that clusters are well reproduced besides being in the proper large scale environment. Comparing random and local-Universe like simulations obtained with differently grouped observational catalogs of peculiar velocities, this paper shows that the grouping scheme used to remove non-linear motions in the catalogs that constrain the simulations affects the quality of the numerical clusters. With a less aggressive grouping scheme - galaxies still falling onto clusters are preserved - combined with a bias minimization scheme, the mass of the dark matter halos, simulacra for 5 local clusters - Virgo, Centaurus, Coma, Hydra and Perseus - is increased by 39% closing the gap with observational mass estimates. Simulacra are found on average in 89% of the simulations, an increase of 5% with respect to the previous grouping scheme. The only exception is Perseus. Since the Perseus-Pisces region is not well covered by the used peculiar velocity catalog, the latest release let us foresee a better simulacrum for Perseus in a near future.

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.

With the advent of wide-field cosmological surveys, we are approaching samples of hundreds of thousands of galaxy clusters. While such large numbers will help reduce statistical uncertainties, the control of systematics in cluster masses becomes ever more crucial. Here we examine the effects of an important source of systematic uncertainty in galaxy-based cluster mass estimation techniques: the presence of significant dynamical substructure. Dynamical substructure manifests as dynamically distinct subgroups in phase-space, indicating an 'unrelaxed' state. This issue affects around a quarter of clusters in a generally selected sample. We employ a set of mock clusters whose masses have been measured homogeneously with commonly-used galaxy-based mass estimation techniques (kinematic, richness, caustic, radial methods). We use these to study how the relation between observationally estimated and true cluster mass depends on the presence of substructure, as identified by various popular diagnostics. We find that the scatter for an ensemble of clusters does not increase dramatically for clusters with dynamical substructure. However, we find a systematic bias for all methods, such that clusters with significant substructure have higher measured masses than their relaxed counterparts. This bias depends on cluster mass: the most massive clusters are largely unaffected by the presence of significant substructure, but masses are significantly overestimated for lower mass clusters, by $\sim10\%$ at $10^{14}$ and $\geq20\%$ for $\leq10^{13.5}$. The use of cluster samples with different levels of substructure can, therefore, bias certain cosmological parameters up to a level comparable to the typical uncertainties in current cosmological studies.

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.

The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.

High quality reconstructions of the three dimensional velocity and density fields of the local Universe are essential to study the local Large Scale Structure. In this paper, the Wiener Filter reconstruction technique is applied to galaxy radial peculiar velocity catalogs to understand how the Hubble constant (H0) value and the grouping scheme affect the reconstructions. While H0 is used to derive radial peculiar velocities from galaxy distance measurements and total velocities, the grouping scheme serves the purpose of removing non linear motions. Two different grouping schemes (based on the literature and a systematic algorithm) as well as five H0 values ranging from 72 to 76 km/s/Mpc are selected. The Wiener Filter is applied to the resulting catalogs. Whatever grouping scheme is used, the larger H0 is, the larger the infall onto the local Volume is. However, this conclusion has to be strongly mitigated: a bias minimization scheme applied to the catalogs after grouping suppresses this effect. At fixed H0, reconstructions obtained with catalogs grouped with the different schemes exhibit structures at the proper location in both cases but the latter are more contrasted in the less aggressive scheme case: having more constraints permits an infall from both sides onto the structures to reinforce their overdensity. Such findings highlight the importance of a balance between grouping to suppress non linear motions and preserving constraints to produce an infall onto structures expected to be large overdensities. Such an observation is promising to perform constrained simulations of the local Universe including its massive clusters.

We have analyzed all available archival XMM-Newton observations of X Comae, a bright X-ray quasar behind the Coma cluster, to study the properties of the warm-hot intergalactic medium in the vicinity of the nearest massive galaxy cluster. The RGS observations confirm the possible presence of a Ne IX K-alpha absorption line at the redshift of Coma, although with a limited statistical significance. This analysis is therefore in line with the earlier analysis by Takei et al. (2007) based on a sub-set of these data. Its large column density and optical depth, however, point to implausible conditions for the absorbing medium, thereby casting serious doubts to its reality. Chandra has never observed X Comae and therefore cannot provide additional information on this source. We combine upper limits to the presence of other X-ray absorption lines (notably from O VII and O VIII) at the redshift of Coma with positive measurements of the soft excess emission from Coma measured by ROSAT (Bonamente et al. 2003). The combination of emission from warm-hot gas at kT~1/4 keV and upper limits from absorption lines provide useful constraints on the density and the sightline length of the putative WHIM towards Coma. We conclude that the putative warm-hot medium towards Coma is consistent with expected properties, with a baryon overdensity >=10 and a sightline extent of order of tens of Mpc.

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.

We study the morphology, luminosity and mass of the superclusters from the BOSS Great Wall (BGW), a recently discovered very rich supercluster complex at the redshift $z = 0.47$. We have employed the Minkowski functionals to quantify supercluster morphology. We calculate supercluster luminosities and masses using two methods. Firstly, we used data about the luminosities and stellar masses of high stellar mass galaxies with $\log(M_*/h^{-1}M_\odot) \geq 11.3$. Secondly, we applied a scaling relation that combines morphological and physical parameters of superclusters to obtain supercluster luminosities, and obtained supercluster masses using the mass-to-light ratios found for local rich superclusters. We find that the BGW superclusters are very elongated systems, with shape parameter values of less than $0.2$. This value is lower than that found for the most elongated local superclusters. The values of the fourth Minkowski functional $V_3$ for the richer BGW superclusters ($V_3 = 7$ and $10$) show that they have a complicated and rich inner structure. We identify two Planck SZ clusters in the BGW superclusters, one in the richest BGW supercluster, and another in one of the poor BGW superclusters. The luminosities of the BGW superclusters are in the range of $1 - 8\times~10^{13}h^{-2}L_\odot$, and masses in the range of $0.4 - 2.1\times~10^{16}h^{-1}M_\odot$. Supercluster luminosities and masses obtained with two methods agree well. We conclude that the BGW is a complex of massive, luminous and large superclusters with very elongated shape. The search and detailed study, including the morphology analysis of the richest superclusters and their complexes from observations and simulations can help us to understand formation and evolution of the cosmic web.

We investigate the impact of filament and void environments on galaxies, looking for residual effects beyond the known relations with environment density. We quantified the host environment of galaxies as the distance to the spine of the nearest filament, and compared various galaxy properties within 12 bins of this distance. We considered galaxies up to 10 $h^{-1}$Mpc from filaments, i.e. deep inside voids. The filaments were defined by a point process (the Bisous model) from the Sloan Digital Sky Survey data release 10. In order to remove the dependence of galaxy properties on the environment density and redshift, we applied weighting to normalise the corresponding distributions of galaxy populations in each bin. After the normalisation with respect to environment density and redshift, several residual dependencies of galaxy properties still remain. Most notable is the trend of morphology transformations, resulting in a higher elliptical-to-spiral ratio while moving from voids towards filament spines, bringing along a corresponding increase in the $g-i$ colour index and a decrease in star formation rate. After separating elliptical and spiral subsamples, some of the colour index and star formation rate evolution still remains. The mentioned trends are characteristic only for galaxies brighter than about $M_{r} = -20$ mag. Unlike some other recent studies, we do not witness an increase in the galaxy stellar mass while approaching filaments. The detected transformations can be explained by an increase in the galaxy-galaxy merger rate and/or the cut-off of extragalactic gas supplies (starvation) near and inside filaments. Unlike voids, large-scale galaxy filaments are not a mere density enhancement, but have their own specific impact on the constituent galaxies, reducing the star formation rate and raising the chances of elliptical morphology also at a fixed environment density level.

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.

Our goal is to study the evolution of the $B-$band luminosity function (LF) since $z=1$ using ALHAMBRA data. We used the photometric redshift and the $I-$band selection magnitude probability distribution functions (PDFs) of those ALHAMBRA galaxies with $I\leq24$ mag to compute the posterior LF. We statistically studied quiescent and star-forming galaxies using the template information encoded in the PDFs. The LF covariance matrix in redshift-magnitude-galaxy type space was computed, including the cosmic variance. That was estimated from the intrinsic dispersion of the LF measurements in the 48 ALHAMBRA sub-fields. The uncertainty due to the photometric redshift prior is also included in our analysis. We modelled the LF with a redshift-dependent Schechter function affected by the same selection effects than the data. The measured ALHAMBRA LF at $0.2\leq z<1$ and the evolving Schechter parameters both for quiescent and star-forming galaxies agree with previous results in the literature. The estimated redshift evolution of $M_B^* \propto Qz$ is $Q_{\rm SF}=-1.03\pm0.08$ and $Q_{\rm Q}=-0.80\pm0.08$, and of $\log \phi^* \propto Pz$ is $P_{\rm SF}=-0.01\pm0.03$ and $P_{\rm Q}=-0.41\pm0.05$. The measured faint-end slopes are $\alpha_{\rm SF}=-1.29\pm0.02$ and $\alpha_{\rm Q}=-0.53\pm0.04$. We find a significant population of faint quiescent galaxies, modelled by a second Schechter function with slope $\beta=-1.31\pm0.11$. We find a factor $2.55\pm0.14$ decrease in the luminosity density $j_B$ of star-forming galaxies, and a factor $1.25\pm0.16$ increase in the $j_B$ of quiescent ones since $z=1$, confirming the continuous build-up of the quiescent population with cosmic time. The contribution of the faint quiescent population to $j_B$ increases from 3% at $z=1$ to 6% at $z=0$. The developed methodology will be applied to future multi-filter surveys such as J-PAS.

We study the alignment of galaxies relative to their local environment in SDSS-DR8 and, using these data, we discuss evolution scenarios for different types of galaxies. We defined a vector field of the direction of anisotropy of the local environment of galaxies. We summed the unit direction vectors of all close neighbours of a given galaxy in a particular way to estimate this field. We found the alignment angles between the spin axes of disc galaxies, or the minor axes of elliptical galaxies, and the direction of anisotropy. The distributions of cosines of these angles are compared to the random distributions to analyse the alignment of galaxies. Sab galaxies show perpendicular alignment relative to the direction of anisotropy in a sparse environment, for single galaxies and galaxies of low luminosity. Most of the parallel alignment of Scd galaxies comes from dense regions, from 2...3 member groups and from galaxies with low luminosity. The perpendicular alignment of S0 galaxies does not depend strongly on environmental density nor luminosity; it is detected for single and 2...3 member group galaxies, and for main galaxies of 4...10 member groups. The perpendicular alignment of elliptical galaxies is clearly detected for single galaxies and for members of < 11 member groups; the alignment increases with environmental density and luminosity. We confirm the existence of fossil tidally induced alignment of Sab galaxies at low z. The alignment of Scd galaxies can be explained via the infall of matter to filaments. S0 galaxies may have encountered relatively massive mergers along the direction of anisotropy. Major mergers along this direction can explain the alignment of elliptical galaxies. Less massive, but repeated mergers are possibly responsible for the formation of elliptical galaxies in sparser areas and for less luminous elliptical galaxies.

The nature versus nurture scenario in galaxy and group evolution is a long-standing problem not yet fully understood on cosmological scales. We study the properties of groups and their central galaxies in different large-scale environments defined by the luminosity density field and the cosmic web filaments. We use the luminosity density field constructed using 8 Mpc/h smoothing to characterize the large-scale environments and the Bisous model to extract the filamentary structures in different large-scale environments. We find differences in the properties of central galaxies and their groups in and outside of filaments at fixed halo and large-scale environments. In high-density environments, the group mass function has higher number densities in filaments compared to that outside of filaments towards the massive end. The relation is opposite in low-density environments. At fixed group mass and large-scale luminosity density, groups in filaments are slightly more luminous and their central galaxies have redder colors, higher stellar masses, and lower specific star formation rates than those outside of filaments. However, the differences in central galaxy and group properties in and outside of filaments are not clear in some group mass bins. We show that the differences in central galaxy properties are due to the higher abundances of elliptical galaxies in filaments. Filamentary structures in the cosmic web are not simply visual associations of galaxies, but rather play an important role in shaping the properties of groups and their central galaxies. The differences in central galaxy and group properties in and outside of cosmic web filaments are not simple effects related to large-scale environmental density. The results point towards an efficient mechanism in cosmic web filaments which quench star formation and transform central galaxy morphology from late to early types.

The radiation of stars heats dust grains in the diffuse interstellar medium and in star-forming regions in galaxies. Modelling this interaction provides information on dust in galaxies, a vital ingredient for their evolution. It is not straightforward to identify the stellar populations heating the dust, and to link attenuation to emission on a sub-galactic scale. Radiative transfer models are able to simulate this dust-starlight interaction in a realistic, three-dimensional setting. We investigate the dust heating mechanisms on a local and global galactic scale, using the Andromeda galaxy (M31) as our laboratory. We perform a series of panchromatic radiative transfer simulations of Andromeda with our code SKIRT. The high inclination angle of M31 complicates the 3D modelling and causes projection effects. However, the observed morphology and flux density are reproduced fairly well from UV to sub-millimeter wavelengths. Our model reveals a realistic attenuation curve, compatible with previous, observational estimates. We find that the dust in M31 is mainly (91 % of the absorbed luminosity) heated by the evolved stellar populations. The bright bulge produces a strong radiation field and induces non-local heating up to the main star-forming ring at 10 kpc. The relative contribution of unevolved stellar populations to the dust heating varies strongly with wavelength and with galactocentric distance.The dust heating fraction of unevolved stellar populations correlates strongly with NUV-r colour and specific star formation rate. These two related parameters are promising probes for the dust heating sources at a local scale.

In the cosmic web, galaxy superclusters or their high-density cores are the largest objects that may collapse at present or during the future evolution. We study the dynamical state and possible future evolution of galaxy superclusters from the Sloan Great Wall (SGW), the richest galaxy system in the nearby Universe. We calculated supercluster masses using dynamical masses of galaxy groups and stellar masses of galaxies. We employed normal mixture modelling to study the structure of rich SGW superclusters and search for components (cores) in superclusters. We analysed the radial mass distribution in the high-density cores of superclusters centred approximately at rich clusters and used the spherical collapse model to study their dynamical state. We found that the lower limit of the total mass of the SGW is approximately $M = 2.5\times~10^{16}h^{-1}M_\odot$. Different mass estimators of superclusters agree well, the main uncertainties in masses of superclusters come from missing groups and clusters. We detected three high-density cores in the richest SGW supercluster (SCl~027) and two in the second richest supercluster (SCl~019). They have masses of $1.2 - 5.9 \times~10^{15}h^{-1}M_\odot$ and sizes of up to $\approx 60 h^{-1}$ Mpc. The high-density cores of superclusters are very elongated, flattened perpendicularly to the line of sight. The comparison of the radial mass distribution in the high-density cores with the predictions of spherical collapse model suggests that their central regions with radii smaller than $8 h^{-1}$Mpc and masses of up to $M = 2\times~10^{15}h^{-1}M_\odot$ may be collapsing. The rich SGW superclusters with their high-density cores represent dynamically evolving environments for studies of the properties of galaxies and galaxy systems.

The distribution of smaller satellite galaxies around large central galaxies has attracted attention because peculiar spatial and kinematic configurations have been detected in some systems. A particularly striking example of such behavior is seen in the satellite system of the Andromeda galaxy, where around 80\% are on the nearside of that galaxy, facing the Milky Way. Motivated by this departure from anisotropy, we examined the spatial distribution of satellites around pairs of galaxies in the SDSS. By stacking tens of thousands of satellites around galaxy pairs we found that satellites tend to bulge towards the other central galaxy, preferably occupying the space between the pair, rather than being spherically or axis-symmetrically distributed around each host. The bulging is a function of the opening angle examined and is fairly strong -- there are up to $\sim$10\% more satellites in the space between the pair, than expected from uniform. Consequently, it is a statistically very strong signal, being inconsistent with a uniform distribution at the 5$\sigma$ level. The possibility that the observed signal is the result of the overlap of two haloes with extended satellite distributions, is ruled out by testing this hypothesis by performing the same tests on isolated galaxies (and their satellites) artificially placed at similar separations. These findings highlight the unrelaxed and interacting nature of galaxies in pairs.

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.

To understand the role of the environment in galaxy formation, evolution, and present-day properties, it is essential to study the multi-frequency behavior of different galaxy populations under various environmental conditions. We crossmatch the SDSS DR10 group catalog with GAMA Data Release 2 and Wide-field Survey Explorer (WISE) data to construct a catalog of 1651 groups and 11436 galaxies containing photometric information in 15 different wavebands ranging from ultraviolet (0.152 {\mu}m) to mid-infrared (22 {\mu}m). We perform the spectral energy distribution (SED) fitting of galaxies using the MAGPHYS code and estimate the rest frame luminosities and stellar masses. We use the 1/Vmax method to estimate the galaxy stellar mass and luminosity functions, and the luminosity density field of galaxies to define the large scale environment of galaxies. The stellar mass functions of both central and satellite galaxies in groups are different in low and high density large scale environments. Satellite galaxies in high density environments have a steeper low mass end slope compared to low density environments, independently of the galaxy morphology. Central galaxies in low density environments have a steeper low mass end slope but the difference disappears for fixed galaxy morphology. The characteristic stellar mass of satellite galaxies is higher in high density environments and the difference exists only for galaxies with elliptical morphologies. Galaxy formation in groups is more efficient in high density large scale environments. Groups in high density environments have higher abundances of satellite galaxies, irrespective of the satellite galaxy morphology. The elliptical satellite galaxies are generally more massive in high density environments. The stellar masses of spiral satellite galaxies show no dependence on the large scale environment.

Superclusters are the largest relatively isolated systems in the cosmic web. Using the SDSS BOSS survey we search for the largest superclusters in the redshift range $0.43<z<0.71$. We generate a luminosity-density field smoothed over $8 h^{-1}\mathrm{Mpc}$ to detect the large-scale over-density regions. Each individual over-density region is defined as single supercluster in the survey. We define the superclusters in the way that they are comparable with the superclusters found in the SDSS main survey. We found a system we call the BOSS Great Wall (BGW), which consists of two walls with diameters 186 and 173 $h^{-1}$Mpc, and two other major superclusters with diameters of 64 and 91 $h^{-1}$Mpc. As a whole, this system consists of 830 galaxies with the mean redshift 0.47. We estimate the total mass to be approximately $2\times10^{17}h^{-1}M_\odot$. The morphology of the superclusters in the BGW system is similar to the morphology of the superclusters in the Sloan Great Wall region. The BGW is one of the most extended and massive system of superclusters yet found in the Universe.

(Abridged) We have analyzed Chandra LETG and XMM-Newton RGS spectra towards the z=0.177 quasar PG 1116+215, a sightline that is rendered particularly interesting by the HST detection of several OVI and HI broad Lyman-alpha absorption lines that may be associated with the warm-hot intergalactic medium. We performed a search for resonance K-alpha absorption lines from OVII and OVIII at the redshifts of the detected far-ultraviolet lines. We detected an absorption line in the Chandra spectra at 5.2 sigma confidence level at wavelengths corresponding to OVIII K-alpha at z=0.0911+-0.0004+-0.0005 (statistical followed by systematic error). This redshift is within 3 sigma of that of a HI broad Lyman-alpha of b=130 km/s at z=0.09279+-0.00005. We have also analyzed the available XMM-Newton RGS data towards PG 1116+215. Unfortunately, the XMM-Newton data are not suitable to investigate this line because of instrumental features at the wavelengths of interest. At the same redshift, the Chandra and XMM-Newton spectra have OVII K-alpha absorption line features of significance 1.5 sigma and 1.8 sigma, respectively. We also analyzed the available SDSS spectroscopic galaxy survey data towards PG 1116+215 in the redshift range of interest. We found evidence for a galaxy filament that intersects the PG 1116+215 sightline and additional galaxy structures that may host WHIM. The combination of HST, Chandra, XMM-Newton and SDSS data indicates that we have likely detected a multi-temperature WHIM at z=0.091-0.093 towards PG 1116+215.

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.

Tidal torque theory suggests that galaxies gain angular momentum in the linear stage of structure formation. Such a theory predicts alignments between the spin of haloes and tidal shear field. However, non-linear evolution and angular momentum acquisition may alter this prediction significantly. In this paper, we use a reconstruction of the cosmic shear field from observed peculiar velocities combined with spin axes extracted from galaxies within $115\, \mathrm{Mpc} $ ($\sim8000 \, {\mathrm {km}}{\mathrm s}^{-1}$) from 2MRS catalog, to test whether or not galaxies appear aligned with principal axes of shear field. Although linear reconstructions of the tidal field have looked at similar issues, this is the first such study to examine galaxy alignments with velocity-shear field. Ellipticals in the 2MRS sample, show a statistically significant alignment with two of the principal axes of the shear field. In general, elliptical galaxies have their short axis aligned with the axis of greatest compression and perpendicular to the axis of slowest compression. Spiral galaxies show no signal. Such an alignment is significantly strengthened when considering only those galaxies that are used in velocity field reconstruction. When examining such a subsample, a weak alignment with the axis of greatest compression emerges for spiral galaxies as well. This result indicates that although velocity field reconstructions still rely on fairly noisy and sparse data, the underlying alignment with shear field is strong enough to be visible even when small numbers of galaxies are considered - especially if those galaxies are used as constraints in the reconstruction.

This paper describes several applications in astronomy and cosmology that are addressed using probabilistic modelling and statistical inference.

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.

We propose a new approach to the missing baryons problem. Building on the common assumption that the missing baryons are in the form of the Warm Hot Intergalactic Medium (WHIM), we further assumed here that the galaxy luminosity density can be used as a tracer of the WHIM. The latter assumption is supported by our finding of a significant correlation between the WHIM density and the galaxy luminosity density in the hydrodynamical simulations of Cui et al. (2012). We further found that the fraction of the gas mass in the WHIM phase is substantially (by a factor of $\sim$1.6) higher within the large scale galactic filaments, i.e. $\sim$70\%, compared to the average in the full simulation volume of $\sim$0.1\,Gpc$^3$. The relation between the WHIM overdensity and the galaxy luminosity overdensity within the galactic filaments is consistent with linear: $\delta_{\rm whim}\,=\,0.7\,\pm\,0.1\,\times\,\delta_\mathrm{LD}^{0.9 \pm 0.2}$. We applied our procedure to the line of sight to the blazar H2356-309 and found evidence for the WHIM in correspondence of the Sculptor Wall (z $\sim$0.03 and $\log{N_H}$ = $19.9^{+0.1}_{-0.3}$) and Pisces-Cetus superclusters (z $\sim$0.06 and $\log{N_H}$ = $19.7^{+0.2}_{-0.3}$), in agreement with the redshifts and column densities of the X-ray absorbers identified and studied by Fang et al. (2010) and Zappacosta et al. (2010). This agreement indicates that the galaxy luminosity density and galactic filaments are reliable signposts for the WHIM and that our method is robust in estimating the WHIM density. The signal that we detected cannot originate from the halos of the nearby galaxies since they cannot account for the large WHIM column densities that our method and X-ray analysis consistently find in the Sculptor Wall and Pisces-Cetus superclusters.

Context. Mass loss is an important property in evolution models of massive stars. As up to 90% of the massive stars have a visual or spectroscopic companion and many of them exhibit mass exchange, mass-loss rates can be acquired through the period study of massive binaries. Aims. Using our own photometric observations as well as archival data, we look for variations in orbital periods of seven massive eclipsing binary systems in the Cygnus OB2 association and estimate their mass-loss rates and stellar parameters. Methods. We use a Bayesian parameter estimation method to simultaneously fit the period and period change to all available data and a stellar modelling tool to model the binary parameters from photometric and radial-velocity data. Results. Four out of the seven selected binaries show non-zero period change values at two-sigma confidence level. We also report for the first time the eclipsing nature of a star MT059.

We study the distribution, masses, and dynamical properties of galaxy groups in the A2142 supercluster. We analyse the global luminosity density distribution in the supercluster and divide the supercluster into the high-density core and the low-density outskirts regions. We find galaxy groups and filaments in the regions of different global density, calculate their masses and mass-to-light ratios and analyse their dynamical state with several 1D and 3D statistics. We use the spherical collapse model to study the dynamical state of the supercluster. We show that in A2142 supercluster groups and clusters with at least ten member galaxies lie along an almost straight line forming a 50 Mpc/h long main body of the supercluster. The A2142 supercluster has a very high density core surrounded by lower-density outskirt regions. The total estimated mass of the supercluster is M_est = 6.2 10^{15}M_sun. More than a half of groups with at least ten member galaxies in the supercluster lie in the high-density core of the supercluster, centered at the rich X-ray cluster A2142. Most of the galaxy groups in the core region are multimodal. In the outskirts of the supercluster, the number of groups is larger than in the core, and groups are poorer. The orientation of the cluster A2142 axis follows the orientations of its X-ray substructures and radio halo, and is aligned along the supercluster axis. The high-density core of the supercluster with the global density D8 > 17 and perhaps with D8 > 13 may have reached the turnaround radius and started to collapse. A2142 supercluster with luminous, collapsing core and straight body is an unusual object among galaxy superclusters. In the course of the future evolution the supercluster may be split into several separate systems.

The cosmic web that characterizes the large-scale structure of the Universe can be quantified by a variety of methods. For example, large redshift surveys can be used in combination with point process algorithms to extract long curvilinear filaments in the galaxy distribution. Alternatively, given a full 3D reconstruction of the velocity field, kinematic techniques can be used to decompose the web into voids, sheets, filaments and knots. In this paper we look at how two such algorithms - the Bisous model and the velocity shear web - compare with each other in the local Universe (within 100 Mpc), finding good agreement. This is both remarkable and comforting, given that the two methods are radically different in ideology and applied to completely independent and different data sets. Unsurprisingly, the methods are in better agreement when applied to unbiased and complete data sets, like cosmological simulations, than when applied to observational samples. We conclude that more observational data is needed to improve on these methods, but that both methods are most likely properly tracing the underlying distribution of matter in the Universe.

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. Gravitational collapse theory and numerical simulations suggest that the velocity field within large-scale galaxy filaments is dominated by motions along the filaments. Aims. Our aim is to check whether observational data reveal any preferred orientation of galaxy pairs with respect to the underlying filaments as a result of the expectedly anisotropic velocity field. Methods. We use galaxy pairs and galaxy filaments identified from the Sloan Digital Sky Survey data. For filament extraction, we use the Bisous model that is based the marked point process technique. During the filament detection, we use the centre point of each pair instead of the positions of galaxies to avoid a built-in influence of pair orientation on the filament construction. For pairs lying within filaments (3012 cases), we calculate the angle between the line connecting galaxies of each pair and their host filament. To avoid redshift-space distortions, the angle is measured in the plain of the sky. Results. The alignment analysis shows that the orientation of galaxy pairs correlates strongly with their host filaments. The alignment signal is stronger for loose pairs, with at least 25% excess of aligned pairs compared to a random distribution. The alignment of galaxy pairs and filaments measured from the observational data is in good concordance with the alignment in the Millennium simulation and thus provides support to the {\Lambda}CDM formalism.

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

We investigate whether the satellite luminosity function (LF) of primary galaxies identified in the Sloan Digital Sky Survey (SDSS) depends on whether the host galaxy is in a filament or not. Isolated primary galaxies are identified in the SDSS spectroscopic sample while potential satellites (that are up to 4 magnitudes fainter than their hosts) are searched for in the much deeper photometric sample. Filaments are constructed from the galaxy distribution by the "Bisous" process. Isolated primary galaxies are divided into two subsamples: those in filaments and those not in filaments. We examine the stacked mean satellite LF of both the filament and non-filament sample and find that, on average, the satellite LFs of galaxies in filaments is significantly higher than those of galaxies not in filaments. The filamentary environment can increases the abundance of the brightest satellites ($M_\mathrm{sat.} < M_\mathrm{prim.} + 2.0$), by a factor of $\sim 2$ compared with non-filament isolated galaxies. This result is independent of primary galaxy magnitude although the satellite LF of galaxies in the faintest magnitude bin, is too noisy to determine if such a dependence exists. Since our filaments are extracted from a spectroscopic flux-limited sample, we consider the possibility that the difference in satellite LF is due to a redshift, colour or environmental bias, finding these to be insufficient to explain our result. The dependence of the satellite LF on the cosmic web suggests that the filamentary environment may have a strong effect on the efficiency of galaxy formation.

We have developed a new method to approach the missing baryons problem. We assume that the missing baryons reside in a form of Warm Hot Intergalactic Medium, i.e. the WHIM. Our method consists of (a) detecting the coherent large scale structure in the spatial distribution of galaxies that traces the Cosmic Web and that in hydrodynamical simulations is associated to the WHIM, (b) map its luminosity into a galaxy luminosity density field, (c) use numerical simulations to relate the luminosity density to the density of the WHIM, (d) apply this relation to real data to trace the WHIM using the observed galaxy luminosities in the Sloan Digital Sky Survey and 2dF redshift surveys. In our application we find evidence for the WHIM along the line of sight to the Sculptor Wall, at redshifts consistent with the recently reported X-ray absorption line detections. Our indirect WHIM detection technique complements the standard method based on the detection of characteristic X-ray absorption lines, showing that the galaxy luminosity density is a reliable signpost for the WHIM. For this reason, our method could be applied to current galaxy surveys to optimise the observational strategies for detecting and studying the WHIM and its properties. Our estimates of the WHIM hydrogen column density in Sculptor agree with those obtained via the X-ray analysis. Due to the additional column density estimate, our method has potential for improving the constrains of the physical parameters of the WHIM as derived with X-ray absorption, and thus for improving the understanding of the missing baryons problem.

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.

The imprints of dust-starlight interactions are visible in scaling relations between stellar characteristics, star formation parameters and dust properties. We aim to examine dust scaling relations on a sub-kpc resolution in the Andromeda galaxy (M31) by comparing the properties on a local and global scale to other galaxies of the local universe. New Herschel observations are combined with available data from GALEX, SDSS, WISE and Spitzer to construct a dataset covering UV to submm wavelengths. We work at the resolution of the SPIRE $500\; \mu$m beam, with pixels corresponding to physical regions of 137 x 608 pc in the galaxy's disk. A panchromatic spectral energy distribution was modelled for each pixel and several dust scaling relations are investigated. We find, on a sub-kpc scale, strong correlations between $M_d/M_\star$ and NUV-r, and between $M_d/M_\star$ and $\mu_\star$ (the stellar mass surface density). Striking similarities with corresponding relations based on integrated galaxies are found. We decompose M31 in four macro-regions based on their FIR morphology; the bulge, inner disk, star forming ring and the outer disk. All regions closely follow the galaxy-scale average trends. The specific star formation characteristics we derive for these macro-regions give strong hints of an inside-out formation of the bulge-disk geometry, as well as an internal downsizing process. However, within each macro-region, a great diversity in individual micro-regions is found. Furthermore, we confirm that dust in the bulge of M31 is heated only by the old stellar populations. In general, the local dust scaling relations indicate that the dust content in M31 is maintained by a subtle interplay of past and present star formation. The similarity with galaxy-based relations strongly suggests that they are in situ correlations, with underlying processes that must be local in nature. (Abriged)

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).

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.

We compare the galaxy populations in superclusters of different morphology in the nearby Universe (180 < d < 270 Mpc) to see whether the inner structure and overall morphology of superclusters are important in shaping galaxy properties in superclusters. Supercluster morphology has been found with Minkowski functionals. We analyse the probability density distributions of colours, morphological types, stellar masses, star formation rates (SFR) of galaxies, and the peculiar velocities of the main galaxies in groups in superclusters of filament and spider types, and in the field. We show that the fraction of red, early-type, low SFR galaxies in filament-type superclusters is higher than in spider-type superclusters; in low-density global environments their fraction is lower than in superclusters. In all environments the fraction of red, high stellar mass, and low SFR galaxies in rich groups is higher than in poor groups. In superclusters of spider morphology red, high SFR galaxies have higher stellar masses than in filament-type superclusters. Groups of equal richness host galaxies with larger stellar masses, a larger fraction of early-type and red galaxies, and a higher fraction of low SFR galaxies, if they are located in superclusters of filament morphology. The peculiar velocities of the main galaxies in groups from superclusters of filament morphology are higher than in those of spider morphology. Groups with higher peculiar velocities of their main galaxies in filament-type superclusters are located in higher density environment than those with low peculiar velocities. There are significant differences between galaxy populations of the individual richest superclusters. Therefore both local (group) and global (supercluster) environments and even supercluster morphology play an important role in the formation and evolution of galaxies.

The main feature of the spatial large-scale galaxy distribution is its intricate network of galaxy filaments. This network is spanned by the galaxy locations that can be interpreted as a three-dimensional point distribution. The global properties of the point process can be measured by different statistical methods, which, however, do not describe directly the structure elements. The morphology of the large scale structure, on the other hand, is an important property of the galaxy distribution. Here we apply an object point process with interactions (the Bisous model) to trace and extract the filamentary network in the presently largest galaxy redshift survey, the Sloan Digital Sky Survey (SDSS). We search for filaments in the galaxy distribution that have a radius of about 0.5 Mpc/h. We divide the detected network into single filaments and present a public catalogue of filaments. We study the filament length distribution and show that the longest filaments reach the length of 60 Mpc/h. The filaments contain 35-40% of the total galaxy luminosity and they cover roughly 5-8% of the total volume, in good agreement with N-body simulations and previous observational results.

The Millennium N-body simulation and the Sloan Digital Sky Survey seventh data release (SDSS DR7) galaxy and galaxy group catalogues are compared to study the properties of galaxy groups and the distribution of galaxies in groups. We construct mock galaxy group catalogues for a Millennium semi-analytical galaxy catalogue by using the same friends-of-friends method, which was used by Tago et al to analyse the SDSS data. We analyse in detail the group luminosities, group richnesses, virial radii, sizes of groups and their rms velocities for four volume-limited samples from observations and simulations. Our results show that the spatial densities of groups agree within one order of magnitude in all samples with a rather good agreement between the mock catalogues and observations. All group property distributions have similar shapes and amplitudes for richer groups. For galaxy pairs and small groups, the group properties for observations and simulations are clearly different. In addition, the spatial distribution of galaxies in small groups is different: at the outskirts of the groups the galaxy number distributions do not agree, although the agreement is relatively good in the inner regions. Differences in the distributions are mainly due to the observational limitations in the SDSS sample and to the problems in the semi-analytical methods that produce too compact and luminous groups.

The large-scale structure of the Universe is characterised by a web-like structure made of voids, sheets, filaments, and knots. The structure of this so-called cosmic web is dictated by the local velocity shear tensor. In particular, the local direction of a filament should be strongly aligned with e3, the eigenvector associated with the smallest eigenvalue of the tensor. That conjecture is tested here on the basis of a cosmological simulation. The cosmic web delineated by the halo distribution is probed by a marked point process with interactions (the Bisous model), detecting filaments directly from the halo distribution (P-web). The detected P-web filaments are found to be strongly aligned with the local e3: the alignment is within 30 degree for 80% of the elements. This indicates that large-scale filaments defined purely from the distribution of haloes carry more than just morphological information, although the Bisous model does not make any prior assumption on the underlying shear tensor. The P-web filaments are also compared to the structure revealed from the velocity shear tensor itself (V-web). In the densest regions, the P- and V-web filaments overlap well (90%), whereas in lower density regions, the P-web filaments preferentially mark sheets in the V-web.

Galaxies properties are known to be affected by their environment. One important question is how their angular momentum reflects the surrounding cosmic web. We use the SDSS to investigate the spin axes of spiral and elliptical galaxies relative to their surrounding filament/sheet orientations. To detect filaments a marked point process with interactions (the "Bisous model") is used. Sheets are found by detecting "flattened" filaments. The minor axes of ellipticals are found to be preferentially perpendicular to hosting filaments. A weak correlation is found with sheets. These findings are consistent with the notion that elliptical galaxies formed via mergers which predominantly occurred along the filaments. The spin axis of spiral galaxies is found to align with the host filament, with no correlation between spiral spin and sheet normal. When examined as a function of distance from the filament axis, a much stronger correlation is found in outer parts, suggesting that the alignment is driven by the laminar infall of gas from sheets to filaments. When compared with numerical simulations our results suggest that the connection between DM halo and galaxy spin is not straightforward. Our results provide an important input to the understanding of how galaxies acquire their angular momentum.

We show that the cosmic-ray positron excess measured by PAMELA and AMS could be induced by Dark Matter annihilations in a local over-density. In such a context leptophilic DM is not needed and good fits to positron data, in agreement with antiproton and gamma-ray measurements, are obtained for DM annihilations to $WW$, $hh$, $ZZ$, $ t \bar t$, $b \bar b$, $q\bar q$ channels. The classic Dark Matter candidates, such as the pure supersymmetric Wino with standard thermal annihilation cross-section, can fit the positron excess, without invoking any additional assumption on Dark Matter properties.

We search for spectral features in Fermi-LAT gamma-rays coming from regions corresponding to eighteen brightest nearby galaxy clusters determined by the magnitude of their signal line-of-site integrals. We observe a double peak-like excess over the diffuse power-law background at photon energies 110 GeV and 130 GeV with the global statistical significance up to 3.6\sigma, confirming independently earlier claims of the same excess from Galactic centre. Interpreting this result as a signal of dark matter annihilations to two monochromatic photon channels in galaxy cluster haloes, and fixing the annihilation cross section from the Galactic centre data, we determine the annihilation boost factor due to dark matter subhaloes from data. Our results contribute to discrimination of the dark matter annihilations from astrophysical processes and from systematic detector effects as the possible explanations to the Fermi-LAT excess.

We analyze publicly available Fermi-LAT high-energy gamma-ray data and confirm the existence of clear spectral feature peaked at E=130GeV. Scanning over the Galaxy we identify several disconnected regions where the observed excess originates from. Our best optimized fit is obtained for the central region of Galaxy with a clear peak at 130GeV with local statistical significance 4.5 sigma. The observed excess is not correlated with Fermi bubbles. We compute the photon spectra induced by dark matter annihilations into two and four standard model particles, the latter via two light intermediate states, and fit the spectra with data. Since our fits indicate sharper and higher signal peak than in the previous works, data favors dark matter direct two-body annihilation channels into photons or other channels giving only line-like spectra. If Einasto halo profile correctly predicts the central cusp of Galaxy, dark matter annihilation cross-section to two photons is of order ten percent of the standard thermal freeze-out cross-section. The large dark matter two-body annihilation cross-section to photons may signal a new resonance that should be searched for at the CERN LHC experiments.

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.

Majority of all galaxies reside in groups of less than 50 member galaxies. These groups are distributed in various large-scale environments from voids to superclusters. Evolution of galaxies is affected by the environment in which they reside. Our aim is to study the effects that the local group scale and the supercluster scale environment have on galaxies. We use a luminosity-density field to determine density of the large-scale environment of galaxies in groups of various richness. We calculate fractions of different types of galaxies in groups with richnesses up to 50 member galaxies and in different large-scale environments from voids to superclusters. The fraction of passive elliptical galaxies rises and the fraction of star-forming spiral galaxies declines when the richness of a group of galaxies rises from two to approximately ten galaxies. On the large scale, the passive elliptical galaxies become more numerous than star-forming spirals when the environmental density grows to the density level typical for superclusters. The large-scale environment affects the level of these fractions in groups: galaxies in equally rich groups are more likely to be elliptical in supercluster environments than in lower densities. The crossing point, where the number of passive and star-forming galaxies is equal, happens in groups with lower richness in superclusters than in voids. Galaxies in low-density areas require richer groups to evolve from star-forming to passive. Groups in superclusters are on average more luminous than groups in large-scale environments with lower density. Our results suggest that the evolution of galaxies is affected by both, by the group in which the galaxy resides, and by its large-scale environment. Galaxies in lower-density regions develop later than galaxies in similar mass groups in high-density environments.

Galaxies are not distributed randomly in the cosmic web but are instead arranged in filaments and sheets surrounding cosmic voids. Observationally there is still no convincing evidence of a link between the properties of galaxies and their host structures. However, by the tidal torque theory (our understanding of the origin of galaxy angular momentum), such a link should exist. Using the presently largest spectroscopic galaxy redshift survey (SDSS) we study the connection between the spin axes of galaxies and the orientation of their host filaments. We use a three dimensional field of orientations to describe cosmic filaments. To restore the inclination angles of galaxies, we use a 3D photometric model of galaxies that gives these angles more accurately than traditional 2D models. We found evidence that the spin axes of bright spiral galaxies have a weak tendency to be aligned parallel to filaments. For elliptical/S0 galaxies, we have a statistically significant result that their spin axes are aligned preferentially perpendicular to the host filaments; we show that this signal practically does not depend on the accuracy of the estimated inclination angles for elliptical/S0 galaxies.

We look for possible spectral features and systematic effects in Fermi-LAT publicly available high-energy gamma-ray data by studying photons from the Galactic centre, nearby galaxy clusters, nearby brightest galaxies, AGNs, unassociated sources, hydrogen clouds and Earth Limb. Apart from already known 130 GeV gamma-ray excesses from the first two sources, we find no new statistically significant signal from others. Much of our effort goes to studying Earth Limb photons. In the energy range 30 GeV to 200 GeV the Earth Limb gamma-ray spectrum follows power-law with spectral index 2.87\pm 0.04 at 95 % CL, in a good agreement with the PAMELA measurement of cosmic ray proton spectral index between 2.82-2.85, confirming the physical origin of the Limb gamma-rays. In small subsets of Earth Limb data with small photon incidence angle it is possible to obtain spectral features at different energies, including at 130 GeV, but determination of background, thus their significances, has large uncertainties in those cases. We observe systematic 2\sigma level differences in the Earth Limb spectra of gamma-rays with small and large incidence angles. The behaviour of those spectral features as well as background indicates that they are likely statistical fluctuations.

Stellar mass distribution in M31 is estimated using optical and infrared imaging data. Combining the derived stellar mass model with various kinematical data, properties of the DM halo of the galaxy are constrained. SDSS observations through the ugriz filters and the Spitzer imaging at 3.6 microns are used to sample the SED of the galaxy at each imaging pixel. Intrinsic dust extinction effects are taken into account by using far-infrared observations. Synthetic SEDs created with different stellar population synthesis models are fitted to the observed SEDs, providing estimates for the stellar mass surface density. The stellar mass distribution of the galaxy is described with a 3D model consisting of a nucleus, a bulge, a disc, a young disc and a halo component, each following the Einasto density distribution (relations between different functional forms of the Einasto density distribution are given in App. B). By comparing the stellar mass distribution to the observed rotation curve and kinematics of outer globular clusters and satellite galaxies, the DM halo parameters are estimated. Stellar population synthesis models suggest that M31 is dominated by old stars throughout the galaxy. The total stellar mass is (10-15)10^10Msun, 30% of which is in the bulge and 56% in the disc. None of the tested DM distribution models can be falsified on the basis of the stellar matter distribution and the rotation curve of the galaxy. The virial mass of the DM halo is (0.8-1.1)10^12Msun and the virial radius is 189-213kpc, depending on the DM distribution. The central density of the DM halo is comparable to that of nearby dwarf galaxies, low-surface-brightness galaxies and distant massive disc galaxies, thus the evolution of central DM halo properties seems to be regulated by similar processes for a broad range of halo masses, environments, and cosmological epochs.

We search for the presence of double gamma-ray line from unassociated Fermi-LAT sources including detailed Monte Carlo simulations to study its global statistical significance. Applying the Su & Finkbeiner selection criteria for high-energy photons we obtain a similar excess over the power-law background from 12 unassociated sources. However, the Fermi-LAT energy resolution and the present low statistics does not allow to distinguish a double peak from a single one with any meaningful statistical significance. We study the statistical significance of the fit to data with Monte Carlo simulations and show that the fit agrees almost perfectly with the expectations from random scan over the sky. We conclude that the claimed high-energy gamma-ray excess over the power-law background from unassociated sources is nothing but an artifact of the applied selection criteria and no preference to any excess can be claimed with the present statistics.

We study the relations between the multimodality of galaxy clusters drawn from the SDSS DR8 and the environment where they reside. As cluster environment we consider the global luminosity density field, supercluster membership, and supercluster morphology. We use 3D normal mixture modelling, the Dressler-Shectman test, and the peculiar velocity of cluster main galaxies as signatures of multimodality of clusters. We calculate the luminosity density field to study the environmental densities around clusters, and to find superclusters where clusters reside. We determine the morphology of superclusters with the Minkowski functionals and compare the properties of clusters in superclusters of different morphology. We apply principal component analysis to study the relations between the multimodality parametres of clusters and their environment simultaneously. We find that multimodal clusters reside in higher density environment than unimodal clusters. Clusters in superclusters have higher probability to have substructure than isolated clusters. The superclusters can be divided into two main morphological types, spiders and filaments. Clusters in superclusters of spider morphology have higher probabilities to have substructure and larger peculiar velocities of their main galaxies than clusters in superclusters of filament morphology. The most luminous clusters are located in the high-density cores of rich superclusters. Five of seven most luminous clusters, and five of seven most multimodal clusters reside in spider-type superclusters; four of seven most unimodal clusters reside in filament-type superclusters. Our study shows the importance of the role of superclusters as high density environment which affects the properties of galaxy systems in them.

Aims. We intend to compile a new galaxy group and cluster sample of the latest available SDSS data, adding several parameter for the purpose of studying the supercluster network, galaxy and group evolution, and their connection to the surrounding environment. Methods. We used a modified friends-of-friends (FoF) method with a variable linking length in the transverse and radial directions to eliminate selection effects and to find reliably as many groups as possible. Using the galaxies as a basis, we calculated the luminosity density field. Results. We create a new catalogue of groups and clusters for the SDSS data release 8 sample. We find and add environmental parameters to our catalogue, together with other galaxy parameters (e.g., morphology), missing from our previous catalogues. We take into account various selection effects caused by a magnitude limited galaxy sample. Our final sample contains 576493 galaxies and 77858 groups. The group catalogue is available at http://www.aai.ee/~elmo/dr8groups/ and from the Strasbourg Astronomical Data Center (CDS).

We search for the presence of substructure, a non-Gaussian, asymmetrical velocity distribution of galaxies, and large peculiar velocities of the main galaxies in galaxy clusters with at least 50 member galaxies, drawn from the SDSS DR8. We employ a number of 3D, 2D, and 1D tests to analyse the distribution of galaxies in clusters: 3D normal mixture modelling, the Dressler-Shectman test, the Anderson-Darling and Shapiro-Wilk tests and others. We find the peculiar velocities of the main galaxies, and use principal component analysis to characterise our results. More than 80% of the clusters in our sample have substructure according to 3D normal mixture modelling, the Dressler-Shectman (DS) test shows substructure in about 70% of the clusters. The median value of the peculiar velocities of the main galaxies in clusters is 206 km/s (41% of the rms velocity). The velocities of galaxies in more than 20% of the clusters show significant non-Gaussianity. While multidimensional normal mixture modelling is more sensitive than the DS test in resolving substructure in the sky distribution of cluster galaxies, the DS test determines better substructure expressed as tails in the velocity distribution of galaxies. Richer, larger, and more luminous clusters have larger amount of substructure and larger (compared to the rms velocity) peculiar velocities of the main galaxies. Principal component analysis of both the substructure indicators and the physical parameters of clusters shows that galaxy clusters are complicated objects, the properties of which cannot be explained with a small number of parameters or delimited by one single test. The presence of substructure, the non-Gaussian velocity distributions, as well as the large peculiar velocities of the main galaxies, shows that most of the clusters in our sample are dynamically young.

We have constructed a set of supercluster catalogues for the galaxies from the SDSS survey main and luminous red galaxy (LRG) flux-limited samples. To delineate superclusters, we calculated luminosity density fields using the B3-spline kernel of the radius of 8 Mpc/h for the main sample and 16 Mpc/h for the LRG sample and define regions with densities over a selected threshold as superclusters, while utilising almost the whole volume of both samples. We created two types of catalogues, one with an adaptive local threshold and a set of catalogues with different global thresholds. We describe the supercluster catalogues and their general properties. Using smoothed bootstrap, we find uncertainty estimates for the density field and use these to attribute confidence levels to the catalogue objects. We have also created a test catalogue for the galaxies from the Millennium simulation to compare the simulated and observed superclusters and to clarify the methods we use. We find that the superclusters are well-defined systems, and the properties of the superclusters of the main and LRG samples are similar. We also show that with adaptive local thresholds we get a sample of superclusters, the properties of which do not depend on their distance from the observer. The Millennium galaxy catalogue superclusters are similar to those observed.

According to the modern cosmological paradigm galaxies and galaxy systems form from tiny density perturbations generated during the very early phase of the evolution of the Universe. Using numerical simulations we study the evolution of phases of density perturbations of different scales to understand the formation and evolution of the cosmic web. We apply the wavelet analysis to follow the evolution of high-density regions (clusters and superclusters) of the cosmic web. We show that the positions of maxima and minima of density waves (their spatial phases) almost do not change during the evolution of the structure. Positions of extrema of density perturbations are the more stable, the larger is the wavelength of perturbations. Combining observational and simulation data we conclude that the skeleton of the cosmic web was present already in an early stage of structure evolution.

Our goal is to see how density waves of different scale combine to form voids between galaxy systems of various scale. We perform numerical simulations of structure formation in cubes of size 100 and 256 Mpc/h, with resolutions 256^3 and 512^3 particles and cells. To understand the role of density perturbations of various scale we cut power spectra at scales from 8 to 128 Mpc/h, using in all series identical initial random realisations. We find that small haloes and short filaments form all over the simulation box, if perturbations only up to scale 8 Mpc/h are present. The phenomenon of large multi-scale voids in the cosmic web requires the presence of an extended spectrum of primordial density perturbations. The void phenomenon is due to the action of two processes: the synchronisation of density perturbations of medium and large scales, and the suppression of galaxy formation in low-density regions by the combined action of negative sections of medium- and large-scale density perturbations, so that their densities are less than the mean density, and thus during the evolution their densities decrease.

We follow the evolution of galaxy systems in numerical simulation. Our goal is to understand the role of density perturbations of various scales in the formation and evolution of the cosmic web. We perform numerical simulations with the full power spectrum of perturbations, and with spectrum cut at long wavelengths. Additionally, we have one model, where we cut the intermediate waves. We analyze the density field and study the void sizes and density field clusters in different models. Our analysis shows that the fine structure (groups and clusters of galaxies) is created by small-scale density perturbations of scale $\leq 8$ \Mpc. Filaments of galaxies and clusters are created by perturbations of intermediate scale from $\sim 8$ to $\sim 32$ \Mpc, superclusters of galaxies by larger perturbations. We conclude that the scale of the pattern of the cosmic web is determined by density perturbations of scale up to $\sim 100$ \Mpc. Larger perturbations do not change the pattern of the web, but modulate the richness of galaxy systems, and make voids emptier. The stop of the increase of the scale of the pattern of the cosmic web with increasing scale of density perturbations can probably be explained as the freezing of the web at redshift $z\simeq 0.7$.

We apply the principal component analysis and Spearman's correlation test to study the properties of superclusters drawn from the SDSS DR7. We analyse possible selection effects in the supercluster catalogue, study the physical and morphological properties of superclusters, find their possible subsets, and determine scaling relations for superclusters. We show that the parameters of superclusters do not correlate with their distance. The correlations between the physical and morphological properties of superclusters are strong. Superclusters can be divided into two populations according to their total luminosity. High-luminosity superclusters form two sets, more elongated systems with the shape parameter K_1/K_2 < 0.5 and less elongated ones with K_1/K_2 > 0.5. The first two principal components account for more than 90% of the variance in the supercluster parameters and define the fundamental plane, which characterises the physical and morphological properties of superclusters. We use principal component analysis to derive scaling relations for superclusters, in which we combine the physical and morphological parameters of superclusters. Structure formation simulations for different cosmologies, and more data about the local and high redshift superclusters are needed to understand better the evolution and the properties of superclusters.

We study the morphology of a set of superclusters drawn from the SDSS DR7. We calculate the luminosity density field to determine superclusters from a flux- limited sample of galaxies from SDSS DR7, and select superclusters with 300 and more galaxies for our study. The morphology of superclusters is described with the fourth Minkowski functional V3, the morphological signature (the curve in the shapefinder's K1-K2 plane) and the shape parameter (the ratio of the shapefinders K1/K2). We investigate the supercluster sample using multidimensional normal mixture modelling, and use Abell clusters to identify our superclusters with known superclusters and to study the large-scale distribution of superclusters. The superclusters in our sample form three chains of superclusters; one of them is the Sloan Great Wall. Most superclusters have filament-like overall shapes. Superclusters can be divided into two sets; more elongated superclusters are more luminous, richer, have larger diameters, and a more complex fine structure than less elongated superclusters. The fine structure of superclusters can be divided into four main morphological types: spiders, multispiders, filaments, and multibranching filaments. We present the 2D and 3D distribution of galaxies and rich groups, the fourth Minkowski functional, and the morphological signature for all superclusters. Widely different morphologies of superclusters show that their evolution has been dissimilar. A study of a larger sample of superclusters from observations and simulations is needed to understand the morphological variety of superclusters and the possible connection between the morphology of superclusters and their large-scale environment.

We present the results of the study of the morphology and galaxy content of the Sloan Great Wall (SGW). We use the luminosity density field to determine superclusters in the SGW, and the fourth Minkowski functional V_3 and the morphological signature (the K_1-K_2 shapefinders curve) to show the different morphologies of the SGW, from a single filament to a multibranching, clumpy planar system. The richest supercluster in the SGW, SCl~126 and especially its core resemble a very rich filament, while another rich supercluster in the SGW, SCl~111, resembles a "multispider" - an assembly of high density regions connected by chains of galaxies. Using Minkowski functionals we study the substructure of individual galaxy populations determined by their color in these superclusters. We assess the statistical significance of the results with the halo model and smoothed bootstrap. We study the galaxy content and the properties of groups of galaxies in two richest superclusters of the SGW, paying special attention to bright red galaxies (BRGs) and to the first ranked galaxies in SGW groups. About 1/3 of BRGs are spirals. The scatter of colors of elliptical BRGs is smaller than that of spiral BRGs. About half of BRGs and of first ranked galaxies in groups have large peculiar velocities. Groups with elliptical BRGs as their first ranked galaxies populate superclusters more uniformly than the groups, which have a spiral BRG as its first ranked galaxy. The galaxy and group content of the core of the supercluster SCl~126 shows several differences in comparison with the outskirts of this supercluster and with the supercluster SCl~111. Our results suggest that the formation history and evolution of individual neighbour superclusters in the SGW has been different.

We study the influence of the environment on the evolution of galaxies by investigating the luminosity function (LF) of galaxies of different morphological types and colours at different environmental density levels. We construct the LFs separately for galaxies of different morphology (spiral and elliptical) and of different colours (red and blue) using data from the Sloan Digital Sky Survey (SDSS), correcting the luminosities for the intrinsic absorption. We use the global luminosity density field to define different environments, and analyse the environmental dependence of galaxy morphology and colour. The smoothed bootstrap method is used to calculate confidence regions of the derived luminosity functions. We find a strong environmental dependency for the LF of elliptical galaxies. The LF of spiral galaxies is almost environment independent, suggesting that spiral galaxy formation mechanisms are similar in different environments. Absorption by the intrinsic dust influences the bright-end of the LF of spiral galaxies. After attenuation correction, the brightest spiral galaxies are still about 0.5 mag less luminous than the brightest elliptical galaxies, except in the least dense environment, where spiral galaxies dominate the LF at every luminosity. Despite the extent of the SDSS survey, the influence of single rich superclusters is present in the galactic LF of the densest environment.

The objective of this work is to obtain an extinction-corrected distribution of optical surface brightness and colour indices of the large nearby galaxy M 31 using homogeneous observational data and a model for intrinsic extinction. We process the Sloan Digital Sky Survey (SDSS) images in ugriz passbands and construct corresponding mosaic images, taking special care of subtracting the varying sky background. We apply the galactic model developed in Tempel et al. (2010) and far-infrared imaging to correct the photometry for intrinsic dust effects. We obtain observed and dust-corrected distributions of the surface brightness of M 31 and a map of line-of-sight extinctions inside the galaxy. Our extinction model suggests that either M 31 is intrinsically non-symmetric along the minor axis or the dust properties differ from those of the Milky Way. Assuming the latter case, we present the surface brightness distributions and integral photometry for the Sloan filters as well as the standard UBVRI system. We find the following intrinsic integral colour indices for M 31: (U-B)_0=0.35; (B-V)_0=0.86; (V-R)_0=0.63; (R-I)_0=0.53; the total intrinsic absorption-corrected luminosities of M 31 in the B and the V filters are 4.10 and 3.24 mag, respectively.

We present the results of the study of the substructure and galaxy content of ten rich clusters of galaxies in three different superclusters of the Sloan Great Wall. We determine the substructure in clusters using the 'Mclust' package from the 'R' statistical environment and analyse their galaxy content. We analyse the distribution of the peculiar velocities of galaxies in clusters and calculate the peculiar velocity of the first ranked galaxy. We show that clusters in our sample have more than one component; in some clusters different components also have different galaxy content. We find that in some clusters with substructure the peculiar velocities of the first ranked galaxies are large. All clusters in our sample host luminous red galaxies. They can be found both in the central areas of clusters as well as in the outskirts, some of them have large peculiar velocities. About 1/3 of red galaxies in clusters are spirals. The scatter of colours of red ellipticals is in most clusters larger than that of red spirals. The presence of substructure in rich clusters, signs of possible mergers and infall, as well as the large peculiar velocities of the first ranked galaxies suggest that the clusters in our sample are not yet virialized. We present merger trees of dark matter haloes in an N-body simulation to demonstrate the formation of present-day dark matter haloes via multiple mergers during their evolution. In simulated dark matter haloes we find a substructure similar to that in observed clusters.

The correlation function of the distribution of matter in the universe shows, at large scales, baryon acoustic oscillations, which were imprinted prior to recombination. This feature was first detected in the correlation function of the luminous red galaxies (LRG) of the Sloan Digital Sky Survey (SDSS). The final release (DR7) of the SDSS has been recently made available, and the useful volume is about two times bigger than in the old sample. We present here, for the first time, the redshift space correlation function of this sample at large scales together with that for one shallower, but denser volume-limited subsample drawn from the 2dF redshift survey. We test the reliability of the detection of the acoustic peak at about 100 Mpc/h and the behaviour of the correlation function at larger scales by means of careful estimation of errors. We confirm the presence of the peak in the latest data although broader than in previous detections.

We create a model for recovering the intrinsic, absorption-corrected surface brightness distribution of a galaxy and apply the model to the M31. We construct a galactic model as a superposition of axially symmetric stellar components and a dust disc to analyse the intrinsic absorption efects. Dust column density is assumed to be proportional to the far-infrared flux of the galaxy. Along each line of sight, the observed far-infrared spectral energy distribution is approximated with modified black body functions considering dust components with different temperatures, allowing to determine the temperatures and relative column densities of the dust components. We apply the model to the nearby galaxy M31 using the Spitzer Space Telescope far-infrared observations for mapping dust distribution and temperature. A warm and a cold dust component are distinguished. The temperature of the warm dust in M31 varies between 56 and 60 K and is highest in the spiral arms; the temperature of the cold component is mostly 15-19 K and rises up to about 25 K at the centre of the galaxy. The intensity-weighted mean temperature of the dust decreases from T ~32 K at the centre to T ~20 K at R ~7 kpc and outwards. We also calculate the intrinsic UBVRIL surface brightness distributions and the spatial luminosity distribution. The intrinsic dust extinction in the V-colour rises from 0.25 mag at the centre to 0.4-0.5 mag at R = 6-13 kpc and decreases smoothly thereafter. The calculated total extinction-corrected luminosity of M31 is L_B = (3.64 pm 0.15) 10^10L_sun, corresponding to an absolute luminosity M_B = (-20.89 pm 0.04) mag. Of the total B-luminosity, 20% (0.24 mag) is obscured from us by the dust inside M31. The intrinsic shape of the bulge is slightly prolate in our best-fit model.

For the first time spectroscopic galaxy redshift surveys are reaching the scales where galaxies can be studied together with the nearest quasars. This gives an opportunity to study the dependence between the activity of a quasar and its environment in a more extensive way than before. We study the spatial distribution of galaxies and groups of galaxies in the environments of low redshift quasars in the Sloan Digital Sky Survey (SDSS). Our aim is to understand how the nearby quasars are embedded in the local and global density field of galaxies and how the environment affects quasar activity. We analyse the environments of nearby quasars using number counts of galaxies. We also study the dependence of group properties to their distance to the nearest quasar. The large scale environments are studied by analysing the locations of quasars in the luminosity density field. Our study of the number counts of galaxies in quasar environments shows an underdensity of bright galaxies at a few Mpc from quasars. Also, the groups of galaxies that have a quasar closer than 2Mpc are poorer and less luminous than in average. Our analysis on the luminosity density field shows that quasars clearly avoid rich superclusters. Nearby quasars seem to be located in outskirts of superclusters or in filaments connecting them. Our results suggest that quasar evolution may be affected by density variations both on supercluster scales and in the local environment.

Aims. We use the 2dF Galaxy Redshift Survey to derive the luminosity function (LF) of the first-ranked (brightest) group/cluster galaxies, the LF of second-ranked, satellite and isolated galaxies, and the LF of groups of galaxies. Methods. We investigate the LFs of different samples in various environments: in voids, filaments, superclusters and supercluster cores. We compare the derived LFs with the Schechter and double-power-law analytical expressions. We also analyze the luminosities of isolated galaxies. Results. We find strong environmental dependency of luminosity functions of all populations. The luminosities of first-ranked galaxies have a lower limit, depending on the global environment (higher in supercluster cores, and absent in voids). The LF of second-ranked galaxies in high-density regions is similar to the LF of first-ranked galaxies in a lower-density environment. The brightest isolated galaxies can be identified with first-ranked galaxies at distances where the remaining galaxies lie outside the observational window used in the survey. Conclusions. The galaxy and cluster LFs can be well approximated by a double-power-law; the widely used Schechter function does not describe well the bright end and the bend of the LFs. Properties of the LFs reflect differences in the evolution of galaxies and their groups in different environments.

We construct a structural model of the Andromeda Galaxy, simultaneously corresponding to observed photometrical and kinematical data and chemical abundances. In this paper we present the observed surface brightness, colour and metallicity distributions, and compare them to the model galaxy. In Paper II (Tempel, Tamm & Tenjes 2007) we present similar data for the kinematics, and derive the mass distribution of the galaxy. On the basis of U, B, V, R, I and L luminosity distributions, we construct the model galaxy as a superposition of four axially symmetric stellar components: a bulge, a disc, an inner halo and an extended diffuse halo. By using far-infrared imaging data of M31 and a thin dust disc assumption, we derive dust-free surface brightness and colour distributions. We find the total absorption corrected luminosity of M31 to be L_B = (3.3+/-0.7)x10^10 L_sun, corresponding to an absolute luminosity M_B = -20.8+/-0.2 mag. Of the total luminosity, 41% (0.57 mag) is obscured from us by the dust inside M31. Using chemical evolution models, we calculate mass-to-light ratios of the components, correspoding to the colour indices and metallicities. We find the total intrinsic mass-to-light ratio of the visible matter to be M/L_B=3.1-5.8 M_sun/L_sun and the total mass of visible matter M_vis =(10-19)x10^10 M_sun. The use of the model parameters for a dynamical analysis and for determining dark matter distribution is presented in Paper II.

In the present paper we derive the density distribution of dark matter (DM) in a well-observed nearby disc galaxy, the Andromeda galaxy. From photometrical and chemical evolution models constructed in the first part of the study (Tamm, Tempel & Tenjes 2007 (arXiv:0707.4375), hereafter Paper I) we can calculate the mass distribution of visible components (the bulge, the disc, the stellar halo, the outer diffuse stellar halo). In the dynamical model we calculate stellar rotation velocities along the major axis and velocity dispersions along the major, minor and intermediate axes of the galaxy assuming triaxial velocity dispersion ellipsoid. Comparing the calculated values with the collected observational data, we find the amount of DM, which must be added to reach an agreement with the observed rotation and dispersion data. We conclude that within the uncertainties, the DM distributions by Moore, Burkert, Navarro, Frenk & White (NFW) and the Einasto fit with observations nearly at all distances. The NFW and Einasto density distributions give the best fit with observations. The total mass of M 31 with the NFW DM distribution is 1.19*10^12 M_sun, the ratio of the DM mass to the visible mass is 10.0. For the Einasto DM distribution, these values are 1.28*10^12 M_sun and 10.8. The ratio of the DM mass to the visible mass inside the Holmberg radius is 1.75 for the NFW and the Einasto distributions. For different cuspy DM distributions, the virial mass is in a range 6.9-7.9*10^11 M_sun and the virial radius is ~150 kpc. The DM mean densities inside 10 pc for cusped models are 33 and 16 M_sun pc^-3 for the NFW and the Einasto profiles, respectively. For the cored Burkert profile, this value is 0.06 M_sun pc^-3.

We extract groups of galaxies from the SDSS Data Release 5 with the purpose of studying the supercluster-void network and environmental properties of groups therein. We use a modified friends-of-friends (FoF) method with adopted variable linking length in transverse and radial direction to eliminate selection effects and to find reliably as many groups as possible to track the supercluster network. We take into account various selection effects due to the use of a magnitude limited sample. To determine linking length scaling we study the luminosity-density relation in observed groups. We follow the changes in group sizes and mean galaxy number densities within groups when shifting nearby groups to larger distances. As a result we show that the linking length should be a slowly growing function with distance. Our final sample contains 17143 groups in the equatorial, and 33219 groups in the northern part of the DR5 survey with membership N_g >= 2. The group catalogue is available at our web-site http://www.obs.ee/~erik/index.html. Due to a narrow magnitude window in the SDSS the group catalogue based on this survey has been obtained by moderately growing linking length scaling law up to redshift z = 0.12. Above this redshift the scaling law turns down. In the redshift range z=0.12 - 0.2 only the cores are detected. Along with applying weights when calculating luminosities it is possible to use groups for determination of the large-scale luminosity-density field.

In the present paper we develop an algorithm allowing to calculate line-of-sight velocity dispersions in an axisymmetric galaxy outside of the galactic plane. When constructing a self-consistent model, we take into account the galactic surface brightness distribution, stellar rotation curve and velocity dispersions. This algorithm is applied to a Sa galaxy NGC 4594 = M 104, for which there exist velocity dispersion measurements outside of the galactic major axis. The mass distribution model is constructed in two stages. In the first stage we construct a luminosity distribution model, where only galactic surface brightness distribution is taken into account. Thereafter, in the second stage we develop on the basis of the Jeans equations a detailed mass distribution model and calculate line-of-sight velocity dispersions and the stellar rotation curve. Here a dark matter halo is added to visible components. Calculated dispersions are compared with observations along different slit positions perpendicular and parallel to the galactic major axis. In the best-fitting model velocity dispersion ellipsoids are radially elongated. Outside the galactic plane velocity dispersion behaviour is more sensitive to the dark matter density distribution and allows to estimate dark halo parameters.