We address the 'plane of satellites problem' by studying planar configurations around two disc galaxies with no late major mergers, formed in zoom-in hydro-simulations. Due to the current lack of good quality kinematic data for M31 satellites, we use only positional information. So far, positional analyses of simulations are unable to find planes as thin and populated as the observed ones. Moreover, they miss systematicity and detail in the plane-searching techniques, as well as in the study of the properties and quality of planes, both in simulations or real data. To fill this gap, i) we extend the 4-galaxy-normal density plot method (Pawlowski et al. 2013) in a way designed to efficiently identify the best quality planes (i.e., thin and populated) without imposing extra constraints on their properties, and ii), we apply it for the first time to simulations. Using zoom-in simulations allows us to mimic MW/M31-like systems regarding the number of satellites involved as well as the galactic disc mass and morphology, in view of possible disc effects. At all timesteps analyzed in both simulations we find satellite planar configurations that are compatible, along given time intervals, with all the spatial characteristics of observed planes identified using the same methodology. However, the fraction of co-orbiting satellites within them is in general low, suggesting time-varying satellite membership. We conclude that high-quality positional planes of satellites are not infrequent in LCDM-formed disc galaxies with a quiet assembly history. Detecting kinematically-coherent, time-persistent planes demands considering the full six-dimensional phase-space information of satellites.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

We compute and study the infrared and sub-mm properties of high redshift ($z \gtrsim 1$) simulated clusters and proto-clusters. The results of a large set of hydro-dynamical zoom-in simulations including active galactic nuclei (AGN) feedback, have been treated with the recently developed radiative transfer code GRASIL-3D, which accounts for the effect of dust reprocessing in an arbitrary geometry. Here, we have slightly generalized the code to adapt it to the present purpose. Then we have post-processed boxes of physical size 2 Mpc encompassing each of the 24 most massive clusters identified at z=0, at several redshifts between 0.5 and 3, producing IR and sub-mm mock images of these regions and SEDs of the radiation coming out from them. While this field is in its infancy from the observational point of view, rapid development is expected in the near future thanks to observations performed in the far IR and sub-mm bands. Notably, we find that in this spectral regime our prediction are little affected by the assumption required by this post-processing, and the emission is mostly powered by star formation rather than accretion onto super massive black hole (SMBH). The comparison with the little observational information currently available, highlights that the simulated cluster regions never attain the impressive star formation rates suggested by these observations. This problem becomes more intriguing taking into account that the brightest cluster galaxies (BCGs) in the same simulations turn out to be too massive. It seems that the interplay between the feedback schemes and the star formation model should be revised, possibly incorporating a positive feedback mode.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

With the ever increasing size and complexity of fully self-consistent simulations of galaxy formation within the framework of the cosmic web, the demands upon object finders for these simulations has simultaneously grown. To this extent we initiated the Halo Finder Comparison Project that gathered together all the experts in the field and has so far led to two comparison papers, one for dark matter field haloes (Knebe et al. 2011), and one for dark matter subhaloes (Onions et al. 2012). However, as state-of-the-art simulation codes are perfectly capable of not only following the formation and evolution of dark matter but also account for baryonic physics (e.g. hydrodynamics, star formation, feedback) object finders should also be capable of taking these additional processes into consideration. Here we report on a comparison of codes as applied to the Constrained Local UniversE Simulation (CLUES) of the formation of the Local Group which incorporates much of the physics relevant for galaxy formation. We compare both the properties of the three main galaxies in the simulation (representing the MW, M31, and M33) as well as their satellite populations for a variety of halo finders ranging from phase-space to velocity-space to spherical overdensity based codes, including also a mere baryonic object finder. We obtain agreement amongst codes comparable to (if not better than) our previous comparisons, at least for the total, dark, and stellar components of the objects. However, the diffuse gas content of the haloes shows great disparity, especially for low-mass satellite galaxies. This is primarily due to differences in the treatment of the thermal energy during the unbinding procedure. We acknowledge that the handling of gas in halo finders is something that needs to be dealt with carefully, and the precise treatment may depend sensitively upon the scientific problem being studied.

We have implemented a chemical evolution model on the parallel AP3M+SPH DEVA code which we use to perform high resolution simulations of spiral galaxy formation. It includes feedback by SNII and SNIa using the Qij matrix formalism. We also include a diffusion mechanism that spreads newly introduced metals. The gas cooling rate depends on its specific composition. We study the stellar populations of the resulting bulges finding a potential scenario where they seem to be composed of two populations: an old, metal poor, $\alpha$-enriched population, formed in a multiclump scenario at the beginning of the simulation and a younger one, formed by slow accretion of satellites or gas, possibly from the disk due to instabilities.

