Researcher profile

S. Zavatarelli

S. Zavatarelli contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
21works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2026arXiv

Detector characterization for a new $^{12}$C+$^{12}$C reaction study at LUNA

The $^{12}$C+$^{12}$C fusion reaction plays a crucial role in stellar evolution, including the occurrence of supernova explosions, and in the synthesis of the chemical elements. However, our understanding of its cross section remains severely deficient, particularly below $E_\textrm{cm}=2.5$\,MeV, the energy range of interest for astrophysics. To address these unresolved issues, the LUNA collaboration will conduct a dedicated study of the $^{12}$C+$^{12}$C reaction at the Bellotti Ion Beam Facility (Bellotti IBF) located deep underground within the Gran Sasso National Laboratory (LNGS) in Italy. Based on the combination of passive and active shields, this campaign aims to achieve unprecedented sensitivity in measuring the cross sections of the two key reaction channels, $^{12}$C($^{12}$C,$α$)$^{20}$Ne and $^{12}$C($^{12}$C,$p$)$^{23}$Na in the low-energy regime via $γ$-ray detection. Here, we report on a sensitivity study for the upcoming campaign with a focus on the characterization of two detectors, namely a HPGe detector and a NaI(Tl) array. Furthermore, their intrinsic contamination is thoroughly investigated since this could potentially influence the overall sensitivity. Assuming typical beam intensities of the Bellotti IBF, we will be able to investigate reaction rates significantly below 100 counts per day. In case of the $^{12}$C+$^{12}$C reaction we therefore expect to acquire experimental data well below the current limit of $E_\textrm{cm}=2.1\,$MeV. The results are supported by simulations to highlight the advantageous low-background environment, essential for high-precision nuclear astrophysics studies.

preprint2026arXiv

Optimizing the potential of KM3NeT in detecting core-collapse supernovae

Core-collapse supernovae mark the end of life of massive stars. However, despite their importance in astrophysics, their underlying mechanisms remain unclear. Neutrinos that emerge from the dense core of the star offer a promising way to study supernova dynamics. A strategy is presented to improve the potential of the KM3NeT neutrino telescope to detect core-collapse supernovae in our Galaxy or the Large Magellanic Cloud by further exploiting the properties of its optical modules equipped with multiple photomultipliers. A supernova burst is expected to produce a sudden hit rate increase in the KM3NeT detectors. New observables have been defined for individual optical modules that exploit the geometry and time distribution of the detected hits, enabling a better discrimination between signal and background signatures. In addition, a thorough investigation of the related systematic uncertainties is presented for the first time. When implemented, this new methodology allowed KM3NeT to probe 46% more Galactic core-collapse supernova candidates than with the previous trigger strategy, reaching the dense Galactic bulge. It is now expected that, once completed, KM3NeT will achieve full Galactic sensitivity to core-collapse supernovae.

preprint2022arXiv

Search for Low-Energy Signals from Fast Radio Bursts with the Borexino Detector

The search for neutrino events in correlation with several of the most intense fast radio bursts (FRBs) has been performed using the Borexino data. We have searched for signals with visible energies above $250$~keV within a time window of $\pm$1000~s corresponding to the detection time of a particular FRB. We also applied an alternative approach based on searching for specific shapes of neutrino-electron scattering spectra in the full exposure spectrum of the Borexino detector. In particular, two incoming neutrino spectra were considered: the monoenergetic line and the spectrum expected from supernovae. The same spectra were considered for electron antineutrinos detected through the inverse beta-decay reaction. No statistically significant excess over the background was observed. As a result, the strongest upper limits on FRB-associated neutrino fluences of all flavors have been obtained in the $0.5 - 50$~MeV neutrino energy range.

preprint2022arXiv

Search for non-standard neutrino interactions with 10 years of ANTARES data

Non-standard interactions of neutrinos arising in many theories beyond the Standard Model can significantly alter matter effects in atmospheric neutrino propagation through the Earth. In this paper, a search for deviations from the prediction of the standard 3-flavour atmospheric neutrino oscillations using the data taken by the ANTARES neutrino telescope is presented. Ten years of atmospheric neutrino data collected from 2007 to 2016, with reconstructed energies in the range from $\sim$16 GeV to $100$ GeV, have been analysed. A log-likelihood ratio test of the dimensionless coefficients $\varepsilon_{μτ}$ and $\varepsilon_{ττ} - \varepsilon_{μμ}$ does not provide clear evidence of deviations from standard interactions. For normal neutrino mass ordering, the combined fit of both coefficients yields a value 1.7$σ$ away from the null result. However, the 68% and 95% confidence level intervals for $\varepsilon_{μτ}$ and $\varepsilon_{ττ} - \varepsilon_{μμ}$, respectively, contain the null value. Best fit values, one standard deviation errors and bounds at the 90% confidence level for these coefficients are given for both normal and inverted mass orderings. The constraint on $\varepsilon_{μτ}$ is among the most stringent to date and it further restrains the strength of possible non-standard interactions in the $μ- τ$ sector.

