Search

Context. The formation of stars has been subject to extensive studies in the past decades from molecular cloud to protoplanetary disk scales. It is still not fully understood how the surrounding material in a protostellar system, that often shows asymmetric structures with complex kinematic properties, feeds the central protostar(s) and their disk(s). Aims. We study the spatial morphology and kinematic properties of the molecular gas surrounding the IRS3A and IRS3B protostellar systems in the L1448N region located in the Perseus molecular cloud. Methods. We present 1 mm NOEMA observations of the PRODIGE large program and analyze the kinematic properties of molecular lines. Given the complexity of the spectral profiles, the lines are fitted with up to three Gaussian velocity components. The clustering algorithm DBSCAN is used to disentangle the velocity components into the underlying physical structure. Results. We discover an extended gas bridge (~3000 au) surrounding both the IRS3A and IRS3B systems in six molecular line tracers (C18O, SO, DCN, H2CO, HC3N, and CH3OH). This gas bridge is oriented along the northeast-southwest direction and shows clear velocity gradients on the order of 100 km/s/pc towards the IRS3A system. We find that the observed velocity profile is consistent with analytical streamline models of gravitational infall towards IRS3A. The high-velocity C18O (2-1) emission towards IRS3A indicates a protostellar mass of ~1.2 Msun. Conclusions. While high angular resolution continuum data often show IRS3A and IRS3B in isolation, molecular gas observations reveal that these systems are still embedded within a large-scale mass reservoir with a complex spatial morphology as well as velocity profiles. The kinematic properties of the extended gas bridge are consistent with gravitational infall toward the IRS3A protostar., Comment: 16 pages, 10 figures; accepted for publication in A&A

{Sulfur Monoxide is known to be a good shock tracer in molecular clouds and protostar environments, but its abundance is difficult to reproduce even using state-of-the-art astrochemical models. } {We investigate the properties of the observed SO emission in the protoplanetary disk of AB Auriga, a Herbig Ae star of 2.4 $\Msun$ located at 156 pc. The AB Aur system is unique because it exhibits a dust trap and at least one young putative planet orbiting at about 30 au from the central star.} {We reduced ALMA archival data (projects 2019.1.00579.S, 2021.1.00690.S and 2021.1.01216.S) and analyzed the three detected SO lines (SO $6_5-5_4$, $6_7-5_6$ and $5_6-4_5$). We also present C$^{17}$O and C$^{18}$O 2-1 data to complement the interpretation of the SO data.} {For the three SO lines, the maximum SO emission in the ring is not located in the dust trap. Moreover, the inner radius of the SO ring is significantly larger than the CO emission inner radius, $\sim 160$ au versus $\sim 90$ au. The SO emission traces gas located in part beyond the dust ring. This emission likely originates from shocks at the interface of the outer spirals, observed in CO and scattered light emission, and the molecular and dust ring. SO is also detected inside the cavity, at a radius $\sim 20-30$ au and with a rotation velocity compatible with the proto-planet P1. We speculate that this SO emission originates from accretion shocks onto the circumplanetary disk of the putative proto-planet P1.} {These observations confirm that SO is a good tracer of shocks in protoplanetary disks and could be a powerful, new tool to detect embedded (proto-)planets.} Conclusions. These observations confirm that SO is a good tracer of shocks in protoplanetary disks and could be a powerful, new tool to detect embedded (proto-)planets.

DarkSide-20k (DS-20k) is a dark matter detection experiment under construction at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy. It utilises ~100 t of low radioactivity argon from an underground source (UAr) in its inner detector, with half serving as target in a dual-phase time projection chamber (TPC). The UAr cryogenics system must maintain stable thermodynamic conditions throughout the experiment's lifetime of >10 years. Continuous removal of impurities and radon from the UAr is essential for maximising signal yield and mitigating background. We are developing an efficient and powerful cryogenics system with a gas purification loop with a target circulation rate of 1000 slpm. Central to its design is a condenser operated with liquid nitrogen which is paired with a gas heat exchanger cascade, delivering a combined cooling power of >8 kW. Here we present the design choices in view of the DS-20k requirements, in particular the condenser's working principle and the cooling control, and we show test results obtained with a dedicated benchmarking platform at CERN and LNGS. We find that the thermal efficiency of the recirculation loop, defined in terms of nitrogen consumption per argon flow rate, is 95 % and the pressure in the test cryostat can be maintained within $\pm$(0.1-0.2) mbar. We further detail a 5-day cool-down procedure of the test cryostat, maintaining a cooling rate typically within -2 K/h, as required for the DS-20k inner detector. Additionally, we assess the circuit's flow resistance, and the heat transfer capabilities of two heat exchanger geometries for argon phase change, used to provide gas for recirculation. We conclude by discussing how our findings influence the finalisation of the system design, including necessary modifications to meet requirements and ongoing testing activities., Comment: 45 pages, 24 figures

