Search

Silicon photonics is a leading platform for realizing the vast numbers of physical qubits needed for useful quantum information processing because it leverages mature complementary metal-oxide-semiconductor (CMOS) manufacturing to integrate on-chip thousands of optical devices for generating and manipulating quantum states of light. A challenge to the practical operation and scale-up of silicon quantum-photonic integrated circuits, however, is the need to control their extreme sensitivity to process and temperature variations, free-carrier and self-heating nonlinearities, and thermal crosstalk. To date these challenges have been partially addressed using bulky off-chip electronics, sacrificing many benefits of a chip-scale platform. Here, we demonstrate the first electronic-photonic quantum system-on-chip (EPQSoC) consisting of quantum-correlated photon-pair sources stabilized via on-chip feedback control circuits, all fabricated in a high-volume 45nm CMOS microelectronics foundry. We use non-invasive photocurrent sensing in a tunable microring cavity photon-pair source to actively lock it to a fixed pump laser while operating in the quantum regime, enabling large scale microring-based quantum systems. In this first demonstration of such a capability, we achieve a high CAR of 134 with an ultra-low g(2)(0) of 0.021 at 2.2 kHz off-chip detected pair rate and 3.3 MHz/mW2 on-chip pair generation efficiency, and over 100 kHz off-chip detected pair rate at higher pump powers (1.5 MHz on-chip). These sources maintain stable quantum properties in the presence of temperature variations, operating reliably in practical settings with many adjacent devices creating thermal disturbances on the same chip. Such dense electronic-photonic integration enables implementation and control of quantum-photonic systems at the scale required for useful quantum information processing with CMOS-fabricated chips.

Highly forbidden $\beta$ decays provide a sensitive test to nuclear models in a regime in which the decay goes through high spin-multipole states, similar to the neutrinoless double-$\beta$ decay process. There are only 3 nuclei ($^{50}$V, $^{113}$Cd, $^{115}$In) which undergo a $4^{\rm th}$ forbidden non-unique $\beta$ decay. In this work, we compare the experimental $^{113}$Cd spectrum to theoretical spectral shapes in the framework of the spectrum-shape method. We measured with high precision, with the lowest energy threshold and the best energy resolution ever, the $\beta$ spectrum of $^{113}$Cd embedded in a 0.43 kg CdWO$_4$ crystal, operated over 26 days as a bolometer at low temperature in the Canfranc underground laboratory (Spain). We performed a Bayesian fit of the experimental data to three nuclear models (IBFM-2, MQPM and NSM) allowing the reconstruction of the spectral shape as well as the half-life. The fit has two free parameters, one of which is the effective weak axial-vector coupling constant, $g_A^{\text{eff}}$, which resulted in $g_A^{\text{eff}}$ between 1.0 and 1.2, compatible with a possible quenching. Based on the fit, we measured the half-life of the $^{113}$Cd $\beta$ decay including systematic uncertainties as $7.73^{+0.60}_{-0.57} \times 10^{15}$ yr, in agreement with the previous experiments. These results represent a significant step towards a better understanding of low-energy nuclear processes., Comment: Accepted for publication by EPJC

In this article, we establish a "Gehring lemma" for a real function satisfying a reverse H\"older inequality on all "kinetic cylinders" contained in a large one: it asserts that the integrability degree of the function improves under such an assumption. The kinetic cylinders are derived from the non-commutative group of invariances of the Kolmogorov equation. Our contributions here are (1) the extension of Gehring's Lemma to this kinetic (hypoelliptic) scaling used to generate the cylinders, (2) the localisation of the lemma in this hypoelliptic context (using ideas from the elliptic theory), (3) the streamlining of a short and quantitative proof. We then use this lemma to establish that the velocity gradient of weak solutions to linear kinetic equations of Fokker-Planck type with rough coefficients have Lebesgue integrability strictly greater than two, while the natural energy estimate merely ensures that it is square integrable. Our argument here is new but relies on Poincar\'e-type inequalities established in previous works.

We prove that the Fisher information is monotone decreasing in time along solutions of the space-homogeneous Boltzmann equation for a large class of collision kernels covering all classical interactions derived from systems of particles. For general collision kernels, a sufficient condition for the monotonicity of the Fisher information along the flow is related to the best constant for an integro-differential inequality for functions on the sphere, which belongs in the family of the Log-Sobolev inequalities. As a consequence, we establish the existence of global smooth solutions to the space-homogeneous Boltzmann equation in the main situation of interest where this was not known, namely the regime of very soft potentials. This is opening the path to the completion of both the classical program of qualitative study of space-homogeneous Boltzmann equation, initiated by Carleman, and the program of using the Fisher information in the study of the Boltzmann equation, initiated by McKean. From the proofs and discussion emerges a strengthened picture of the links between kinetic theory, information theory and log-Sobolev inequalities., Comment: 35 pages

Trefftz schemes are high-order Galerkin methods whose discrete spaces are made of elementwise exact solutions of the underlying PDE. Trefftz basis functions can be easily computed for many PDEs that are linear, homogeneous, and have piecewise-constant coefficients. However, if the equation has variable coefficients, exact solutions are generally unavailable. Quasi-Trefftz methods overcome this limitation relying on elementwise "approximate solutions" of the PDE, in the sense of Taylor polynomials. We define polynomial quasi-Trefftz spaces for general linear PDEs with smooth coefficients and source term, describe their approximation properties and, under a non-degeneracy condition, provide a simple algorithm to compute a basis. We then focus on a quasi-Trefftz DG method for variable-coefficient elliptic diffusion-advection-reaction problems, showing stability and high-order convergence of the scheme. The main advantage over standard DG schemes is the higher accuracy for comparable numbers of degrees of freedom. For non-homogeneous problems with piecewise-smooth source term we propose to construct a local quasi-Trefftz particular solution and then solve for the difference. Numerical experiments in 2 and 3 space dimensions show the excellent properties of the method both in diffusion-dominated and advection-dominated problems., Comment: 26 pages, 6 figures, 2 tables, added some remarks and one figure