We present a series of cosmological N-body simulations which make use of the hydrodynamic approach to the evolution of structures (Dominguez 2000). This approach addresses explicitly the existence of a finite spatial resolution and the dynamical effect of subresolution degrees of freedom. We adapt this method to cosmological simulations of the standard LCDM structure formation scenario and study the effects induced at redshift z=0 by this novel approach on the large-scale clustering patterns as well as (individual) dark matter halos. Comparing these simulations to usual N-body simulations, we find that (i) the new (hydrodynamic) model entails a proliferation of low--mass halos, and (ii) dark matter halos have a higher degree of rotational support. These results agree with the theoretical expectation about the qualitative behaviour of the "correction terms" introduced by the hydrodynamic approach: these terms act as a drain of inflow kinetic energy and a source of vorticity by the small-scale tidal torques and shear stresses.

We present results of a study of the stellar age distributions in the sample of elliptical-like objects (ELOs) identified at z=0 in four simulations operating in the context of a concordance cosmological model. The simulations show that the formation of most stars in each ELO of the sample is a consequence of violent dynamical events, either fast multiclump collapse at high z, or mergers at lower z. This second way can explain the age spread as well as the dynamical peculiarities observed in some ellipticals, but its relative weight is never dominant and decreases as the ELO mass at the halo scale, $M_{vir}$, increases, to such an extent that some recent mergers contributing an important fraction to the total ELO mass can possibly contribute only a small fraction of new born stars. More massive objects have older means and narrower spreads in their stellar age distributions than less massive ones. The ELO sample shows also a tight correlation between $M_{vir}$ and the central stellar l.o.s. velocity dispersion, $\sigma_{los}$. This gives a trend of the means and spreads of ELO stellar populations with $\sigma_{los}$ that is consistent, even quantitatively, with the age effects observationally detected in the stellar populations of elliptical galaxies. Therefore, these effects can be explained as the observational manifestation of the intrinsic correlations found in the ELO sample between $M_{vir}$ and the properties of the stellar age distribution, on the one hand, and $M_{vir}$ and $\sigma_{los}$, on the other hand. These correlations hint, for the first time, at a possible way to reconcile age effects in ellipticals, and, particularly, the increase of $\alpha / <Fe>$ ratios with $\sigma_{los}$, with the hierarchical clustering paradigm.

We describe DEVA, a multistep AP3M-like-SPH code particularly designed to study galaxy formation and evolution in connection with the global cosmological model. This code uses a formulation of SPH equations which ensures both energy and entropy conservation by including the so-called $\bn h$ terms. Particular attention has also been paid to angular momentum conservation and to the accuracy of our code. We find that, in order to avoid unphysical solutions, our code requires that cooling processes must be implemented in a non-multistep way. We detail various cosmological simulations which have been performed to test our code and also to study the influence of the $\bn h$ terms. Our results indicate that such correction terms have a non-negligible effect on some cosmological simulations, especially on high density regions associated either to shock fronts or central cores of collapsed objects. Moreover, they suggest that codes paying a particular attention to the implementation of conservation laws of physics at the scales of interest, can attain good accuracy levels in conservation laws with limited computational resources.

We use cosmological SPH simulations to study the effects of mergers in the star formation history of galactic objects in hierarchical clustering scenarios. We find that during some merger events, gaseous discs can experience two starbursts: the first one during the orbital decay phase, due to gas inflows driven as the satellite approaches, and the second one, when the two baryonic clumps collide. A trend for these first induced starbursts to be more efficient at transforming the gas into stars is also found. We detect that systems which do not experience early gas inflows have well-formed stellar bulges and more concentrated potential wells, which seem to be responsible for preventing further gas inward transport triggered by tidal forces. Our results constitute the first proof that bulges can form as the product of collapse, collisions and secular evolution in a cosmological framework, and they are consistent with a rejuvenation of the stellar population in bulges at intermediate z with, at least, 50% of the stars (in SCDM) being formed at high z. (Abridged)

We present results from a careful and detailed analysis of the structural and dynamical properties of a sample of 29 disc-like objects identified at z=0 in three AP3M-SPH fully consistent cosmological simulations. These simulations are realizations of a CDM hierarchical model, where an inefficient Schmidt law-like algorithm to model the stellar formation process has been implemented. We focus on properties that can be constrained with available data from observations of spiral galaxies, namely, the bulge and disc structural parameters and the rotation curves. Comparisons with data from Broeils (1992), de Jong (1996) and Courteau (1996, 1997) give satisfactory agreement, in contrast with previous findings using other codes. This suggests that the stellar formation implementation we have used has succeded in forming compact bulges that stabilize disc-like structures in the violent phases of their assembly, while in the quiescent phases the gas has cooled and collapsed according with the Fall & Efstathiou standard model of disc formation.

We analyze the structural and dynamical properties of disk-like objects formed in fully consistent cosmological simulations with an inefficient star formation algorithm. Comparison with data of similar observable properties of spiral galaxies gives satisfactory agreement.

We analyze the structural and dynamical properties of disk-like objects formed in fully consistent cosmological simulations which include inefficient star formation. Comparison with data of similar observable properties of spiral galaxies gives satisfactory agreement, in contrast with previous findings using other codes. This suggests that the stellar formation implementation used has allowed the formation of disks as well as guaranteed their stability.