The dual-phase liquid argon time projection chamber is presently one of the leading technologies to search for dark matter particles with masses below 10 GeV/c$^2$. This was demonstrated by the DarkSide-50 experiment with approximately 50 kg of low-radioactivity liquid argon as target material. The next generation experiment DarkSide-20k, currently under construction, will use 1,000 times more argon and is expected to start operation in 2027. Based on the DarkSide-50 experience, here we assess the DarkSide-20k sensitivity to models predicting light dark matter particles, including Weakly Interacting Massive Particles (WIMPs) and sub-GeV/c$^2$ particles interacting with electrons in argon atoms. With one year of data, a sensitivity improvement to dark matter interaction cross-sections by at least one order of magnitude with respect to DarkSide-50 is expected for all these models. A sensitivity to WIMP--nucleon interaction cross-sections below $1\times10^{-42}$ cm$^2$ is achievable for WIMP masses above 800 MeV/c$^2$. With 10 years exposure, the neutrino fog can be reached for WIMP masses around 5 GeV/c$^2$., Comment: submitted to Nature Communications

We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2-1) isotopologue emission observed at 0.9 with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution PIs: P. Caselli & Th. Henning). Our sample consists of Class II disks with no evidence of strong radial substructures. We use thesedata to constrain the thermal and chemical structure of these disks through theoretical models for gas emission. To fit the combined optically thick and thin CO line data in Fourier space, we developed the DiskCheF code, which includes the parameterized disk physical structure, machine-learning (ML) accelerated chemistry, and the RADMC-3D line radiative transfer module. A key novelty of DiskCheF is the fast and feasible ML-based chemistry trained on the extended grid of the disk physical-chemical models precomputed with the ANDES2 code. This ML approach allows complex chemical kinetics models to be included in a time-consuming disk fitting without the need to run a chemical code. We present a novel approach to incorporate chemistry into disk modeling without the need to explicitly calculate a chemical network every time. Using this new disk modeling tool, we successfully fit the 12CO, 13CO, and C18O (2-1) data from the CI, CY, DL, DM, DN, and IQ Tau disks. The combination of optically thin and optically thick CO lines allows us to simultaneously constrain the disk temperature and mass distribution, and derive the CO-based gas masses. These values are in reasonable agreement with the disk dust masses rescaled by a factor of 100 as well as with other indirect gas measurements., Comment: 25 pages, 19 figures

Deuterated molecules and their molecular D/H-ratios (RD(D)) are important diagnostic tools to study the physical conditions of star-forming regions. The degree of deuteration, RD(D), can be significantly enhanced over the elemental D/H-ratio depending on physical parameters. Within the Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE), we aim to explore the large-scale distribution of deuterated molecules in the nearby Cygnus-X region. We focus on the analysis of large-scale structures of deuterated molecules in the filamentary region hosting the prominent Hii region DR21 and DR21(OH). Here we discuss the HCO+, HNC and HCN molecules and their deuterated isotopologues DCO+, DNC and DCN. The spatial distributions of integrated line emissions from DCO+, DNC, and DCN reveal morphological differences. DCO+ displays the most extended emission, characterized by several prominent peaks. Likewise, DNC exhibits multiple peaks, although its emission appears less extended compared to DCO+. In contrast to the extended emission of DCO+ and DNC, DCN appears the least extended, with distinct peaks. Focusing only on the regions where all three molecules are observed, the mean deuteration ratios for each species are 0.01 for both DNC and DCN, and = 0.005 for DCO+. Anti-correlations are found with deuterated molecules and dust temperature or N(H2). The strongest anti-correlation is found with RD(DCO+) and N(H2). The anti-correlation of RD(DCO+) and N(H2) is suggested to be a result of a combination of an increased photodissociation degree and shocks. A strong positive correlation between the ratio of integrated intensities of DCN and DNC with their 13C-isotopologues, are found in high column density regions. The positive relationship between the ratios implies that the D-isotopologue of the isomers could potentially serve as a tracer for the kinetic gas temperature., Comment: 24 pages, 21 figures, accepted to A&A

Experiments aimed at direct searches for WIMP dark matter require highly effective reduction of backgrounds and control of any residual radioactive contamination. In particular, neutrons interacting with atomic nuclei represent an important class of backgrounds due to the expected similarity of a WIMP-nucleon interaction, so that such experiments often feature a dedicated neutron detector surrounding the active target volume. In the context of the development of DarkSide-20k detector at INFN Gran Sasso National Laboratory (LNGS), several R&D projects were conceived and developed for the creation of a new hybrid material rich in both hydrogen and gadolinium nuclei to be employed as an essential element of the neutron detector. Thanks to its very high cross-section for neutron capture, gadolinium is one of the most widely used elements in neutron detectors, while the hydrogen-rich material is instrumental in efficiently moderating the neutrons. In this paper results from one of the R&Ds are presented. In this effort the new hybrid material was obtained as a poly(methyl methacrylate) (PMMA) matrix, loaded with gadolinium oxide in the form of nanoparticles. We describe its realization, including all phases of design, purification, construction, characterization, and determination of mechanical properties of the new material.

