We present a detailed characterization of planes of satellites in the Milky Way (MW) and M31 systems. To this end we introduce an extension to the '4-galaxy-normal density plot' method \citep{Pawlowski13}, by which plot over-densities signal the normal direction to predominant planes of satellites within a given sample. For a given over-density, the extension provides a \textit{collection} of planes, each including a different number of objects $N_{\rm sat}$. We apply this method to the position data of confirmed MW and M31 satellites and quantify the quality of planes through the outputs of a Tensor of Inertia plane-fitting technique. Plane quality is quantified in terms of population ($N_{\rm sat}$) and flattening (the short-to-long axis ratio $c/a$ or the rms thickness normal to the plane). Therefore, planes with the same population or flattening can be compared with each other allowing us to single-out best-quality planes. For the first time, we study the second-most predominant planar configuration of satellites in M31, singling out a plane with 18 satellite members that shows a quality comparable to the Great Plane of Andromeda (GPoA, with $N_{\rm sat}=19$) despite it being more affected by distance uncertainties. This structure is viewed nearly face-on from the MW and is approximately normal to the GPoA. Overall, we find planes of satellites around the MW and M31 with higher qualities than those previously reported with a given $N_{\rm sat}$. We also show that mass plays no role in determining a satellite's membership or not to the respective best-quality planes.

We study the morphological and structural properties of the host galaxies associated with 57 optically-selected luminous type 2 AGN at $z\sim$0.3-0.4: 16 high-luminosity Seyfert 2 (HLSy2, 8.0$\le$log($L_{\rm [OIII]}/L_{\odot})<$8.3) and 41 obscured quasars (QSO2, log($L_{\rm [OIII]}/L_{\odot})\ge$8.3). With this work, the total number of QSO2 at $z<1$ with parametrized galaxies increases from $\sim$35 to 76. Our analysis is based on HST WFPC2 and ACS images that we fit with {\sc GALFIT}. HLSy2 and QSO2 show a wide diversity of galaxy hosts. The main difference lies in the higher incidence of highly-disturbed systems among QSO2. This is consistent with a scenario in which galaxy interactions are the dominant mechanism triggering nuclear activity at the highest AGN power. There is a strong dependence of galaxy properties with AGN power (assuming $L_ {\rm [OIII]}$ is an adequate proxy). The relative contribution of the spheroidal component to the total galaxy light (B/T) increases with $L_ {\rm [OIII]}$. While systems dominated by the spheoridal component spread across the total range of $L_ {\rm [OIII]}$, most disk-dominated galaxies concentrate at log($L_{\rm [OIII]}/L_{\odot})<$8.6. This is expected if more powerful AGN are powered by more massive black holes which are hosted by more massive bulges or spheroids. The average galaxy sizes ($\langle r_{\rm e} \rangle$) are 5.0$\pm$1.5 kpc for HLSy2 and 3.9$\pm$0.6 kpc for HLSy2 and QSO2 respectively. These are significantly smaller than those found for QSO1 and narrow line radio galaxies at similar $z$. We put the results of our work in context of related studies of AGN with quasar-like luminosities.

The significant diversity of rotation curve (RC) shapes in dwarf galaxies has recently emerged as a challenge to LCDM: in dark matter (DM) only simulations, DM halos have a universal cuspy density profile that results in self-similar RC shapes. We compare RC shapes of simulated galaxies from the NIHAO project with observed galaxies from the homogeneous SPARC dataset. The DM halos of the NIHAO galaxies can expand to form cores, with the degree of expansion depending on their stellar-to-halo mass ratio. By means of the $V_{\rm 2kpc}-V_{\rm Rlast}$ relation (where $V_{\rm Rlast}$ is the outermost measured rotation velocity), we show that both the average trend and the scatter in RC shapes of NIHAO galaxies are in reasonable agreement with SPARC: this represents a significant improvement compared to simulations that do not result in DM core formation, suggesting that halo expansion is a key process in matching the diversity of dwarf galaxy RCs. Note that NIHAO galaxies can reproduce even the extremely slowly rising RCs of IC 2574 and UGC 5750. Revealingly, the range where observed galaxies show the highest diversity corresponds to the range where core formation is most efficient in NIHAO simulations, 50$<V_{\rm Rlast}$\km s$^{-1}<$100. A few observed galaxies in this range cannot be matched by any NIHAO RC nor by simulations that predict a universal halo profile. Interestingly, the majority of these are starbursts or emission-line galaxies, with steep RCs and small effective radii. Such galaxies represent an interesting observational target providing new clues to the process/viability of cusp-core transformation, the relationship between starburst and inner potential well, and the nature of DM.

