We perform a set of non-radiative hydro-dynamical (NHD) simulations of a rich cluster sized dark matter halo from the Phoenix project with 3 different numerical resolutions, to investigate the effect of hydrodynamics alone on the subhalo population in the halo. Compared to dark matter only (DMO) simulations of the same halo, subhaloes are less abundant for relatively massive subhaloes ($M_{sub} > 2.5 \times 10^9h^{-1}M_{\odot}$, or $V_{max} > 70 kms^{-1}$) but more abundant for less massive subhaloes in the NHD simulations. This results in different shapes in the subhalo mass/$V_{max}$ function in two different sets of simulations. At given subhalo mass, the subhaloes less massive than $10^{10} h^{-1}M_{\odot}$ have larger $V_{max}$ in the NHD than DMO simulations, while $V_{max}$ is similar for the subhaloes more massive than the mass value. This is mainly because the progenitors of present day low mass subhaloes have larger concentration parameters in the NHD than DMO simulations. The survival number fraction of the accreted low mass progenitors of the main halo at redshift 2 is about 50 percent higher in the NHD than DMO simulations.

Using a significantly enlarged \textit{Fermi}-LAT BL Lac objects (BL Lacs) sample, we construct the gamma-ray luminosity function (GLF) of BL Lacs, by the joint use of the space density distribution and source counts distribution. We use three well-studied forms of the GLF, i.e., the forms of pure density evolution (PDE), pure luminosity evolution (PLE) and luminosity-dependent density evolution (LDDE). The Markov Chain Monte Carlo (MCMC) technique is used to constrain model parameters. Our results suggest that LDDE model can give the best description for the BL Lac GLF. And the model shows that the BL Lacs with a harder GeV spectrum and a less luminosity evolve as strongly as FSRQs, and the evolution decreases as increasing luminosity. We also model the average photon spectra of BL Lacs with a double power laws model. Using this modeled spectra, BL Lacs contribute $ \sim 20\%$ of the total extragalactic gamma-ray background (EGB) at $E>100$ MeV, $\sim 100\%$ of the EGB at $E>50$ GeV, and the unresolved BL Lacs contribute $\sim 20 \%$ of the isotropic diffuse gamma-ray background (IGRB) at $E>100$ MeV. A prediction of the TeV EGB spectra are given, which may be tested by the future detectors.

We compare galaxy morphology predicted by the Illustris hydrodynamical simulation and a semi-analytic model (SAM) grafted in the halo merger trees from the Illustris-Dark matter simulation. Morphology is classified according to the luminous profile and the kinematic bulge-to-total ratio for Illustris galaxy, and the bulge-to-total stellar mass ratio for SAM galaxy. For late-type galaxies in the Illustris catalogue, most of their counterparts in the SAM model have the same type, and the consistency between two models is higher for lower mass galaxies. For early-type (ET) galaxies in Illustris, the consistency is quite low for the counterparts except for most massive galaxies. By comparing in detail the growth histories of some matched galaxy pairs of Milky Way mass in Illustris and the SAM model, we notice two aspects of differences in determining galaxy morphology between the two galaxy formation implementations. First, in the SAM, major merger and frequent minor mergers result in the growth of bulges and turn the galaxy into ET, while bulge formation is not connected to mergers as tightly as in SAM for the Illustris galaxies. In addition, the satellite stellar mass can decrease significantly due to tidal stripping before merging into the central galaxy in Illustris, while it does not decrease in the SAM model. This results in less mergers with large (stellar) mass ratios in the Illustris simulation, and less effect of mergers on shaping galaxy morphology.