Complex Organic Molecules (COMs) have been found toward low-mass protostars but the origins of the COM emission are still unclear. It can be associated with, for example, hot corinos, outflows, and/or accretion shock/disk atmosphere. We have conducted NOEMA observations toward SVS13A from the PROtostars & DIsks: Global Evolution (PRODIGE) program. Our previous \ce{DCN} observations reveal a possible infalling streamer, which may affect the chemistry of the central protobinary by inducing accretion outbursts and/or shocked gas. Here, we further analyze six O-bearing COMs: CH3OH, aGg'-(CH2OH)2, C2H5OH, CH2(OH)CHO, CH3CHO, and CH3OCHO. Although the COM emission is not spatially resolved, we constrain the source sizes to $\lesssim0.3-0.4$ arcsec (90$-$120 au) by conducting uv-domain Gaussian fitting. Interestingly, the high-spectral resolution data reveal complex line profiles with multiple peaks showing differences between these six O-bearing COMs. The LTE fitting unveils differences in excitation temperatures and emitting areas among these COMs. We further conduct multiple-velocity-component LTE fitting to decompose the line emission into different kinematic components. Up to 6 velocity components are found from the LTE modeling. The temperature, column density, and source size of these components from each COM are obtained. We find a variety in excitation temperatures ($100-500$ K) and source sizes (D$\sim10-70$ au) from these kinematic components from different COMs. The emission of each COM can trace several components and different COMs most likely trace different regions. Given this complex structure, we suggest that the central region is inhomogeneous and unlikely to be heated by only protostellar radiation. We conclude that accretion shocks induced by the large-scale infalling streamer likely exist and contribute to the complexity of the COM emission., Comment: 29 pages, 18 figures, accepted to A&A

We are performing a line survey of 8 planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from >20 C-, O,- N-, and S-bearing species. The observations in four spectral settings at 210-280 GHz with $1\sigma$ rms sensitivity of ~ 8-12 mJy/beam at 0.9" and 0.3 km/s resolution will be completed in 2024. The uv-visibilities are fitted with the DiskFit model to obtain key stellar and disk properties. In this paper, the combined $^{12}$CO, $^{13}$CO and C$^{18}$O $J = 2-1$ data are presented. We find that the CO fluxes and disk masses inferred from dust continuum tentatively correlate with the CO emission sizes. We constrain dynamical stellar masses, geometries, temperatures, the CO column densities and gas masses for each disk. The best-fit temperatures at 100 au are ~ 17-37 K, and decrease radially with the power-law exponent q ~ 0.05-0.76. The inferred CO column densities decrease radially with the power-law exponent p ~ 0.2-3.1. The gas masses estimated from $^{13}$CO (2-1) are ~ $0.001-0.2 M_\textrm{Sun}$. The best-fit CO column densities point to severe CO freeze-out in the disks. The DL Tau disk is an outlier, and has either stronger CO depletion or lower gas mass than the rest of the sample. The CO isotopologue ratios are roughly consistent with the observed values in disks and the low-mass star-forming regions., Comment: 47 pages, 10 figures, 6 tables, accepted for publication by A&A on 22/02/2024

The stability of a dark matter detector on the timescale of a few years is a key requirement due to the large exposure needed to achieve a competitive sensitivity. It is especially crucial to enable the detector to potentially detect any annual event rate modulation, an expected dark matter signature. In this work, we present the performance history of the DarkSide-50 dual-phase argon time projection chamber over its almost three-year low-radioactivity argon run. In particular, we focus on the electroluminescence signal that enables sensitivity to sub-keV energy depositions. The stability of the electroluminescence yield is found to be better than 0.5%. Finally, we show the temporal evolution of the observed event rate around the sub-keV region being consistent to the background prediction., Comment: 13 pages, 5 figures

Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions., Comment: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

This work is devoted to the problem of synchronization of two Morris-Lecar neuron models. The Morris-Lecar model is a second-order system of differential equations, which describes an uneasy relationship between the membrane potential and the activation of ion channels inside the membrane. Synchronization, which is a state of a network when models begin to act similarly in some sense, of such models is interesting not only from a mathematical point of view, but also from a biological one, because synchronous activity plays a very important role in brain functioning. The speed gradient algorithm, which is a continuous version of gradient algorithms, was applied to solve this problem. The algorithm of coupling strength control was obtained. It ensures the achievement of the control goal. The MATLAB modelling demonstrated the correctness of the obtained result and fast convergence rate of corresponding models' variables to each other., Comment: 5 pages, 6 figures, in Russian, supported by grant SPbU ID 84912397

Borexino could efficiently distinguish between alpha and beta radiation in its liquid scintillator by the characteristic time profile of their scintillation pulse. This alpha/beta discrimination, first demonstrated at the tonne scale in the Counting Test Facility prototype, was used throughout the lifetime of the experiment between 2007 and 2021. With this method, alpha events are identified and subtracted from the beta-like solar neutrino events. This is particularly important in liquid scintillator as alpha scintillation is quenched many-fold. In Borexino, the prominent Po-210 decay peak was a background in the energy range of electrons scattered from Be-7 solar neutrinos. Optimal alpha-beta discrimination was achieved with a "multi-layer perceptron neural network", which its higher ability to leverage the timing information of the scintillation photons detected by the photomultiplier tubes. An event-by-event, high efficiency, stable, and uniform pulse shape discrimination was essential in characterising the spatial distribution of background in the detector. This benefited most Borexino measurements, including solar neutrinos in the \pp chain and the first direct observation of the CNO cycle in the Sun. This paper presents the key milestones in alpha/beta discrimination in Borexino as a term of comparison for current and future large liquid scintillator detectors, Comment: 13 pages, 14 figures

