Search

Motivated by the growing use of Artificial Intelligence (AI) tools in control design, this paper takes the first steps towards bridging the gap between results from Direct Gradient methods for the Linear Quadratic Regulator (LQR), and neural networks. More specifically, it looks into the case where one wants to find a Linear Feed-Forward Neural Network (LFFNN) feedback that minimizes a LQR cost. This paper starts by computing the gradient formulas for the parameters of each layer, which are used to derive a key conservation law of the system. This conservation law is then leveraged to prove boundedness and global convergence of solutions to critical points, and invariance of the set of stabilizing networks under the training dynamics. This is followed by an analysis of the case where the LFFNN has a single hidden layer. For this case, the paper proves that the training converges not only to critical points but to the optimal feedback control law for all but a set of measure-zero of the initializations. These theoretical results are followed by an extensive analysis of a simple version of the problem (the ``vector case''), proving the theoretical properties of accelerated convergence and robustness for this simpler example. Finally, the paper presents numerical evidence of faster convergence of the training of general LFFNNs when compared to traditional direct gradient methods, showing that the acceleration of the solution is observable even when the gradient is not explicitly computed but estimated from evaluations of the cost function.

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 DArk Matter In CCDs at Modane (DAMIC-M) experiment is designed to search for light dark matter (m$_{\chi}$<10\,GeV/c$^2$) at the Laboratoire Souterrain de Modane (LSM) in France. DAMIC-M will use skipper charge-coupled devices (CCDs) as a kg-scale active detector target. Its single-electron resolution will enable eV-scale energy thresholds and thus world-leading sensitivity to a range of hidden sector dark matter candidates. A DAMIC-M prototype, the Low Background Chamber (LBC), has been taking data at LSM since 2022. The LBC provides a low-background environment, which has been used to characterize skipper CCDs, study dark current, and measure radiopurity of materials planned for DAMIC-M. It also allows testing of various subsystems like readout electronics, data acquisition software, and slow control. This paper describes the technical design and performance of the LBC.

The novelty of the current work is precisely to propose a statistical procedure to combine estimates of the modal parameters provided by any set of Operational Modal Analysis (OMA) algorithms so as to avoid preference for a particular one and also to derive an approximate joint probability distribution of the modal parameters, from which engineering statistics of interest such as mean value and variance are readily provided. The effectiveness of the proposed strategy is assessed considering measured data from an actual centrifugal compressor. The statistics obtained for both forward and backward modal parameters are finally compared against modal parameters identified during standard stability verification testing (SVT) of centrifugal compressors prior to shipment, using classical Experimental Modal Analysis (EMA) algorithms. The current work demonstrates that combination of OMA algorithms can provide quite accurate estimates for both the modal parameters and the associated uncertainties with low computational costs., Comment: 6 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

Numerical ``direct'' approaches to time-optimal control often fail to find solutions that are singular in the sense of the Pontryagin Maximum Principle, performing better when searching for saturated (bang-bang) solutions. In previous work by one of the authors, singular solutions were shown to exist for the time-optimal control problem for fully actuated mechanical systems under hard torque constraints. Explicit formulas, based on a Lie theoretic analysis of the problem, were given for singular segments of trajectories, but the global structure of solutions remains unknown. In this work, we review the aforementioned framework, and show how to effectively combine these formulas with the use of general-purpose optimal control software packages. By using the explicit formula given by the theory in the intervals where the numerical solution enters a singular arc, we not only obtain an algebraic expression for the control in that interval but we are also able to remove artifacts present in the numerical solution. In this way, the best features of numerical algorithms and theory complement each other and provide a better picture of the global optimal structure. We illustrate the technique on a two degree of freedom robotic arm example, using two distinct optimal control numerical software packages running on different programming languages.