Recent determinations of the radial distributions of mono-metallicity populations (MMPs, i.e., stars in narrow bins in [Fe/H] within wider [$\alpha$/Fe] ranges) by the SDSS-III/APOGEE DR12 survey cast doubts on the classical thin - thick disk dichotomy. The analysis of these observations lead to the non-$[\alpha$/Fe] enhanced populations splitting into MMPs with different surface densities according to their [Fe/H]. By contrast, $[\alpha$/Fe] enhanced (i.e., old) populations show an homogeneous behaviour. We analyze these results in the wider context of disk formation within non-isolated halos embedded in the Cosmic Web, resulting in a two-phase mass assembly. By performing hydrodynamical simulations in the context of the $\rm \Lambda CDM$ model, we have found that the two phases of halo mass assembly (an early, fast phase, followed by a slow one, with low mass assembly rates) are very relevant to determine the radial structure of MMP distributions, while radial mixing has only a secondary role, depending on the coeval dynamical and/or destabilizing events. Indeed, while the frequent dynamical violent events occuring at high redshift remove metallicity gradients, and imply efficient stellar mixing, the relatively quiescent dynamics after the transition keeps [Fe/H] gaseous gradients and prevents newly formed stars to suffer from strong radial mixing. By linking the two-component disk concept with the two-phase halo mass assembly scenario, our results set halo virialization (the event marking the transition from the fast to the slow phases) as the separating event marking periods characterized by different physical conditions under which thick and thin disk stars were born.

Many massive star forming disc galaxies in the redshift range 3 to 0.5 are observed to have a clumpy morphology showing giant clumps of size $\sim$1 kpc and masses of about $10^7M_{\odot}$ to $10^{10} M_{\odot}$. The nature and fate of these giant clumps is still under debate. In this work we use 19 high-resolution simulations of disc galaxies from the NIHAO sample to study the formation and the evolution of clumps in the discs of high redshift galaxies. We use mock HST - CANDELS observations created with the radiative transfer code GRASIL-3D to carry out, for the first time, a quantitative comparison of the observed fraction of clumpy galaxies and its evolution with redshift with simulations. We find a good agreement between the observed clumpy fraction and the one of the NIHAO galaxies. We find that dust attenuation can suppress intrinsically bright clumps and enhance less luminous ones. In our galaxy sample we only find clumps in light (u-band) from young stars but not in stellar mass surface density maps. This means that the NIHAO sample does not show clumpy stellar discs but rather a clumpy light distribution originating from clumpy star formation events. The clumps found in the NIHAO sample match observed age/color gradients as a function of distance from the galaxy center but they show no sign of inward migration. Clumps in our simulations disperse on timescales of a about a hundred Myr and their contribution to bulge growth is negligible.