We use spectra and maps of the $J=1-0$ and $J=2-1$ DCO$^{+}$, DCN, DNC, $\rm N_2D^+$ lines and $1_{11}-1_{01}$ ortho- and para-NH$_{2}$D lines, obtained with the IRAM-30m telescope, as well as observations of their hydrogenated isotopologues to study deuteration processes in five high-mass star-forming regions. The temperature was estimated from CH$_3$CCH lines, also observed with the IRAM-30m telescope, and from NH$_3$ lines, observed with the 100-m radio telescope in Effelsberg, as well as using the integrated intensity ratios of the $J=1-0$ H$^{13}$CN and HN$^{13}$C lines and their main isotopologues. Applying a non-LTE radiative transfer model with RADEX, the gas density and the molecular column densities were estimated. D/H ratios are $0.001-0.05$ for DCO$^{+}$, $0.001-0.02$ for DCN, $0.001-0.05$ for DNC and $0.02-0.4$ for NH$_{2}$D. The D/H ratios decrease with increasing temperature in the range of $\rm 20-40 \,K$ and slightly vary at densities $n(\rm H_2) \sim 10^4-10^6\, cm^{-3}$. The deuterium fraction of $\rm N_2H^{+}$ is $0.008-0.1$ at temperatures in the range of $\rm 20-25\, K$ and at a density of $\sim 10^5\, \rm cm^{-3}$. We also estimate relative abundances and find $ \sim 10^{-11}-10^{-9}$ for DCO$^{+}$ and DNC, $ \sim 10^{-11}-10^{-10}$ for $\rm N_2D^+$ and $ \sim 10^{-10}-10^{-8}$ for NH$_{2}$D. The relative abundances of these species decrease with increasing temperature. However, the DCN/H$_2$ ratio is almost constant ($\sim 10^{-10}$). The observational results agree with the predictions of chemical models (although in some cases there are significant differences)., Comment: 13 pages + appendix (25 pages total), 11 figures. Accepted for publication in MNRAS

Molecular outflows are believed to be a key ingredient in the process of star formation. The molecular outflow associated with DR21 Main in Cygnus-X is one of the most extreme, in mass and size, molecular outflows in the Milky Way. The outflow is suggested to belong to a rare class of explosive outflows which are formed by the disintegration of protostellar systems.We aim to explore the morphology, kinematics,and energetics of the DR21 Main outflow, and compare those properties to confirmed explosive outflows to unravel the underlying driving mechanism behind DR21. Line and continuum emission are studied at a wavelength of 3.6\,mm with IRAM 30 m and NOEMA telescopes as part of the Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE) program. The spectra include ($J= 1-0$) transitions of HCO$^+$, HCN, HNC, N$_2$H$^+$, H$_2$CO, CCH tracing different temperature and density regimes of the outflowing gas at high-velocity resolution ($\sim$ 0.8 km s$^{-1}$). The map encompasses the entire DR21 Main outflow and covers all spatial scales down to a resolution of ~3" ($\sim$ 0.02 pc). Integrated intensity maps of the HCO$^+$ emission reveal a strongly collimated bipolar outflow with significant overlap of the blue- and red-shifted emission. The opening angles of both outflow lobes decrease with velocity, from $\sim80$ to 20$^{\circ}$ for the velocity range from 5 to 45 km s$^{-1}$ relative to the source velocity. No evidence is found for the presence of elongated, "filament-like" structures expected in explosive outflows. N$_2$H$^+$ emission near the western outflow lobe reveals the presence of a dense molecular structure which appears to be interacting with the DR21 Main outflow. The overall morphology as well as the detailed kinematics of the DR21 Main outflow is more consistent with that of a typical bipolar outflow instead of an explosive counterpart., Comment: 26 pages, 37 figures. Accepted for publication in A&A

The disk-outflow connection plays a key role in extracting excess angular momentum from a forming protostar. We have previously reported the discovery of a small molecular outflow from the edge-on T Tauri star in the Bok globule CB26 that shows a peculiar velocity pattern, reminiscent of an outflow that corotates with the disk. We report new, high-resolution mm-interferometric observations of CB26 with the aim of revealing the morphology and kinematics of the outflow at the disk-outflow interface. The IRAM PdBI was used to observe CO(2-1) at 1.3mm with a resolution of 0.5". Using a physical model of the disk, which was derived from the dust emission, we employed chemo-dynamical modeling combined with line radiative transfer to constrain kinematic parameters and to construct a model of the CO emission from the disk that allowed us to separate the emission of the disk from that of the outflow. Our observations confirm the disk-wind nature of the rotating molecular outflow from CB26. The new high-resolution data reveal an X-shaped morphology of the CO emission close to the disk, and vertical streaks extending from the disk surface with a small half-opening angle of ~7deg, which can be traced out to vertical heights of ~500au. We interpret this emission as the combination of the disk atmosphere and a well-collimated disk wind, which we trace down to vertical heights of 40au, where it is launched from the surface of the flared disk at radii of 20-45au. The observed CO outflow has a total momentum flux of 1e-5 Msun km/s/yr, which is nearly three orders of magnitude larger than the maximum thrust that can be provided by the luminosity of the central star. We conclude that photoevaporation cannot be the main driving mechanism for this outflow, but it must be predominantly an MHD disk wind. It is thus far the best-resolved rotating disk wind observed to be launched from a circumstellar disk., Comment: Accepted by A&A, 25 pages, 19 figures

