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.

We study the optical morphology of galaxies in a large-scale hydrodynamic cosmological simulation, the EAGLE simulation. Galaxy morphologies were characterized using non-parametric statistics (Gini, $M_{20}$, Concentration and Asymmetry) derived from mock images computed using a 3D radiative transfer technique and post-processed to approximate observational surveys. The resulting morphologies were contrasted to observational results from a sample of $\log_{10}(M_{*}/M_\odot) > 10$ galaxies at $z \sim 0.05$ in the GAMA survey. We find that the morphologies of EAGLE galaxies reproduce observations, except for asymmetry values which are larger in the simulated galaxies. Additionally, we study the effect of spatial resolution in the computation of non-parametric morphologies, finding that Gini and Asymmetry values are systematically reduced with decreasing spatial resolution. Gini values for lower mass galaxies are especially affected. Comparing against other large scale simulations, the non-parametric statistics of EAGLE galaxies largely agree with those found in IllustrisTNG. Additionally, EAGLE galaxies mostly reproduce observed trends between morphology and star formation rate and galaxy size. Finally, We also find a significant correlation between optical and kinematic estimators of morphologies, although galaxy classification based on an optical or a kinematic criteria results in different galaxy subsets. The correlation between optical and kinematic morphologies is stronger in central galaxies than in satellites, indicating differences in morphological evolution.

We identify a simulated Milky Way analog in the EAGLE suite of cosmological hydrodynamical simulations. This galaxy not only shares similar global properties as the Milky Way, but was specifically selected because its merger history resembles that currently known for the Milky Way. In particular we find that this Milky Way analog has experienced its last significant merger (with a stellar mass ratio $\sim 0.2$) at $z\sim 1.2$. We show that this merger affected both the dynamical properties of the stars present at the time, contributing to the formation of a thick disk, and also leading to a significant increase in the star formation rate of the host. This object is thus particularly suitable for understanding the early evolutionary history of the Milky Way. It is also an ideal candidate for re-simulation with much higher resolution as this would allow addressing a plethora of interesting questions such as, for example, the specific distribution of dark matter near the Sun.

We study the properties of long gamma-ray bursts (LGRBs) using a large scale hydrodynamical cosmological simulation, the Illustris simulation. We determine the LGRB host populations under different thresholds for the LGRB progenitor metallicities, according to the collapsar model. We compare the simulated sample of LGRBs hosts with recent, largely unbiased, host samples: BAT6 and SHOALS. We find that at $z<1$ simulated hosts follow the mass-metallicity relation and the fundamental metallicity relation simultaneously, but with a paucity of high-metallicity hosts, in accordance with observations. We also find a clear increment in the mean stellar mass of LGRB hosts and their SFR with redshift up to $z<3$ on account of the metallicity dependence of progenitors. We explore the possible origin of LGRBs in metal rich galaxies, and find that the intrinsic metallicity dispersion in galaxies could explain their presence. LGRB hosts present a tighter correlation between galaxy metallicity and internal metallicity dispersion compared to normal star forming galaxies. We find that the Illustris simulations favours the existence of a metallicity threshold for LGRB progenitors in the range 0.3 - 0.6 Z$_\odot$

Understanding the formation and evolution of early-type, spheroid-dominated galaxies is an open question within the context of the hierarchical clustering scenario, particularly, in low-density environments. Our goal is to study the main structural, dynamical, and stellar population properties and assembly histories of field spheroid-dominated galaxies formed in a LCDM scenario to assess to what extend they are consistent with observations. We selected spheroid-dominated systems from a LCDM simulation that includes star formation, chemical evolution and Supernova feedback. A sample of 18 field systems with Mstar <= 6x10^10 Msun that are dominated by the spheroid component. For this sample we estimate the fundamental relations of ellipticals and then compared with current observations. The simulated spheroid galaxies have sizes in good agreement with observations. The bulges follow a Sersic law with Sersic indexes that correlate with the bulge-to-total mass ratios. The structural-dynamical properties of the simulated galaxies are consistent with observed Faber-Jackson, Fundamental Plane, and Tully-Fisher relations. However, the simulated galaxies are bluer and with higher star formation rates than observed isolated early-type galaxies. The archaeological mass growth histories show a slightly delayed formation and more prominent inside-out growth mode than observational inferences based on the fossil record method. The main structural and dynamical properties of the simulated spheroid-dominated galaxies are consistent with observations. This is remarkable since none of them has been tuned to be reproduced. However, the simulated galaxies are blue and star-forming, and with later stellar mass growth histories as compared to observational inferences. This is mainly due to the persistence of extended discs in the simulations. Abridged