We prove that the entropy solution to a scalar conservation law posed on the real line with a flux that is discontinuous at one point (in the space variable) coincides with the derivative of the solution to a Hamilton-Jacobi (HJ) equation whose Hamiltonian is discontinuous. Flux functions (Hamiltonians) are not assumed to be convex in the state (gradient) variable. The proof consists in proving the convergence of two numerical schemes. We rely on the theory developed by B.~Andreianov, K.~H.~Karlsen and N.~H.~Risebro (\textit{Arch. Ration. Mech. Anal.}, 2011) for such scalar conservation laws and on the viscosity solution theory developed by the authors (\textit{arxiv}, 2023) for the corresponding HJ equation. This study allows us to characterise certain germs introduced in the AKR theory (namely maximal and complete ones) and relaxation operators introduced in the viscosity solution framework.

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

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

We explore the bound neutrons decay into invisible particles (e.g., $n\rightarrow 3 \nu$ or $nn \rightarrow 2 \nu$) in the JUNO liquid scintillator detector. The invisible decay includes two decay modes: $ n \rightarrow { inv} $ and $ nn \rightarrow { inv} $. The invisible decays of $s$-shell neutrons in $^{12}{\rm C}$ will leave a highly excited residual nucleus. Subsequently, some de-excitation modes of the excited residual nuclei can produce a time- and space-correlated triple coincidence signal in the JUNO detector. Based on a full Monte Carlo simulation informed with the latest available data, we estimate all backgrounds, including inverse beta decay events of the reactor antineutrino $\bar{\nu}_e$, natural radioactivity, cosmogenic isotopes and neutral current interactions of atmospheric neutrinos. Pulse shape discrimination and multivariate analysis techniques are employed to further suppress backgrounds. With two years of exposure, JUNO is expected to give an order of magnitude improvement compared to the current best limits. After 10 years of data taking, the JUNO expected sensitivities at a 90% confidence level are $\tau/B( n \rightarrow { inv} ) > 5.0 \times 10^{31} \, {\rm yr}$ and $\tau/B( nn \rightarrow { inv} ) > 1.4 \times 10^{32} \, {\rm yr}$., Comment: 28 pages, 7 figures, 4 tables

The current experiments searching for neutrinoless double-$\beta$ ($0\nu\beta\beta$) decay also collect large statistics of Standard Model allowed two-neutrino double-$\beta$ ($2\nu\beta\beta$) decay events. These can be used to search for Beyond Standard Model (BSM) physics via $2\nu\beta\beta$ decay spectral distortions. $^{100}$Mo has a natural advantage due to its relatively short half-life, allowing higher $2\nu\beta\beta$ decay statistics at equal exposures compared to the other isotopes. We demonstrate the potential of the dual read-out bolometric technique exploiting a $^{100}$Mo exposure of 1.47 kg $\times$ y, acquired in the CUPID-Mo experiment at the Modane underground laboratory (France). We set limits on $0\nu\beta\beta$ decays with the emission of one or more Majorons, on $2\nu\beta\beta$ decay with Lorentz violation, and $2\nu\beta\beta$ decay with a sterile neutrino emission. In this analysis, we investigate the systematic uncertainty induced by modeling the $2\nu\beta\beta$ decay spectral shape parameterized through an improved model, an effect never considered before. This work motivates searches for BSM processes in the upcoming CUPID experiment, which will collect the largest amount of $2\nu\beta\beta$ decay events among the next-generation experiments.

We show the existence and uniqueness of fundamental solution operators to Kolmo\-gorov-Fokker-Planck equations with rough (measurable) coefficients and local or integral diffusion on finite and infinite time strips. In the local case, that is to say when the diffusion operator is of differential type, we prove $\L^2$ decay using Davies' method and the conservation property. We also prove that the existence of a generalized fundamental solution with the expected pointwise Gaussian upper bound is equivalent to Moser's $\L^2-\L^\infty$ estimates for local weak solutions to the equation and its adjoint. When coefficients are real, this gives the existence and uniqueness of such a generalized fundamental solution and a new and natural way to obtain pointwise decay., Comment: Revision following the referee's suggestion. Edition of the presentation of the homogeneous Sobolev norms for clarity. Most changes are in Section 2