We analyse the mass assembly of central galaxies in the EAGLE hydrodynamical simulations. We build merger trees to connect galaxies to their progenitors at different redshifts and characterize their assembly histories by focusing on the time when half of the galaxy stellar mass was assembled into the main progenitor. We show that galaxies with stellar mass $M_*<10^{10.5}M_{\odot}$ assemble most of their stellar mass through star formation in the main progenitor (`in-situ' star formation). This can be understood as a consequence of the steep rise in star formation efficiency with halo mass for these galaxies. For more massive galaxies, however, an increasing fraction of their stellar mass is formed outside the main progenitor and subsequently accreted. Consequently, while for low-mass galaxies the assembly time is close to the stellar formation time, the stars in high-mass galaxies typically formed long before half of the present-day stellar mass was assembled into a single object, giving rise to the observed anti-hierarchical downsizing trend. In a typical present-day $M_*\geq10^{11}M_{\odot}$ galaxy, around $20\%$ of the stellar mass has an external origin. This fraction decreases with increasing redshift. Bearing in mind that mergers only make an important contribution to the stellar mass growth of massive galaxies, we find that the dominant contribution comes from mergers with galaxies of mass greater than one tenth of the main progenitor's mass. The galaxy merger fraction derived from our simulations agrees with recent observational estimates.

We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged)

We present the evolution of galaxy sizes, from redshift 2 to 0, for actively star forming and passive galaxies in the cosmological hydrodynamical 1003 cMpc3 simulation of the EAGLE project. We find that the sizes increase with stellar mass , but that the relation weakens with increasing redshift. Separating galaxies by their star formation activity, we find that passive galaxies are typically smaller than active galaxies at fixed stellar mass. These trends are consistent with those found in observations and the level of agreement between the predicted and observed size - mass relation is of order 0.1 dex for z < 1 and 0.2-0.3 dex from redshift 1 to 2. We use the simulation to compare the evolution of individual galaxies to that of the population as a whole. While the evolution of the size-stellar mass relation for active galaxies provides a good proxy for the evolution of individual galaxies, the evolution of individual passive galaxies is not well represented by the observed size - mass relation due to the evolving number density of passive galaxies. Observations of z \approx 2 galaxies have revealed an abundance of massive red compact galaxies, that depletes below z \approx 1. We find that a similar population forms naturally in the simulation. Comparing these galaxies to their z = 0 descendants, we find that all compact galaxies grow in size due to the high-redshift stars migrating outwards. Approximately 60% of the compact galaxies increase in size further due to renewed star formation and/or mergers.

We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.

The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation.

By means of N-body simulations we investigate the impact of minor mergers on the angular momentum and dynamical properties of the merger remnant. Our simulations cover a range of initial orbital characteristics and gas-to-stellar mass fractions (from 0 to 20%), and include star formation and supernova feedback. We confirm and extend previous results by showing that the specific angular momentum of the stellar component always decreases independently of the orbital parameters or morphology of the satellite, and that the decrease in the rotation velocity of the primary galaxy is accompanied by a change in the anisotropy of the orbits. However, the decrease affects only the old stellar population, and not the new population formed from gas during the merging process. This means that the merging process induces an increasing difference in the rotational support of the old and young stellar components, with the old one lagging with respect to the new. Even if our models are not intended specifically to reproduce the Milky Way and its accretion history, we find that, under certain conditions, the modeled rotational lag found is compatible with that observed in the Milky Way disk, thus indicating that minor mergers can be a viable way to produce it. The lag can increase with the vertical distance from the disk midplane, but only if the satellite is accreted along a direct orbit, and in all cases the main contribution to the lag comes from stars originally in the primary disk rather than from stars in the satellite galaxy. We also discuss the possibility of creating counter-rotating stars in the remnant disk, their fraction as a function of the vertical distance from the galaxy midplane, and the cumulative effect of multiple mergers on their creation.

By means of N-body/SPH simulations we investigate the morphological properties of thick disks formed through minor mergers. We show that the vertical surface density profile of the post-merger thick disk follows a sech function and has an excess in the regions far from the disk mid-plane (z>2kpc). This stellar excess also follows a sech function with a larger scale height than the main thick disk component, and it is usually dominated by stars from the primary galaxy. Stars in the excess have a rotational velocity lower than that of stars in the thick disk, and they may thus be confused with stars in the inner galactic halo. The thick disk scale height increases with radius and the rate of its increase is smaller for more gas rich primary galaxies. On the contrary, the scale height of the stellar excess is independent of both radius and gas fraction. We also find that the post-merger thick disk has a radial scale length which is 10-50% larger than that of the thin disk. Two consecutive mergers have basically the same effect on heating the stellar disk as a single merger of the same total mass. To investigate how thick disks produced through secular processes may differ from those produced by minor mergers, we also simulated gravitationally unstable gas-rich disks. These disks do not produce either a stellar excess or a ratio of thick to thin disk scale lengths greater than one. Our simulation results are consistent with observations of the ratio of thick to thin disk scale lengths of the Milky Way and nearby galaxies, and with the Toomre diagram of the Milky Way. We conclude that minor mergers are a viable mechanism for the creation of galactic thick disks and investigating stars at several kpc above the mid-plane of the Milky Way and other galaxies may provide a quantitative method for studying the minor merger history of galaxies.

We study the distribution of orbital eccentricities of stars in thick disks generated by the heating of a pre-existing thin stellar disk through a minor merger (mass ratio 1:10), using N-body/SPH numerical simulations of interactions that span a range of gas fractions in the primary disk and initial orbital configurations. The resulting eccentricity distributions have an approximately triangular shape, with a peak at 0.2-0.35, and a relatively smooth decline towards higher values. Stars originally in the satellite galaxy tend to have higher eccentricities (on average from e = 0.45 to e = 0.75), which is in general agreement with the models of Sales and collaborators, although in detail we find fewer stars with extreme values and no evidence of their secondary peak around e = 0.8. The absence of this high-eccentricity feature results in a distribution that qualitatively matches the observations. Moreover, the increase in the orbital eccentricities of stars in the solar neighborhood with vertical distance from the Galactic mid-plane recently found by Diericxk and collaborators can be qualitatively reproduced by our models, but only if the satellite is accreted onto a direct orbit. We thus speculate that if minor mergers were the dominant means of formating the Milky Way thick disk, the primary mechanism should be merging with satellite(s) on direct orbits.

We have investigated the impact of dissipationless minor galaxy mergers on the angular momentum of the remnant. Our simulations cover a range of initial orbital characteristics and the system consists of a massive galaxy with a bulge and disk merging with a much less massive (one-tenth or one-twentieth) gasless companion which has a variety of morphologies (disk- or elliptical-like) and central baryonic mass concentrations. During the process of merging, the orbital angular momentum is redistributed into the internal angular momentum of the final system; the internal angular momentum of the primary galaxy can increase or decrease depending on the relative orientation of the orbital spin vectors (direct or retrograde), while the initially non-rotating dark matter halo always gains angular momentum. The specific angular momentum of the stellar component always decreases independent of the orbital parameters or morphology of the satellite, the decrease in the rotation velocity of the primary galaxy is accompanied by a change in the anisotropy of the orbits, and the ratio of rotation speed to velocity dispersion of the merger remnant is lower than the initial value, not only due to an increase in the dispersion but also to the slowing -down of the disk rotation. We briefly discuss several astrophysical implications of these results, suggesting that minor mergers do not cause a "random walk" process of the angular momentum of the stellar disk component of galaxies, but rather a steady decrease. Minor mergers may play a role in producing the large scatter observed in the Tully-Fisher relation for S0 galaxies, as well as in the increase of the velocity dispersion and the decrease in $v/\sigma$ at large radii as observed in S0 galaxies.

Using a purely analytic approach to gaseous and dark matter halos, we study the cross-correlation between the Sunyaev-Zel'dovich (SZ) sky and galaxy survey under flat sky approximation, in an attempt to acquire the redshift information of the SZ map. The problem can be greatly simplified when it is noticed that the signals of the SZ-galaxy correlation arise only from hot gas and galaxies inside the same massive halos (i.e. clusters), and field galaxies make almost no contribution to the cross-correlation. Under the assumption that both the hot gas and galaxies trace the common gravitational potential of dark halos, we calculate the expected cross SZ-galaxy power spectra for the WMAP/Planck SZ maps and the SDSS galaxy sample at small scales $100<l<1000$. It turns out, however, that it is not presently feasible to measure such small angular cross power spectra because of the high noise levels at $l>400$ with the WMAP/Planck experiments. Future SZ observations with better angular resolutions and sufficiently wide sky coverages will be needed if this technique is applied for the statistical measurement of redshift distribution of the SZ sources.

We explore the application of artificial neural networks (ANNs) for the estimation of atmospheric parameters (Teff, logg, and [Fe/H]) for Galactic F- and G-type stars. The ANNs are fed with medium-resolution (~ 1-2 A) non flux-calibrated spectroscopic observations. From a sample of 279 stars with previous high-resolution determinations of metallicity, and a set of (external) estimates of temperature and surface gravity, our ANNs are able to predict Teff with an accuracy of ~ 135-150 K over the range 4250 <= Teff <= 6500 K, logg with an accuracy of ~ 0.25-0.30 dex over the range 1.0 <= logg <= 5.0 dex, and [Fe/H] with an accuracy ~ 0.15-0.20 dex over the range -4.0 <= [Fe/H] <= +0.3. Such accuracies are competitive with the results obtained by fine analysis of high-resolution spectra. It is noteworthy that the ANNs are able to obtain these results without consideration of photometric information for these stars. We have also explored the impact of the signal-to-noise ratio (S/N) on the behavior of ANNs, and conclude that, when analyzed with ANNs trained on spectra of commensurate S/N, it is possible to extract physical parameter estimates of similar accuracy with stellar spectra having S/N as low as 13. Taken together, these results indicate that the ANN approach should be of primary importance for use in present and future large-scale spectroscopic surveys.

New generation large-aperture telescopes, multi-object spectrographs, and large format detectors are making it possible to acquire very large samples of stellar spectra rapidly. In this context, traditional star-by-star spectroscopic analysis are no longer practical. New tools are required that are capable of extracting quickly and with reasonable accuracy important basic stellar parameters coded in the spectra. Recent analyses of Artificial Neural Networks (ANNs) applied to the classification of astronomical spectra have demonstrated the ability of this concept to derive estimates of temperature and luminosity. We have adapted the back-propagation ANN technique developed by von Hippel et al. (1994) to predict effective temperatures, gravities and overall metallicities from spectra with resolving power ~ 2000 and low signal-to-noise ratio. We show that ANN techniques are very effective in executing a three-parameter (Teff,log g,[Fe/H]) stellar classification. The preliminary results show that the technique is even capable of identifying outliers from the training sample.