We report the first measurement of CNO solar neutrinos by Borexino that uses the Correlated Integrated Directionality (CID) method, exploiting the sub-dominant Cherenkov light in the liquid scintillator detector. The directional information of the solar origin of the neutrinos is preserved by the fast Cherenkov photons from the neutrino scattered electrons, and is used to discriminate between signal and background. The directional information is independent from the spectral information on which the previous CNO solar neutrino measurements by Borexino were based. While the CNO spectral analysis could only be applied on the Phase-III dataset, the directional analysis can use the complete Borexino data taking period from 2007 to 2021. The absence of CNO neutrinos has been rejected with >5{\sigma} credible level using the Bayesian statistics. The directional CNO measurement is obtained without an external constraint on the $^{210}$Bi contamination of the liquid scintillator, which was applied in the spectral analysis approach. The final and the most precise CNO measurement of Borexino is then obtained by combining the new CID-based CNO result with an improved spectral fit of the Phase-III dataset. Including the statistical and the systematic errors, the extracted CNO interaction rate is $R(\mathrm{CNO})=6.7^{+1.2}_{-0.8} \, \mathrm{cpd/100 \, tonnes}$. Taking into account the neutrino flavor conversion, the resulting CNO neutrino flux at Earth is $\Phi_\mathrm{CNO}=6.7 ^{+1.2}_{-0.8} \times 10^8 \, \mathrm{cm^{-2} s^{-1}}$, in agreement with the high metallicity Standard Solar Models. The results described in this work reinforce the role of the event directional information in large-scale liquid scintillator detectors and open up new avenues for the next-generation liquid scintillator or hybrid neutrino experiments.

Dark matter induced event rate in an Earth-based detector is predicted to show an annual modulation as a result of the Earth's orbital motion around the Sun. We searched for this modulation signature using the ionization signal of the DarkSide-50 liquid argon time projection chamber. No significant signature compatible with dark matter is observed in the electron recoil equivalent energy range above $40~{\rm eV_{ee}}$, the lowest threshold ever achieved in such a search., Comment: 8 pages, 4 figures

The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars., Comment: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

The search for neutrino events in correlation with gravitational wave (GW) events for three observing runs (O1, O2 and O3) from 09/2015 to 03/2020 has been performed using the Borexino data-set of the same period. We have searched for signals of neutrino-electron scattering with visible energies above 250 keV within a time window of 1000 s centered at the detection moment of a particular GW event. The search was done with three visible energy thresholds of 0.25, 0.8 and 3.0 MeV.Two types of incoming neutrino spectra were considered: the mono-energetic line and the spectrum expected from supernovae. The same spectra were considered for electron antineutrinos detected through inverse beta-decay (IBD) reaction. GW candidates originated by merging binaries of black holes (BHBH), neutron stars (NSNS) and neutron star and black hole (NSBH) were analysed separately. Additionally, the subset of most intensive BHBH mergers at closer distances and with larger radiative mass than the rest was considered. In total, follow-ups of 74 out of 93 gravitational waves reported in the GWTC-3 catalog were analyzed and no statistically significant excess over the background was observed. As a result, the strongest upper limits on GW-associated neutrino and antineutrino fluences for all flavors (\nu_e, \nu_\mu, \nu_\tau) have been obtained in the (0.5 - 5.0) MeV neutrino energy range., Comment: 13 pages, 8 figures

We obtained deep images of G 196-3 B and VHS J1256-1257 b with the NOrthern Extended Millimeter Array (NOEMA) at 1.3 mm. These data were combined with recently published Atacama Large Millimeter Array (ALMA) and Very Large Array (VLA) data of VHS J1256-1257 b at 0.87 mm and 0.9 cm, respectively. Neither G 196-3 B nor VHS J1256-1257 b were detected in the NOEMA, ALMA and VLA data. At 1.3 mm, we imposed flux upper limits of 0.108 mJy (G 196-3 B) and 0.153 mJy (VHS J1256-1257 b) with a 3-sigma confidence. Using the flux upper limits at the millimeter and radio wavelength regimes, we derived maximum values of 0.016 M$_{\rm Earth}$ and 0.004 M$_{\rm Earth}$ for the mass of any cold dust that might be surrounding G 196-3 B and VHS J1256-1257 b, respectively. We put our results in the context of other deep millimeter observations of free-floating and companion objects with substellar masses smaller than 20 M$_{\rm Jupiter}$ and ages between 1 and a few hundred million years. Only two very young objects are detected out of a few tens concluding, as other groups did before, that the disks around these very low-mass objects must have small masses and possibly reduced sizes. If debris disks around substellar objects scale down in a similar manner as protoplanetary disks do, millimeter observations of moderately young brown dwarfs and planets must be at least two orders of magnitude deeper for being able to detect and characterize their surrounding debris disks., Comment: Accepted fro publication in A&A

