We have developed a new tool to analyse galaxies in the EAGLE simulations as close as possible to observations.We investigated the evolution of their kinematic properties by means of the angular momentum proxy parameter,$ \lambda_{Re} $for galaxies with $M_{*} \ge 5 \times 10^{9} M_{\odot}$ in the RefL0100N1504 simulation up to redshift two (z = 2). Galaxies in the simulation show a wide variety of kinematic features, similiar to those found in integral-field spectroscopic studies. At z=0 the distribution of galaxies in the ${\lambda}_{Re}-{\epsilon}$ plane is also in good agreement with results from observations. Scaling relations at z = 0 indicate that there is critical mass, $M_{crit} = 10^{10.3} M_{\odot}$, that divides two different regimes when we include the ${\lambda}_{Re}$ parameter. The simulation shows that the distribution of galaxies in the ${\lambda}_{Re}-{\epsilon}$ plane evolves with time until z = 2 when galaxies are equally distributed both in ${\lambda}_{Re}$ and ${\epsilon}$. We studied the evolution of ${\lambda}_{Re}$ with time and found that there is no connection between the angular momentum at z = 2 and z = 0. All systems reach their maximum ${\lambda}_{Re}$ at z = 1 and then steadily lose angular momentum regardless of their merger history, except for the high star-forming systems that sustain that maximum value over time. The evolution of the Re in galaxies that have not experienced any merger in the last 10 Gyr can be explained by their level of gas accretion.

Early-type galaxies are considered to be the end-products of massive galaxy formation. Optical spectroscopic studies reveal that massive early-type galaxies formed the bulk of their stars over short timescales (1<Gyr) and at high redshift (z>2), followed by passive evolution to the present. However, their optical spectra are insensitive to constrain small episodes of recent star formation, since they are dominated by old stars. Fortunately, this problem can be tackled in the ultraviolet range. While recent studies that make use of ultraviolet absorption lines have suggested the presence of young stars in a few early-type galaxies, the age and mass fractions of young stars and their dependence on galaxy mass, is unknown. Here we report a detailed study of these young stellar populations, from high-quality stacked spectra of 28,663 galaxies from the BOSS survey, analysing optical and ultraviolet absorption lines simultaneously. We find that residual star formation is ubiquitous in massive early-type galaxies, measuring average mass fractions of ~0.5% in young stars in the last 2 Gyr of their evolution. This fraction shows a decreasing trend with galaxy stellar mass, consistent with a down-sizing scenario. We also find that synthetic galaxies from state-of-the-art cosmological numerical simulations significantly overproduce both intermediate and young stellar populations. Therefore, our results pose stringent constraints on numerical simulations of galaxy formation.

A recent determination of the relationships between the X-ray luminosity of the ISM (Lx) and the stellar and total mass, for a sample of nearby early-type galaxies (ETGs), is used to investigate the origin of the hot gas, via a comparison with the results of hydrodynamical simulations of the ISM evolution for a large set of isolated ETGs. After the epoch of major galaxy formation (after z~2), the ISM is replenished by stellar mass losses and SN ejecta, at the rate predicted by stellar evolution, and is depleted by star formation; it is heated by the thermalization of stellar motions, SNe explosions and the mechanical (from winds) and radiative AGN feedback. The models agree well with the observed relations, even for the largely different Lx values at the same mass, thanks to the sensitivity of the gas flow to many galaxy properties; this holds for models including AGN feedback, and those without. Therefore, the mass input from the stellar population is able to account for a major part of the observed Lx; and AGN feedback, while very important to maintain massive ETGs in a time-averaged quasi-steady state, keeping low star formation and the black hole mass, does not dramatically alter the gas content originating in stellar recycled material. These conclusions are based on theoretical predictions for the stellar population contributions in mass and energy, and on a self-consistent modeling of AGN feedback.