The spectral energy distributions (SEDs) of low-mass low-metallicity (dwarf) galaxies are a challenging piece of the puzzle of galaxy formation in the near Universe. These SEDs show some particular features in the submillimeter to far-infrared wavelength range compared to normal, larger, and metal-richer galaxies that cannot be explained by the current models. These can be summarized as: a broadening of the IR peak, which implies a warmer dust component; an excess of emission in the submm ($\sim$500 $\mu$m), that causes a flattening of the submm/FIR slope; and a very low intensity of PAH emission features. With the aim of explaining these emission patterns, the SEDs of a sample of 27 simulated dwarf galaxies were calculated using the GRASIL-3D radiation transfer code. This code has the particularity that it separately treats the radiative transfer through dust grains within molecular clouds and within the cirrus, the dense and diffuse components of the gas phase, respectively. The simulated galaxies have stellar masses ranging from 10$^6$-10$^9$ M$_\odot$, and were obtained from a single simulation run within a Local Group environment with initial conditions from the CLUES project. We report a study of the IRAS, Spitzer and Herschel bands luminosities, and of the SFRs, dust, and gas (HI and H$_2$) mass contents. We find a satisfactory agreement with observational data, with GRASIL-3D naturally reproducing the specific spectral features mentioned above. We conclude that the GRASIL-3D two-component dust model gives a physical interpretation of the emission of dwarf galaxies: molecular clouds and cirrus represent the warm and cold dust components, respectively, needed to reproduce observational data.

We present predictions of Spectral Energy Distributions (SEDs), from the UV to the FIR, of simulated galaxies at $z=0$. These were obtained by post-processing the results of an N-body+hydro simulation of a small cosmological volume, that uses the Multi-Phase Particle Integrator (MUPPI) for star formation and stellar feedback, with the GRASIL-3D radiative transfer code, that includes reprocessing of UV light by dust. Physical properties of galaxies resemble observed ones, though with some tension at small and large stellar masses. Comparing predicted SEDs of simulated galaxies with different samples of local galaxies, we find that these resemble observed ones, when normalised at 3.6 $\mu$m. A comparison with the Herschel Reference Survey shows that, when binning galaxies in Star Formation Rate (SFR), average SEDs are reproduced to within a factor of $\sim2$ even in normalization, while binning in stellar mass highlights the same tension that is present in the stellar mass -- SFR plane. We use our sample to investigate the correlation of IR luminosity in Spitzer and Herschel bands with several galaxy properties. SFR is the quantity that best correlates with IR light up to $160\ \mu$m, while at longer wavelengths better correlations are found with molecular mass and, at $500\ \mu$m, with dust mass. However, using the position of the FIR peak as a proxy for cold dust temperature, we assess that heating of cold dust is mostly determined by SFR, with stellar mass giving only a minor contribution. We finally show how our sample of simulated galaxies can be used as a guide to understand the physical properties and selection biases of observed samples.

Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulated spiral galaxies with masses comparable to the Milky Way. The simulations span a range of star formation and energy feedback strengths and prescriptions, particle- and grid-based hydrodynamical implementations, as well as initial conditions/assembly history. Results. Consistently, we find that the steeper, negative, radial metallicity gradients seen in the mid-plane flatten with increasing height away from the plane. In simulations with stronger (and/or more spatially-extended) feedback, the negative radial gradients invert, becoming positive for heights in excess of 1 kpc. Such behaviour is consistent with that inferred from recent observations. Our measurements of the vertical metallicity gradients show no clear correlation with galactocentric radius, and are in good agreement with those observed in the Milky Way's thick disc (locally). Conclusions. Simulations employing stronger/more extended feedback prescriptions possess radial and vertical metallicity and kinematic gradients more in line with recent observations. The inverted, positive, radial metallicity gradients seen in the simulated thick stellar discs originate from a population of younger, more metal-rich, stars formed in-situ, superimposed upon a background population of older migrators from the inner disc; the contrast provided by the former increases radially, due to the inside-out growth of the disc. A similar behaviour may be responsible for the same flattening seen in the radial gradients with scaleheight in the Milky Way.

Using 22 hydrodynamical simulated galaxies in a LCDM cosmological context we recover not only the observed baryonic Tully-Fisher relation, but also the observed "mass discrepancy--acceleration" relation, which reflects the distribution of the main components of the galaxies throughout their disks. This implies that the simulations, which span the range 52 < V$_{\rm flat}$ < 222 km/s where V$_{\rm flat}$ is the circular velocity at the flat part of the rotation curve, and match galaxy scaling relations, are able to recover the observed relations between the distributions of stars, gas and dark matter over the radial range for which we have observational rotation curve data. Furthermore, we explicitly match the observed baryonic to halo mass relation for the first time with simulated galaxies. We discuss our results in the context of the baryon cycle that is inherent in these simulations, and with regards to the effect of baryonic processes on the distribution of dark matter.