We present a novel approach for the search of dark matter in the DarkSide-50 experiment, relying on Bayesian Networks. This method incorporates the detector response model into the likelihood function, explicitly maintaining the connection with the quantity of interest. No assumptions about the linearity of the problem or the shape of the probability distribution functions are required, and there is no need to morph signal and background spectra as a function of nuisance parameters. By expressing the problem in terms of Bayesian Networks, we have developed an inference algorithm based on a Markov Chain Monte Carlo to calculate the posterior probability. A clever description of the detector response model in terms of parametric matrices allows us to study the impact of systematic variations of any parameter on the final results. Our approach not only provides the desired information on the parameter of interest, but also potential constraints on the response model. Our results are consistent with recent published analyses and further refine the parameters of the detector response model., Comment: 24 pages, 12 figures, 1 table

Influence of the prolonged irradiation by fission products of 252Cf radionuclide on the operational parameters of silicon-lithium Si(Li) p-i-n detectors, Si surface barrier detectors and Si planar p+n detector was investigated. The obtained results revealed a linear shift of the fission fragment peaks positions towards the lower energies with increase of the irradiation dose for all investigated detectors. The rate of the peaks shift was found to depend strongly on the detector type and the strength of the electric field in the detectors active region, but not on the temperature of irradiation (room or liquid nitrogen temperature). Based on the obtained results, the possibility of integration of the investigated types of Si semiconductor detectors in a radionuclide neutron calibration source is considered., Comment: 6 pages, 2 figures, 1 table

Expansion of wind-blown bubbles or HII regions lead to formation of shocks in the interstellar medium, which compress surrounding gas into dense layers. We made spatially and velocity-resolved observations of the RCW~120 PDR and nearby molecular gas with CO(6-5) and 13CO(6-5) lines and distinguished a bright CO-emitting layer, which we related with the dense shocked molecular gas moving away from the ionizing star due to expansion of HII region. Simulating gas density and temperature, as well as brightness of several CO and C+ emission lines from the PDR, we found reasonable agreement with the observed values. Analysing gas kinematics, we revealed the large-scale shocked PDR and also several dense environments of embedded protostars and outflows. We observe the shocked layer as the most regular structure in the CO(6-5) map and in the velocity space, when the gas around YSOs is dispersed by the outflows., Comment: accepted for publication in MNRAS

Aims. We present high sensitivity and high-spectral resolution NOEMA observations of the Class 0/I binary system SVS13A, composed of the low-mass protostars VLA4A and VLA4B with a separation of ~90 au. VLA4A is undergoing an accretion burst that enriches the chemistry of the surrounding gas. This gives us an excellent opportunity to probe the chemical and physical conditions as well as the accretion process. Methods. We observe the (12K-11K) lines of CH3CN and CH313CN, the DCN (3-2) line, and the C18O (2-1) line toward SVS13A using NOEMA. Results. We find complex line profiles at disk scales which cannot be explained by a single component or pure Keplerian motion. By adopting two velocity components to model the complex line profiles, we find that the temperatures and densities are significantly different between these two components. This suggests that the physical conditions of the emitting gas traced via CH3CN can change dramatically within the circumbinary disk. In addition, combining our observations of DCN (3-2) with previous ALMA high-angular-resolution observations, we find that the binary system (or VLA4A) might be fed by an infalling streamer from envelope scales (~700 au). If this is the case, this streamer contributes to the accretion of material onto the system with a rate of at least 1.4x10-6 Msun yr-1. Conclusions. We conclude that the CH3CN emission in SVS13A traces hot gas from a complex structure. This complexity might be affected by a streamer that is possibly infalling and funneling material into the central region., Comment: 20 pages, 19 figures, accepted to A&A

We present a description of the originally developed $\beta$-spectrometer consisting of two Si(Li)-detectors with sensitive area thickness above 8~mm and $4\pi$-geometry. The full absorption spectrometer allows for direct measurements of $\beta$-spectra, disregarding the corrections to the response function induced by the electron backscattering from the crystal surface. In case of $\beta$-spectra of transitions to the excited state of the daughter isotope additional 3" BGO-detector is used in order to detect the $\gamma$-quanta in coincidence with the pair of Si(Li)-spectrometers., Comment: 9 pages, 7 figures

Dark matter lighter than 10 GeV/c$^2$ encompasses a promising range of candidates. A conceptual design for a new detector, DarkSide-LowMass, is presented, based on the DarkSide-50 detector and progress toward DarkSide-20k, optimized for a low-threshold electron-counting measurement. Sensitivity to light dark matter is explored for various potential energy thresholds and background rates. These studies show that DarkSide-LowMass can achieve sensitivity to light dark matter down to the solar neutrino floor for GeV-scale masses and significant sensitivity down to 10 MeV/c$^2$ considering the Migdal effect or interactions with electrons. Requirements for optimizing the detector's sensitivity are explored, as are potential sensitivity gains from modeling and mitigating spurious electron backgrounds that may dominate the signal at the lowest energies.

