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.

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.

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)

This paper considers a new method for the binary asteroid orbit determination problem. The method is based on the Bayesian approach with a global optimisation algorithm. The orbital parameters to be determined are modelled through an a posteriori distribution made of a priori and likelihood terms. The first term constrains the parameters space and it allows the introduction of available knowledge about the orbit. The second term is based on given observations and it allows us to use and compare different observational error models. Once the a posteriori model is built, the estimator of the orbital parameters is computed using a global optimisation procedure: the simulated annealing algorithm. The maximum a posteriori (MAP) techniques are verified using simulated and real data. The obtained results validate the proposed method. The new approach guarantees independence of the initial parameters estimation and theoretical convergence towards the global optimisation solution. It is particularly useful in these situations, whenever a good initial orbit estimation is difficult to get, whenever observations are not well-sampled, and whenever the statistical behaviour of the observational errors cannot be stated Gaussian like.

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.

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.

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

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.

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.

Context. The main feature of the spatial large-scale galaxy distribution is an intricate network of galaxy filaments. Although many attempts have been made to quantify this network, there is no unique and satisfactory recipe for that yet. Aims. The present paper compares the filaments in the real data and in the numerical models, to see if our best models reproduce statistically the filamentary network of galaxies. Methods. We apply an object point process with interactions (the Bisous process) to trace and describe the filamentary network both in the observed samples (the 2dFGRS catalogue) and in the numerical models that have been prepared to mimic the data.We compare the networks. Results. We find that the properties of filaments in numerical models (mock samples) have a large variance. A few mock samples display filaments that resemble the observed filaments, but usually the model filaments are much shorter and do not form an extended network. Conclusions. We conclude that although we can build numerical models that are similar to observations in many respects, they may fail yet to explain the filamentary structure seen in the data. The Bisous-built filaments are a good test for such a structure.

This paper tackles important aspects of comets dynamics from a statistical point of view. Existing methodology uses numerical integration for computing planetary perturbations for simulating such dynamics. This operation is highly computational. It is reasonable to wonder whenever statistical simulation of the perturbations can be much more easy to handle. The first step for answering such a question is to provide a statistical study of these perturbations in order to catch their main features. The statistical tools used are order statistics and heavy tail distributions. The study carried out indicated a general pattern exhibited by the perturbations around the orbits of the important planet. These characteristics were validated through statistical testing and a theoretical study based on Opik theory.

We propose to apply an object point process to automatically delineate filaments of the large-scale structure in redshift catalogues. We illustrate the feasibility of the idea on an example of the recent 2dF Galaxy Redshift Survey, describe the procedure, and characterise the results.

We propose to apply a marked point process to automatically delineate filaments of the large-scale structure in redshift catalogues. We illustrate the feasibility of the idea on an example of simulated catalogues, describe the procedure, and characterize the results. We find the distribution of the length of the filaments, and suggest how to use this approach to obtain other statistical characteristics of filamentary networks.