Imaging devices made from an ionization target layer of amorphous selenium (aSe) coupled to a silicon complementary metal-oxide-semiconductor (CMOS) active pixel array for charge readout are a promising technology for neutrino physics. The high spatial resolution in a solid-state target provides unparalleled rejection of backgrounds from natural radioactivity in the search for neutrinoless $\beta\beta$ decay and for electron neutrino ($\nu_e$) spectroscopy with $^{82}$Se. In this white paper, we summarize the broad scientific program of a large imaging detector with a 10-ton target of $^{82}$Se. We review the detector technology, and outline the ongoing research program to realize this experiment.

The NEXT collaboration is dedicated to the study of double beta decays of $^{136}$Xe using a high-pressure gas electroluminescent time projection chamber. This advanced technology combines exceptional energy resolution ($\leq 1\%$FWHM at the$Q_{\beta\beta}$ value of the neutrinoless double beta decay) and powerful topological event discrimination. Building on the achievements of the NEXT-White detector, the NEXT-100 detector started taking data at the Laboratorio Subterr\'aneo de Canfranc (LSC) in May of 2024. Designed to operate with xenon gas at 13.5 bar, NEXT-100 consists of a time projection chamber where the energy and the spatial pattern of the ionising particles in the detector are precisely retrieved using two sensor planes (one with photo-multiplier tubes and the other with silicon photo-multipliers). In this paper, we provide a detailed description of the NEXT-100 detector, describe its assembly, present the current estimation of the radiopurity budget, and report the results of the commissioning run, including an assessment of the detector stability.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

The ANAIS experiment is intended to search for dark matter annual modulation with ultrapure NaI(Tl) scintillators in order to provide a model independent confirmation or refutation of the long-standing DAMA/LIBRA positive annual modulation signal in the low energy detection rate, using the same target and technique. Other experiments exclude the region of parameters singled out by DAMA/LIBRA. However, these experiments use different target materials, so the comparison of their results depends on the models assumed for the dark matter particle and its distribution in the galactic halo. ANAIS-112, consisting of nine 12.5 kg NaI(Tl) modules produced by Alpha Spectra Inc., disposed in a 3$\times$3 matrix configuration, is taking data smoothly with excellent performance at the Canfranc Underground Laboratory, Spain, since August, 2017. Last published results corresponding to three-year exposure were compatible with the absence of modulation and incompatible with DAMA/LIBRA for a sensitivity above 2.5$\sigma$ C.L. Present status of the experiment and a reanalysis of the first 3 years data using new filtering protocols based on machine-learning techniques are reported. This reanalysis allows to improve the sensitivity previously achieved for the DAMA/LIBRA signal. Updated sensitivity prospects are also presented: with the improved filtering, testing the DAMA/LIBRA signal at 5$\sigma$ will be within reach in 2025.

Neutrinoless double-beta decay is a key process in particle physics. Its experimental investigation is the only viable method that can establish the Majorana nature of neutrinos, providing at the same time a sensitive inclusive test of lepton number violation. CROSS (Cryogenic Rare-event Observatory with Surface Sensitivity) aims at developing and testing a new bolometric technology to be applied to future large-scale experiments searching for neutrinoless double-beta decay of the promising nuclei $^{100}$Mo and $^{130}$Te. The limiting factor in large-scale bolometric searches for this rare process is the background induced by surface radioactive contamination, as shown by the results of the CUORE experiment. The basic concept of CROSS consists of rejecting this challenging background component by pulse-shape discrimination, assisted by a proper coating of the faces of the crystal containing the isotope of interest and serving as energy absorber of the bolometric detector. In this paper, we demonstrate that ultra-pure superconductive Al films deposited on the crystal surfaces act successfully as pulse-shape modifiers, both with fast and slow phonon sensors. Rejection factors higher than 99.9% of $\alpha$ surface radioactivity have been demonstrated in a series of prototypes based on crystals of Li$_2$MoO$_4$ and TeO$_2$. We have also shown that point-like energy depositions can be identified up to a distance of $\sim 1$ mm from the coated surface. The present program envisions an intermediate experiment to be installed underground in the Canfranc laboratory (Spain) in a CROSS-dedicated facility. This experiment, comprising $\sim 3\times 10^{25}$ nuclei of $^{100}$Mo, will be a general test of the CROSS technology as well as a worldwide competitive search for neutrinoless double-beta decay, with sensitivity to the effective Majorana mass down to 70 meV in the most favorable conditions.