For a linear Hamming metric code of length n over a finite field, the number of distinct weights of its codewords is at most n. The codes achieving the equality in the above bound were called full weight spectrum codes. In this paper we will focus on the analogous class of codes within the framework of cyclic subspace codes. Cyclic subspace codes have garnered significant attention, particularly for their applications in random network coding to correct errors and erasures. We investigate one-orbit cyclic subspace codes that are full weight spectrum in this context. Utilizing number theoretical results and combinatorial arguments, we provide a complete classification of full weight spectrum one-orbit cyclic subspace codes.

Subspace codes have recently been used for error correction in random network coding. In this work, we focus on one-orbit cyclic subspace codes. If $S$ is an $\mathbb{F}_q$-subspace of $\mathbb{F}_{q^n}$, then the one-orbit cyclic subspace code defined by $S$ is \[ \mathrm{Orb}(S)=\{\alpha S \colon \alpha \in \mathbb{F}_{q^n}^*\}, \]where $\alpha S=\lbrace \alpha s \colon s\in S\rbrace$ for any $\alpha\in \mathbb{F}_{q^n}^*$. Few classification results of subspace codes are known, therefore it is quite natural to initiate a classification of cyclic subspace codes, especially in the light of the recent classification of the isometries for cyclic subspace codes. We consider three-dimensional one-orbit cyclic subspace codes, which are divided into three families: the first one containing only $\mathrm{Orb}(\mathbb{F}_{q^3})$; the second one containing the optimum-distance codes; and the third one whose elements are codes with minimum distance $2$. We study inequivalent codes in the latter two families., Comment: Accepted to the 2024 IEEE International Symposium on Information Theory (ISIT 2024)

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.

The MEMENTO (MEtallic Melt Evolution in Next-step TOkamaks) code is a new numerical implementation of the physics model originally developed for the MEMOS-U code with the objective to self-consistently describe the generation of melt and its subsequent large scale dynamics in fusion devices and to assess the damage of metallic reactor armor under powerful normal and off-normal plasma events. The model has been validated in multiple dedicated EUROfusion experiments. MEMENTO solves the heat and phase transfer problem coupled with the incompressible Navier-Stokes equations in the shallow water approximation for the thin liquid film over the solid metal and with the current propagation equations on a domain that features a time-evolving deforming metal-plasma interface. The code utilizes non-uniform and adaptive meshing along with sub-cycling in time facilitated by the AMReX open-source framework as well as AMReX's built-in parallelization capabilities., Comment: 15 pages, 12 figures

The standard approach to test for deviations from General Relativity on cosmological scales is to combine measurements of the growth rate of structure with gravitational lensing. In this study, we show that this method suffers from an important limitation with regard to these two probes: models of dark matter with additional interactions can lead to the very same observational signatures found in modified gravity (and viceversa). Using synthetic data of redshift-space distortions, weak lensing, and cosmic microwave background, we demonstrate that this degeneracy is inevitable between modifications of gravity and a dark fifth force. We then show that the coming generation of surveys, in particular the Square Kilometer Array, will allow us to break the degeneracy between such models through measurements of gravitational redshift. Performing a Markov Chain Monte Carlo analysis of the synthetic data set, we quantify the extent to which gravitational redshift can distinguish between two representative classes of models, Generalized Brans-Dicke (modified gravity) and Coupled Quintessence (fifth force)., Comment: 10 pages, 3 figures

Virtual Reality (VR) has emerged as a promising tool for enhancing social skills and emotional well-being in individuals with Autism Spectrum Disorder (ASD). Through a technical exploration, this study employs a multiplayer serious gaming environment within VR, engaging 34 individuals diagnosed with ASD and employing high-precision biosensors for a comprehensive view of the participants' arousal and responses during the VR sessions. Participants were subjected to a series of 3 virtual scenarios designed in collaboration with stakeholders and clinical experts to promote socio-cognitive skills and emotional regulation in a controlled and structured virtual environment. We combined the framework with wearable non-invasive sensors for bio-signal acquisition, focusing on the collection of heart rate variability, and respiratory patterns to monitor participants behaviors. Further, behavioral assessments were conducted using observation and semi-structured interviews, with the data analyzed in conjunction with physiological measures to identify correlations and explore digital-intervention efficacy. Preliminary analysis revealed significant correlations between physiological responses and behavioral outcomes, indicating the potential of physiological feedback to enhance VR-based interventions for ASD. The study demonstrated the feasibility of using real-time data to adapt virtual scenarios, suggesting a promising avenue to support personalized therapy. The integration of quantitative physiological feedback into digital platforms represents a forward step in the personalized intervention for ASD. By leveraging real-time data to adjust therapeutic content, this approach promises to enhance the efficacy and engagement of digital-based therapies., Comment: 19 pages, 6 figures