We study non-parametric morphologies of mergers events in a cosmological context, using the Illustris project. We produce mock g-band images comparable to observational surveys from the publicly available Illustris simulation idealized mock images at $z=0$. We then measure non parametric indicators: asymmetry, Gini, $M_{20}$, clumpiness and concentration for a set of galaxies with $M_* >10^{10}$ M$_\odot$. We correlate these automatic statistics with the recent merger history of galaxies and with the presence of close companions. Our main contribution is to assess in a cosmological framework, the empirically derived non-parametric demarcation line and average time-scales used to determine the merger rate observationally. We found that 98 per cent of galaxies above the demarcation line have a close companion or have experienced a recent merger event. On average, merger signatures obtained from the $G-M_{20}$ criteria anticorrelate clearly with the elapsing time to the last merger event. We also find that the asymmetry correlates with galaxy pair separation and relative velocity, exhibiting the larger enhancements for those systems with pair separations $d < 50$ h$^{-1}$ kpc and relative velocities $V < 350$ km s$^{-1}$. We find that the $G-M_{20}$ is most sensitive to recent mergers ($\sim0.14$ Gyr) and to ongoing mergers with stellar mass ratios greater than 0.1. For this indicator, we compute a merger average observability time-scale of $\sim0.2$ Gyr, in agreement with previous results and demonstrate that the morphologically derived merger rate recovers the intrinsic total merger rate of the simulation and the merger rate as a function of stellar mass.

We study the nature of long gamma ray burst (LGRB) progenitors using cosmological simulations of structure formation and galactic evolution. LGRBs are potentially excellent tracers of stellar evolution in the early universe. We developed a Monte Carlo numerical code which generates LGRBs coupled to cosmological simulations. The simulations allows us to follow the ormation of galaxies self-consistently. We model the detectability of LGRBs and their host galaxies in order to compare results with observational data obtained by high-energy satellites. Our code also includes stochastic effects in the observed rate of LGRBs.

X-ray observations unveiled various types of radio-silent Isolated Neutron Stars (INSs), phenomenologically very diverse, e.g. the Myr old X-ray Dim INS (XDINSs) and the kyr old magnetars. Although their phenomenology is much diverse, the similar periods (P=2--10 s) and magnetic fields (~10^{14} G) suggest that XDINSs are evolved magnetars, possibly born from similar populations of supermassive stars. One way to test this hypothesis is to identify their parental star clusters by extrapolating backward the neutron star velocity vector in the Galactic potential. By using the information on the age and space velocity of the XDINS RX J1856.5-3754, we computed backwards its orbit in the Galactic potential and searched for its parental stellar cluster by means of a closest approach criterion. We found a very likely association with the Upper Scorpius OB association, for a neutron star age of 0.42+/-0.08 Myr, a radial velocity V_r^NS =67+/- 13$ km s^{-1}, and a present-time parallactic distance d_\pi^NS = 123^{+11}_{-15} pc. Our result confirms that the "true" neutron star age is much lower than the spin-down age (tau_{sd}=3.8 Myrs), and is in good agreement with the cooling age, as computed within standard cooling scenarios. The mismatch between the spin-down and the dynamical/cooling age would require either an anomalously large breaking index (n~20) or a decaying magnetic field with initial value B_0 ~ 10^{14} G. Unfortunately, owing to the uncertainty on the age of the Upper Scorpius OB association and the masses of its members we cannot yet draw firm conclusions on the estimated mass of the RX J1856.5-3754 progenitor.