The Neutrino Experiment with a Xenon Time-Projection Chamber (NEXT) is intended to investigate the neutrinoless double beta decay of 136Xe, which requires a severe suppression of potential backgrounds; therefore, an extensive screening and selection process is underway to control the radiopurity levels of the materials to be used in the experimental set-up of NEXT. The detector design combines the measurement of the topological signature of the event for background discrimination with the energy resolution optimization. Separate energy and tracking readout planes are based on different sensors: photomultiplier tubes for calorimetry and silicon multi-pixel photon counters for tracking. The design of a radiopure tracking plane, in direct contact with the gas detector medium, was specially challenging since the needed components like printed circuit boards, connectors, sensors or capacitors have typically, according to available information in databases and in the literature, activities too large for experiments requiring ultra-low background conditions. Here, the radiopurity assessment of tracking readout components based on gamma-ray spectroscopy using ultra-low background germanium detectors at the Laboratorio Subterraneo de Canfranc (Spain) is described. According to the obtained results, radiopure enough printed circuit boards made of kapton and copper, silicon photomultipliers and other required components, fulfilling the requirement of an overall background level in the region of interest of at most 8 10-4 counts keV-1 kg-1 y-1, have been identified.

The Neutrino Experiment with a Xenon TPC (NEXT), intended to investigate the neutrinoless double beta decay using a high-pressure xenon gas TPC filled with Xe enriched in 136Xe at the Canfranc Underground Laboratory in Spain, requires ultra-low background conditions demanding an exhaustive control of material radiopurity and environmental radon levels. An extensive material screening process is underway for several years based mainly on gamma-ray spectroscopy using ultra-low background germanium detectors in Canfranc but also on mass spectrometry techniques like GDMS and ICPMS. Components from shielding, pressure vessel, electroluminescence and high voltage elements and energy and tracking readout planes have been analyzed, helping in the final design of the experiment and in the construction of the background model. The latest measurements carried out will be presented and the implication on NEXT of their results will be discussed. The commissioning of the NEW detector, as a first step towards NEXT, has started in Canfranc; in-situ measurements of airborne radon levels were taken there to optimize the system for radon mitigation and will be shown too.

The DAMA experiment using ultra low background NaI(Tl) crystal scintillators has measured an annual modulation effect in the keV region which satisfies all the peculiarities of an effect induced by Dark Matter particles. In this paper we analyze this annual modulation effect in terms of mirror Dark Matter, an exact duplicate of ordinary matter from parallel hidden sector, which chemical composition is dominated by mirror helium while it can also contain significant fractions of heavier elements as Carbon and Oxygen. Dark mirror atoms are considered to interact with the target nuclei in the detector via Rutherford-like scattering induced by kinetic mixing between mirror and ordinary photons, both being massless. In the present analysis we consider various possible scenarios for the mirror matter chemical composition. For all the scenarios, the relevant ranges for the kinetic mixing parameter have been obtained taking also into account various existing uncertainties in nuclear and particle physics quantities.

ANAIS is a direct detection dark matter experiment aiming at the study of the annual modulation expected in the interaction rate. It uses same target and technique than DAMA/LIBRA experiment, which reported a highly significant positive modulation compatible with that expected for dark matter particles distributed in the galactic halo. However, other very sensitive experiments do not find any hint of particles with the required properties, although comparison is model dependent. In 2017, ANAIS-112 experiment was installed at the Canfranc Underground Laboratory (LSC), in Spain, and after the commissioning run for calibration and general assessment, ANAIS-112 started data taking in dark matter mode on August 3$^{rd}$, 2017. It consists of nine NaI(Tl) modules, amounting 112.5 kg of mass in total. ANAIS-112 will be able to test the DAMA/LIBRA result with the achieved background and threshold at three sigma level in five years of data taking. Here we report on the experimental apparatus and detector performance after the first year of data taking. Total live time available amounts to 341.72 days, being the corresponding exposure 105.32 kg x yr.