The devices designed for the Internet-of-Things encompass a large variety of distinct processor architectures, forming a highly heterogeneous zoo. In order to tackle this, we employ a simulator to estimate the performance of the matrix-matrix multiplication (GEMM) kernel on processors designed to operate at the edge. Our simulator adheres to the modern implementations of GEMM, advocated by GotoBLAS2, BLIS, OpenBLAS, etc., to carefully account for the amount of data transfers across the memory hierarchy of different algorithmic variants of the kernel. %Armed with this tool, A small collection of experiments provide the necessary data to calibrate the simulator and deliver highly accurate estimations of the execution time for a given processor architecture., Comment: 12 pages, 2 Tables, 6 Figures

In many particle physics experiments the data processing is based on the analysis of the digitized waveforms provided by the detector. While the waveform amplitude is usually correlated to the event energy, the shape may carry useful information for event discrimination. Thanks to the high signal to noise ratio they provide, matched filters are often applied. Their original design is however intended for the estimation of the waveform amplitude only. In this work we introduce an analytical extension of the original matched filter for the estimation of a possible shape deformation with respect to a reference template. The new filter is validated on simulations and, with respect to shape parameters calculated on unfiltered waveforms or derived from the original matched filter, it improves the discrimination capability by at least a factor 2 both at low and high signal to noise ratios, making it applicable to the data of several experiments., Comment: 10 pages, 5 figures

The finite temperature version of the Vashishta-Singwi (VS) dielectric scheme for the paramagnetic warm dense uniform electron fluid is revisited correcting for an earlier thermodynamic derivative error. The VS scheme handles quantum mechanical effects at the level of the random phase approximation and treats correlations via the density expansion of a generalized Singwi-Tosi-Land-Sj\"olander (STLS) closure that inserts a parameter determined by enforcing the compressibility sum rule. Systematic comparison with quasi-exact results, based on quantum Monte Carlo simulations, reveals a structural superiority of the VS scheme towards strong coupling and a thermodynamic superiority of the STLS scheme courtesy of a favorable cancellation of errors. Guidelines are provided for the construction of dielectric schemes that are expected to be more accurate but computationally costly., Comment: 24 pages, 16 figures, 2 tables

Frequency-selective or even frequency-tunable Terahertz (THz) photodevices are critical components for many technological applications that require nanoscale manipulation, control and confinement of light. Within this context, gate-tunable phototransistors based on plasmonic resonances are often regarded as the most promising devices for frequency-selective detection of THz fields. The exploitation of constructive interference of plasma waves in such detectors not only promises frequency selectivity, but also a pronounced sensitivity enhancement at the target frequencies. However, clear signatures of plasmon-assisted resonances in THz detectors have been only revealed at cryogenic temperatures so far, and remain unobserved at application-relevant room-temperature conditions. In this work, we demonstrate the sought-after room-temperature resonant detection of THz radiation in short-channel gated photodetectors made from high-quality single-layer graphene. The survival of this intriguing resonant regime at room-temperature ultimately relies on the weak intrinsic electron-phonon scattering in graphene, which avoids the damping of the plasma oscillations.