preprint2022arXiv

Search for secluded dark matter towards the Galactic Centre with the ANTARES neutrino telescope

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.

preprint2022arXiv

Search for solar atmospheric neutrinos with the ANTARES neutrino telescope

Solar Atmospheric Neutrinos (SA$ν$s) are produced by the interaction of cosmic rays with the solar medium. The detection of SA$ν$s would provide useful information on the composition of primary cosmic rays as well as the solar density. These neutrinos represent an irreducible source of background for indirect searches for dark matter towards the Sun and the measurement of their flux would allow for a better assessment of the uncertainties related to these searches. In this paper we report on the analysis performed, based on an unbinned likelihood maximisation, to search for SA$ν$s with the ANTARES neutrino telescope. After analysing the data collected over 11 years, no evidence for a solar atmospheric neutrino signal has been found. An upper limit at 90\% confidence level on the flux of solar atmospheric neutrinos has been obtained, equal to 7$\times$$10^{-11}$ [TeV$^{-1}$cm$^{-2}$s$^{-1}$] at E$_ν=$ 1 TeV for the reference cosmic ray model assumed.

preprint2022arXiv

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

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.

preprint2022arXiv

The study of key reactions shaping the post main-sequence evolution of massive stars in underground facilities

The chemical evolution of the Universe and several phases of the stellar life are regulated by minute nuclear reactions. The key point for each of these reactions is the value of cross sections at the energies at which they take place in stellar environments. Direct cross-section measurements are mainly hampered by the very low counting rate and by cosmic background, nevertheless they have become possible by combining the best experimental techniques with the cosmic silence of an underground laboratory. In the nineties the LUNA (Laboratory for Underground Nuclear Astrophysics) collaboration opened the era of underground nuclear astrophysics installing first a home-made 50 kV and later on, a second 400 kV accelerator under the Gran Sasso mountain in Italy: in 25 years of experimental activity, important reactions responsible for hydrogen burning could have been studied down to the relevant energies thanks to the high current proton and helium beams provided by the machines. The interest to the next and warmer stages of star evolution (i.e. post main sequence, helium and carbon burning) drove a new project based on an ion accelerator in the MV range called LUNA-MV, able to deliver proton, helium and carbon beams. The present contribution is aimed to discuss the state-of-the-art for some selected key processes of post main sequence stellar phases: the 12Cag and the 12C+12C fundamental for helium and carbon burning phases, and 13Can and 22Neag that are relevant to the synthesis of heavy elements in AGB stars. The perspectives opened by an underground MV-facility will be highlighted.

preprint2021arXiv

Correlated and Integrated Directionality for sub-MeV solar neutrinos in Borexino

Liquid scintillator detectors play a central role in the detection of neutrinos from various sources. In particular, it is the only technique used so far for the precision spectroscopy of sub-MeV solar neutrinos, as demonstrated by the Borexino experiment at the Gran Sasso National Laboratory in Italy. The benefit of a high light yield, and thus a low energy threshold and a good energy resolution, comes at the cost of the directional information featured by water Cherenkov detectors, measuring $^8$B solar neutrinos above a few MeV. In this paper we provide the first directionality measurement of sub-MeV solar neutrinos which exploits the correlation between the first few detected photons in each event and the known position of the Sun for each event. This is also the first signature of directionality in neutrinos elastically scattering off electrons in a liquid scintillator target. This measurement exploits the sub-dominant, fast Cherenkov light emission that precedes the dominant yet slower scintillation light signal. Through this measurement, we have also been able to extract the rate of $^{7}$Be solar neutrinos in Borexino. The demonstration of directional sensitivity in a traditional liquid scintillator target paves the way for the possible exploitation of the Cherenkov light signal in future kton-scale experiments using liquid scintillator targets. Directionality is important for background suppression as well as the disentanglement of signals from various sources.

preprint2021arXiv

Determining the Neutrino Mass Ordering and Oscillation Parameters with KM3NeT/ORCA