We present a detailed analysis of the local evolution of 206 Lagrangian Volumes (LVs) selected at high redshift around galaxy seeds, identified in a large-volume $\Lambda$ cold dark matter ($\Lambda$CDM) hydrodynamical simulation. The LVs have a mass range of $1 - 1500 \times 10^{10} M_\odot$. We follow the dynamical evolution of the density field inside these initially spherical LVs from $z=10$ up to $z_{\rm low} = 0.05$, witnessing highly non-linear, anisotropic mass rearrangements within them, leading to the emergence of the local cosmic web (CW). These mass arrangements have been analysed in terms of the reduced inertia tensor $I_{ij}^r$, focusing on the evolution of the principal axes of inertia and their corresponding eigendirections, and paying particular attention to the times when the evolution of these two structural elements declines. In addition, mass and component effects along this process have also been investigated. We have found that deformations are led by dark matter dynamics and they transform most of the initially spherical LVs into prolate shapes, i.e. filamentary structures. An analysis of the individual freezing-out time distributions for shapes and eigendirections shows that first most of the LVs fix their three axes of symmetry (like a skeleton) early on, while accretion flows towards them still continue. Very remarkably, we have found that more massive LVs fix their skeleton earlier on than less massive ones. We briefly discuss the astrophysical implications our findings could have, including the galaxy mass-morphology relation and the effects on the galaxy-galaxy merger parameter space, among others.

Using cosmological galaxy simulations from the MaGICC project, we study the evolution of the stellar masses, star formation rates and gas phase abundances of star forming galaxies. We derive the stellar masses and star formation rates using observational relations based on spectral energy distributions by applying the new radiative transfer code GRASIL-3D to our simulated galaxies. The simulations match well the evolution of the stellar mass-halo mass relation, have a star forming main sequence that maintains a constant slope out to redshift z $\sim$ 2, and populate projections of the stellar mass - star formation - metallicity plane, similar to observed star forming disc galaxies. We discuss small differences between these projections in observational data and in simulations, and the possible causes for the discrepancies. The light-weighted stellar masses are in good agreement with the simulation values, the differences between the two varying between 0.06 dex and 0.20 dex. We also find a good agreement between the star formation rate tracer and the true (time-averaged) simulation star formation rates. Regardless if we use mass- or light-weighted quantities, our simulations indicate that bursty star formation cycles can account for the scatter in the star forming main sequence.

We introduce a new model for the spectral energy distribution of galaxies, GRASIL-3D, which includes a careful modelling of the dust component of the interstellar medium. GRASIL-3D is an entirely new model based on the formalism of an existing and widely applied spectrophotometric model, GRASIL, but specifically designed to be interfaced with galaxies with any arbitrarily given geometry, such as galaxies calculated by theoretical hydrodynamical galaxy formation codes. GRASIL-3D is designed to separately treat radiative transfer in molecular clouds and in the diffuse cirrus component. The code has a general applicability to the outputs of simulated galaxies, either from Lagrangian or Eulerian hydrodynamic codes. As an application, the new model has been interfaced to the P-DEVA and GASOLINE smoothed-particle hydrodynamic codes, and has been used to calculate the spectral energy distribution for a variety of simulated galaxies from UV to sub-millimeter wavelengths, whose comparison with observational data gives encouraging results. In addition, GRASIL-3D allows 2D images of such galaxies to be obtained, at several angles and in different bands.

