The measurement of the flux of muons produced in cosmic ray air showers is essential for the study of primary cosmic rays. Such measurements are important in extensive air shower detectors to assess the energy spectrum and the chemical composition of the cosmic ray flux, complementary to the information provided by fluorescence detectors. Detailed simulations of the cosmic ray air showers are carried out, using codes such as CORSIKA, to estimate the muon flux at sea level. These simulations are based on the choice of hadronic interaction models, for which improvements have been implemented in the post-LHC era. In this work, a deficit in simulations that use state-of-the-art QCD models with respect to the measurement deep underwater with the KM3NeT neutrino detectors is reported. The KM3NeT/ARCA and KM3NeT/ORCA neutrino telescopes are sensitive to TeV muons originating mostly from primary cosmic rays with energies around 10 TeV. The predictions of state-of-the-art QCD models show that the deficit with respect to the data is constant in zenith angle; no dependency on the water overburden is observed. The observed deficit at a depth of several kilometres is compatible with the deficit seen in the comparison of the simulations and measurements at sea level.

Position calibration in the deep sea is typically done by means of acoustic multilateration using three or more acoustic emitters installed at known positions. Rather than using hydrophones as receivers that are exposed to the ambient pressure, the sound signals can be coupled to piezo ceramics glued to the inside of existing containers for electronics or measuring instruments of a deep sea infrastructure. The ANTARES neutrino telescope operated from 2006 until 2022 in the Mediterranean Sea at a depth exceeding 2000m. It comprised nearly 900 glass spheres with 432mm diameter and 15mm thickness, equipped with photomultiplier tubes to detect Cherenkov light from tracks of charged elementary particles. In an experimental setup within ANTARES, piezo sensors have been glued to the inside of such - otherwise empty - glass spheres. These sensors recorded signals from acoustic emitters with frequencies from 46545 to 60235Hz. Two waves propagating through the glass sphere are found as a result of the excitation by the waves in the water. These can be qualitatively associated with symmetric and asymmetric Lamb-like waves of zeroth order: a fast (early) one with $v_e \approx 5$mm/$\mu$s and a slow (late) one with $v_\ell \approx 2$mm/$\mu$s. Taking these findings into account improves the accuracy of the position calibration. The results can be transferred to the KM3NeT neutrino telescope, currently under construction at multiple sites in the Mediterranean Sea, for which the concept of piezo sensors glued to the inside of glass spheres has been adapted for monitoring the positions of the photomultiplier tubes.

Interest for studying cosmic neutrinos using deep-sea detectors has increase after the discovery of a diffuse flux of cosmic neutrinos by the IceCube collaboration and the possibility of wider multi-messenger studies with the observations of gravitational waves. The ANTARES detector was the first neutrino telescope in seawater, operating successfully in the Mediterranean Sea for more than a decade and a half. All challenges related to the operation in the deep sea were accurately addressed by the collaboration. Deployment and connection operations became smoother over time; data taking and constant re-calibration of the detector due to the variable environmental conditions were fully automated. A wealth of results on the subject of astroparticle physics, particle physics and multi-messenger astronomy have been obtained, despite the relative modest size of the detector, paving the way to a new generation of larger undersea detectors. This review summarizes the efforts by the ANTARES collaboration that made the possibility to operate neutrino telescopes in seawater a reality and the results obtained in this endeavor.

A search for a neutrino signal from WIMP pair annihilations in the centre of the Earth has been performed with the data collected with the ANTARES neutrino telescope from 2007 to 2012. The event selection criteria have been developed and tuned to maximise the sensitivity of the experiment to such a neutrino signal. No significant excess of neutrinos over the expected background has been observed. Upper limits at $90\%$ C.L. on the WIMP annihilation rate in the Earth and the spin independent scattering cross-section of WIMPs to nucleons $\sigma^{SI}_p$ were calculated for WIMP pair annihilations into either $\tau^{+}\tau^{-}$, $W^+W^-$, $b\overline{b}$ or the non-SUSY $\nu_{\mu}\bar{\nu}_{\mu}$ as a function of the WIMP mass (between $25\,\mathrm{GeV/c^2}$ and $1000\,\mathrm{GeV/c^2}$) and as a function of the thermally averaged annihilation cross section times velocity $\langle\sigma_{A} v\rangle_{Earth}$ of the WIMPs in the centre of the Earth. For masses of the WIMP close to the mass of iron nuclei ($50\,\mathrm{GeV/c^2}$), the obtained limits on $\sigma^{SI}_p$ are more stringent than those obtained by other indirect searches.

