Dark matter models involving a very light bosonic particle, generally known as Fuzzy Dark Matter (FDM), have been recently attracting great interest in the cosmology community, as their wave-like phenomenology would simultaneously explain the longstanding mis-detection of a dark matter particle and help easing the small-scale issues related to the standard Cold Dark Matter (CDM) scenario. With the present work, we initiate a series of papers aiming at investigating the evolution of FDM structures in a cosmological framework performed with our N-body code AX-GADGET, detailing for the first time in the literature how the actual scaling relations between solitonic cores and host haloes properties are significantly affected by the dynamical state, morphology and merger history of the individual systems. In particular, in this first paper we confirm the ability of AX-GADGET to correctly reproduce the typical FDM solitonic core and we employ it to study the non-linear evolution of eight FDM haloes in their cosmological context through the zoom-in simulation approach. We find that the scaling relations identified in previous works for isolated systems are generally modified for haloes evolving in a realistic cosmological environment, and appear to be valid only as a limit for the most relaxed and spherically symmetric systems.

The advent of a new generation of large-scale galaxy surveys is pushing cosmological numerical simulations in an uncharted territory. The simultaneous requirements of high resolution and very large volume pose serious technical challenges, due to their computational and data storage demand. In this paper, we present a novel approach dubbed Dynamic Zoom Simulations -- or DZS -- developed to tackle these issues. Our method is tailored to the production of lightcone outputs from N-body numerical simulations, which allow for a more efficient storage and post-processing compared to standard comoving snapshots, and more directly mimic the format of survey data. In DZS, the resolution of the simulation is dynamically decreased outside the lightcone surface, reducing the computational work load, while simultaneously preserving the accuracy inside the lightcone and the large-scale gravitational field. We show that our approach can achieve virtually identical results to traditional simulations at half of the computational cost for our largest box. We also forecast this speedup to increase up to a factor of 5 for larger and/or higher-resolution simulations. We assess the accuracy of the numerical integration by comparing pairs of identical simulations run with and without DZS. Deviations in the lightcone halo mass function, in the sky-projected lightcone, and in the 3D matter lightcone always remain below 0.1%. In summary, our results indicate that the DZS technique may provide a highly-valuable tool to address the technical challenges that will characterise the next generation of large-scale cosmological simulations.

Fuzzy Dark Matter (FDM) represents an alternative and intriguing description of the standard Cold Dark Matter (CDM) fluid, able to explain the lack of direct detection of dark matter particles in the GeV sector and to alleviate small scales tensions in the cosmic large-scale structure formation. Cosmological simulations of FDM models in the literature were performed either with very expensive high-resolution grid-based simulations of individual haloes or through N-body simulations encompassing larger cosmic volumes but resorting on significant approximations in the FDM non-linear dynamics to reduce their computational cost. With the use of the new N-body cosmological hydrodynamical code AX-GADGET, we are now able not only to overcome such numerical problems, but also to combine a fully consistent treatment of FDM dynamics with the presence of gas particles and baryonic physical processes, in order to quantify the FDM impact on specific astrophysical observables. In particular, in this paper we perform and analyse several hydrodynamical simulations in order to constrain the FDM mass by quantifying the impact of FDM on Lyman-$\alpha$ forest observations, as obtained for the first time in the literature in a N-body setup without approximating the FDM dynamics. We also study the statistical properties of haloes, exploiting the large available sample, to extract information on how FDM affects the abundance, the shape, and density profiles of dark matter haloes.

We present a new module of the parallel N-Body code P-GADGET3 for cosmological simulations of light bosonic non-thermal dark matter, often referred as Fuzzy Dark Matter (FDM). The dynamics of the FDM features a highly non-linear Quantum Potential (QP) that suppresses the growth of structures at small scales. Most of the previous attempts of FDM simulations either evolved suppressed initial conditions, completely neglecting the dynamical effects of QP throughout cosmic evolution, or resorted to numerically challenging full-wave solvers. The code provides an interesting alternative, following the FDM evolution without impairing the overall performance. This is done by computing the QP acceleration through the Smoothed Particle Hydrodynamics (SPH) routines, with improved schemes to ensure precise and stable derivatives. As an extension of the P-GADGET3 code, it inherits all the additional physics modules implemented up to date, opening a wide range of possibilities to constrain FDM models and explore its degeneracies with other physical phenomena. Simulations are compared with analytical predictions and results of other codes, validating the QP as a crucial player in structure formation at small scales.

According to the second Fermi LAT Catalog (2FGL), about one third of the gamma-ray sources listed have no assigned counterparts at lower energies. Many statistical methods have been developed to find proper counterparts for these sources. We explore the sky area covered at low radio frequency by Westerbork in the Southern Hemisphere (WISH) survey to search for blazar-like associations among the unidentified gamma-ray sources listed in the 2FGL (UGSs). Searching the WISH and NRAO VLA Sky Survey (NVSS) radio surveys within the positional uncertainty regions of the 2FGL UGSs, we select as gamma-ray blazar candidates the radio sources characterized by flat radio spectra between 352 MHz and 1400 MHz. We propose new gamma-ray blazar associations for eight UGSs and we also discuss their spectral properties at low radio frequencies. We compare the radio flux density distribution of the low radio frequency gamma-ray blazar candidates with that of gamma-ray blazars associated with other methods. We find significant differences between these distributions. Finally, we discuss the results of this association method and its possible applicability to other regions of the sky and future radio surveys.

About one third of the gamma-ray sources listed in the second Fermi LAT catalog (2FGL) have no firmly established counterpart at lower energies so being classified as unidentified gamma-ray sources (UGSs). Here we propose a new approach to find candidate counterparts for the UGSs based on the 325 MHz radio survey performed with Westerbork Synthesis Radio Telescope (WSRT) in the northern hemisphere. First we investigate the low-frequency radio properties of blazars, the largest known population of gamma-ray sources; then we search for sources with similar radio properties combining the information derived from the Westerbork Northern Sky Survey (WENSS) with those of the NRAO VLA Sky survey (NVSS). We present a list of candidate counterparts for 32 UGSs with at least one counterpart in the WENSS. We also performed an extensive research in literature to look for infrared and optical counterparts of the gamma-ray blazar candidates selected with the low-frequency radio observations to confirm their nature. On the basis of our multifrequency research we identify 23 new gamma-ray blazar candidates out of 32 UGSs investigated. Comparison with previous results on the UGSs are also presented. Finally, we speculate on the advantages on the use of the low-frequency radio observations to associate UGSs and to search for gamma-ray pulsar candidates.