We analyze and compare the bulges of a sample of L* spiral galaxies in hydrodynamical simulations in a cosmological context, using two different codes, P-DEVA and GASOLINE. The codes regulate star formation in very different ways, with P-DEVA simulations inputing low star formation efficiency under the assumption that feedback occurs on subgrid scales, while the GASOLINE simulations have feedback which drives large scale outflows. In all cases, the marked knee-shape in mass aggregation tracks, corresponding to the transition from an early phase of rapid mass assembly to a later slower one, separates the properties of two populations within the simulated bulges. The bulges analyzed show an important early starburst resulting from the collapse-like fast phase of mass assembly, followed by a second phase with lower star formation, driven by a variety of processes such as disk instabilities and/or mergers. Classifying bulge stellar particles identified at z=0 into old and young according to these two phases, we found bulge stellar sub-populations with distinct kinematics, shapes, stellar ages and metal contents. The young components are more oblate, generally smaller, more rotationally supported, with higher metallicity and less alpha-element enhanced than the old ones. These results are consistent with the current observational status of bulges, and provide an explanation for some apparently paradoxical observations, such as bulge rejuvenation and metal-content gradients observed. Our results suggest that bulges of L* galaxies will generically have two bulge populations which can be likened to classical and pseudo-bulges, with differences being in the relative proportions of the two, which may vary due to galaxy mass and specific mass accretion and merger histories.

We present a detailed analysis of the global and fine structure of four middle-mass disc galaxies obtained from simulations in a $\Lambda$CDM scenario. These objects have photometric D/T ratios in good agreement with those observed for late-type spirals, as well as kinematic properties in agreement with the observational Tully-Fisher relation. We identify the different dynamical components at z=0 on the basis of both orbital parameters and the binding energy of stars in the galaxy. In this way, we recognize a slowly rotating centrally concentrated spheroid, and two disc components supported by rotation: a thin disc with stars in nearly circular orbits, and a thick disc with orbital parameters transitional between the thin disc and the spheroid. The spheroidal component is composed mainly by old, metal-poor and {\alpha}-enhanced stars. The distribution of metals in this component shows, however, a clear bimodality with a low-metallicity peak, which could be related to a classical bulge, and a high-metallicity peak, which could be related to a pseudo-bulge. The thin disc appears in our simulations as the youngest and most metal-rich component. The radial distribution of ages and colours in this component are U-shaped: the new stars are forming in the inner regions, and then migrate through secular processes. Finally, we also find a thick disc containing about 16% of the total stellar mass and with properties that are intermediate between those of the thin disc and the spheroid. Its low-metallicity stars are {\alpha}-enhanced when compared to thin disc stars of the same metallicity. The structural parameters (e.g., the scale height) of the simulated thick discs suggest that such a component could result from the combination of different thickening mechanisms that include merger-driven processes, but also long-lived internal perturbations of the thin disc. [Abridged]

The mass assembly and star formation histories of massive galaxies identified at low redshift z in different cosmological hydrodynamical simulations, have been studied through a detailed follow-up backwards in time of their constituent mass elements (sampled by particles) of different types. Then, the configurations they depict at progressively higher zs have been analysed. The analyses show that these histories share common generic patterns, irrespective of particular circumstances. In any case, the results we have found are different depending on the particle type. The most outstanding differences follow. We have found that by z ~ 3.5 - 6, mass elements identified as stellar particles at z=0 exhibit a gaseous cosmic-web-like morphology with scales of ~ 1 physical Mpc, where the densest mass elements have already turned into stars by z ~ 6. These settings are in fact the densest pieces of the cosmic web, where no hot particles show up, and dynamically organized as a hierarchy of flow convergence regions, that is, attraction basins for mass flows. On the other hand, mass elements identified at the diffuse hot coronae surrounding massive galaxies at z = 0, do not display a clear web-like morphology at any z. Diffuse gas is heated when flow convergence regions go through contractive deformations, and most of it keeps hot and with low density along the evolution. To shed light on the physical foundations of the behaviour our analyses show up, as well as on their possible observational implications, these patterns have been confronted with some generic properties of singular flows as described by the adhesion model. We have found that these common patterns simulations show can be interpreted as a consequence of flow properties, that, moreover, could explain different generic observational results on massive galaxies or their samples. We briefly discuss some of them.[Abridged]