The ANTARES detector is at present the most sensitive neutrino telescope in the Northern Hemisphere. The highly significant cosmic neutrino excess observed by the Antarctic IceCube detector can be studied with ANTARES, exploiting its complementing field of view, exposure, and lower energy threshold. Searches for an all-flavor diffuse neutrino signal, covering 9 years of ANTARES data taking, are presented in this letter. Upward-going events are used to reduce the atmospheric muon background. This work includes for the first time in ANTARES both track-like (mainly $\nu_\mu)$ and shower-like (mainly $\nu_e$) events in this kind of analysis. Track-like events allow for an increase of the effective volume of the detector thanks to the long path traveled by muons in rock and/or sea water. Shower-like events are well reconstructed only when the neutrino interaction vertex is close to, or inside, the instrumented volume. A mild excess of high-energy events over the expected background is observed in 9 years of ANTARES data in both samples. The best fit for a single power-law cosmic neutrino spectrum, in terms of per-flavor flux at 100 TeV, is $\Phi_0^{1f}(100\ \textrm{TeV}) = \left(1.7\pm 1.0\right) \times$10$^{-18}$\,GeV$^{-1}$\,cm$^{-2}$\,s$^{-1}$\,sr$^{-1}$ with spectral index $\Gamma = 2.4^{+0.5}_{-0.4}$. The null cosmic flux assumption is rejected with a significance of 1.6$\sigma$.

A search for Secluded Dark Matter annihilation in the Sun using 2007-2012 data of the ANTARES neutrino telescope is presented. Three different cases are considered: a) detection of dimuons that result from the decay of the mediator, or neutrino detection from: b) mediator that decays into a dimuon and, in turn, into neutrinos, and c) mediator that decays directly into neutrinos. As no significant excess over background is observed, constraints are derived on the dark matter mass and the lifetime of the mediator.

High-energy neutrinos could be produced in the interaction of charged cosmic rays with matter or radiation surrounding astrophysical sources. To look for transient sources associated with neutrino emission, a follow-up program of neutrino alerts has been operating within the ANTARES Collaboration since 2009. This program, named TAToO, has triggered robotic optical telescopes (MASTER, TAROT, ROTSE and the SVOM ground based telescopes) immediately after the detection of any relevant neutrino candidate and scheduled several observations in the weeks following the detection. A subset of ANTARES events with highest probabilities of being of cosmic origin has also been followed by the Swift and the INTEGRAL satellites, the Murchison Widefield Array radio telescope and the H.E.S.S. high-energy gamma-ray telescope. The results of twelve years of observations are reported. No optical counterpart has been significantly associated with an ANTARES candidate neutrino signal during image analysis. Constraints on transient neutrino emission have been set. In September 2015, ANTARES issued a neutrino alert and during the follow-up, a potential transient counterpart was identified by Swift and MASTER. A multi-wavelength follow-up campaign has allowed to identify the nature of this source and has proven its fortuitous association with the neutrino. The return of experience is particularly important for the design of the alert system of KM3NeT, the next generation neutrino telescope in the Mediterranean Sea.

This study presents a novel search for magnetic monopoles using data collected over a 14 year period (2008-2022) by the ANTARES neutrino telescope. The interaction of magnetic monopoles with matter was modeled according to Kazama, Yang, and Goldhaber cross-section. Upper limits on the flux of magnetic monopoles are obtained for velocities both above and below the Cherenkov threshold. No events consistent with the passage of magnetic monopoles were detected, enabling the setting of an upper flux limit for relativistic magnetic monopoles of the order of $10^{-18} \mathrm{cm}^{-2} \mathrm{s}^{-1} \mathrm{sr}^{-1}$.

Gamma-ray bursts are promising candidate sources of high-energy astrophysical neutrinos. The recent GRB 221009A event, identified as the brightest gamma-ray burst ever detected, provides a unique opportunity to investigate hadronic emissions involving neutrinos. The KM3NeT undersea neutrino detectors participated in the worldwide follow-up effort triggered by the event, searching for neutrino events. In this letter, we summarize subsequent searches, in a wide energy range from MeV up to a few PeVs. No neutrino events are found in any of the searches performed. Upper limits on the neutrino emission associated with GRB 221009A are computed.

The flux of very high-energy neutrinos produced in our Galaxy by the interaction of accelerated cosmic rays with the interstellar medium is not yet determined. The characterization of this flux will shed light on Galactic accelerator features, gas distribution morphology and Galactic cosmic ray transport. The central Galactic plane can be the site of an enhanced neutrino production, thus leading to anisotropies in the extraterrestrial neutrino signal as measured by the IceCube Collaboration. The ANTARES neutrino telescope, located in the Mediterranean Sea, offers a favourable view on this part of the sky, thereby allowing for a contribution to the determination of this flux. The expected diffuse Galactic neutrino emission can be obtained linking a model of generation and propagation of cosmic rays with the morphology of the gas distribution in the Milky Way. In this paper, the so-called "Gamma model" introduced recently to explain the high-energy gamma ray diffuse Galactic emission, is assumed as reference. The neutrino flux predicted by the "Gamma model" depends of the assumed primary cosmic ray spectrum cut-off. Considering a radially-dependent diffusion coefficient, this proposed scenario is able to account for the local cosmic ray measurements, as well as for the Galactic gamma ray observations. Nine years of ANTARES data are used in this work to search for a possible Galactic contribution according to this scenario. All flavour neutrino interactions are considered. No excess of events is observed and an upper limit is set on the neutrino flux of $1.1$ ($1.2$) times the prediction of the "Gamma model" assuming the primary cosmic ray spectrum cut-off at 5 (50) PeV. This limit excludes the diffuse Galactic neutrino emission as the major cause of the "spectral anomaly" between the two hemispheres measured by IceCube.