In this article, we establish embeddings {\`a} la Lions and transfer of regularity {\`a} la Bouchut for a large scale of kinetic spaces. We use them to identify a notion of weak solutions to Kolmogorov-Fokker-Planck equations with (local or integral) diffusion and rough (measurable) coefficients under minimal requirements. We prove their existence and uniqueness for a large class of source terms, first in full space for the time, position and velocity variables and then for the kinetic Cauchy problem on infinite and finite time intervals., Comment: 40 pages, submitted. A correction to the argument in Theorem 6.7 and the corresponding argument in Section 7. Statements unchanged. Some typos eliminated

Neutrino detectors are among the largest photon detection instruments, built to capture scarce photons upon energy deposition. Many discoveries in neutrino physics, including the neutrino itself, are inseparable from the advances in photon detection technology, particularly in photo-sensors and readout electronics, to yield ever higher precision and richer detection information. The measurement of the energy of neutrinos, referred to as calorimetry, can be achieved in two distinct approaches: photon-counting, where single-photon can be counted digitally, and photon-integration, where multi-photons are aggregated and estimated via analogue signals. The energy is pursued today to reach permille level systematics control precision in ever-vast volumes, exemplified by experiments like JUNO. The unprecedented precision brings to the foreground the systematics due to calorimetric response entanglements in energy, position and time that were negligible in the past, thus driving further innovation in calorimetry. This publication describes a novel articulation that detectors can be endowed with multiple photon detection systems. This multi-calorimetry approach opens the notion of dual-calorimetry detector, consisting of both photon-counting and photon-integration systems, as a cost-effective evolution from the single calorimetry setups used over several decades for most experiments so far. The dual-calorimetry design exploits unique response synergies between photon-counting and photon-integration systems, including correlations and cancellations in calorimetric responses, to maximise the mitigation of response entanglements, thereby yielding permille-level high-precision calorimetry.

We prove that the set of singular times for weak solutions of the homogeneous Boltzmann equation with very soft potentials constructed as in Villani (1998) has Hausdorff dimension at most $\frac{|\gamma+2s|}{2s}$ with $\gamma \in [-4s,-2s)$ and $s \in (0,1)$.

In this note, we consider an evolution coercive Hamilton-Jacobi equation posed in a domain and supplemented with a boundary condition. We are interested in proving a comparison principle in the case where the time and the (normal) gradient variables are strongly coupled at the boundary. We elaborate on a method introduced by P.-L. Lions and P. Souganidis (Atti Accad. Naz. Lincei, 2017) to extend their comparison principle to more general boundary conditions and to Hamiltonians that are not globally Lipschitz continuous in the time variable. Their argument relies on a single blow-up procedure after rescaling the semi-solutions to be compared. We refer to our technique as the twin blow-ups method since two blow-ups are performed simultaneously, one for each variable of the doubling variable method. The Lipschitz regularity of the regularized subsolution provides a key Lipschitz inequality satisfied by the two blow-up limits, that are a priori allowed to be infinite. For expository reasons, the result is presented here in the framework of space dimension one and the general case is treated in a companion paper.

This work is concerned with Hamilton-Jacobi equations of evolution type posed in domains and supplemented with boundary conditions. Hamiltonians are coercive but are neither convex nor quasiconvex. We analyse boundary conditions when understood in the sense of viscosity solutions. This analysis is based on the study of boundary conditions of evolution type. More precisely, we give a new formula for the relaxed boundary conditions derived by J. Guerand (J. Differ. Equations, 2017). This new point of view unveils a connection between the relaxation operator and the classical Godunov flux from the theory of conservation laws. We apply our methods to two classical boundary value problems. It is shown that the relaxed Neumann boundary condition is expressed in terms of Godunov's flux while the relaxed Dirichlet boundary condition reduces to an obstacle problem at the boundary associated with the lower non-increasing envelope of the Hamiltonian.

The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neutrinos at the Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton liquid scintillator detector currently under construction in South China. The real-time monitoring system is designed to ensure both prompt alert speed and comprehensive coverage of progenitor stars. It incorporates prompt monitors on the electronic board as well as online monitors at the data acquisition stage. Assuming a false alert rate of 1 per year, this monitoring system exhibits sensitivity to pre-SN neutrinos up to a distance of approximately 1.6 (0.9) kiloparsecs and SN neutrinos up to about 370 (360) kiloparsecs for a progenitor mass of 30 solar masses, considering both normal and inverted mass ordering scenarios. The pointing ability of the CCSN is evaluated by analyzing the accumulated event anisotropy of inverse beta decay interactions from pre-SN or SN neutrinos. This, along with the early alert, can play a crucial role in facilitating follow-up multi-messenger observations of the next galactic or nearby extragalactic CCSN., Comment: 24 pages, 9 figures, accepted for the publication at JCAP

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

Neutrinoless double beta decay ($0\nu\beta\beta$) is a yet unobserved nuclear process which would demonstrate Lepton Number violation, a clear evidence of beyond Standard Model physics. The process two neutrino double beta decay ($2\nu\beta\beta)$ is allowed by the Standard Model and has been measured in numerous experiments. In this letter, we report a measurement of $2\nu\beta\beta$ decay half-life of $^{100}$Mo to the ground state of $^{100}$Ru of $(7.07~\pm~0.02~\text{(stat.)}~\pm~0.11~\text{(syst.)})~\times~10^{18}$~yr by the CUPID-Mo experiment. With a relative precision of $\pm~1.6$ \% this is the most precise measurement to date of a $2\nu\beta\beta$ decay rate in $^{100}$Mo. In addition, we constrain higher-order corrections to the spectral shape which provides complementary nuclear structure information. We report a novel measurement of the shape factor $\xi_{3,1}=0.45~\pm 0.03~\text{(stat.)} \ \pm 0.05 \ \text{(syst.)}$, which is compared to theoretical predictions for different nuclear models. We also extract the first value for the effective axial vector coupling constant obtained from a spectral shape study of $2\nu\beta\beta$ decay., Comment: 8 pages, 5 figures

CUPID-Mo, located in the Laboratoire Souterrain de Modane (France), was a demonstrator for the next generation $0\nu\beta\beta$ decay experiment, CUPID. It consisted of an array of 20 enriched Li$_{2}$$ ^{100}$MoO$_4$ bolometers and 20 Ge light detectors and has demonstrated that the technology of scintillating bolometers with particle identification capabilities is mature. Furthermore, CUPID-Mo can inform and validate the background prediction for CUPID. In this paper, we present a detailed model of the CUPID-Mo backgrounds. This model is able to describe well the features of the experimental data and enables studies of the $2\nu\beta\beta$ decay and other processes with high precision. We also measure the radio-purity of the Li$_{2}$$^{100}$MoO$_4$ crystals which are found to be sufficient for the CUPID goals. Finally, we also obtain a background index in the region of interest of 3.7$^{+0.9}_{-0.8}$(stat)$^{+1.5}_{-0.7}$(syst)$\times10^{-3}$counts/$\Delta$E$_{FWHM}$/mol$_{iso}$/yr, the lowest in a bolometric $0\nu\beta\beta$ decay experiment.

CUPID is a next-generation bolometric experiment aiming at searching for neutrinoless double-beta decay with ~250 kg of isotopic mass of $^{100}$Mo. It will operate at $\sim$10 mK in a cryostat currently hosting a similar-scale bolometric array for the CUORE experiment at the Gran Sasso National Laboratory (Italy). CUPID will be based on large-volume scintillating bolometers consisting of $^{100}$Mo-enriched Li$_2$MoO$_4$ crystals, facing thin Ge-wafer-based bolometric light detectors. In the CUPID design, the detector structure is novel and needs to be validated. In particular, the CUORE cryostat presents a high level of mechanical vibrations due to the use of pulse tubes and the effect of vibrations on the detector performance must be investigated. In this paper we report the first test of the CUPID-design bolometric light detectors with NTD-Ge sensors in a dilution refrigerator equipped with a pulse tube in an above-ground lab. Light detectors are characterized in terms of sensitivity, energy resolution, pulse time constants, and noise power spectrum. Despite the challenging noisy environment due to pulse-tube-induced vibrations, we demonstrate that all the four tested light detectors comply with the CUPID goal in terms of intrinsic energy resolution of 100 eV RMS baseline noise. Indeed, we have measured 70--90 eV RMS for the four devices, which show an excellent reproducibility. We have also obtained outstanding energy resolutions at the 356 keV line from a $^{133}$Ba source with one light detector achieving 0.71(5) keV FWHM, which is -- to our knowledge -- the best ever obtained when compared to $\gamma$ detectors of any technology in this energy range., Comment: Prepared for submission to JINST; 16 pages, 7 figures, and 1 table

An array of twelve 0.28 kg lithium molybdate (LMO) low-temperature bolometers equipped with 16 bolometric Ge light detectors, aiming at optimization of detector structure for CROSS and CUPID double-beta decay experiments, was constructed and tested in a low-background pulse-tube-based cryostat at the Canfranc underground laboratory in Spain. Performance of the scintillating bolometers was studied depending on the size of phonon NTD-Ge sensors glued to both LMO and Ge absorbers, shape of the Ge light detectors (circular vs. square, from two suppliers), in different light collection conditions (with and without reflector, with aluminum coated LMO crystal surface). The scintillating bolometer array was operated over 8 months in the low-background conditions that allowed to probe a very low, $\mu$Bq/kg, level of the LMO crystals radioactive contamination by $^{228}$Th and $^{226}$Ra., Comment: Prepared for submission to JINST; 23 pages, 9 figures, and 4 tables

The Jiangmen Underground Neutrino Observatory (JUNO), the first multi-kton liquid scintillator detector, which is under construction in China, will have a unique potential to perform a real-time measurement of solar neutrinos well below the few MeV threshold typical for Water Cherenkov detectors. JUNO's large target mass and excellent energy resolution are prerequisites for reaching unprecedented levels of precision. In this paper, we provide estimation of the JUNO sensitivity to 7Be, pep, and CNO solar neutrinos that can be obtained via a spectral analysis above the 0.45 MeV threshold. This study is performed assuming different scenarios of the liquid scintillator radiopurity, ranging from the most opti mistic one corresponding to the radiopurity levels obtained by the Borexino experiment, up to the minimum requirements needed to perform the neutrino mass ordering determination with reactor antineutrinos - the main goal of JUNO. Our study shows that in most scenarios, JUNO will be able to improve the current best measurements on 7Be, pep, and CNO solar neutrino fluxes. We also perform a study on the JUNO capability to detect periodical time variations in the solar neutrino flux, such as the day-night modulation induced by neutrino flavor regeneration in Earth, and the modulations induced by temperature changes driven by helioseismic waves.

Random coincidences of events could be one of the main sources of background in the search for neutrino-less double-beta decay of $^{100}$Mo with macro-bolometers, due to their modest time resolution. Scintillating bolometers as those based on Li$_2$MoO$_4$ crystals and employed in the CROSS and CUPID experiments can eventually exploit the coincident fast signal detected in a light detector to reduce this background. However, the scintillation provides a modest signal-to-noise ratio, making difficult a pile-up pulse-shape recognition and rejection at timescales shorter than a few ms. Neganov-Trofimov-Luke assisted light detectors (NTL-LDs) offer the possibility to effectively increase the signal-to-noise ratio, preserving a fast time-response, and enhance the capability of pile-up rejection via pulse shape analysis. In this article we present: a) an experimental work performed with a Li$_2$MoO$_4$ scintillating bolometer, studied in the framework of the CROSS experiment, and utilizing a NTL-LD; b) a simulation method to reproduce, synthetically, randomly coincident two-neutrino double-beta decay events; c) a new analysis method based on a pulse-shape discrimination algorithm capable of providing high pile-up rejection efficiencies. We finally show how the NTL-LDs offer a balanced solution between performance and complexity to reach background index $\sim$$10^{-4}$ counts/keV/kg/year with 280~g Li$_2$MoO$_4$ ($^{100}$Mo enriched) bolometers at 3034 keV, the Q-value of the double-beta decay, and target the goal of a next generation experiment like CUPID.

Subcellular RNA localization is a critical mechanism for the spatial control of gene expression. Its mechanism and precise functional role is not yet very well understood. Single Molecule Fluorescence in Situ Hybridization (smFISH) images allow for the detection of individual RNA molecules with subcellular accuracy. In return, smFISH requires robust methods to quantify and classify RNA spatial distribution. Here, we present PointFISH, a novel computational approach for the recognition of RNA localization patterns. PointFISH is an attention-based network for computing continuous vector representations of RNA point clouds. Trained on simulations only, it can directly process extracted coordinates from experimental smFISH images. The resulting embedding allows scalable and flexible spatial transcriptomics analysis and matches performance of hand-crafted pipelines.

Deep learning provides us with powerful methods to perform nucleus or cell segmentation with unprecedented quality. However, these methods usually require large training sets of manually annotated images, which are tedious and expensive to generate. In this paper we propose to use In Silico Labeling (ISL) as a pretraining scheme for segmentation tasks. The strategy is to acquire label-free microscopy images (such as bright-field or phase contrast) along fluorescently labeled images (such as DAPI or CellMask). We then train a model to predict the fluorescently labeled images from the label-free microscopy images. By comparing segmentation performance across several training set sizes, we show that such a scheme can dramatically reduce the number of required annotations.

The CUPID-Mo experiment, located at Laboratoire Souterrain de Modane (France), was a demonstrator experiment for CUPID. It consisted of an array of 20 Li$_2^{100}$MoO$_4$ (LMO) calorimeters each equipped with a Ge light detector (LD) for particle identification. In this work, we present the result of a search for two-neutrino and neutrinoless double beta decays of $^{100}$Mo to the first 0$^+$ and $2^+$ excited states of $^{100}$Ru using the full CUPID-Mo exposure (2.71 kg$\times$yr of LMO). We measure the half-life of $2\nu\beta\beta$ decay to the $0^{+}_1$ state as $T_{1/2}^{2\nu \rightarrow 0_1^+}=7.5\pm 0.8 \ \text{(stat.)} \ ^{+ 0.4}_{-0.3} \ \text{(syst.)} )\times 10^{20} \ \mathrm{yr}$. The bolometric technique enables measurement of the electron energies as well as the gamma rays from nuclear de-excitation and this allows us to set new limits on the two-neutrino decay to the $2_1^+$ state of $T^{2\nu \rightarrow 2_1^+}_{1/2}>4.4\times 10^{21} \ \mathrm{yr} \ \text{(90 % c.i.)}$ and on the neutrinoless modes of $T_{1/2}^{0\nu\rightarrow 2_1^+}>2.1\times10^{23} \ \mathrm{yr}\ \text{(90 % c.i.)}$, $T_{1/2}^{0\nu\rightarrow 0_1^+}>1.2\times10^{23} \ \mathrm{yr}\ \text{(90 % c.i.)}$. Information on the electrons spectral shape is obtained which allows us to make the first comparison of the single state (SSD) and higher state (HSD) $2\nu\beta\beta$ decay models for the $0_1^+$ excited state of $^{100}$Ru., Comment: 22 pages, 15 figures

This paper deals with the space-homogenous Landau equation with very soft potentials, including the Coulomb case. This nonlinear equation is of parabolic type with diffusion matrix given by the convolution product of the solution with the matrix $a_{ij} (z)=|z|^\gamma (|z|^2 \delta_{ij} - z_iz_j)$ for $\gamma \in [-3,-2)$. We derive local truncated entropy estimates and use them to establish two facts. Firstly, we prove that the set of singular points (in time and velocity) for the weak solutions constructed as in [C. Villani, Arch. Rational Mech. Anal. 143 (1998), 273-307] has zero $\mathscr{P}^{m_\ast}$ parabolic Hausdorff measure with $m_\ast:= \frac72 |2+\gamma|$. Secondly, we prove that if such a weak solution is axisymmetric, then it is smooth away from the symmetry axis. In particular, radially symmetric weak solutions are smooth away from the origin., Comment: 58 pages. Hausdorff dimension improved

Adipose tissue (AT) transcriptome studies provide holistic pictures of adaptation to weight and related bioclinical settings changes. Objective To implement AT gene expression profiling and investigate the link between changes in bioclinical parameters and AT gene expression during 3 steps of a 2-phase dietary intervention (DI). Methods AT transcriptome profiling was obtained from sequencing 1051 samples, corresponding to 556 distinct individuals enrolled in a weight loss intervention (8-week low-calorie diet (LCD) at 800 kcal/day) followed with a 6-month ad libitum randomized DI. Transcriptome profiles obtained with QuantSeq sequencing were benchmarked against Illumina RNAseq. Reverse transcription quantitative polymerase chain reaction was used to further confirm associations. Cell specificity was assessed using freshly isolated cells and THP-1 cell line. Results During LCD, 5 modules were found, of which 3 included at least 1 bioclinical variable. Change in body mass index (BMI) connected with changes in mRNA level of genes with inflammatory response signature. In this module, change in BMI was negatively associated with changes in expression of genes encoding secreted protein (GDF15, CCL3, and SPP1). Through all phases of the DI, change in GDF15 was connected to changes in SPP1, CCL3, LIPA and CD68. Further characterization showed that these genes were specific to macrophages (with LIPA, CD68 and GDF15 expressed in anti-inflammatory macrophages) and GDF15 also expressed in preadipocytes. Conclusion Network analyses identified a novel AT feature with GDF15 upregulated with calorie restriction induced weight loss, concomitantly to macrophage markers. In AT, GDF15 was expressed in preadipocytes and macrophages where it was a hallmark of anti-inflammatory cells.

Current experiments to search for broken lepton-number symmetry through the observation of neutrinoless double-beta decay ($0\mathrm{\nu\beta\beta}$) provide the most stringent limits on the Majorana nature of neutrinos and the effective Majorana neutrino mass ($m_{\beta\beta}$). The next-generation experiments will focus on the sensitivity to the $0\mathrm{\nu\beta\beta}$ half-life of $\mathcal{O}(10^{27}$--$10^{28}$~years$)$ and $m_{\beta\beta}\lesssim15$~meV, which would provide complete coverage of the so-called Inverted Ordering region of the neutrino mass parameter space. By taking advantage of recent technological breakthroughs, new, future calorimetric experiments at the 1-ton scale can increase the sensitivity by at least another order of magnitude, exploring the large fraction of the parameter space that corresponds to the Normal neutrino mass ordering. In case of a discovery, such experiments could provide important insights toward a new understanding of the mechanism of $0\mathrm{\nu\beta\beta}$. We present here a series of projects underway that will provide advancements in background reduction, cryogenic readout, and physics searches beyond $0\mathrm{\nu\beta\beta}$, all moving toward the next-to-next generation CUPID-1T detector., Comment: contribution to Snowmass 2021

The CUPID-Mo experiment to search for 0$\nu\beta\beta$ decay in $^{100}$Mo has been recently completed after about 1.5 years of operation at Laboratoire Souterrain de Modane (France). It served as a demonstrator for CUPID, a next generation 0$\nu\beta\beta$ decay experiment. CUPID-Mo was comprised of 20 enriched Li$_2$$^{100}$MoO$_4$ scintillating calorimeters, each with a mass of $\sim$ 0.2 kg, operated at $\sim$20 mK. We present here the final analysis with the full exposure of CUPID-Mo ($^{100}$Mo exposure of 1.47 kg$\times$yr) used to search for lepton number violation via 0$\nu\beta\beta$ decay. We report on various analysis improvements since the previous result on a subset of data, reprocessing all data with these new techniques. We observe zero events in the region of interest and set a new limit on the $^{100}$Mo 0$\nu\beta\beta$ decay half-life of $T^{0\nu}_{1/2} > 1.8 \times 10^{24}$ year (stat.+syst.) at 90% CI. Under the light Majorana neutrino exchange mechanism this corresponds to an effective Majorana neutrino mass of $\left < (0.28$--$0.49)$ eV, dependent upon the nuclear matrix element utilized.

CUPID will be a next generation experiment searching for the neutrinoless double $\beta$ decay, whose discovery would establish the Majorana nature of the neutrino. Based on the experience achieved with the CUORE experiment, presently taking data at LNGS, CUPID aims to reach a background free environment by means of scintillating Li$_{2}$$^{100}$MoO$_4$ crystals coupled to light detectors. Indeed, the simultaneous heat and light detection allows us to reject the dominant background of $\alpha$ particles, as proven by the CUPID-0 and CUPID-Mo demonstrators. In this work we present the results of the first test of the CUPID baseline module. In particular, we propose a new optimized detector structure and light sensors design to enhance the engineering and the light collection, respectively. We characterized the heat detectors, achieving an energy resolution of (5.9 $\pm$ 0.2) keV FWHM at the $Q$-value of $^{100}$Mo (about 3034 keV). We studied the light collection of the baseline CUPID design with respect to an alternative configuration which features gravity-assisted light detectors' mounting. In both cases we obtained an improvement in the light collection with respect to past measures and we validated the particle identification capability of the detector, which ensures an $\alpha$ particle rejection higher than 99.9%, fully satisfying the requirements for CUPID., Comment: 10 pages, 5 figures

The earth observation industry provides satellite imagery with high spatial resolution and short revisit time. To allow efficient operational employment of these images, automating certain tasks has become necessary. In the defense domain, aircraft detection on satellite imagery is a valuable tool for analysts. Obtaining high performance detectors on such a task can only be achieved by leveraging deep learning and thus us-ing a large amount of labeled data. To obtain labels of a high enough quality, the knowledge of military experts is needed.We propose a hybrid clustering active learning method to select the most relevant data to label, thus limiting the amount of data required and further improving the performances. It combines diversity- and uncertainty-based active learning selection methods. For aircraft detection by segmentation, we show that this method can provide better or competitive results compared to other active learning methods.

Trefftz methods are numerical methods for the approximation of solutions to boundary and/or initial value problems. They are Galerkin methods with particular test and trial functions, which solve locally the governing partial differential equation (PDE). This property is called the Trefftz property. Quasi-Trefftz methods were introduced to leverage the advantages of Trefftz methods for problems governed by variable coefficient PDEs, by relaxing the Trefftz property into a so-called quasi-Trefftz property: test and trial functions are not exact solutions but rather local approximate solutions to the governing PDE. In order to develop quassi-Trefftz methods for aero-acoustics problems governed by the convected Helmholtz equation, the present work tackles the question of the definition, construction and approximation properties of three families of quasi-Trefftz functions: two based on generalizations on plane wave solutions, and one polynomial. The polynomial basis shows significant promise as it does not suffer from the ill-conditioning issue inherent to wave-like bases.

This article deals with kinetic Fokker-Planck equations with essentially bounded coefficients. A weak Harnack inequality for non-negative super-solutions is derived by considering their Log-transform and following S. N. Kruzhkov (1963). Such a result rests on a new weak Poincar{\'e} inequality sharing similarities with the one introduced by W. Wang and L. Zhang in a series of works about ultraparabolic equations (2009, 2011, 2017). This functional inequality is combined with a classical covering argument recently adapted by L. Silvestre and the second author (2020) to kinetic equations.

A scintillating bolometer based on a large cubic Li$_{2}$$^{100}$MoO$_4$ crystal (45 mm side) and a Ge wafer (scintillation detector) has been operated in the CROSS cryogenic facility at the Canfranc underground laboratory in Spain. The dual-readout detector is a prototype of the technology that will be used in the next-generation $0\nu2\beta$ experiment CUPID. The measurements were performed at 18 and 12 mK temperature in a pulse tube dilution refrigerator. This setup utilizes the same technology as the CUORE cryostat that will host CUPID and so represents an accurate estimation of the expected performance. The Li$_{2}$$^{100}$MoO$_4$ bolometer shows a high energy resolution of 6 keV FWHM at the 2615 keV $\gamma$ line. The detection of scintillation light for each event triggered by the Li$_{2}$$^{100}$MoO$_4$ bolometer allowed for a full separation ($\sim$8$\sigma$) between $\gamma$($\beta$) and $\alpha$ events above 2 MeV. The Li$_{2}$$^{100}$MoO$_4$ crystal also shows a high internal radiopurity with $^{228}$Th and $^{226}$Ra activities of less than 3 and 8 $\mu$Bq/kg, respectively. Taking also into account the advantage of a more compact and massive detector array, which can be made of cubic-shaped crystals (compared to the cylindrical ones), this test demonstrates the great potential of cubic Li$_{2}$$^{100}$MoO$_4$ scintillating bolometers for high-sensitivity searches for the $^{100}$Mo $0\nu2\beta$ decay in CROSS and CUPID projects., Comment: 19 pages, 7 figures, 1 table

The CUPID Collaboration is designing a tonne-scale, background-free detector to search for double beta decay with sufficient sensitivity to fully explore the parameter space corresponding to the inverted neutrino mass hierarchy scenario. One of the CUPID demonstrators, CUPID-Mo, has proved the potential of enriched Li$_{2}$$^{100}$MoO$_4$ crystals as suitable detectors for neutrinoless double beta decay search. In this work, we characterised cubic crystals that, compared to the cylindrical crystals used by CUPID-Mo, are more appealing for the construction of tightly packed arrays. We measured an average energy resolution of (6.7$\pm$0.6) keV FWHM in the region of interest, approaching the CUPID target of 5 keV FWHM. We assessed the identification of $\alpha$ particles with and without a reflecting foil that enhances the scintillation light collection efficiency, proving that the baseline design of CUPID already ensures a complete suppression of this $\alpha$-induced background contribution. We also used the collected data to validate a Monte Carlo simulation modelling the light collection efficiency, which will enable further optimisations of the detector.

Precise characterization of detector time resolution is of crucial importance for next-generation cryogenic-bolometer experiments searching for neutrinoless double-beta decay, such as CUPID, in order to reject background due to pile-up of two-neutrino double-beta decay events. In this paper, we describe a technique developed to study the pile-up rejection capability of cryogenic bolometers. Our approach, which consists of producing controlled pile-up events with a programmable waveform generator, has the benefit that we can reliably and reproducibly control the time separation and relative energy of the individual components of the generated pile-up events. The resulting data allow us to optimize and benchmark analysis strategies to discriminate between individual and pile-up pulses. We describe a test of this technique performed with a small array of detectors at the Laboratori Nazionali del Gran Sasso, in Italy; we obtain a 90% rejection efficiency against pulser-generated pile-up events with rise time of ~15ms down to time separation between the individual events of about 2ms.

Trefftz methods are high-order Galerkin schemes in which all discrete functions are elementwise solution of the PDE to be approximated. They are viable only when the PDE is linear and its coefficients are piecewise constant. We introduce a 'quasi-Trefftz' discontinuous Galerkin method for the discretisation of the acoustic wave equation with piecewise-smooth wavespeed: the discrete functions are elementwise approximate PDE solutions. We show that the new discretisation enjoys the same excellent approximation properties as the classical Trefftz one, and prove stability and high-order convergence of the DG scheme. We introduce polynomial basis functions for the new discrete spaces and describe a simple algorithm to compute them. The technique we propose is inspired by the generalised plane waves previously developed for time-harmonic problems with variable coefficients; it turns out that in the case of the time-domain wave equation under consideration the quasi-Trefftz approach allows for polynomial basis functions., Comment: 29 pages, 11 figures, 4 tables

A Bacillus subtilis derivative was obtained from strain ATCC 6633 by replacement of the native promoter of the mycosubtilin operon by a constitutive promoter originating from the replication gene repU of the Staphylococcus aureus plasmid pUB110. The recombinant strain, designated BBG100, produced up to 15-fold more mycosubtilin than the wild type produced. The overproducing phenotype was related to enhancement of the antagonistic activities against several yeasts and pathogenic fungi. Hemolytic activities were also clearly increased in the modified strain. Mass spectrometry analyses of enriched mycosubtilin extracts showed similar patterns of lipopeptides for BBG100 and the wild type. Interestingly, these analyses also revealed a new form of mycosubtilin which was more easily detected in the BBG100 sample. When tested for its biocontrol potential, wild-type strain ATCC 6633 was almost ineffective for reducing a Pythium infection of tomato seedlings. However, treatment of seeds with the BBG100 overproducing strain resulted in a marked increase in the germination rate of seeds. This protective effect afforded by mycosubtilin overproduction was also visualized by the significantly greater fresh weight of emerging seedlings treated with BBG100 compared to controls or seedlings inoculated with the wild-type strain.

Generalized Plane Waves (GPWs) were introduced to take advantage of Trefftz methods for problems modeled by variable coefficient equations. Despite the fact that GPWs do not satisfy the Trefftz property, i.e. they are not exact solutions to the governing equation, they instead satisfy a quasi-Trefftz property: they are only approximate solutions. They lead to high order numerical methods, and this quasi-Trefftz property is critical for their numerical analysis. The present work introduces a new family of GPWs, amplitude-based. The motivation lies in the poor behavior of the phase-based GPW approximations in the pre-asymptotic regime, which will be tamed by avoiding high degree polynomials within an exponential. The new ansatz is introduces higher order terms in the amplitude rather than the phase of a plane wave as was initially proposed. The new functions' construction and the study of their interpolation properties are guided by the roadmap proposed in [16]. For the sake of clarity, the first focus is on the two-dimensional Helmholtz equation with spatially-varying wavenumber. The extension to a range of operators allowing for anisotropy in the first and second order terms follows. Numerical simulations illustrate the theoretical study of the new quasi-Trefftz functions.

The Boltzmann equation is a nonlinear partial differential equation that plays a central role in statistical mechanics. From the mathematical point of view, the existence of global smooth solutions for arbitrary initial data is an outstanding open problem. In the present article, we review a program focused on the type of particle interactions known as non-cutoff. It is dedicated to the derivation of a priori estimates in $C^\infty$, depending only on physically meaningful conditions. We prove that the solution will stay uniformly smooth provided that its mass, energy and entropy densities remain bounded, and away from vacuum.

This article is concerned with a porous medium equation whose pressure law is both nonlinear and nonlocal, namely $\partial_t u = { \nabla \cdot} \left(u \nabla(-\Delta)^{\frac{\alpha}{2}-1}u^{m-1} \right)$ where $u:\mathbb{R}_+\times \mathbb{R}^N \to \mathbb{R}_+$, for $0<\alpha<2$ and $m\geq2$. We prove that the $L^1\cap L^\infty$ weak solutions constructed by Biler, Imbert and Karch (2015) are locally H{\"o}lder-continuous in time and space. In this article, the classical parabolic De Giorgi techniques for the regularity of PDEs are tailored to fit this particular variant of the PME equation. In the spirit of the work of Caffarelli, Chan and Vasseur (2011), the two main ingredients are the derivation of local energy estimates and a so-called ''intermediate value lemma''. For $\alpha\leq1$, we adapt the proof of Caffarelli, Soria and V{\'a}zquez (2013), who treated the case of a linear pressure law. We then use a non-linear drift to cancel out the singular terms that would otherwise appear in the energy estimates.

We derive $C^\infty$ a priori estimates for solutions of the inhomogeneous Boltzmann equation without cut-off, conditional to point-wise bounds on their mass, energy and entropy densities. We also establish decay estimates for large velocities, for all derivatives of the solution.

In this self-contained document, we aim to present the basic theoretical building blocks to understand modeling of stellarator magnetic fields, some of the challenges associated with modeling, and optimization for designing stellarators. As often as possible, the ideas will be presented using equations and pictures, and references to other relevant introductory material will be included. This document is accessible to those who may not have a physics background but are interested in applications of mathematical and computational tools to stellarator research.

This work focuses on the study of partial differential equation (PDE) based basis function for Discontinuous Galerkin methods to solve numerically wave-related boundary value problems with variable coefficients. To tackle problems with constant coefficients, wave-based methods have been widely studied in the literature: they rely on the concept of Trefftz functions, i.e. local solutions to the governing PDE, using oscillating basis functions rather than polynomial functions to represent the numerical solution. Generalized Plane Waves (GPWs) are an alternative developed to tackle problems with variable coefficients, in which case Trefftz functions are not available. In a similar way, they incorporate information on the PDE, however they are only approximate Trefftz functions since they don't solve the governing PDE exactly, but only an approximated PDE. Considering a new set of PDEs beyond the Helmholtz equation, we propose to set a roadmap for the construction and study of local interpolation properties of GPWs. Identifying carefully the various steps of the process, we provide an algorithm to summarize the construction of these functions, and establish necessary conditions to obtain high order interpolation properties of the corresponding basis.

We prove that the set of singular times for weak solutions of the space homogeneous Landau equation with Coulomb potential constructed as in [C. Villani, Arch. Rational Mech. Anal. 143 (1998), 273-307] has Hausdorff dimension at most 1/2.

Numerical simulations of fluids in astrophysics and computational fluid dynamics (CFD) are among the most computationally-demanding calculations, in terms of sustained floating-point operations per second, or FLOP/s. It is expected that these numerical simulations will significantly benefit from the future Exascale computing infrastructures, that will perform 10^18 FLOP/s. The performance of the SPH codes is, in general, adversely impacted by several factors, such as multiple time-stepping, long-range interactions, and/or boundary conditions. In this work an extensive study of three SPH implementations SPHYNX, ChaNGa, and XXX is performed, to gain insights and to expose any limitations and characteristics of the codes. These codes are the starting point of an interdisciplinary co-design project, SPH-EXA, for the development of an Exascale-ready SPH mini-app. We implemented a rotating square patch as a joint test simulation for the three SPH codes and analyzed their performance on a modern HPC system, Piz Daint. The performance profiling and scalability analysis conducted on the three parent codes allowed to expose their performance issues, such as load imbalance, both in MPI and OpenMP. Two-level load balancing has been successfully applied to SPHYNX to overcome its load imbalance. The performance analysis shapes and drives the design of the SPH-EXA mini-app towards the use of efficient parallelization methods, fault-tolerance mechanisms, and load balancing approaches., Comment: arXiv admin note: substantial text overlap with arXiv:1809.08013

The study of positivity of solutions to the Boltzmann equation goes back to Carleman (1933), and the initial argument of Carleman was developed byPulvirenti-Wennberg (1997), the second author and Briant (2015). The appearance of a lower bound with Gaussian decay had however remained an open question for long-range interactions (the so-called non-cutoff collision kernels). We answer this question and establish such Gaussian lower bound for solutions to the Boltzmann equation without cutoff, in the case of hard and moderately soft potentials, with spatial periodic conditions, and under the sole assumption that hydrodynamic quantities (local mass, local energy and local entropy density) remain bounded. The paper is mostly self-contained, apart from the uniform upper bound on the solution established by the third author (2016)., Comment: 14 pages

We present a boundary integral equation solver for computing Taylor relaxed states in non-axisymmetric solid and shell-like toroidal geometries. The computation of Taylor states in these geometries is a key element for the calculation of stepped pressure stellarator equilibria. The integral representation of the magnetic field in this work is based on the generalized Debye source formulation, and results in a well-conditioned second-kind boundary integral equation. The integral equation solver is based on a spectral discretization of the geometry and unknowns, and the computation of the associated weakly-singular integrals is performed with high-order quadrature based on a partition of unity. The resulting scheme for applying the integral operator is then coupled with an iterative solver and suitable preconditioners. Several numerical examples are provided to demonstrate the accuracy and efficiency of our method, and a direct comparison with the leading code in the field is reported.