Relaxed, massive galactic objects have been identified at redshifts z = 4;5; and 6 in hydrodynamical simulations run in a large cosmological volume. This allowed us to analyze the assembly patterns of the high mass end of the galaxy distribution at these high zs, by focusing on their structural and dynamical properties. Our simulations indicate that massive objects at high redshift already follow certain scaling relations. These relations define virial planes at the halo scale, whereas at the galactic scale they define intrinsic dynamical planes that are, however, tilted relative to the virial plane. Therefore, we predict that massive galaxies must lie on fundamental planes from their formation. We briefly discuss the physical origin of the tilt in terms the physical processes underlying massive galaxy formation at high z, in the context of a two-phase galaxy formation scenario. Specifically, we have found that it lies on the different behavior of the gravitationally heated gas as compared with cold gas previously involved in caustic formation, and the mass dependence of the energy available to heat the gas.

We present SPH cosmological simulations of the formation of three disk galaxies with a detailed treatment of chemical evolution and cooling. The resulting galaxies have properties compatible with observations: relatively high disk-to-total ratios, thin stellar disks and good agreement with the Tully-Fisher and the luminosity-size relations. They present a break in the luminosity profile at 3.0 +- 0.5 disk scale lengths, while showing an exponential mass profile without any apparent breaks, in line with recent observational results. Since the stellar mass profile is exponential, only differences in the stellar populations can be the cause of the luminosity break. Although we find a cutoff for the star formation rate imposed by a density threshold in our star formation model, it does not coincide with the luminosity break and is located at 4.3 +- 0.4 disk scale lengths, with star formation going on between both radii. The color profiles and the age profiles are "U-shaped", with the minimum for both profiles located approximately at the break radius. The SFR to stellar mass ratio increases until the break, explaining the coincidence of the break with the minimum of the age profile. Beyond the break we find a steep decline in the gas density and, consequently, a decline in the SFR and redder colors. We show that most stars (64-78%) in the outer disk originate in the inner disk and afterwards migrate there. Such stellar migrations are likely the main origin of the U-shaped age profile and, therefore, of the luminosity break.

Aims:In this study we investigate the evolution of shape and kinematics of elliptical galaxies in a cosmological framework. Methods: We use a set of hydrodynamic, self-consistent simulations operating in the context of a concordance cosmological model where relaxed elliptical-like objects (ELOs) were identified at redshifts z=0, z=0.5, z=1 and z=1.5. Results: The population of elliptical systems analysed here present a systematic change through time, i.e. evolution, by becoming rounder in general at z=0 and, at the same time more velocity dispersion supported. This is found to be primarily due to major dry mergers where only a modest amount of angular momentum is involved into the merger event. Despite the general trend, in a significant amount of cases the merger event involves a higher specific angular momentum, which in general causes the system to acquire a higher rotational support and/or a more oblate shape. These evolutionary patterns are still present when we study our systems in projection, mimicking real observations, and thus they should become apparent in future observations.

We describe an SPH model for chemical enrichment and radiative cooling in cosmological simulations of structure formation. This model includes: i) the delayed gas restitution from stars by means of a probabilistic approach designed to reduce the statistical noise and, hence, to allow for the study of the inner chemical structure of objects with moderately high numbers of particles; ii) the full dependence of metal production on the detailed chemical composition of stellar particles by using, for the first time in SPH codes, the Qij matrix formalism that relates each nucleosynthetic product to its sources; and iii) the full dependence of radiative cooling on the detailed chemical composition of gas particles, achieved through a fast algorithm using a new metallicity parameter zeta(T) that gives the weight of each element on the total cooling function. The resolution effects and the results obtained from this SPH chemical model have been tested by comparing its predictions in different problems with known theoretical solutions. We also present some preliminary results on the chemical properties of elliptical galaxies found in self-consistent cosmological simulations. Such simulations show that the above zeta-cooling method is important to prevent an overestimation of the metallicity-dependent cooling rate, whereas the Qij formalism is important to prevent a significant underestimation of the [alpha/Fe] ratio in simulated galaxy-like objects.

