We present a computational framework for "painting" galaxies on top of the Dark Matter Halo/Sub-Halo hierarchy obtained from N-body simulations. The method we use is based on the sub-halo clustering and abundance matching (SCAM) scheme which requires observations of the 1- and 2-point statistics of the target (observed) population we want to reproduce. This method is particularly tailored for high redshift studies and thereby relies on the observed high-redshift galaxy luminosity functions and correlation properties. The core functionalities are written in c++ and exploit Object Oriented Programming, with a wide use of polymorphism, to achieve flexibility and high computational efficiency. In order to have an easily accessible interface, all the libraries are wrapped in python and provided with an extensive documentation. We validate our results and provide a simple and quantitative application to reionization, with an investigation of physical quantities related to the galaxy population, ionization fraction and bubble size distribution.

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

Using a new XMM-Newton observation, we have characterized the X-ray properties of the middle-aged radio-quiet gamma-ray pulsar J0357+3205 (named Morla) and its tail. The X-ray emission from the pulsar is consistent with a magnetospheric non-thermal origin plus a thermal emission from a hot spot (or hot spots). The lack of a thermal component from the whole surface makes Morla the coldest neutron star in its age range. We found marginal evidence for a double-peaked modulation of the X-ray emission. The study of the 9'-long tail confirmed the lack of extended emission near the pulsar itself. The tail shows a very asymmetric brightness profile and its spectrum lacks any spatial variation. We found the nebular emission to be inconsistent with a classical bow-shock, ram-pressure dominated pulsar wind nebula. We propose thermal bremsstrahlung as an alternative mechanism for Morla's tail emission. In this scenario, the tail emission comes from the shocked interstellar medium (ISM) material heated up to X-ray temperatures. This can fully explain the peculiar features of the tail, assuming a hot, moderately dense interstellar medium around the pulsar. For a bremsstrahlung-emitting tail, we can estimate the pulsar distance to be between 300 and 900 pc. A pulsar velocity of about 1900 km/s is required - which would make Morla the pulsar with the largest velocity - and high inclination angles (>70 degrees) are preferred. We propose Morla's nebula as the first example of a new "turtle's tail" class of thermally-emitting nebulae associated to high velocity pulsars.

RX J0822-4300 is the Central Compact Object associated with the Puppis A supernova remnant. Previous X-ray observations suggested RX J0822-4300 to be a young neutron star with a weak dipole field and a peculiar surface temperature distribution dominated by two antipodal spots with different temperatures and sizes. An emission line at 0.8 keV was also detected. We performed a very deep (130 ks) observation with XMM-Newton, which allowed us to study in detail the phase-resolved properties of RX J0822-4300. Our new data confirm the existence of a narrow spectral feature, best modelled as an emission line, only seen in the `Soft' phase interval - when the cooler region is best aligned to the line of sight. Surprisingly, comparison of our recent observations to the older ones yields evidence for a variation in the emission line component, which can be modelled as a decrease in the central energy from ~0.80 keV in 2001 to ~0.73 keV in 2009--2010. The line could be generated via cyclotron scattering of thermal photons in an optically thin layer of gas, or - alternatively - it could originate in low-rate accretion by a debris disk. In any case, a variation in energy, pointing to a variation of the magnetic field in the line emitting region, cannot be easily accounted for.

1E 1207.4-5209, the peculiar Central Compact object in the G296.5+10.0 supernova remnant, has been proposed to be an "anti-magnetar" - a young neutron star born with a weak dipole field. Accretion, possibly of supernova fallback material, has also been invoked to explain a large surface temperature anisotropy as well as the generation of peculiar cyclotron absorption features superimposed to its thermal spectrum. Interestingly enough, a faint optical/infrared source was proposed as a possible counterpart to 1E 1207.4-5209, but later questioned, based on coarse positional coincidence. Considering the large offset of 1E 1207.4-5209 with respect to the center of its host supernova remnant, the source should move at ~70 mas/yr. Thus, we tested the association by measuring the proper motion of the proposed optical counterpart. Using HST observations spanning 3.75 years, we computed a 3 sigma upper limit of 7 mas/yr. Absolute astrometry on the same HST data set also places the optical source significantly off the 99% confidence Chandra position. This allows us to safely rule out the association. Using the HST data set, coupled to ground-based observations collected at the ESO/VLT, we set the deepest limits ever obtained to the optical/infrared emission from 1E 1207.4-5209. By combining such limits to the constraints derived from X-ray timing, we rule out accretion as the source of the thermal anisotropy of the neutron star.

The study of the younger, and brighter, pulsars is important to understand the optical emission properties of isolated neutron stars. PSRB0540-69, the second brightest (V~22) optical pulsar, is obviously a very interesting target for these investigations. The aim of this work is threefold: constraining the pulsar proper motion and its velocity on the plane of the sky through optical astrometry, obtaining a more precise characterisation of the pulsar optical spectral energy distribution (SED) through a consistent set of multi-band, high-resolution, imaging photometry observations, measuring the pulsar optical phase-averaged linear polarisation, for which only a preliminary and uncertain measurement was obtained so far from ground-based observations. We performed high-resolution observations of PSRB0540-69 with the WFPC2 aboard the HST, in both direct imaging and polarimetry modes. From multi-epoch astrometry we set a 3sigma upper limit of 1 mas/yr on the pulsar proper motion, implying a transverse velocity <250 km/s at the 50 kpc LMC distance. Moreover, we determined the pulsar absolute position with an unprecedented accuracy of 70 mas. From multi-band photometry we characterised the pulsar power-law spectrum and we derived the most accurate measurement of the spectral index (0.70+/-0.07) which indicates a spectral turnover between the optical and X-ray bands. Finally, from polarimetry we obtained a new measurement of the pulsar phase-averaged polarisation degree (16+/-4%),consistent with magnetosphere models depending on the actual intrinsic polarisation degree and depolarisation factor, and we found that the polarisation vector (22+/-12deg position angle) is possibly aligned with the semi-major axis of the pulsar-wind nebula and with the apparent proper motion direction of its bright emission knot.