The next generation of water Cherenkov neutrino telescopes in the Mediterranean Sea are under construction offshore France (KM3NeT/ORCA) and Sicily (KM3NeT/ARCA). The KM3NeT/ORCA detector features an energy detection threshold which allows to collect atmospheric neutrinos to study flavour oscillation. This paper reports the KM3NeT/ORCA sensitivity to this phenomenon. The event reconstruction, selection and classification are described. The sensitivity to determine the neutrino mass ordering was evaluated and found to be 4.4 $σ$ if the true ordering is normal and 2.3 $σ$ if inverted, after three years of data taking. The precision to measure $Δm^2_{32}$ and $θ_{23}$ were also estimated and found to be $85\cdot10^{-6}$ eV$^2$ and $(^{+1.9}_{-3.1})^{\circ}$ for normal neutrino mass ordering and, $75\cdot10^{-6}$ eV$^2$ and $(^{+2.0}_{-7.0})^{\circ}$ for inverted ordering. Finally, a unitarity test of the leptonic mixing matrix by measuring the rate of tau neutrinos is described. Three years of data taking were found to be sufficient to exclude $ν_τ$ and $\barν_τ$ event rate variations larger than 20% at $3σ$ level.

preprint2021arXiv

First Directional Measurement of sub-MeV Solar Neutrinos with Borexino

We report the measurement of sub-MeV solar neutrinos through the use of their associated Cherenkov radiation, performed with the Borexino detector at the Laboratori Nazionali del Gran Sasso. The measurement is achieved using a novel technique that correlates individual photon hits of events to the known position of the Sun. In an energy window between 0.54 MeV to 0.74 MeV, selected using the dominant scintillation light, we have measured 10887$^{+2386}_{-2103} (\mathrm{stat.})\pm 947 (\mathrm{syst.})$ ($68\%$ confidence interval) solar neutrinos out of 19904 total events. This corresponds to a $^{7}$Be neutrino interaction rate of 51.6$^{+13.9}_{-12.5}$ counts/(day$\cdot$ 100 ton), which is in agreement with the Standard Solar Model predictions and the previous spectroscopic results of Borexino. The no-neutrino hypothesis can be excluded with $>$5$σ$ confidence level. For the first time, we have demonstrated the possibility of utilizing the directional Cherenkov information for sub-MeV solar neutrinos, in a large-scale, high light yield liquid scintillator detector. This measurement provides an experimental proof of principle for future hybrid event reconstruction using both Cherenkov and scintillation signatures simultaneously.

preprint2021arXiv

Implementation and first results of the KM3NeT real-time core-collapse supernova neutrino search

The KM3NeT research infrastructure is under construction in the Mediterranean Sea. KM3NeT will study atmospheric and astrophysical neutrinos with two multi-purpose neutrino detectors, ARCA and ORCA, primarily aimed at GeV-PeV neutrinos. Thanks to the multi-photomultiplier tube design of the digital optical modules, KM3NeT is capable of detecting the neutrino burst from a Galactic or near-Galactic core-collapse supernova. This potential is already exploitable with the first detection units deployed in the sea. This paper describes the real-time implementation of the supernova neutrino search, operating on the two KM3NeT detectors since the first months of 2019. A quasi-online astronomy analysis is introduced to study the time profile of the detected neutrinos for especially significant events. The mechanism of generation and distribution of alerts, as well as the integration into the SNEWS and SNEWS 2.0 global alert systems are described. The approach for the follow-up of external alerts with a search for a neutrino excess in the archival data is defined. Finally, an overview of the current detector capabilities and a report after the first two years of operation are given.

preprint2021arXiv

Nanobeacon: A time calibration device for the KM3NeT neutrino telescope

The KM3NeT Collaboration is currently constructing a multi-site high-energy neutrino telescope in the Mediterranean Sea consisting of matrices of pressure-resistant glass spheres, each holding a set of 31 small-area photomultipliers. The main goals of the telescope are the observation of neutrino sources in the Universe and the measurement of the neutrino oscillation parameters with atmospheric neutrinos. Both extraterrestrial and atmospheric neutrinos are detected through the Cherenkov light induced in seawater by charged particles produced in neutrino interactions in the surrounding medium. A relative time synchronization between photomultipliers of the order of 1 ns is needed to guarantee the required angular resolution of the detector. Due to the large detector volumes to be instrumented by KM3NeT, a cost reduction of the different systems is a priority. To this end, the inexpensive Nanobeacon has been designed and developed by the KM3NeT Collaboration to be used for detector time-calibration studies. At present, more than 600 Nanobeacons have been already produced. The characterization of the optical pulse and the wavelength emission profile of the devices are critical for the time calibration. In this paper, the main features of the Nanobeacon design, production and operation, together with the main properties of the light pulse generated are described.