Context. In the past few years, there has been a rise in the detection of streamers, asymmetric flows of material directed toward the protostellar disk with material from outside the star's natal core. It is unclear how they affect the process of mass accretion, in particular beyond the Class 0 phase. Aims. We investigate the gas kinematics around Per-emb-50, a Class I source in the crowded star-forming region NGC 1333. Our goal is to study how the mass infall proceeds from envelope to disk scales in this source. Results. We discover a streamer delivering material toward Per-emb-50 in H$_2$CO and C$^{18}$O emission. The streamer's emission can be well described by the analytic solutions for an infalling parcel of gas along a streamline with conserved angular momentum, both in the image plane and along the line of sight velocities. The streamer has a mean infall rate of $1.3 \times 10^{ -6}$ M$_{ \odot}$ yr$^{ -1}$, $5 -10$ times higher than the current accretion rate of the protostar. SO and SO$_2$ emission reveal asymmetric infall motions in the inner envelope, additional to the streamer around Per-emb-50. Furthermore, the presence of SO$_2$ could mark the impact zone of the infalling material. Conclusions. The streamer delivers sufficient mass to sustain the protostellar accretion rate and might produce an accretion burst, which would explain the protostar's high luminosity with respect to other Class I sources. Our results highlight the importance of late infall for protostellar evolution: streamers might provide a significant amount of mass for stellar accretion after the Class 0 phase., Comment: 20 pages, 14 figures, accepted for publication in A&A

Dark matter elastic scattering off nuclei can result in the excitation and ionization of the recoiling atom through the so-called Migdal effect. The energy deposition from the ionization electron adds to the energy deposited by the recoiling nuclear system and allows for the detection of interactions of sub-GeV/c$^2$ mass dark matter. We present new constraints for sub-GeV/c$^2$ dark matter using the dual-phase liquid argon time projection chamber of the DarkSide-50 experiment with an exposure of (12306 $\pm$ 184) kg d. The analysis is based on the ionization signal alone and significantly enhances the sensitivity of DarkSide-50, enabling sensitivity to dark matter with masses down to 40 MeV/c$^2$. Furthermore, it sets the most stringent upper limit on the spin independent dark matter nucleon cross section for masses below $3.6$ GeV/c$^2$., Comment: 7 pages, 3 figures

We present a search for dark matter particles with sub-GeV/$c^2$ masses whose interactions have final state electrons using the DarkSide-50 experiment's (12306 $\pm$ 184) kg d low-radioactivity liquid argon exposure. By analyzing the ionization signals, we exclude new parameter space for the dark matter-electron cross section $\bar{\sigma}_e$, the axioelectric coupling constant $g_{Ae}$, and the dark photon kinetic mixing parameter $\kappa$. We also set the first dark matter direct-detection constraints on the mixing angle $\left|U_{e4}\right|^2$ for keV sterile neutrinos., Comment: 6 pages, 2 figures

We report on the search for dark matter WIMPs in the mass range below 10 GeV/c$^2$, from the analysis of the entire dataset acquired with a low-radioactivity argon target by the DarkSide-50 experiment at LNGS. The new analysis benefits from more accurate calibration of the detector response, improved background model, and better determination of systematic uncertainties, allowing us to accurately model the background rate and spectra down to 0.06 keV$_{er}$. A 90% C.L. exclusion limit for the spin-independent cross section of 3 GeV/c$^2$ mass WIMP on nucleons is set at 6$\times$10$^{-43}$ cm$^2$, about a factor 10 better than the previous DarkSide-50 limit. This analysis extends the exclusion region for spin-independent dark matter interactions below the current experimental constraints in the $[1.2, 3.6]$ GeV/c$^2$ WIMP mass range., Comment: 11 pages, 12 figures

We present an improved measurement of the CNO solar neutrino interaction rate at Earth obtained with the complete Borexino Phase-III dataset. The measured rate R$_{\rm CNO}$ = $6.7^{+2.0}_{-0.8}$ counts/(day$ \cdot$ 100 tonnes), allows us to exclude the absence of the CNO signal with about 7$\sigma$ C.L. The correspondent CNO neutrino flux is $6.6^{+2.0}_{-0.9} \times 10^8$ cm$^{-2}$ s$^{-1}$, taking into account the neutrino flavor conversion. We use the new CNO measurement to evaluate the C and N abundances in the Sun with respect to the H abundance for the first time with solar neutrinos. Our result of $N_{\rm CN}$ = $(5.78^{+1.86}_{-1.00})\times10^{-4}$ displays a $\sim$2$\sigma$ tension with the "low metallicity" spectroscopic photospheric measurements. On the other hand, our result used together with the $^7$Be and $^8$B solar neutrino fluxes, also measured by Borexino, permits to disfavour at 3.1$\sigma$ C.L. the "low metallicity" SSM B16-AGSS09met as an alternative to the "high metallicity" SSM B16-GS98., Comment: 8 pages, 6 figures

Since the beginning of 2012, the Borexino collaboration has been reporting precision measurements of the solar neutrino fluxes, emitted in the proton-proton chain and in the Carbon-Nitrogen-Oxygen cycle. The experimental sensitivity achieved in Phase-II and Phase-III of the Borexino data taking made it possible to detect the annual modulation of the solar neutrino interaction rate due to the eccentricity of Earth's orbit, with a statistical significance greater than 5$\sigma$. This is the first precise measurement of the Earth's orbital parameters based solely on solar neutrinos and an additional signature of the solar origin of the Borexino signal. The complete periodogram of the time series of the Borexino solar neutrino detection rate is also reported, exploring frequencies between one cycle/year and one cycle/day. No other significant modulation frequencies are found. The present results were uniquely made possible by Borexino's decade-long high-precision solar neutrino detection., Comment: 12 pages, 10 figures

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$\sigma$ 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., Comment: 6 pages, 4 figures, short letter of arXiv:2109.04770

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., Comment: 14 pages, 10 figures