Angiogenesis is the process wherein endothelial cells (ECs) form sprouts that elongate from the pre-existing vasculature to create new vascular networks. In addition to its essential role in normal development, angiogenesis plays a vital role in pathologies such as cancer, diabetes and atherosclerosis. Mathematical and computational modelling has contributed to unravelling its complexity. Many existing theoretical models of angiogenic sprouting are based on the 'snail-trail' hypothesis. This framework assumes that leading ECs positioned at sprout tips migrate towards low-oxygen regions while other ECs in the sprout passively follow the leaders' trails and proliferate to maintain sprout integrity. However, experimental results indicate that, contrary to the snail-trail assumption, ECs exchange positions within developing vessels, and the elongation of sprouts is primarily driven by directed migration of ECs. The functional role of cell rearrangements remains unclear. This review of the theoretical modelling of angiogenesis is the first to focus on the phenomenon of cell mixing during early sprouting. We start by describing the biological processes that occur during early angiogenesis, such as phenotype specification, cell rearrangements and cell interactions with the microenvironment. Next, we provide an overview of various theoretical approaches that have been employed to model angiogenesis, with particular emphasis on recent in silico models that account for the phenomenon of cell mixing. Finally, we discuss when cell mixing should be incorporated into theoretical models and what essential modelling components such models should include in order to investigate its functional role.

We use the high-resolution HR-DEMNUni simulations to compute cross-correlations of the Cosmic Neutrino Background C$_{\nu}$B with other observables, namely: cold dark matter density, effective weak lensing convergence, neutrino velocity and neutrino deflection angle. We provide, for the first time, this high-fidelity cross-correlation maps in order to illustrate how much can be learned from them once the cosmic neutrino background is measured. As cross-correlations will have more signal than auto-correlations of the C$_{\nu}$B itself, these might be the ones to be measured first by Ptolemy-like experiments. Our predictions thus provide the imprint of what massive neutrinos should look like from cosmological observations.

The electromagnetic coupling axion-photon in a microwave cavity is revisited with the Boundary Integral - Resonant Mode Expansion (BI-RME) 3D technique. Such full-wave modal technique has been applied for the rigorous analysis of the excitation of a microwave cavity with an axion field. In this scenario, the electromagnetic field generated by the axion-photon coupling can be assumed to be driven by equivalent electrical charge and current densities. These densities have been inserted in the general BI-RME 3D equations, which express the RF electromagnetic field existing within a cavity as an integral involving the Dyadic Green functions of the cavity (under Coulomb gauge) as well as such densities. This method is able to take into account any arbitrary spatial and temporal variation of both magnitude and phase of the axion field. Next, we have obtained a simple network driven by the axion current source, which represents the coupling between the axion field and the resonant modes of the cavity. With this approach, it is possible to calculate the extracted and dissipated RF power as a function of frequency along a broad band and without Cauchy-Lorentz approximations, obtaining the spectrum of the electromagnetic field generated in the cavity, and dealing with modes relatively close to the axion resonant mode. Moreover, with this technique we have a complete knowledge of the signal extracted from the cavity, not only in magnitude but also in phase. This can be an interesting issue for future analysis where the axion phase is an important parameter.

In preparation to the CROSS $2\beta$ decay experiment, we installed a new detector suspension with magnetic dumping inside a pulse-tube cryostat of a dedicated low-background facility at the LSC (Spain). The suspension was tested with two scintillating bolometers based on large-volume 116CdWO4 (CWO-enr) and Li2MoO4 (LMO) crystals. The former, a reference device, was used for testing new noise conditions and for comparing bolometric performance of an advanced Li2MoO4 crystal developed in the framework of the CLYMENE project, in view of next-generation double-beta decay experiments like CUPID. We cooled down detectors to 15 mK and achieved high performance for all tested devices. In particular both CWO-enr and LMO bolometers demonstrated the energy resolution of 6 keV FWHM for the 2.6 MeV gamma quanta, among the best for thermal detectors based on such compounds. The baseline noise resolution (FWHM) of the CWO-enr detector was improved by 2 keV, compared to the best previous measurement of this detector in the CROSS facility, while the noise of the Ge-based optical bolometer was improved by a factor 2, to 100 eV FWHM. Despite of the evident progress in the improving of noise conditions of the set-up, we see high-frequency harmonics of a pulse-tube induced noise, suggesting a noise pick-up by cabling. Another Ge light detector was assisted with the signal amplification exploiting the Neganov-Trofimov-Luke effect, which allowed to reach 20 eV FWHM noise resolution by applying 60 V electrode bias. Highly-efficient particle identification was achieved with both detectors, despite a low scintillation efficiency of the LMO material. The radiopurity level of the LMO crystal is rather high; only traces of 210Po and 226Ra were detected (0.1 mBq/kg each), while the 228Th activity is expected to be at least an order of magnitude lower, as well as a 40K activity is found to be < 6 mBq/kg.