preprint2021arXiv

Search for neutrinos from the tidal disruption events AT2019dsg and AT2019fdr with the ANTARES telescope

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.

preprint2020arXiv

A new approach to monitor 13C-targets degradation in situ for 13C(alpha,n)16O cross-section measurements at LUNA

Direct measurements of reaction cross-sections at astrophysical energies often require the use of solid targets able to withstand high ion beam currents for extended periods of time. Thus, monitoring target thickness, isotopic composition, and target stoichiometry during data taking is critical to account for possible target modifications and to reduce uncertainties in the final cross-section results. A common technique used for these purposes is the Nuclear Resonant Reaction Analysis (NRRA), which however requires that a narrow resonance be available inside the dynamic range of the accelerator used. In cases when this is not possible, as for example the 13C(alpha,n)16O reaction recently studied at low energies at the Laboratory for Underground Nuclear Astrophysics (LUNA) in Italy, alternative approaches must be found. Here, we present a new application of the shape analysis of primary gamma rays emitted by the 13C(p,g)14N radiative capture reaction. This approach was used to monitor 13C target degradation {\em in situ} during the 13C(alpha,n)16O data taking campaign. The results obtained are in agreement with evaluations subsequently performed at Atomki (Hungary) using the NRRA method. A preliminary application for the extraction of the 13C(alpha,n)16O reaction cross-section at one beam energy is also reported.

preprint2020arXiv

Comprehensive geoneutrino analysis with Borexino

This paper presents a geoneutrino measurement using 3262.74 days of data taken with the Borexino detector at LNGS in Italy. By observing $52.6 ^{+9.4}_{-8.6} ({\rm stat}) ^{+2.7}_{-2.1}({\rm sys})$ geoneutrinos (68% interval) from $^{238}$U and $^{232}$Th, a signal of $47.0^{+8.4}_{-7.7}\,({\rm stat)}^{+2.4}_{-1.9}\,({\rm sys})$ TNU with $^{+18.3}_{-17.2}$% total precision was obtained. This result assumes the same Th/U mass ratio found in chondritic CI meteorites but compatible results were found when contributions from $^{238}$U and $^{232}$Th were fit as free parameters. Antineutrino background from reactors is fit unconstrained and found compatible with the expectations. The null-hypothesis of observing a signal from the mantle is excluded at a 99.0% C.L. when exploiting the knowledge of the local crust. Measured mantle signal of $21.2 ^{+9.6}_{-9.0} ({\rm stat})^{+1.1}_{-0.9} ({\rm sys})$ TNU corresponds to the production of a radiogenic heat of $24.6 ^{+11.1}_{-10.4}$ TW (68% interval) from $^{238}$U and $^{232}$Th in the mantle. Assuming 18% contribution of $^{40}$K in the mantle and $8.1^{+1.9}_{-1.4}$ TW of radiogenic heat of the lithosphere, the Borexino estimate of the total Earth radiogenic heat is $38.2 ^{+13.6}_{-12.7}$ TW, corresponding to a convective Urey ratio of 0.78$^{+0.41}_{-0.28}$. These values are compatible with different geological models, however there is a 2.4$σ$ tension with those which predict the lowest concentration of heat-producing elements. By fitting the data with a constraint on the reactor antineutrino background, the existence of a hypothetical georeactor at the center of the Earth having power greater than 2.4 TW at 95% C.L. is excluded. Particular attention is given to all analysis details, which should be of interest for the next generation geoneutrino measurements.

preprint2020arXiv

Constraints on Flavor-Diagonal Non-Standard Neutrino Interactions from Borexino Phase-II