We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups., Comment: 34 pages, 17 figures, accepted for publication in A&A

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., Comment: 19 pages, 15 Figures, revised version after comments from PRD Referees, shorter letter submitted with the title: "First Directional Measurement of sub-MeV Solar Neutrinos with Borexino"

Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources., Comment: 16 pages, 11 figures, Accepted for publication in A&A

DarkSide-50 has demonstrated the high potential of dual-phase liquid argon time projection chambers in exploring interactions of WIMPs in the GeV/c$^2$ mass range. The technique, based on the detection of the ionization signal amplified via electroluminescence in the gas phase, allows to explore recoil energies down to the sub-keV range. We report here on the DarkSide-50 measurement of the ionization yield of electronic recoils down to $\sim$180~eV$_{er}$, exploiting $^{37}$Ar and $^{39}$Ar decays, and extrapolated to a few ionization electrons with the Thomas-Imel box model. Moreover, we present a model-dependent determination of the ionization response to nuclear recoils down to $\sim$500~eV$_{nr}$, the lowest ever achieved in liquid argon, using \textit{in situ} neutron calibration sources and external datasets from neutron beam experiments., Comment: 11 pages, 12 figures, 1 table

Finding unequivocal evidence of dark matter interactions in a particle detector is a major objective of physics research. Liquid argon time projection chambers offer a path to probe Weakly Interacting Massive Particles scattering cross sections on nucleus down to the so-called neutrino floor, in a mass range from few GeV's to hundredths of TeV's. Based on the successful operation of the DarkSide-50 detector at LNGS, a new and more sensitive experiment, DarkSide-20k, has been designed and is now under construction. A thorough understanding of the DarkSide-50 detector response and, therefore, of all kind of observed events, is essential for an optimal design of the new experiment. In this paper, we report on a particular set of events, which were not used for dark matter searches. Namely, standard two-pulse scintillation-ionization signals accompanied by a small amplitude third pulse, originating from single or few electrons, in a time window of less than a maximum drift time. We compare our findings to those of a recent paper of the LUX Collaboration (D.S.Akerib et al. Phys.Rev.D 102, 092004). Indeed, both experiments observe events related to photoionization of the cathode. From the measured rate of these events, we estimate for the first time the quantum efficiency of the tetraphenyl butadiene deposited on the DarkSide-50 cathode at wavelengths around 128 nm, in liquid argon. Also, both experiments observe events likely related to photoionization of impurities in the liquid. The probability of photoelectron emission per unit length turns out to be one order of magnitude smaller in DarkSide-50 than in LUX. This result, together with the much larger measured electron lifetime, coherently hints toward a lower concentration of contaminants in DarkSide-50 than in LUX., Comment: 10 pages, 8 figures, 1 table

Cosmogenic radio-nuclei are an important source of background for low-energy neutrino experiments. In Borexino, cosmogenic $^{11}$C decays outnumber solar $pep$ and CNO neutrino events by about ten to one. Highly efficient identification of this background is mandatory for these neutrino analyses. We present here the details of the most consolidated strategy, used throughout Borexino solar neutrino measurements. It hinges upon finding the space-time correlations between $^{11}$C decays, the preceding parent muons and the accompanying neutrons. This article describes the working principles and evaluates the performance of this Three-Fold Coincidence (TFC) technique in its two current implementations: a hard-cut and a likelihood-based approach. Both show stable performances throughout Borexino Phases II (2012-2016) and III (2016-2020) data sets, with a $^{11}$C tagging efficiency of $\sim$90 % and $\sim$63-66 % of the exposure surviving the tagging. We present also a novel technique that targets specifically $^{11}$C produced in high-multiplicity during major spallation events. Such $^{11}$C appear as a burst of events, whose space-time correlation can be exploited. Burst identification can be combined with the TFC to obtain about the same tagging efficiency of $\sim$90 % but with a higher fraction of the exposure surviving, in the range of $\sim$66-68 %., Comment: 21 pages, 14 figures (but 15 files, one figure being made of 2 images), 3 tables

Borexino is a liquid scintillator detector located at the Laboratori Nazionale del Gran Sasso, Italy with the main goal to measure solar neutrinos. The experiment recently provided the first direct experimental evidence of CNO-cycle neutrinos in the Sun, rejecting the no-CNO signal hypothesis with a significance greater than 5$\sigma$ at 99\%C.L. The intrinsic $^{210}$Bi is an important background for this analysis due to its similar spectral shape to that of CNO neutrinos. $^{210}$Bi can be measured through its daughter $^{210}$Po which can be distinguished through an event-by-event basis via pulse shape discrimination. However, this required reducing the convective motions in the scintillator that brought additional $^{210}$Po from peripheral sources. This was made possible through the thermal insulation and stabilization campaign performed between 2015 and 2016. This article will explain the strategy and the different methods performed to extract the $^{210}$Bi upper limit in Phase-III (Jul 2016- Feb 2020) of the experiment through the analysis of $^{210}$Po in the cleanest region of the detector called the Low Polonium Field., Comment: Contribution to the 2021 Neutrinos session of the 55th Rencontres de Moriond