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.

The spatial and temporal relationships between stellar age, kinematics, and chemistry are a fundamental tool for uncovering the physics driving galaxy formation and evolution. Observationally, these trends are derived using carefully selected samples isolated via the application of appropriate magnitude, colour, and gravity selection functions of individual stars; conversely, the analysis of chemodynamical simulations of galaxies has traditionally been restricted to the age, metallicity, and kinematics of `composite' stellar particles comprised of open cluster-mass simple stellar populations. As we enter the Gaia era, it is crucial that this approach changes, with simulations confronting data in a manner which better mimics the methodology employed by observers. Here, we use the \textsc{SynCMD} synthetic stellar populations tool to analyse the metallicity distribution function of a Milky Way-like simulated galaxy, employing an apparent magnitude plus gravity selection function similar to that employed by the RAdial Velocity Experiment (RAVE); we compare such an observationally-motivated approach with that traditionally adopted - i.e., spatial cuts alone - in order to illustrate the point that how one analyses a simulation can be, in some cases, just as important as the underlying sub-grid physics employed.

Using a cosmological hydrodynamical simulation of a galaxy of similar mass to the Large Magellanic Cloud (LMC), we examine the predicted characteristics of its lowest metallicity populations. In particular, we emphasise the spatial distributions of first (Pop III) and second (polluted by only immediate Pop III ancestors) generation stars. We find that primordial composition stars form not only in the central galaxy's progenitor, but also in locally collapsed sub-halos during the early phases of galaxy formation. The lowest metallicity stars in these sub-halos end up in a relatively extended distribution around the host, with these accreted stars possessing present-day galactocentric distances as great as ~40kpc. By contrast, the earliest stars formed within the central galaxy remain in the inner region, where the vast majority of star formation occurs, for the entirety of the simulation. Consequently, the fraction of stars that are from the earliest generation increases strongly with radius.

Within a cosmological hydrodynamical simulation, we form a disc galaxy with sub- components which can be assigned to a thin stellar disc, thick disk, and a low mass stellar halo via a chemical decomposition. The thin and thick disc populations so selected are distinct in their ages, kinematics, and metallicities. Thin disc stars are young (<6.6 Gyr), possess low velocity dispersion ({\sigma}U,V,W = 41, 31, 25 km/s), high [Fe/H], and low [O/Fe]. The thick disc stars are old (6.6<age<9.8 Gyrs), lag the thin disc by \sim21 km/s, possess higher velocity dispersion ({\sigma}U,V,W = 49, 44, 35 km/s), relatively low [Fe/H] and high [O/Fe]. The halo component comprises less than 4% of stars in the "solar annulus" of the simulation, has low metallicity, a velocity ellipsoid defined by ({\sigma}U,V,W = 62, 46, 45 km/s) and is formed primarily in-situ during an early merger epoch. Gas-rich mergers during this epoch play a major role in fuelling the formation of the old disc stars (the thick disc). This is consistent with studies which show that cold accretion is the main source of a disc galaxy's baryons. Our simulation initially forms a relatively short (scalelength \sim1.7 kpc at z=1) and kinematically hot disc, primarily from gas accreted during the galaxy's merger epoch. Far from being a competing formation scenario, migration is crucial for reconciling the short, hot, discs which form at high redshift in {\Lambda}CDM, with the properties of the thick disc at z=0. The thick disc, as defined by its abundances maintains its relatively short scale-length at z = 0 (2.31 kpc) compared with the total disc scale-length of 2.73 kpc. The inside-out nature of disc growth is imprinted the evolution of abundances such that the metal poor {\alpha}-young population has a larger scale-length (4.07 kpc) than the more chemically evolved metal rich {\alpha}-young population (2.74 kpc).

We perform a combined X-ray and strong lensing analysis of RX J1347.5-1145, one of the most luminous galaxy clusters at X-ray wavelengths. We show that evidence from strong lensing alone, based on published VLT and new HST data, strongly argues in favor of a complex structure. The analysis takes into account arc positions, shapes and orientations and is done thoroughly in the image plane. The cluster inner regions are well fitted by a bimodal mass distribution, with a total projected mass of $M_{tot} = (9.9 \pm 0.3)\times 10^{14} M_\odot/h$ within a radius of $360 \mathrm{kpc}/h$ ($1.5'$). Such a complex structure could be a signature of a recent major merger as further supported by X-ray data. A temperature map of the cluster, based on deep Chandra observations, reveals a hot front located between the first main component and an X-ray emitting South Eastern sub-clump. The map also unveils a filament of cold gas in the innermost regions of the cluster, most probably a cooling wake caused by the motion of the cD inside the cool core region. A merger scenario in the plane of the sky between two dark matter sub-clumps is consistent with both our lensing and X-ray analyses, and can explain previous discrepancies with mass estimates based on the virial theorem.

We study the anomalous flux ratio which is observed in some four-image lens systems, where the source lies close to a fold caustic. In this case two of the images are close to the critical curve and their flux ratio should be equal to unity, instead in several cases the observed value differs significantly. The most plausible solution is to invoke the presence of substructures, as for instance predicted by the Cold Dark Matter scenario, located near the two images. In particular, we analyze the two fold lens systems PG1115+080 and B1555+375, for which there are not yet satisfactory models which explain the observed anomalous flux ratios. We add to a smooth lens model, which reproduces well the positions of the images but not the anomalous fluxes, one or two substructures described as singular isothermal spheres. For PG1115+080 we consider a smooth model with the influence of the group of galaxies described by a SIS and a substructure with mass $\sim 10^{5} M_{\odot}$ as well as a smooth model with an external shear and one substructure with mass $\sim 10^{8} M_{\odot}$ . For B1555+375 either a strong external shear or two substructures with mass $\sim 10^{7} M_{\odot}$ reproduce the data quite well.

Warm Dark Matter (WDM) has been invoked to resolve apparent conflicts of Cold Dark Matter (CDM) models with observations on subgalactic scales. In this work we provide a new and independent lower limit for the WDM particle mass (e.g. sterile neutrino) through the analysis of image fluxes in gravitationally lensed QSOs. Starting from a theoretical unperturbed cusp configuration we analyze the effects of intergalactic haloes in modifying the fluxes of QSO multiple images, giving rise to the so-called anomalous flux ratio. We found that the global effect of such haloes strongly depends on their mass/abundance ratio and it is maximized for haloes in the mass range $10^6-10^8 \Msun$. This result opens up a new possibility to constrain CDM predictions on small scales and test different warm candidates, since free streaming of warm dark matter particles can considerably dampen the matter power spectrum in this mass range. As a consequence, while a ($\Lambda$)CDM model is able to produce flux anomalies at a level similar to those observed, a WDM model, with an insufficiently massive particle, fails to reproduce the observational evidences. Our analysis suggests a lower limit of a few keV ($m_{\nu} \sim 10$) for the mass of warm dark matter candidates in the form of a sterile neutrino. This result makes sterile neutrino Warm Dark Matter less attractive as an alternative to Cold Dark Matter, in good agreement with previous findings from Lyman-$\alpha$ forest and Cosmic Microwave Background analysis.

Aims. Falomo et al. (2005) discovered a disk-like galaxy at ~ 1.2 arcsec from the QSO Q0045-3337 by means of ESO VLT adaptive optics. They estimated a galaxy Einstein radius (for a point mass) of comparable size, thus pointing up the existence of a new, rare, spiral lens candidate, despite no evident image splitting. Here we discuss the possible lensing effect of the galaxy in some more detail. Methods. We performed two dimensional surface photometry on the VLT image of the galaxy, confirming its spiral nature. We then verified if simple mass models, partially constrained by observational data, require unrealistic parameters to produce a still hidden second quasar image. We also evaluated the respective viability of an instrumental or a lensing origin of the observed QSO deformation. Results. After galaxy model subtraction, we found a residual image, likely not related to gravitational lensing. Existing data are not sufficient to assess the presence of image splitting, nor to constrain the number of images or discriminate between galaxy mass models. Conclusion. Further observations are mandatory to progress in the study of this remarkable system, that could shed more light on the lensing behaviour of spiral galaxies.

It has been shown that the dissipative gas infall during galaxy formation has the capability to modify the shape of dark matter halos. In this paper we perform the first detailed analysis of particle orbits in a cosmological dark matter halos to understand {\it how} and {\it why} baryons alter its shape. We perform a series of numerical experiments where we grow a baryonic core inside a live dark matter halo extracted from a cosmological simulation. We follow the orbits of more than 300 particles with more that 50000 timesteps. Our results clearly show that the dissipational component is responsible for repeatedly deflecting orbits which visit often the center of the system. Moreover the gravitational potential time dependence associated with the growth of the baryonic mass, shifts the frequencies of the orbits, making them extremely chaotic. This randomization makes the orbits explore a large phase space. When this effect takes place for a significant number of orbits it will be manifested in the density distribution as an approach to a rounder configuration. As a consequence, the influence of the central mass on the shape of the phase space decreases with increasing distance from the center. We discuss the importance of future analysis of controlled experiments (i.e. using analytic potentials instead of live DM halos) to better decipher the dynamics of this phenomenon.

It has been argued that the flux anomalies detected in gravitationally lensed QSOs are evidence for substructures in the foreground lensing haloes. In this paper we investigate this issue in greater detail focusing on the Cusp relation which corresponds to images of a source located to the cusp of the inner caustic curve. We use numerical simulations combined with a Monte Carlo approach to study the effects of the expected power law distribution of substructures within LCDM haloes on the multiple images. Generally, the high number of anomalous flux ratios in the cusp configurations is unlikely explained by 'simple' perturbers (subhaloes) inside the lensing galaxy, either modeled by point masses or extended NFW subhaloes. We considered in our analysis a mass range of 10^5-10^7 Msun for the subhaloes. We also demonstrate that including the effects of the surrounding mass distribution, such as other galaxies close to the primary lens, does not change the results. We conclude that triple images of lensed QSOs do not show any direct evidence for dark dwarf galaxies such as cold dark matter substructure.