We present a new, multi-dimensional implementation of the Advanced Spectral Leakage (ASL) scheme with the purpose of modelling neutrino-matter interactions in neutron star mergers. A major challenge is the neutrino absorption in the semi-transparent regime, which is responsible for driving winds from the merger remnant. The composition of such winds is crucial in the understanding of the electromagnetic emission in the recently observed macronova following GW170817. Compared to the original version, we introduce an optical-depth-dependent flux factor to model the average angle of neutrino propagation, and a modulation that accounts for flux anisotropies in non-spherical geometries. We scrutinise our approach by first comparing the new scheme against the original one for a spherically symmetric core-collapse supernova snapshot, both in 1D and in 3D, and additionally against a two-moment (M1) scheme as implemented in 1D into the code GR1D. The luminosities and mean energies agree to a few percents in most tests. Finally, for the case of a binary merger remnant snapshot we compare the new ASL scheme with the M1 scheme that is implemented in the Eulerian adaptive mesh refinement code FLASH. We find that the neutrino absorption distribution in the semi-transparent regime is overall well reproduced. Both approaches agree to within $\lesssim 15\%$ for the average energies and to better than $\sim 35 \%$ in the total luminosities.

On the ground of the large number of gamma-ray bursts (GRBs) detected with cosmological redshift, we classified GRBs in seven subclasses, all with binary progenitors originating gravitational waves (GWs). Each binary is composed by combinations of carbon-oxygen cores (CO$_{\rm core}$), neutron stars (NSs), black holes (BHs) and white dwarfs (WDs). The long bursts, traditionally assumed to originate from a BH with an ultra-relativistic jetted emission, not emitting GWs, have been subclassified as (I) X-ray flashes (XRFs), (II) binary-driven hypernovae (BdHNe), and (III) BH-supernovae (BH-SNe). They are framed within the induced gravitational collapse (IGC) paradigm with progenitor a CO$_{\rm core}$-NS/BH binary. The supernova (SN) explosion of the CO$_{\rm core}$ triggers an accretion process onto the NS/BH. If the accretion does not lead the NS to its critical mass, an XRF occurs, while when the BH is present or formed by accretion, a BdHN occurs. When the binaries are not disrupted, XRFs lead to NS-NS and BdHNe lead to NS-BH. The short bursts, originating in NS-NS, are subclassified as (IV) short gamma-ray flashes (S-GRFs) and (V) short GRBs (S-GRBs), the latter when a BH is formed. There are (VI) ultra-short GRBs (U-GRBs) and (VII) gamma-ray flashes (GRFs), respectively formed in NS-BH and NS-WD. We use the occurrence rate and GW emission of these subclasses to assess their detectability by Advanced LIGO-Virgo, eLISA, and resonant bars. We discuss the consequences of our results in view of the announcement of the LIGO-Virgo Collaboration of the source GW 170817 as being originated by a NS-NS.