Black holes (BHs) are believed to be a key ingredient of galaxy formation. However, the galaxy-BH interplay is challenging to study due to the large dynamical range and complex physics involved. As a consequence, hydrodynamical cosmological simulations normally adopt sub-grid models to track the unresolved physical processes, in particular BH accretion; usually the spatial scale where the BH dominates the hydrodynamical processes (the Bondi radius) is unresolved, and an approximate Bondi-Hoyle accretion rate is used to estimate the growth of the BH. By comparing hydrodynamical simulations at different resolutions (300, 30, 3 pc) using a Bondi-Hoyle approximation to sub-parsec runs with non-parameterized accretion, our aim is to probe how well an approximated Bondi accretion is able to capture the BH accretion physics and the subsequent feedback on the galaxy. We analyse an isolated galaxy simulation that includes cooling, star formation, Type Ia and Type II supernovae, BH accretion and AGN feedback (radiation pressure, Compton heating/cooling) where mass, momentum, and energy are deposited in the interstellar medium through conical winds. We find that on average the approximated Bondi formalism can lead to both over- and under-estimations of the BH growth, depending on resolution and on how the variables entering into the Bondi-Hoyle formalism are calculated.

We present 2D hydrodynamical simulations for the evolution of early-type galaxies containing central massive black holes (MBHs), starting at age 2 Gyr. The code contains accurate and physically consistent radiative and mechanical AGN wind feedback, with parsec-scale central resolution. Mass input comes from stellar evolution; energy input includes Type Ia and II supernova and stellar heating; star-formation is included. Realistic, axisymmetric dynamical models for the galaxies are built solving the Jeans' equations. The lowest mass models (Mstar = 8 10^{10}Msun) develop global outflows sustained by SNIa's heating, ending with a significantly lower amount of hot gas and new stars. In more massive models, nuclear outbursts last to the present epoch, with large and frequent fluctuations in nuclear emission and from the gas (Lx). Each burst last ~ 10^{7.5} yr, during which (for r < 2-3 kpc) cold, inflowing, and hot, outflowing gas phases coexist. The Lx-T relation for the gas matches that of local galaxies. AGN activity causes positive feedback for star formation. Roughly half of the total mass loss is recycled into new stars (DeltaMstar), just ~ 3% of it is accreted on the MBH, the remainder being ejected from the galaxy. The ratio between the mass of gas expelled to that in to new stars, the load factor, is ~0.6. Rounder galaxies shapes lead to larger final MBH masses, DeltaMstar, and Lx. Almost all the time is spent at very low nuclear luminosities, yet one quarter of the total energy is emitted at an Eddington ratio > 0.1. The duty-cycle of AGN activity approximates 4% (Abridged).

AGN feedback from supermassive black holes (SMBHs) at the center of early type galaxies is commonly invoked as the explanation for the quenching of star formation in these systems. The situation is complicated by the significant amount of mass injected in the galaxy by the evolving stellar population over cosmological times. In absence of feedback, this mass would lead to unobserved galactic cooling flows, and to SMBHs two orders of magnitude more massive than observed. By using high-resolution 2D hydrodynamical simulations with radiative transport and star formation in state-of-the-art galaxy models, we show how the intermittent AGN feedback is highly structured on spatial and temporal scales, and how its effects are not only negative (shutting down the recurrent cooling episodes of the ISM), but also positive, inducing star formation in the inner regions of the host galaxy.

Early-type galaxies (ETGs) host a hot ISM produced mainly by stellar winds, and heated by Type Ia supernovae and the thermalization of stellar motions. High resolution 2D hydrodynamical simulations showed that ordered rotation in the stellar component results in the formation of a centrifugally supported cold equatorial disc. In a recent numerical investigation we found that subsequent generations of stars are formed in this cold disc; this process consumes most of the cold gas, leaving at the present epoch cold masses comparable to those observed. Most of the new stellar mass formed a few Gyrs ago, and resides in a disc.