Sidon spaces have been introduced by Bachoc, Serra and Z\'emor as the $q$-analogue of Sidon sets, classical combinatorial objects introduced by Simon Szidon. In 2018 Roth, Raviv and Tamo introduced the notion of $r$-Sidon spaces, as an extension of Sidon spaces, which may be seen as the $q$-analogue of $B_r$-sets, a generalization of classical Sidon sets. Thanks to their work, the interest on Sidon spaces has increased quickly because of their connection with cyclic subspace codes they pointed out. This class of codes turned out to be of interest since they can be used in random linear network coding. In this work we focus on a particular class of them, the one-orbit cyclic subspace codes, through the investigation of some properties of Sidon spaces and $r$-Sidon spaces, providing some upper and lower bounds on the possible dimension of their \textit{r-span} and showing explicit constructions in the case in which the upper bound is achieved. Moreover, we provide further constructions of $r$-Sidon spaces, arising from algebraic and combinatorial objects, and we show examples of $B_r$-sets constructed by means of them.

Upcoming galaxy surveys provide the necessary sensitivity to measure gravitational redshift, a general relativistic effect that generates a dipole in galaxy clustering data when correlating two distinct populations of galaxies. Here, we study the constraining power of gravitational redshift within the framework of the effective theory of interacting dark energy. This formalism describes linear cosmological perturbations in scalar-tensor theories of gravity with a limited number of free functions, and allows each particle species to be coupled differently to the gravitational sector. In this work, we focus on Horndeski theories with a non-minimal coupling of dark matter to the scalar degree of freedom, yielding a breaking of the weak equivalence principle for this cosmic component, a scenario that is yet untested. We show that the dipole generated by gravitational redshift significantly breaks degeneracies and tightens the constraints on the parameters of the effective theory compared to the standard redshift-space distortion analysis solely based on the even multipoles in the galaxy correlation function, with an improvement of up to $\sim 50\%$ for populations with a galaxy bias difference equal to 1. We make the Python package EF-TIGRE (Effective Field Theory of Interacting dark energy with Gravitational REdshift) developed for this work publicly available., Comment: 24+6 pages, 5+4 figures, code available at https://github.com/Mik3M4n/EF-TIGRE. v2 matches version in press on JCAP

We explore the utilization of the Apache TVM open source framework to automatically generate a family of algorithms that follow the approach taken by popular linear algebra libraries, such as GotoBLAS2, BLIS and OpenBLAS, in order to obtain high-performance blocked formulations of the general matrix multiplication (GEMM). % In addition, we fully automatize the generation process, by also leveraging the Apache TVM framework to derive a complete variety of the processor-specific micro-kernels for GEMM. This is in contrast with the convention in high performance libraries, which hand-encode a single micro-kernel per architecture using Assembly code. % In global, the combination of our TVM-generated blocked algorithms and micro-kernels for GEMM 1)~improves portability, maintainability and, globally, streamlines the software life cycle; 2)~provides high flexibility to easily tailor and optimize the solution to different data types, processor architectures, and matrix operand shapes, yielding performance on a par (or even superior for specific matrix shapes) with that of hand-tuned libraries; and 3)~features a small memory footprint., Comment: 35 pages, 22 figures. Submitted to ACM TOMS

The optimization of the matrix multiplication (or GEMM) has been a need during the last decades. This operation is considered the flagship of current linear algebra libraries such as BLIS, OpenBLAS, or Intel OneAPI because of its widespread use in a large variety of scientific applications. The GEMM is usually implemented following the GotoBLAS philosophy, which tiles the GEMM operands and uses a series of nested loops for performance improvement. These approaches extract the maximum computational power of the architectures through small pieces of hardware-oriented, high-performance code called micro-kernel. However, this approach forces developers to generate, with a non-negligible effort, a dedicated micro-kernel for each new hardware. In this work, we present a step-by-step procedure for generating micro-kernels with the Exo compiler that performs close to (or even better than) manually developed microkernels written with intrinsic functions or assembly language. Our solution also improves the portability of the generated code, since a hardware target is fully specified by a concise library-based description of its instructions., Comment: 11 pages, 18 figures. Presented at CGO 2024. It includes a software artifact step-by-step execution