[Abridged] We have analysed the parameters characterising the mass, size and velocity dispersion both at the baryonic scale and at the halo scales of two samples of relaxed elliptical-like-objects (ELOs) identified, at z=0, in a set of self-consistent hydrodynamical simulations operating in the context of a concordance cosmological model. At the halo scale they have been found to satisfy virial relations; at the scale of the baryonic object the (logarithms of the) ELO stellar masses, projected stellar half-mass radii, and stellar central l.o.s. velocity dispersions define a flattened ellipsoid close to a plane (the intrinsic dynamical plane, IDP), tilted relative to the virial one, whose observational manifestation is the observed FP. The ELO samples have been found to show systematic trends with the mass scale in both, the relative content and the relative distributions of the baryonic and the dark mass ELO components, so that homology is broken in the spatial mass distribution (resulting in the IDP tilt), but ELOs are still a two-parameter family where the two parameters are correlated. The physical origin of these trends presumably lies in the systematic decrease, with increasing ELO mass, of the relative amount of dissipation experienced by the baryonic mass component along ELO stellar mass assembly. ELOs also show kinematical segregation, but it does not appreciably change with the mass scale. The non-homogeneous population of IDPs explains the role played by the virial mass to determine the correlations among intrinsic parameters. In this paper we also show that the central stellar line-of-sight velocity dispersion of ELOs, is a fair empirical estimator of the virial mass, and this explains the central role played by this quantity at determining the observational correlations.

[Abridged] We have analysed the mass and velocity distributions of two samples of relaxed elliptical-like-objects (ELOs) identified, at z=0, in a set of self-consistent hydrodynamical simulations operating in the context of a concordance cosmological model. Our analysis shows that they are embedded in extended, massive dark matter haloes, and they also have an extended corona of hot diffuse gas. Dark matter haloes have experienced adiabatic contraction along their assembly process. The relative ELO dark- to bright-mass content and space distributions show broken homology, and they are consistent with observational results on the dark matter fraction at the central regions, as well as on the gradients of the mass-to-light ratio profiles for boxy ellipticals, as a function of their stellar masses. These results indicate that massive ellipticals miss stars at their central regions, as compared with less massive ones. Our simulations indicate that these missing baryons could be found beyond the virial radii as a hot, diffuse plasma. The projected stellar mass profiles of our ELOs can be well fit by the S\'ersic (1968) formula. Parameters of these fits reproduce observational correlations of the shape parameter with size, luminosity and velocity dispersion. The total mass density profiles show a power-law behaviour over a large r/r_{vir} interval, consistent with data on massive lens ellipticals at shorter radii. The velocity dispersion profiles show kinematical segregation, with no systematic mass dependence and a positive anisotropy, roughly independent of the radial distance outside the central regions. The consistency with observations strongly suggests that they could also describe important intrinsic characteristics of real ellipticals, as well as some of their properties recently inferred from observational data.

In a series of recent papers, a new formalism has been developed that explains the inner structure of dark matter halos as collisionless, dissipationless systems assembled through mergers and accretion at the typical cosmological rate. Nearby ellipticals are also collisionless, dissipationless systems assembling their mass through mergers, but contrarily to the former structures they do not continuously accrete external matter because they are shielded by their host halos. Here we explore the idea that the infall of their own matter ejected within the halo on the occasion of a violent merger can play a role similar to external accretion in halos. The predicted stellar mass density profile fits the observed one, and the empirical total mass density profile is also recovered.

We present results of a study on the evolution of the parameters characterizing the structure and dynamics of the relaxed elliptical-like objects (ELOs) identified at z=0, z=1 and z=1.5 in a set of hydrodynamical, self-consistent simulations operating in the context of a concordance cosmological model. The values of the stellar mass, the stellar half-mass radius and the stellar mean-square velocity have been measured in each ELO and found to populate, at any z, a flattened ellipsoid close to a plane (the dynamical plane, DP). Our simulations indicate that, at the intermediate zs considered, individual ELOs evolve, increasing the values of these parameters as a consequence of on-going mass assembly, but, nevertheless, their DP is roughly preserved within its scatter, in agreement with observations of the Fundamental Plane of ellipticals at different zs. We briefly discuss how this lack of significant dynamical and structural evolution in ELO samples arises, in terms of the two different phases operating in the mass aggregation history of their dark matter halos. According with our simulations, most dissipation involved in ELO formation takes place at the early violent phase, causing the stellar mass, the stellar half-mass radius and the stellar mean-square velocity parameters to settle down to the DP, and, moreover, the transformation of most of the available gas into stars. In the subsequent slow phase, ELO stellar mass growth preferentially occurs through non-dissipative processes, so that the DP is preserved and the ELO star formation rate considerably decreases. These results hint, for the first time, to a possible way of explaining, in the context of cosmological simulations, different apparently paradoxical observational results on ellipticals.