The CROSS experiment will search for neutrinoless double-beta decay using a specific mechanical structure to hold thermal detectors. The design of the structure was tuned to minimize the background contribution, keeping an optimal detector performance. A single module of the structure holds two scintillating bolometers (with a crystal size of 45x45x45 mm and a Ge slab facing the crystal's upper side) in the Cu frame, allowing for a modular construction of a large-scale array. Two designs are released: the initial $Thick$ version contains around 15% of Cu over the crystal mass (lithium molybdate, LMO), while this ratio is reduced to ~6% in a finer ($Slim$) design. Both designs were tested extensively at aboveground (IJCLab, France) and underground (LSC, Spain) laboratories. In particular, at LSC we used a pulse-tube-based CROSS facility to operate a 6-crystal array of LMOs enriched/depleted in $^{100}$Mo. The tested LMOs show high spectrometric performance in both designs; notably, the measured energy resolution is 5--7 keV FWHM at 2615 keV $\gamma$s, nearby the Q-value of $^{100}$Mo (3034 keV). Due to the absence of a reflective cavity around LMOs, a low scintillation signal is detected by Ge bolometers: ~0.3 keV (150 photons) for 1-MeV $\gamma$($\beta$) LMO-event. Despite that, an acceptable separation between $\alpha$ and $\gamma$($\beta$) events is achieved with most devices. The highest efficiency is reached with light detectors in the $Thick$ design thanks to a lower baseline noise width (0.05--0.09 keV RMS) when compared to that obtained in the $Slim$ version (0.10--0.35 keV RMS). Given the pivotal role of bolometric photodetectors for particle identification and random coincidences rejection, we will use the structure here described with upgraded light detectors, featuring thermal signal amplification via the Neganov-Trofimov-Luke effect, as also demonstrated in the present work.

We report on the development of thermal detectors based on large-size tellurium dioxide crystals (45x45x45 mm), containing tellurium enriched in $^{130}$Te to about 91%, for the CROSS double-beta decay experiment. A powder used for the crystals growth was additionally purified by the directional solidification method, resulting in the reduction of the concentration of impurities by a factor 10, to a few ppm of the total concentration of residual elements (the main impurity is Fe). The purest part of the ingot (the first ~200 mm, about 80% of the total length of the cylindrical part of the ingot) was determined by scanning segregation profiles of impurities and used for the $^{130}$TeO$_2$ powder production with no evidence of re-contamination. The crystal growth was verified with precursors produced from powder with natural Te isotopic composition, and two small-size (20x20x10 mm) samples were tested at a sea-level laboratory showing high bolometric and spectrometric performance together with acceptable $^{210}$Po content (below 10 mBq/kg). This growth method was then applied for the production of six large cubic $^{130}$TeO$_2$ crystals and 4 of them were taken randomly to be characterized at the Canfranc underground laboratory, in the CROSS-dedicated low-background cryogenic facility. Two $^{130}$TeO$_2$ samples were coated with a thin, $O$(100 nm), metal film in form of Al layer (on 4 sides) or AlPd grid (on a single side) to investigate the possibility to tag surface events by pulse-shape discrimination. Similarly to the small natural precursors, large-volume $^{130}$TeO$_2$ bolometers show high performance and even better internal purity ($^{210}$Po activity $\sim$ 1 mBq/kg, while activities of $^{228}$Th and $^{226}$Ra are below 0.01 mBq/kg), satisfying requirements for the CROSS and, potentially, next-generation experiments.

We report on the development of scintillating bolometers based on lithium molybdate crystals containing molybdenum depleted in the double-$\beta$ active isotope $^{100}$Mo (Li$_2$$^{100\textrm{depl}}$MoO$_4$). We used two Li$_2$$^{100\textrm{depl}}$MoO$_4$ cubic samples, 45 mm side and 0.28 kg each, produced following purification and crystallization protocols developed for double-$\beta$ search experiments with $^{100}$Mo-enriched Li$_2$MoO$_4$ crystals. Bolometric Ge detectors were utilized to register scintillation photons emitted by the Li$_2$$^{100\textrm{depl}}$MoO$_4$ crystal scintillators. The measurements were performed in the CROSS cryogenic set-up at the Canfranc underground laboratory (Spain). We observed that the Li$_2$$^{100\textrm{depl}}$MoO$_4$ scintillating bolometers are characterized by excellent spectrometric performance ($\sim$3--6 keV FWHM at 0.24--2.6 MeV $\gamma$'s), moderate scintillation signal ($\sim$0.3--0.6 keV/MeV depending on light collection conditions) and high radiopurity ($^{228}$Th and $^{226}$Ra activities are below a few $\mu$Bq/kg), comparable to the best reported results of low-temperature detectors based on Li$_2$MoO$_4$ with natural or $^{100}$Mo-enriched molybdenum content. Prospects of Li$_2$$^{100\textrm{depl}}$MoO$_4$ bolometers for use in rare-event search experiments are briefly discussed.

If neutrinoless double beta decay is discovered, the next natural step would be understanding the lepton number violating physics responsible for it. Several alternatives exist beyond the exchange of light neutrinos. Some of these mechanisms can be distinguished by measuring phase-space observables, namely the opening angle $\cos\theta$ among the two decay electrons, and the electron energy spectra, $T_1$ and $T_2$. In this work, we study the statistical accuracy and precision in measuring these kinematic observables in a future xenon gas detector with the added capability to precisely locate the decay vertex. For realistic detector conditions (a gas pressure of 10 bar and spatial resolution of 4 mm), we find that the average $\overline{\cos\theta}$ and $\overline{T_1}$ values can be reconstructed with a precision of 0.19 and 110 keV, respectively, assuming that only 10 neutrinoless double beta decay events are detected.

Low-energy neutrinos from the cosmic background are captured by objects in the sky that contain material susceptible of single beta decay. Neutrons, which compose most of a neutron star, capture low-energy neutrinos from the cosmic neutrino background and release a high-energy electron in the MeV range. Also, planets contain unstable isotopes that capture the cosmic neutrinos. We show that this process is feasible and results in a non-negligible flux of electrons in the MeV range in neutron stars. We present a novel observable, the redshift evolution of the temperature of neutron stars due to neutrino capture, that could provide a route for detection of the cosmic neutrino background from future gravitational waves observatories. For planets the flux is significantly smaller and a measurement is not possible with currently envisioned technology. While the signature from neutron stars is small and challenging, it could result in a novel way to detect the cosmic neutrino background.

A critical element in the realization of large liquid and gas time projection chambers (TPCs) is the delivery and distribution of high voltages into and around the detector. Such experiments require of order tens of kilovolts to enable electron drift over meter-scale distances. This paper describes the design and operation of the cathode feedthrough and high voltage distribution through the field cage of the NEXT-100 experiment, an underground TPC that will search for neutrinoless double beta decay $0\nu\beta\beta$. The feedthrough has been demonstrated to hold pressures up to 20~bar and sustain voltages as high as -65~kV, and the TPC is operating stably at its design high voltages. The system has been realized within the constraints of a stringent radiopurity budget and is now being used to execute a suite of sensitive double beta decay analyses.

ANAIS (Annual modulation with NaI Scintillators) is a dark matter direct detection experiment consisting of 112.5 kg of NaI(Tl) detectors in operation at the Canfranc Underground Laboratory (LSC), in Spain, since August 2017. ANAIS' goal is to confirm or refute in a model independent way the DAMA/LIBRA positive result: an annual modulation in the low-energy detection rate having all the features expected for the signal induced by dark matter particles in a standard galactic halo. This modulation, observed for about 20 years, is in strong tension with the negative results of other very sensitive experiments, but a model-independent comparison is still lacking. By using the same target material, NaI(Tl), such comparison is more direct and almost independent on dark matter particle and halo models. Here, we present the annual modulation analysis corresponding to three years of ANAIS data (for an effective exposure of 313.95 kg$\times$y), applying a blind procedure which updates that developed for the 1.5 years analysis, and later applied to 2 years. The analysis also improves the background modelling in the fitting of the region of interest rates. We obtain for the best fit in the [1-6] keV ([2-6] keV) energy region a modulation amplitude of -0.0034$\pm$0.0042 cpd/kg/keV (0.0003$\pm$0.0037 cpd/kg/keV), supporting the absence of modulation in our data, and incompatible with DAMA/LIBRA result at 3.3 (2.6) $\sigma$, for a sensitivity of 2.5 (2.7) $\sigma$. Moreover, we include two complementary analyses: a phase-free annual modulation search and the exploration of the possible presence of a periodic signal at other frequencies. Finally, we carry out several consistency checks of our result, and we update the ANAIS-112 projected sensitivity for the scheduled 5 years of operation.