Rare-event search experiments located on-surface, such as short-baseline reactor neutrino experiments, are often limited by muon-induced background events. Highly efficient muon vetos are essential to reduce the detector background and to reach the sensitivity goals. We demonstrate the feasibility of deploying organic plastic scintillators at sub-Kelvin temperatures. For the NUCLEUS experiment, we developed a cryogenic muon veto equipped with wavelength shifting fibers and a silicon photo multiplier operating inside a dilution refrigerator. The achievable compactness of cryostat-internal integration is a key factor in keeping the muon rate to a minimum while maximizing coverage. The thermal and light output properties of a plastic scintillation detector were examined. We report first data on the thermal conductivity and heat capacity of the polystyrene-based scintillator UPS-923A over a wide range of temperatures extending below one Kelvin. The light output was measured down to 0.8K and observed to increase by a factor of 1.61$\pm$0.05 compared to 300K. The development of an organic plastic scintillation muon veto operating in sub-Kelvin temperature environments opens new perspectives for rare-event searches with cryogenic detectors at sites lacking substantial overburden.

Charge-coupled devices (CCDs) are a leading technology in direct dark matter searches because of their eV-scale energy threshold and high spatial resolution. The sensitivity of future CCD experiments could be enhanced by distinguishing nuclear recoil signals from electronic recoil backgrounds in the CCD silicon target. We present a technique for event-by-event identification of nuclear recoils based on the spatial correlation between the primary ionization event and the lattice defect left behind by the recoiling atom, later identified as a localized excess of leakage current under thermal stimulation. By irradiating a CCD with an $^{241}$Am$^{9}$Be neutron source, we demonstrate $>93\%$ identification efficiency for nuclear recoils with energies $>150$ keV, where the ionization events were confirmed to be nuclear recoils from topology. The technique remains fully efficient down to 90 keV, decreasing to 50$\%$ at 8 keV, and reaching ($6\pm2$)$\%$ at 1.5--3.5 keV. Irradiation with a $^{24}$Na $\gamma$-ray source shows no evidence of defect generation by electronic recoils, with the fraction of electronic recoils with energies $<85$ keV that are spatially correlated with defects $<0.1$$\%$., Comment: 9 pages, 7 figures

We present the first continuous operation in a surface lab of BULLKID, a detector for searches of light Dark Matter and precision measurements of the coherent and elastic neutrino-nucleus scattering. The detector consists of an array of 60 cubic silicon particle absorbers of 0.34 g each, sensed by cryogenic kinetic inductance detectors. The data presented focusses on one of the central elements of the array and on its surrounding elements used as veto. The energy spectrum resulting from an exposure of 39 hours to ambient backgrounds, obtained without radiation shields, is flat at the level of $(2.0\pm0.1\,{\rm stat.}\pm0.2\,{\rm syst.})\times10^6$ counts / keV kg days down to the energy threshold of $160\pm13$ eV. The data analysis demonstrates the unique capability of rejecting backgrounds generated from interactions in other sites of the array, stemming from the segmented and monolithic structure of the detector.

The direct search for dark matter in the form of weakly interacting massive particles (WIMP) is performed by detecting nuclear recoils (NR) produced in a target material from the WIMP elastic scattering. A promising experimental strategy for direct dark matter search employs argon dual-phase time projection chambers (TPC). One of the advantages of the TPC is the capability to detect both the scintillation and charge signals produced by NRs. Furthermore, the existence of a drift electric field in the TPC breaks the rotational symmetry: the angle between the drift field and the momentum of the recoiling nucleus can potentially affect the charge recombination probability in liquid argon and then the relative balance between the two signal channels. This fact could make the detector sensitive to the directionality of the WIMP-induced signal, enabling unmistakable annual and daily modulation signatures for future searches aiming for discovery. The Recoil Directionality (ReD) experiment was designed to probe for such directional sensitivity. The TPC of ReD was irradiated with neutrons at the INFN Laboratori Nazionali del Sud, and data were taken with 72 keV NRs of known recoil directions. The direction-dependent liquid argon charge recombination model by Cataudella et al. was adopted and a likelihood statistical analysis was performed, which gave no indications of significant dependence of the detector response to the recoil direction. The aspect ratio R of the initial ionization cloud is estimated to be 1.037 +/- 0.027 and the upper limit is R < 1.072 with 90% confidence level, Comment: 20 pages, 10 figures, submitted to Eur. Phys. J. C

