The Extreme Energy Events (EEE) Project employs Multi-gap Resistive Plate Chamber (MRPC) for studying the secondary cosmic ray muons in Extensive Air Showers. The array consists of about 60 tracking detectors, sparse on Italian territory and at CERN. The MRPCs are flowed with a gas mixture based on $C_2H_2F_4$ and $SF_6$, both of which are fluorinated greenhouse gases with a high environmental impact on the atmosphere. Due to the restrictions imposed by the European Union, these gases are being phased out of production and their cost is largely increasing. The EEE Collaboration started a campaign to reduce the gas emission from its array with the aim of containing costs and decreasing the experiment global warming impact. One method is to reduce the gas rate in each EEE detector. Another is to develop a gas recirculation system, whose a first prototype has been installed at one of the EEE stations located at CERN. Jointly a parallel strategy is focused on searching for environmental friendly gas mixtures which are able to substitute the standard mixture without affecting the MRPC performance. An overview and first results are presented here.

We continue our study of the binary neutron star parameter space by investigating the effect of the spin orientation on the dynamics, gravitational wave emission, and mass ejection during the binary neutron star coalescence. We simulate seven different configurations using multiple resolutions to allow a reasonable error assessment. Due to the particular choice of the setups, five configurations show precession effects, from which two show a precession (`wobbling') of the orbital plane, while three show a `bobbing' motion, i.e., the orbital angular momentum does not precess, while the orbital plane moves along the orbital angular momentum axis. Considering the ejection of mass, we find that precessing systems can have an anisotropic mass ejection, which could lead to a final remnant kick of about $\sim 40 \rm km/s$ for the studied systems. Furthermore, for the chosen configurations, anti-aligned spins lead to larger mass ejecta than aligned spins, so that brighter electromagnetic counterparts could be expected for these configurations. Finally, we compare our simulations with the precessing, tidal waveform approximant IMRPhenomPv2_NRTidalv2 and find good agreement between the approximant and our numerical relativity waveforms with phase differences below 1.2 rad accumulated over the last $\sim$ 16 gravitational wave cycles.

Neutrino-matter interactions play a key role in binary neutron star mergers. Thermodynamics conditions at the surfaces where neutrinos decouple from matter influence neutrino spectra, ultimately affecting the evolution of the remnant and the properties of the ejecta. In this work, we post-process results of general relativistic merger simulations employing microphysical equations of state and approximate neutrino transport to investigate the thermodynamics conditions at which weak and thermal equilibrium freezes out (equilibrium surfaces), as well as conditions at which the transition between diffusion and free-streaming regime occurs (diffusion surfaces). We find that the rest mass density and the neutrino energy are the most relevant quantities in determining the location of the decoupling surfaces. For mean energy neutrinos ($\langle E_{{\nu}_e} \rangle \approx 9~{\rm MeV}$, $\langle E_{{\bar{\nu}}_e} \rangle \approx 15~{\rm MeV}$, $\langle E_{{\nu}_{\mu,\tau}} \rangle \approx 25~{\rm MeV}$), diffusion surfaces are located around $10^{11}{\rm g~cm^{-3}}$ for all neutrino species, while equilibrium surfaces for heavy flavor neutrinos are significantly deeper (several $10^{12}{\rm g~cm^{-3}}$) than the ones of $\bar{\nu}_e$ and $\nu_e$ ($\gtrsim 10^{11}{\rm g~cm^{-3}}$). The resulting decoupling temperatures are in good agreement with the average neutrino energies ($\langle E_{\nu} \rangle \sim 3.15~T$), with the softer equation of state characterized by systematically larger decoupling temperatures ($\Delta T \lesssim 1~{\rm MeV}$). Neutrinos streaming at infinity with different energies come from very different regions of the remnant. The presence of a massive NS or of a BH in the remnant influences the neutrino thermalization process.

The Extreme Energy Events Project is a synchronous sparse array of 52 tracking detectors for studying High Energy Cosmic Rays (HECR) and Cosmic Rays-related phenomena. The observatory is also meant to address Long Distance Correlation (LDC) phenomena: the network is deployed over a broad area covering 10 degrees in latitude and 11 in longitude. An overview of a set of preliminary results is given, extending from the study of local muon flux dependance on solar activity to the investigation of the upward-going component of muon flux traversing the EEE stations; from the search for anisotropies at the sub-TeV scale to the hints for observations of km-scale Extensive Air Shower (EAS).