On October 1, 2019, the IceCube Collaboration detected a muon track neutrino with high probability of being of astrophysical origin, IC191001A. After a few hours, the tidal disruption event (TDE) AT2019dsg, observed by the Zwicky Transient Facility (ZTF), was indicated as the most likely counterpart of the IceCube track. More recently, the follow-up campaign of the IceCube alerts by ZTF suggested a second TDE, AT2019fdr, as a promising counterpart of another IceCube muon track candidate, IC200530A, detected on May 30, 2020. These are the second and third associations between astrophysical sources and high-energy neutrinos after the compelling identification of the blazar TXS 0506+056. Here, the search for ANTARES neutrinos from the directions of AT2019dsg and AT2019fdr using a time-integrated approach is presented. As no significant evidence for space clustering is found in the ANTARES data, upper limits on the one-flavour neutrino flux and fluence are set.

For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

Searches for dark matter (DM) have not provided any solid evidence for the existence of weakly interacting massive particles in the GeV-TeV mass range. Coincidentally, the scale of new physics is being pushed by collider searches well beyond the TeV domain. This situation strongly motivates the exploration of DM masses much larger than a TeV. Secluded scenarios contain a natural way around the unitarity bound on the DM mass, via the early matter domination induced by the mediator of its interactions with the Standard Model. High-energy neutrinos constitute one of the very few direct accesses to energy scales above a few TeV. An indirect search for secluded DM signals has been performed with the ANTARES neutrino telescope using data from 2007 to 2015. Upper limits on the DM annihilation cross section for DM masses up to 6 PeV are presented and discussed.

By constantly monitoring at least one complete hemisphere of the sky, neutrino telescopes are well designed to detect neutrinos emitted by transient astrophysical events. Real-time searches with the ANTARES telescope have been performed to look for neutrino candidates coincident with gamma-ray bursts detected by the Swift and Fermi satellites, highenergy neutrino events registered by IceCube, transient events from blazars monitored by HAWC, photon-neutrino coincidences by AMON notices and gravitational wave candidates observed by LIGO/Virgo. By requiring temporal coincidence, this approach increases the sensitivity and the significance of a potential discovery. Thanks to the good angular accuracy of neutrino candidates reconstructed with the ANTARES telescope, a coincident detection can also improve the positioning area of non-well localised triggers such as those detected by gravitational wave interferometers. This paper summarises the results of the follow-up performed by the ANTARES telescope between 01/2014 and 02/2022, which corresponds to the end of the data taking period.

Since 2015 the LIGO and Virgo interferometers have detected gravitational waves from almost one hundred coalescences of compact objects (black holes and neutron stars). This article presents the results of a search performed with data from the ANTARES telescope to identify neutrino counterparts to the gravitational wave sources detected during the third LIGO/Virgo observing run and reported in the catalogues GWTC-2, GWTC-2.1, and GWTC-3. This search is sensitive to all-sky neutrinos of all flavours and of energies >100 GeV, thanks to the inclusion of both track-like events (mainly induced by $\nu_\mu$ charged-current interactions) and shower-like events (induced by other interaction types). Neutrinos are selected if they are detected within $\pm 500$ s from the GW merger and with a reconstructed direction compatible with its sky localisation. No significant excess is found for any of the 80 analysed GW events, and upper limits on the neutrino emission are derived. Using the information from the GW catalogues and assuming isotropic emission, upper limits on the total energy $E_{\rm tot, \nu}$ emitted as neutrinos of all flavours and on the ratio $f_\nu = E_{\rm tot, \nu}/E_{\rm GW}$ between neutrino and GW emissions are also computed. Finally, a stacked analysis of all the 72 binary black hole mergers (respectively the 7 neutron star - black hole merger candidates) has been performed to constrain the typical neutrino emission within this population, leading to the limits: E_{\rm tot, \nu} < 4.0 \times 10^{53} erg and f_\nu < 0.15 (respectively, E_{\rm tot, \nu} < 3.2 \times 10^{53} erg and f_\nu < 0.88) for $E^{-2}$ spectrum and isotropic emission. Other assumptions including softer spectra and non-isotropic scenarios have also been tested.