We present a novel approach to implement compressive sensing in laser scanning microscopes (LSM), specifically in image scanning microscopy (ISM), using a single-photon avalanche diode (SPAD) array detector. Our method addresses two significant limitations in applying compressive sensing to LSM: the time to compute the sampling matrix and the quality of reconstructed images. We employ a fixed sampling strategy, skipping alternate rows and columns during data acquisition, which reduces the number of points scanned by a factor of four and eliminates the need to compute different sampling matrices. By exploiting the parallel images generated by the SPAD array, we improve the quality of the reconstructed compressive-ISM images compared to standard compressive confocal LSM images. Our results demonstrate the effectiveness of our approach in producing higher-quality images with reduced data acquisition time and potential benefits in reducing photobleaching., Comment: Presented in ISCS23

The Lamport diagram is a pervasive and intuitive tool for informal reasoning about "happens-before" relationships in a concurrent system. However, traditional axiomatic formalizations of Lamport diagrams can be painful to work with in a mechanized setting like Agda. We propose an alternative, inductive formalization -- the causal separation diagram (CSD) -- that takes inspiration from string diagrams and concurrent separation logic, but enjoys a graphical syntax similar to Lamport diagrams. Critically, CSDs are based on the idea that causal relationships between events are witnessed by the paths that information follows between them. To that end, we model happens-before as a dependent type of paths between events. The inductive formulation of CSDs enables their interpretation into a variety of semantic domains. We demonstrate the interpretability of CSDs with a case study on properties of logical clocks, widely-used mechanisms for reifying causal relationships as data. We carry out this study by implementing a series of interpreters for CSDs, culminating in a generic proof of Lamport's clock condition that is parametric in a choice of clock. We instantiate this proof on Lamport's scalar clock, on Mattern's vector clock, and on the matrix clocks of Raynal et al. and of Wuu and Bernstein, yielding verified implementations of each. The CSD formalism and our case study are mechanized in the Agda proof assistant., Comment: This revision is as published in PACMPL through OOPSLA, but with [authorversion] set. Compared to the previous version, the introduction has been almost entirely rewritten

Dark Matter (DM) particles with sufficiently large cross sections may scatter as they travel through Earth's bulk. The corresponding changes in the DM flux give rise to a characteristic daily modulation signal in detectors sensitive to DM-electron interactions. Here, we report results obtained from the first underground operation of the DAMIC-M prototype detector searching for such a signal from DM with MeV-scale mass. A model-independent analysis finds no modulation in the rate of 1$e^-$ events with periods in the range 1-48 h. We then use these data to place exclusion limits on DM in the mass range [0.53, 2.7] MeV/c$^2$ interacting with electrons via a dark photon mediator. Taking advantage of the time-dependent signal we improve by $\sim$2 orders of magnitude on our previous limit obtained from the total rate of 1$e^-$ events, using the same data set. This daily modulation search represents the current strongest limit on DM-electron scattering via ultralight mediators for DM masses around 1 MeV/c$^2$.