We report on results of N-body simulations aimed at testing the hypothesis that galaxies in X-ray emitting (i.e., virialized) Compact Groups are not tidally stripped when they are embedded in a common, massive, quiescent dark matter halo. To disentangle the effects of interactions from spurious effects due to an incorrect choice of the initial galaxy model configurations, these have been chosen to be tidally-limited King spheres, representing systems in quasi-equilibrium within the tidal field of the halo. The potential of the halo has been assumed to be frozen and the braking due to dynamical friction neglected. Our results confirm the hypothesis of low rates of tidal stripping and suggest a scenario for virialized Compact Group evolution in their quiescent phases with only very moderate tidally induced galaxy evolution can be generally expected. This implies the group stability, provided that the dynamical friction timescales in these systems are not much shorter than the Hubble time. We discuss briefly this possibility, in particular taking account of the similarity between the velocity dispersions of a typical virialized Compact Groups and the internal velocity dispersion of typical member galaxies. A number of puzzling observational data on Compact Groups can be easily explained in this framework. Other observations would be better understood as the result of enhanced merging activity in the proto-group environment, leading to virialized Compact Group formation through mergers of lower mass halos, as predicted by hierarchical scenarios of structure formation.

The possibility that quasi-equilibrium self-gravitating galaxy-like configurations exist in a tidal field is analyzed in this paper. More specifically, we address the question of how to predict initial configurations modeling galaxies that are able to survive environmental effects in a dense environment for a Hubble time or so, provided thay dynamical friction is neglected. The initial quasi-equilibrium configurations have been built up taking into account the external tidal field produced by the halo. It modifies the escape velocity field of the configuration, compared with isolated configurations. The survival of the configurations as they orbit inside the halos has been studied through N-body simulations. As a general result, it has been found out that the bulk of the models is conserved along 12.5 Gyears of evolution, and that the low rates of mass losses they experience are consistent with those expected when the adiabatic protection hypothesis is at work. So, solutions for galaxy configurations in tidal quasi-equilibrium have been found, showing that tidal stripping in quiescent phases does not seem to be very important, unless that the density of the galaxy environment at its formation had been much lower than that of the galaxy environment at the point of its orbit where the tidal perturbation is maximum. (Abridged)

We report results on the formation of disk-like structures in two cosmological hydrodynamical simulations in a hierarchical clustering scenario, sharing the same initial conditions. In the first one, a simple and generic implementation of star formation has allowed galaxy-like objects with stellar bulges and extended, populated disks to form. Gas in the disk comes from both, particles that survive mergers keeping in part their angular momentum content, and new gas supply by infall, once the merger process is over, with global specific angular momentum conservation. The stellar bulge forms from gas that has lost most of its angular momentum. In the second simulation, no star formation has been included. In this case, objects consist of an overpopulated central gas concentration, and an extended, underpopulated disk. The central concentration forms from particles that suffer an important angular momentum loss in violent events, and it often contains more than 70% of the object's baryonic mass. The external disk forms by late infall of gas, that roughly conserves its specific angular momentum. The difference between these two simulations is likely to be due to the stabilizing character of the stellar bulge-like cores that form in the first simulation, which diminishes the inflow of gas triggered by mergers and interactions.

We present results that suggest that disc angular momentum catastrophes (DAMCs) plaguing hierarchical hydrodynamical simulations of galaxy formation, not including star formation processes, are due to the inward gas transport that follows bar disk instabilities triggered by interactions and mergers. They also show that DAMCs can be easily avoided by including star-forming processes, as they lead unavoidably to the formation of compact stellar bulges that stabilize disks against bars. The formation of disks similar to those observed demands, in addition, that not all the gas is depleted into stars at high z, so that they can be formed at lower z.

A statistical method is developed to infer the typical density profiles of poor galaxy systems without resort to binning of data or assuming a given center to each system. The method is applied to the accordant redshift quartets in Hickson's compact groups (HCGs). The distribution of separations in these groups is consistent with a unique generalized modified Hubble surface density profile, with best-fit asymptotic slope beta = 1.4 and core radius R_c = 18 h^{-1} kpc, although a King profile (beta = 1) is also consistent with the data (with R_c = 6 h^{-1} kpc). These distributions are more concentrated than what has been previously determined for these groups. HCGs are unlikely to act as strong gravitational lenses, but analogous systems 5 to 10 times more distant should produce a non-negligible fraction of giant arcs.

A detailed analysis of the galaxy distribution in the Southern Sky Redshift Survey (SSRS) by means of the multifractal or scaling formalism is presented. It is shown that galaxies cluster in different ways according to their morphological type as well as their size. Ellipticals are more clustered than spirals, even at scales up to 15 h$^{-1}$ Mpc, whereas no clear segregation between early and late spirals is found. It is also shown that smaller galaxies distribute more homogeneously than larger galaxies.

We present the results of an exhaustive analysis of the morphological segregation of galaxies in the CfA and SSRS catalogs through the scaling formalism. Morphological segregation between ellipticals and spirals has been detected at scales up to 15-20 h$^{-1}$ Mpc in the CfA catalog, and up to 20-30 h$^{-1}$ Mpc in the SSRS catalog. Moreover, it is present not only in the densest areas of the galaxy distribution, but also in zones of moderate density.