High resolution 2D hydrodynamical simulations describing the evolution of the hot ISM in axisymmetric two-component models of early-type galaxies well reproduced the observed trends of the X-ray luminosity ($L_\mathrm{x}$) and temperature ($T_\mathrm{x}$) with galaxy shape and rotation, however they also revealed the formation of an exceedingly massive cooled gas disc in rotating systems. In a follow-up of this study, here we investigate the effects of star formation in the disc, including the consequent injection of mass, momentum and energy in the pre-existing interstellar medium. It is found that subsequent generations of stars originate one after the other in the equatorial region; the mean age of the new stars is $> 5$ Gyr, and the adopted recipe for star formation can reproduce the empirical Kennicutt-Schmidt relation. The results of the previous investigation without star formation, concerning $L_\mathrm{x}$ and $T_\mathrm{x}$ of the hot gas, and their trends with galactic shape and rotation, are confirmed. At the same time, the consumption of most of the cold gas disc into new stars leads to more realistic final systems, whose cold gas mass and star formation rate agree well with those observed in the local universe. In particular, our models could explain the observation of kinematically aligned gas in massive, fast-rotating early-type galaxies.

By means of high resolution 2D hydrodynamical simulations, we study the evolution of the hot ISM for a large set of early-type galaxy models, characterized by various degrees of flattening and internal rotation. The galaxies are described by state-of-the-art axisymmetric two-component models, tailored to reproduce real systems; the dark matter haloes follow the Navarro-Frenk-White or the Einasto profile. The gas is produced by the evolving stars, and heated by Type Ia SNe. We find that, in general, the rotation field of the ISM in rotating galaxies is very similar to that of the stars, with a consequent negligible heating contribution from thermalization of the ordered motions. The relative importance of flattening and rotation in determining the final X-ray luminosity $L_x$ and temperature $T_x$ of the hot haloes is a function of the galactic mass. Flattening and rotation in low mass galaxies favour the establishment of global winds, with the consequent reduction of $L_x$. In medium-to-high mass galaxies, flattening and rotation are not sufficient to induce global winds, however, in the rotating models the nature of the gas flows is deeply affected by conservation of angular momentum, resulting in a reduction of both $L_x$ and $T_x$.

Observational and numerical studies gave hints that the hot gaseous haloes of ETGs may be sensitive to the galaxy internal kinematics. By using high resolution 2D hydro simulations, and realistic two-component (stars plus dark matter) axisymmetric galaxy models, we study the evolution of the hot haloes in a suite of flat ETGs of fixed mass distribution, but with variable amounts of azimuthal velocity dispersion and rotational support, including the possibility of a counter-rotating inner stellar disc. The hot halo is fed by stellar mass losses and heated by SNIa explosions and thermalization of stellar motions. We measure the value of the thermalization parameter gamma (the ratio between the heating due to the relative velocity between the stellar streaming and the ISM bulk flow, and the heating attainable by complete thermalization of the stellar streaming motions). We find that 1) the X-ray emission and the average temperature are larger in fully velocity dispersion supported systems; 2) 0.1<gamma<0.2 for isotropic rotators (with a trend for being larger for lower dark mass models); 3) systems that are isotropic rotators at large radii with an inner counter-rotating disc, or fully velocity dispersion supported systems with an inner rotating disc, have gamma=1, again with a trend to increase for lower dark mass contents. We also find that the lower X-ray luminosities of isotropic rotators cannot be explained just by their low gamma, but are due to the complicated flow structure, consequence of the angular momentum stored at large radii. X-ray emission weighted temperatures and luminosities nicely match observed values; the X-ray isophotes are boxy in case of significant galaxy rotation. Overall, it is found that rotation has an important role to explain the observational result that more rotationally supported ETGs on average show a lower X-ray emission [abridged].

With 2D hydrodynamical simulations, we study the evolution of the hot gas flows in early-type galaxies, focussing on the effects of galaxy rotation on the thermal and dynamical status of the ISM. The galaxy is modelled as a two-component axisymmetric system (stars and dark matter), with a variable amount of azimuthal velocity dispersion and rotational support; the presence of a counter rotating stellar disk is also considered. It is found that the ISM of the rotationally supported (isotropic) model is more prone to thermal instabilities than the fully velocity dispersion counterpart, while its ISM temperature and X-ray luminosity are lower. The model with counter rotation shows an intermediate behaviour.