We present results from a 3.25 kg-day target exposure of two silicon charge-coupled devices (CCDs), each with 24 megapixels and skipper readout, deployed in the DAMIC setup at SNOLAB. With a reduction in pixel readout noise of a factor of 10 relative to the previous detector, we investigate the excess population of low-energy events in the CCD bulk previously observed above expected backgrounds. We address the dominant systematic uncertainty of the previous analysis through a depth fiducialization designed to reject surface backgrounds on the CCDs. The measured bulk ionization spectrum confirms the presence of an excess population of low-energy events in the CCD target with characteristic rate of ${\sim}7$ events per kg-day and electron-equivalent energies of ${\sim}80~$eV, whose origin remains unknown., Comment: 8 pages, 4 figures

The emission of $\gamma$-rays after a neutron capture in a cryogenic detector can generate mono-energetic nuclear recoils in the sub-keV regime, of direct interest for the calibration of Dark Matter and Coherent Elastic Neutrino Nucleus Scattering experiments. Here we show that accurate predictions of the nuclear recoil spectra induced by neutron captures require taking into account the interplay between the development in time of the de-excitation $\gamma$-cascade of the target nucleus and that of the associated atomic collisions in matter. We present detailed simulations coupling the FIFRELIN code for the description of the $\gamma$-cascades and the IRADINA code for the modelling of the fast atomic movements in matter. Nuclear recoil spectra are predicted, and made available to the community, for concrete cases of Al$_2$O$_3$, Si, Ge and CaWO$_4$ crystals exposed to a low intensity beam of thermal neutrons. We find that timing effects cause new calibration peaks to emerge in the recoil spectra and also impact the shape of the continuous recoil distribution. We discuss how they could give access to a rich physics program, spanning the accurate study of the response of cryogenic detectors in the sub-keV range, tests of solid state physics simulations and tests of nuclear models.

Recent research in neural networks and machine learning suggests that using many more parameters than strictly required by the initial complexity of a regression problem can result in more accurate or faster-converging models -- contrary to classical statistical belief. This phenomenon, sometimes known as ``benign overfitting'', raises questions regarding in what other ways might overparameterization affect the properties of a learning problem. In this work, we investigate the effects of overfitting on the robustness of gradient-descent training when subject to uncertainty on the gradient estimation. This uncertainty arises naturally if the gradient is estimated from noisy data or directly measured. Our object of study is a linear neural network with a single, arbitrarily wide, hidden layer and an arbitrary number of inputs and outputs. In this paper we solve the problem for the case where the input and output of our neural-network are one-dimensional, deriving sufficient conditions for robustness of our system based on necessary and sufficient conditions for convergence in the undisturbed case. We then show that the general overparametrized formulation introduces a set of spurious equilibria which lay outside the set where the loss function is minimized, and discuss directions of future work that might extend our current results for more general formulations., Comment: 10 pages, 1 figure, extended conference version

Oscura is a planned light-dark matter search experiment using Skipper-CCDs with a total active mass of 10 kg. As part of the detector development, the collaboration plans to build the Oscura Integration Test (OIT), an engineering test with 10% of the total mass. Here we discuss the early science opportunities with the OIT to search for millicharged particles (mCPs) using the NuMI beam at Fermilab. mCPs would be produced at low energies through photon-mediated processes from decays of scalar, pseudoscalar, and vector mesons, or direct Drell-Yan productions. Estimates show that the OIT would be a world-leading probe for mCPs in the MeV mass range., Comment: 13 pages, 10 figures

Oscura is a proposed multi-kg skipper-CCD experiment designed for a dark matter (DM) direct detection search that will reach unprecedented sensitivity to sub-GeV DM-electron interactions with its 10 kg detector array. Oscura is planning to operate at SNOLAB with 2070 m overburden, and aims to reach a background goal of less than one event in each electron bin in the 2-10 electron ionization-signal region for the full 30 kg-year exposure, with a radiation background rate of 0.01 dru. In order to achieve this goal, Oscura must address each potential source of background events, including instrumental backgrounds. In this work, we discuss the main instrumental background sources and the strategy to control them, establishing a set of constraints on the sensors' performance parameters. We present results from the tests of the first fabricated Oscura prototype sensors, evaluate their performance in the context of the established constraints and estimate the Oscura instrumental background based on these results.