The Borexino detector measures solar neutrino fluxes via neutrino-electron elastic scattering. Observed spectra are determined by the solar-$ν_{e}$ survival probability $P_{ee}(E)$, and the chiral couplings of the neutrino and electron. Some theories of physics beyond the Standard Model postulate the existence of Non-Standard Interactions (NSI's) which modify the chiral couplings and $P_{ee}(E)$. In this paper, we search for such NSI's, in particular, flavor-diagonal neutral current interactions that modify the $ν_e e$ and $ν_τe$ couplings using Borexino Phase II data. Standard Solar Model predictions of the solar neutrino fluxes for both high- and low-metallicity assumptions are considered. No indication of new physics is found at the level of sensitivity of the detector and constraints on the parameters of the NSI's are placed. In addition, with the same dataset the value of $\sin^2θ_W$ is obtained with a precision comparable to that achieved in reactor antineutrino experiments.

preprint2020arXiv

Deep-sea deployment of the KM3NeT neutrino telescope detection units by self-unrolling

KM3NeT is a research infrastructure being installed in the deep Mediterranean Sea. It will house a neutrino telescope comprising hundreds of networked moorings - detection units or strings equipped with optical instrumentation to detect the Cherenkov radiation generated by charged particles from neutrino-induced collisions in its vicinity. In comparison to moorings typically used for oceanography, several key features of the KM3NeT string are different: the instrumentation is contained in transparent and thus unprotected glass spheres; two thin Dyneema ropes are used as strength members; and a thin delicate backbone tube with fibre-optics and copper wires for data and power transmission, respectively, runs along the full length of the mooring. Also, compared to other neutrino telescopes such as ANTARES in the Mediterranean Sea and GVD in Lake Baikal, the KM3NeT strings are more slender to minimise the amount of material used for support of the optical sensors. Moreover, the rate of deploying a large number of strings in a period of a few years is unprecedented. For all these reasons, for the installation of the KM3NeT strings, a custom-made, fast deployment method was designed. Despite the length of several hundreds of metres, the slim design of the string allows it to be compacted into a small, re-usable spherical launching vehicle instead of deploying the mooring weight down from a surface vessel. After being lowered to the seafloor, the string unfurls to its full length with the buoyant launching vehicle rolling along the two ropes.The design of the vehicle, the loading with a string, and its underwater self-unrolling are detailed in this paper.

preprint2020arXiv

Improved measurement of $^8$B solar neutrinos with 1.5 kt y of Borexino exposure

We report on an improved measurement of the $^8$B solar neutrino interaction rate with the Borexino experiment at the Laboratori Nazionali del Gran Sasso. Neutrinos are detected via their elastic scattering on electrons in a large volume of liquid scintillator. The measured rate of scattered electrons above 3 MeV of energy is $0.223\substack{+0.015 \\ -0.016}\,(stat)\,\substack{+0.006 \\ -0.006}\,(syst)$ cpd/100 t, which corresponds to an observed solar neutrino flux assuming no neutrino flavor conversion of $Φ\substack{\rm ES \\ ^8\rm B}=2.57\substack{+0.17 \\ -0.18}(stat)\substack{+0.07\\ -0.07}(syst)\times$10$^6$ cm$^{-2}\,$s$^{-1}$. This measurement exploits the active volume of the detector in almost its entirety for the first time, and takes advantage of a reduced radioactive background following the 2011 scintillator purification campaign and of novel analysis tools providing a more precise modeling of the background. Additionally, we set a new limit on the interaction rate of solar $hep$ neutrinos, searched via their elastic scattering on electrons as well as their neutral current-mediated inelastic scattering on carbon, $^{12}$C($ν,ν'$)$^{12}$C* ($E_γ$= 15.1 MeV).

preprint2020arXiv

Setup commissioning for an improved measurement of the D(p,gamma)3He cross section at Big Bang Nucleosynthesis energies

Among the reactions involved in the production and destruction of deuterium during Big Bang Nucleosynthesis, the deuterium-burning D(p,gamma)3He reaction has the largest uncertainty and limits the precision of theoretical estimates of primordial deuterium abundance. Here we report the results of a careful commissioning of the experimental setup used to measure the cross-section of the D(p,gamma)3He reaction at the Laboratory for Underground Nuclear Astrophysics of the Gran Sasso Laboratory (Italy). The commissioning was aimed at minimising all sources of systematic uncertainty in the measured cross sections. The overall systematic error achieved (< 3 %) will enable improved predictions of BBN deuterium abundance.

preprint2018arXiv

Speeding up complex multivariate data analysis in Borexino with parallel computing based on Graphics Processing Unit

A spectral fitter based on the graphics processor unit (GPU) has been developed for Borexino solar neutrino analysis. It is able to shorten the fitting time to a superior level compared to the CPU fitting procedure. In Borexino solar neutrino spectral analysis, fitting usually requires around one hour to converge since it includes time-consuming convolutions in order to account for the detector response and pile-up effects. Moreover, the convergence time increases to more than two days when including extra computations for the discrimination of $^{11}$C and external $γ$s. In sharp contrast, with the GPU-based fitter it takes less than 10 seconds and less than four minutes, respectively. This fitter is developed utilizing the GooFit project with customized likelihoods, pdfs and infrastructures supporting certain analysis methods. In this proceeding the design of the package, developed features and the comparison with the original CPU fitter are presented.