We report on a study of the parameters characterizing the mass and velocity distributions of two samples of relaxed elliptical-like objects (ELOs) identified, at z=0, in a set of self-consistent hydrodynamical simulations operating in the context of a concordance cosmological model. Star formation (SF) has been phenomenologically implemented in the simulations in the framework of the turbulent sequential scenario through a threshold gas density and an efficiency parameter. Each ELO sample is characterized by the values these parameters take. We have found that the (logarithms of the) ELO stellar masses, projected stellar half-mass radii, and stellar central line-of-sight (LOS) velocity dispersions define dynamical fundamental planes (FPs). Zero points depend on the particular values that the SF parameters take, while slopes do not change. The ELO samples have been found to show systematic trends with the mass scale in both the relative content and the relative distributions of the baryonic and the dark mass ELO components. The physical origin of these trends lies in the systematic decrease, with increasing ELO mass, of the relative dissipation experienced by the baryonic mass component along ELO mass assembly, resulting in a tilt of the dynamical FP relative to the virial plane. The dynamical FPs shown by the two ELO samples are consistent with that shown by the SDSS elliptical sample in the same variables, with no further need for any relevant contribution from stellar population effects to explain the observed tilt. These effects could, however, have contributed to the scatter of the observed FP, as the dynamical FPs have been found to be thinner than the observed one. The results we report on hint, for the first time, at a possible way to understand the tilt of the observed FP in a cosmological context.

We present results of an analysis of the structural and kinematical properties of a sample of elliptical-like objects (ELOs) identified in four hydrodynamical, self-consistent simulations run with the DEVA code (Serna et al. 2003). Star formation has been implemented in the code through a simple phenomenological parameterization, that takes into account stellar physics processes only implicitly through the values of a threshold gas density, $\rho_{\rm g,thres}$, and an efficiency parameter, $c_*$. The four simulations operate in the context of a $\Lambda$CDM cosmological model consistent with observations, resolve ELO mass assembly at scales up to $\simeq$ 2 kpc, and differ in the values of their star formation parameters. Stellar masses, projected half-mass radii and central l.o.s. velocity dispersions, $\sigma_{\rm los, 0}$, have been measured on the ELO sample and their values compared with data from the Sloan digital sky survey. For the first time in self-consistent simulations, a good degree of agreement has been shown, including the Faber-Jackson and the $D_n - \sigma_{\rm los, 0}$ relations, among others, but only when particular values of the $\rho_{\rm g,thres}$ and $c_*$ parameters are used. This demostrates the effect that the star formation parameterization has on the ELO mass distribution. Additionally, our results suggest that it is not strictly necessary, at the scales resolved in this work, to appeal to energy sources other than gravitational (as for example supernovae feedback effects) to account for the structure and kinematics of large ellipticals.

In the framework of the fluctuation-dissipation approach to dynamical friction, we derive an expression giving the orbital energy exchange experienced by a compound body as it moves interacting with a non homogeneous discrete background. The body is assumed to be composed of particles endowed with a velocity spectrum and with a non homogeneous spatial distribution. The Chandrasekhar formula is recovered in the limit of a point-like satellite with zero velocity dispersion and infinite temperature moving through an homogeneous infinite medium. In this same limit, but dropping the zero satellite velocity dispersion ($\sigma_S$) condition, the orbital energy loss is found to be smaller than in the $\sigma_S=0$ case by a factor of up to an order of magnitude in some situations.