Search

In the pursuit of understanding the population of stellar remnants within the Milky Way, we analyze the sample of $\sim 950$ microlensing events observed by the Spitzer Space Telescope between 2014 and 2019. In this study we focus on a sub-sample of nine microlensing events, selected based on their long timescales, small microlensing parallaxes and joint observations by the Gaia mission, to increase the probability that the chosen lenses are massive and the mass is measurable. Among the selected events we identify lensing black holes and neutron star candidates, with potential confirmation through forthcoming release of the Gaia time-series astrometry in 2026. Utilizing Bayesian analysis and Galactic models, along with the Gaia Data Release 3 proper motion data, four good candidates for dark remnants were identified: OGLE-2016-BLG-0293, OGLE-2018-BLG-0483, OGLE-2018-BLG-0662, and OGLE-2015-BLG-0149, with lens masses of $2.98^{+1.75}_{-1.28}~M_{\odot}$, $4.65^{+3.12}_{-2.08}~M_{\odot}$, $3.15^{+0.66}_{-0.64}~M_{\odot}$ and $1.4^{+0.75}_{-0.55}~M_{\odot}$, respectively. Notably, the first two candidates are expected to exhibit astrometric microlensing signals detectable by Gaia, offering the prospect of validating the lens masses. The methodologies developed in this work will be applied to the full Spitzer microlensing sample, populating and analyzing the time-scale ($t_{\rm E}$) vs. parallax ($\pi_{\rm E}$) diagram to derive constraints on the population of lenses in general and massive remnants in particular., Comment: submitted to ApJ

This article presents constraints on dark-matter-electron interactions obtained from the first underground data-taking campaign with multiple SuperCDMS HVeV detectors operated in the same housing. An exposure of 7.63 g-days is used to set upper limits on the dark-matter-electron scattering cross section for dark matter masses between 0.5 and 1000 MeV/$c^2$, as well as upper limits on dark photon kinetic mixing and axion-like particle axioelectric coupling for masses between 1.2 and 23.3 eV/$c^2$. Compared to an earlier HVeV search, sensitivity was improved as a result of an increased overburden of 225 meters of water equivalent, an anticoincidence event selection, and better pile-up rejection. In the case of dark-matter-electron scattering via a heavy mediator, an improvement by up to a factor of 25 in cross-section sensitivity was achieved., Comment: 7 pages + title and references, 4 figures, and 1 table

Cryogenic calorimetric experiments to search for neutrinoless double-beta decay ($0\nu\beta\beta$) are highly competitive, scalable and versatile in isotope. The largest planned detector array, CUPID, is comprised of about 1500 individual Li$_2^{100}$MoO$_{4}$ detector modules with a further scale up envisioned for a follow up experiment (CUPID-1T). In this article, we present a novel detector concept targeting this second stage with a low impedance TES based readout for the Li$_2$MoO$_{4}$ absorber that is easily mass-produced and lends itself to a multiplexed readout. We present the detector design and results from a first prototype detector operated at the NEXUS shallow underground facility at Fermilab. The detector is a 2-cm-side cube with 21$\,$g mass that is strongly thermally coupled to its readout chip to allow rise-times of $\sim$0.5$\,$ms. This design is more than one order of magnitude faster than present NTD based detectors and is hence expected to effectively mitigate backgrounds generated through the pile-up of two independent two neutrino decay events coinciding close in time. Together with a baseline resolution of 1.95$\,$keV (FWHM) these performance parameters extrapolate to a background index from pile-up as low as $5\cdot 10^{-6}\,$counts/keV/kg/yr in CUPID size crystals. The detector was calibrated up to the MeV region showing sufficient dynamic range for $0\nu\beta\beta$ searches. In combination with a SuperCDMS HVeV detector this setup also allowed us to perform a precision measurement of the scintillation time constants of Li$_2$MoO$_{4}$. The crystal showed a significant fast scintillation emission with O(10$\,\mu$s) time-scale, more than an order below the detector response of presently considered light detectors suggesting the possibility of further progress in pile-up rejection through better light detectors in the future.

We measure space- and time-correlated charge jumps on a four-qubit device, operating 107 meters below the Earth's surface in a low-radiation, cryogenic facility designed for the characterization of low-threshold particle detectors. The rock overburden of this facility reduces the cosmic ray muon flux by over 99% compared to laboratories at sea level. Combined with 4$\pi$ coverage of a movable lead shield, this facility enables quantifiable control over the flux of ionizing radiation on the qubit device. Long-time-series charge tomography measurements on these weakly charge-sensitive qubits capture discontinuous jumps in the induced charge on the qubit islands, corresponding to the interaction of ionizing radiation with the qubit substrate. The rate of these charge jumps scales with the flux of ionizing radiation on the qubit package, as characterized by a series of independent measurements on another energy-resolving detector operating simultaneously in the same cryostat with the qubits. Using lead shielding, we achieve a minimum charge jump rate of 0.19$^{+0.04}_{-0.03}$ mHz, almost an order of magnitude lower than that measured in surface tests, but a factor of roughly eight higher than expected based on reduction of ambient gammas alone. We operate four qubits for over 22 consecutive hours with zero correlated charge jumps at length scales above three millimeters., Comment: 12 pages, 6 figures, 4 tables. Minor update to the measured gamma flux ratio (Page 4 and Supplemental Section F) in the LMO detector, from 23 to 20. Typos corrected, references added. Extraneous .tex files have been removed that were causing errors with the "HTML (experimental)" arxiv feature

We present pyLIMASS, a novel algorithm for estimating the physical properties of the lensing system in microlensing events. The main idea of pyLIMASS is to combine all available information regarding the microlensing event, defined as observables, and to estimate the parameter distributions of the system, such as the lens mass and distance. The algorithm is based on isochrones for the stars model and combine the observables using a Gaussian Mixtures approach. After describing the mathematical formalism and its implementation, we discuss the algorithm's performance on simulated and published events. Generally, the pyLIMASS estimations are in good agreement (i.e., within 1-$\sigma$) with the results of the selected published events, making it an effective tool to estimate the lens properties and their distribution. The applicability of the method was tested by using a catalog of realistically simulated events that could be observed by the future Galactic Bulge Time Domain Survey of the Nancy Grace Roman Space Telescope. By solely using constraints from the Roman lightcurves and images, pyLIMASS estimates the masses of the lens of the Roman catalog with a median precision of 20% with almost no bias., Comment: Accepted to AJ

We present the analysis of the microlensing event OGLE-2015-BLG-0845, which was affected by both the microlensing parallax and xallarap effects. The former was detected via the simultaneous observations from the ground and Spitzer, and the latter was caused by the orbital motion of the source star in a relatively close binary. The combination of these two effects led to a mass measurement of the lens object, revealing a low-mass ($0.14 \pm 0.05 M_{\odot}$) M-dwarf at the bulge distance ($7.6 \pm 1.0$ kpc). The source binary consists of a late F-type subgiant and a K-type dwarf of $\sim1.2 M_{\odot}$ and $\sim 0.9 M_{\odot}$, respectively, and the orbital period is $70 \pm 10$ days. OGLE-2015-BLG-0845 is the first single-lens event in which the lens mass is measured via the binarity of the source. Given the abundance of binary systems as potential microlensing sources, the xallarap effect may not be a rare phenomenon. Our work thus highlights the application of the xallarap effect in the mass determination of microlenses, and the same method can be used to identify isolated dark lenses., Comment: New version after the review process. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Recently Ahmadi et al. (2021) and Tagliaferro (2022) proposed some iterative methods for the numerical solution of linear systems which, under the classical hypothesis of strict diagonal dominance, typically converge faster than the Jacobi method, but slower than the forward/backward Gauss-Seidel one. In this paper we introduce a general class of iterative methods, based on suitable splittings of the matrix that defines the system, which include all of the methods mentioned above and have the same cost per iteration in a sequential computation environment. We also introduce a partial order relation in the set of the splittings and, partly theoretically and partly on the basis of a certain number of examples, we show that such partial order is typically connected to the speed of convergence of the corresponding methods. We pay particular attention to the case of linear systems for which the Jacobi iteration matrix is nonnegative, in which case we give a rigorous proof of the correspondence between the partial order relation and the magnitude of the spectral radius of the iteration matrices. Within the considered general class, some new specific promising methods are proposed as well.

Various dark matter search experiments employ phonon-based crystal detectors operated at cryogenic temperatures. Some of these detectors, including certain silicon detectors used by the SuperCDMS Collaboration, are able to achieve single-charge sensitivity when a voltage bias is applied across the detector. The total amount of phonon energy measured by such a detector is proportional to the number of electron-hole pairs created by the interaction. However, crystal impurities and surface effects can cause propagating charges to either become trapped inside the crystal or create additional unpaired charges, producing non-quantized measured energy as a result. A new analytical model for describing these detector response effects in phonon-based crystal detectors is presented. This model improves upon previous versions by demonstrating how the detector response, and thus the measured energy spectrum, is expected to differ depending on the source of events. We use this model to extract detector response parameters for SuperCDMS HVeV detectors, and illustrate how this robust modelling can help statistically discriminate between sources of events in order to improve the sensitivity of dark matter search experiments., Comment: 19 pages, 7 figures

Microcalorimeters that leverage microwave kinetic inductance detectors to read out phonon signals in the particle-absorbing target, referred to as kinetic inductance phonon-mediated (KIPM) detectors, offer an attractive detector architecture to probe dark matter (DM) down to the fermionic thermal relic mass limit. A prototype KIPM detector featuring a single aluminum resonator patterned onto a 1-gram silicon substrate was operated in the NEXUS low-background facility at Fermilab for characterization and evaluation of this detector architecture's efficacy for a dark matter search. An energy calibration was performed by exposing the bare substrate to a pulsed source of 470 nm photons, resulting in a baseline resolution on the energy absorbed by the phonon sensor of $2.1\pm0.2$ eV, a factor of two better than the current state-of-the-art, enabled by millisecond-scale quasiparticle lifetimes. However, due to the sub-percent phonon collection efficiency, the resolution on energy deposited in the substrate is limited to $\sigma_E=318 \pm 28$ eV. We further model the signal pulse shape as a function of device temperature to extract quasiparticle lifetimes, as well as the observed noise spectra, both of which impact the baseline resolution of the sensor.

Coherent elastic neutrino-nucleus scattering (CE$\nu$NS) offers a valuable approach in searching for physics beyond the Standard Model. The Ricochet experiment aims to perform a precision measurement of the CE$\nu$NS spectrum at the Institut Laue-Langevin nuclear reactor with cryogenic solid-state detectors. The experiment plans to employ an array of cryogenic thermal detectors, each with a mass around 30 g and an energy threshold of sub-100 eV. The array includes nine detectors read out by Transition-Edge Sensors (TES). These TES based detectors will also serve as demonstrators for future neutrino experiments with thousands of detectors. In this article we present an update in the characterization and modeling of a prototype TES detector., Comment: Submitted to LTD20 proceeding

This paper deals with the error analysis of the trapezoidal rule for the computation of Fourier type integrals, based on two double exponential transformations. The theory allows to construct algorithms in which the steplength and the number of nodes can be a priori selected. The analysis is also used to design an automatic integrator that can be employed without any knowledge of the function involved in the problem. Several numerical examples, which confirm the reliability of this strategy, are reported.

We complete the analysis of planetary candidates found by the KMT AnomalyFinder for the 2017 prime fields that cover $\sim 13\,{\rm deg}^2$. We report 3 unambiguous planets: OGLE-2017-BLG-0640, OGLE-2017-BLG-1275, and OGLE-2017-BLG-1237. The first two of these were not previously identified, while the last was not previously published due to technical complications induced by a nearby variable. We further report that a fourth anomalous event, the previously recognized OGLE-2017-BLG-1777, is very likely to be planetary, although its light curve requires unusually complex modeling because the lens and source both have orbiting companions. One of the 3 unambiguous planets, OGLE-2017-BLG-1275 is the first AnomalyFinder discovery that has a {\it Spitzer} microlens parallax measurement, $\pi_E \sim 0.045\pm0.015$, implying that this planetary system almost certainly lies in the Galactic bulge. In the order listed, the four planetary events have planet-host mass ratios $q$, and normalized projected separations $s$, of $(\log q,s)$ = $(-2.31,0.61)$, $(-2.06,0.63/1.09)$, $(-2.10,1.04)$, and $(-2.86,0.72)$. Combined with previously published events, the 2017 AnomalyFinder prime fields contain 11 unambiguous planets with well-measured $q$ and one very likely candidate, of which 3 are AnomalyFinder discoveries. In addition to these 12, there are three other unambiguous planets with large uncertainties in $q$., Comment: 67 pages, 13 figures, 16 tables

The future Ricochet experiment aims to search for new physics in the electroweak sector by measuring the Coherent Elastic Neutrino-Nucleus Scattering process from reactor antineutrinos with high precision down to the sub-100 eV nuclear recoil energy range. While the Ricochet collaboration is currently building the experimental setup at the reactor site, it is also finalizing the cryogenic detector arrays that will be integrated into the cryostat at the Institut Laue Langevin in early 2024. In this paper, we report on recent progress from the Ge cryogenic detector technology, called the CryoCube. More specifically, we present the first demonstration of a 30~eVee (electron equivalent) baseline ionization resolution (RMS) achieved with an early design of the detector assembly and its dedicated High Electron Mobility Transistor (HEMT) based front-end electronics. This represents an order of magnitude improvement over the best ionization resolutions obtained on similar heat-and-ionization germanium cryogenic detectors from the EDELWEISS and SuperCDMS dark matter experiments, and a factor of three improvement compared to the first fully-cryogenic HEMT-based preamplifier coupled to a CDMS-II germanium detector. Additionally, we discuss the implications of these results in the context of the future Ricochet experiment and its expected background mitigation performance., Comment: 10 pages, 5 figures, 1 table

Coherent elastic neutrino-nucleus scattering (CE$\nu$NS) offers valuable sensitivity to physics beyond the Standard Model. The Ricochet experiment will use cryogenic solid-state detectors to perform a precision measurement of the CE$\nu$NS spectrum induced by the high neutrino flux from the Institut Laue-Langevin nuclear reactor. The experiment will employ an array of detectors, each with a mass of $\sim$30 g and a targeted energy threshold of 50 eV. Nine of these detectors (the "Q-Array") will be based on a novel Transition-Edge Sensor (TES) readout style, in which the TES devices are thermally coupled to the absorber using a gold wire bond. We present initial characterization of a Q-Array-style detector using a 1 gram silicon absorber, obtaining a baseline root-mean-square resolution of less than 40 eV., Comment: 18 pages, 8 figures

We systematically investigate the claim by Vandorou et al. (2023) to have detected the host star of the low mass-ratio ($q<10^{-4}$) microlensing planet OGLE-2016-BLG-1195Lb, via Keck adaptive optics (AO) measurements $\Delta t=4.12\,$yr after the peak of the event ($t_0$). If correct, this measurement would contradict the microlens parallax measurement derived from Spitzer observations in solar orbit taken near $t_0$. We show that this host identification would be in $4\,\sigma$ conflict with the original ground-based lens-source relative proper-motion measurements. By contrast, Gould (2022) estimated a probability $p=10\%$ that the ``other star'' resolved by single-epoch late-time AO would be a companion to the host or the microlensed source, which is much more probable than a 4$\,\sigma$ statistical fluctuation. In addition, independent of this proper-motion discrepancy, the kinematics of this host-identification are substantially less probable than those of the Spitzer solution. Hence, this identification should not be accepted, pending additional observations that would either confirm or contradict it, which could be taken in 2023. Motivated by this tension, we present two additional investigations. We explore the possibility that Vandorou et al. (2023) identified the wrong ``star'' (or stellar asterism) on which to conduct their analysis. We find that astrometry of KMT and Keck images favors a star (or asterism) lying about 175 mas northwest of the one that they chose. We also present event parameters from a combined fit to all survey data, which yields, in particular, a more precise mass ratio, $q=(4.6\pm 0.4)\times 10^{-5}$. Finally, we discuss the broader implications of minimizing such false positives for the first measurement of the planet mass function, which will become possible when AO on next-generation telescopes are applied to microlensing planets., Comment: 31 pages, 7 figures, 3 tables, 2 appendices, submitted to AAS Journals

We measured the nuclear--recoil ionization yield in silicon with a cryogenic phonon-sensitive gram-scale detector. Neutrons from a mono-energetic beam scatter off of the silicon nuclei at angles corresponding to energy depositions from 4\,keV down to 100\,eV, the lowest energy probed so far. The results show no sign of an ionization production threshold above 100\,eV. These results call for further investigation of the ionization yield theory and a comprehensive determination of the detector response function at energies below the keV scale.

This paper deals with the computation of the Lerch transcendent by means of the Gauss-Laguerre formula. An a priori estimate of the quadrature error, that allows to compute the number of quadrature nodes necessary to achieve an arbitrary precision, is derived. Exploiting the properties of the Gauss-Laguerre rule and the error estimate, a truncated approach is also considered. The algorithm used and its Matlab implementation are reported. The numerical examples confirm the reliability of this approach.

We present a new analysis of previously published of SuperCDMS data using a profile likelihood framework to search for sub-GeV dark matter (DM) particles through two inelastic scattering channels: bremsstrahlung radiation and the Migdal effect. By considering these possible inelastic scattering channels, experimental sensitivity can be extended to DM masses that are undetectable through the DM-nucleon elastic scattering channel, given the energy threshold of current experiments. We exclude DM masses down to $220~\textrm{MeV}/c^2$ at $2.7 \times 10^{-30}~\textrm{cm}^2$ via the bremsstrahlung channel. The Migdal channel search provides overall considerably more stringent limits and excludes DM masses down to $30~\textrm{MeV}/c^2$ at $5.0 \times 10^{-30}~\textrm{cm}^2$., Comment: Submitted to PRD

We report the analysis of microlensing event OGLE-2017-BLG-1038, observed by the Optical Gravitational Lensing Experiment, Korean Microlensing Telescope Network, and Spitzer telescopes. The event is caused by a giant source star in the Galactic Bulge passing over a large resonant binary lens caustic. The availability of space-based data allows the full set of physical parameters to be calculated. However, there exists an eightfold degeneracy in the parallax measurement. The four best solutions correspond to very-low-mass binaries near ($M_1 = 170^{+40}_{-50} M_J$ and $M_2 = 110^{+20}_{-30} M_J$), or well below ($M_1 = 22.5^{+0.7}_{-0.4} M_J$ and $M_2 = 13.3^{+0.4}_{-0.3} M_J$) the boundary between stars and brown dwarfs. A conventional analysis, with scaled uncertainties for Spitzer data, implies a very-low-mass brown dwarf binary lens at a distance of 2 kpc. Compensating for systematic Spitzer errors using a Gaussian process model suggests that a higher mass M-dwarf binary at 6 kpc is equally likely. A Bayesian comparison based on a galactic model favors the larger-mass solutions. We demonstrate how this degeneracy can be resolved within the next ten years through infrared adaptive-optics imaging with a 40 m class telescope., Comment: 20 pages, 11 figures, 4 tables

G4CMP simulates phonon and charge transport in cryogenic semiconductor crystals using the Geant4 toolkit. The transport code is capable of simulating the propagation of acoustic phonons as well as electron and hole charge carriers. Processes for anisotropic phonon propagation, oblique charge-carrier propagation, and phonon emission by accelerated charge carriers are included. The simulation reproduces theoretical predictions and experimental observations such as phonon caustics, heat-pulse propagation times, and mean charge-carrier drift velocities. In addition to presenting the physics and features supported by G4CMP, this report outlines example applications from the dark matter and quantum information science communities. These communities are applying G4CMP to model and design devices for which the energy transported by phonons and charge carriers is germane to the performance of superconducting instruments and circuits placed on silicon and germanium substrates. The G4CMP package is available to download from GitHub: github.com/kelseymh/G4CMP., Comment: 21 pages, 11 figures, 10 tables

This paper deals with the solution of Maxwell's equations to model the electromagnetic fields in the case of a layered earth. The integrals involved in the solution are approximated by means of a novel approach based on the splitting of the reflection term. The inverse problem, consisting in the computation of the unknown underground conductivity distribution from a set of modeled magnetic field components, is also considered. Two optimization algorithms are applied, based on line- and global-search methods, and a new minimization approach is presented. Several EM surveys from the ground surface are simulated, considering the horizontal coplanar (HCP) and perpendicular (PRP) magnetic dipolar configurations. The numerical experiments, carried out for the study of river-levees integrity, allowed to estimate the errors associated to these kind of investigations, and confirm the reliability of the technique.

This paper introduces a very fast method for the computation of the resolvent of fractional powers of operators. The analysis is kept in the continuous setting of (potentially unbounded) self adjoint positive operators in Hilbert spaces. The method is based on the Gauss-Laguerre rule, exploiting a particular integral representation of the resolvent. We provide sharp error estimates that can be used to a priori select the number of nodes to achieve a prescribed tolerance.

Human anatomy, morphology, and associated diseases can be studied using medical imaging data. However, access to medical imaging data is restricted by governance and privacy concerns, data ownership, and the cost of acquisition, thus limiting our ability to understand the human body. A possible solution to this issue is the creation of a model able to learn and then generate synthetic images of the human body conditioned on specific characteristics of relevance (e.g., age, sex, and disease status). Deep generative models, in the form of neural networks, have been recently used to create synthetic 2D images of natural scenes. Still, the ability to produce high-resolution 3D volumetric imaging data with correct anatomical morphology has been hampered by data scarcity and algorithmic and computational limitations. This work proposes a generative model that can be scaled to produce anatomically correct, high-resolution, and realistic images of the human brain, with the necessary quality to allow further downstream analyses. The ability to generate a potentially unlimited amount of data not only enables large-scale studies of human anatomy and pathology without jeopardizing patient privacy, but also significantly advances research in the field of anomaly detection, modality synthesis, learning under limited data, and fair and ethical AI. Code and trained models are available at: https://github.com/AmigoLab/SynthAnatomy., Comment: 13 pages, 3 figures, 2 tables, accepted at SASHIMI MICCAI 2022

Neutrinoless Double-Beta decay (0$\nu\beta\beta$) processes sample a wide range of intermediate forbidden nuclear transitions, which may be impacted by quenching of the axial vector coupling constant ($g_A/g_V$), the uncertainty of which plays a pivotal role in determining the sensitivity reach of 0$\nu\beta\beta$ experiments. In this Letter, we present measurements performed on a high-resolution LiInSe$_{2}$~ bolometer in a ''source=detector'' configuration to measure the spectral shape of the 4-fold forbidden $\beta$-decay of $^{115}$In. The value of $g_A/g_V$ is determined by comparing the spectral shape of theoretical predictions to the experimental $\beta$ spectrum taking into account various simulated background components as well as a variety of detector effects. We find evidence of quenching of $g_A/g_V$ at $>5\sigma$ with a model-dependent quenching factor of $0.655\pm0.002$ as compared to the free-nucleon value for the Interacting Shell Model. We also measured the $^{115}$In half-life to be [$5.18\pm0.06(\text{stat.})^{+0.005}_{-0.015}(\text{sys.})]\times{10}^{14}$ yr within the Interacting Shell Model framework. This work demonstrates the power of the bolometeric technique to perform precision nuclear physics single-$\beta$ decay measurements, which can help reduce the uncertainties in the calculation of $0\nu\beta\beta$ nuclear matrix elements., Comment: 7 Pages, 4 Figures, 2 Tables, Submitted to PRL

CDMSlite Run 2 was a search for weakly interacting massive particles (WIMPs) with a cryogenic 600 g Ge detector operated in a high-voltage mode to optimize sensitivity to WIMPs of relatively low mass from 2 - 20 GeV/$c^2$. In this article, we present an effective field theory (EFT) analysis of the CDMSlite Run 2 data using an extended energy range and a comprehensive treatment of the expected background. A binned likelihood Bayesian analysis was performed on the recoil energy data, taking into account the parameters of the EFT interactions and optimizing the data selection with respect to the dominant background components. Energy regions within 5$\sigma$ of known activation peaks were removed from the analysis. The Bayesian evidences resulting from the different operator hypotheses show that the CDMSlite Run 2 data are consistent with the background-only models and do not allow for a signal interpretation assuming any additional EFT interaction. Consequently, upper limits on the WIMP mass and coupling-coefficient amplitudes and phases are presented for each EFT operator. These limits improve previous CDMSlite Run 2 bounds for WIMP masses above 5 GeV/$c^2$., Comment: 16 pages, 8 figures

Recent experiments searching for sub-GeV/$c^2$ dark matter have observed event excesses close to their respective energy thresholds. Although specific to the individual technologies, the measured excess event rates have been consistently reported at or below event energies of a few-hundred eV, or with charges of a few electron-hole pairs. In the present work, we operated a 1-gram silicon SuperCDMS-HVeV detector at three voltages across the crystal (0 V, 60 V and 100 V). The 0 V data show an excess of events in the tens of eV region. Despite this event excess, we demonstrate the ability to set a competitive exclusion limit on the spin-independent dark matter--nucleon elastic scattering cross section for dark matter masses of $\mathcal{O}(100)$ MeV/$c^2$, enabled by operation of the detector at 0 V potential and achievement of a very low $\mathcal{O}(10)$ eV threshold for nuclear recoils. Comparing the data acquired at 0 V, 60 V and 100 V potentials across the crystal, we investigated possible sources of the unexpected events observed at low energy. The data indicate that the dominant contribution to the excess is consistent with a hypothesized luminescence from the printed circuit boards used in the detector holder.

We complete the analysis of all 2018 prime-field microlensing planets identified by the KMTNet AnomalyFinder. Among the 10 previously unpublished events with clear planetary solutions, 8 are either unambiguously planetary or are very likely to be planetary in nature: OGLE-2018-BLG-1126, KMT-2018-BLG-2004, OGLE-2018-BLG-1647, OGLE-2018-BLG-1367, OGLE-2018-BLG-1544, OGLE-2018-BLG-0932, OGLE-2018-BLG-1212, and KMT-2018-BLG-2718. Combined with the 4 previously published new AnomalyFinder events and 12 previously published (or in preparation) planets that were discovered by eye, thismakes a total of 24 2018 prime-field planets discovered or recovered by AnomalyFinder. Together with a paper in preparation on 2018 sub-prime planets, this work lays the basis for the first statistical analysis of the planet mass-ratio function based on planets identified in KMTNet data. By systematically applying the heuristic analysis of Hwang et al. (2022) to each event, we identify the small modification in their formalism that is needed to unify the so-called close/wide and inner/outer degeneracies, as conjectured by, Comment: 22 pages, 14 tables, 15 figures

We present K2-2016-BLG-0005Lb, a densely sampled, planetary binary caustic-crossing microlensing event found from a blind search of data gathered from Campaign 9 of the Kepler K2 mission (K2C9). K2-2016-BLG-0005Lb is the first bound microlensing exoplanet discovered from space-based data. The event has caustic entry and exit points that are resolved in the K2C9 data, enabling the lens--source relative proper motion to be measured. We have fitted a binary microlens model to the Kepler data, and to simultaneous observations from multiple ground-based surveys. Whilst the ground-based data only sparsely sample the binary caustic, they provide a clear detection of parallax that allows us to break completely the microlensing mass--position--velocity degeneracy and measure the planet's mass directly. We find a host mass of $0.58\pm0.04 ~{\rm M}_\odot$ and a planetary mass of $1.1\pm0.1 ~{\rm M_J}$. The system lies at a distance of $5.2\pm0.2~$kpc from Earth towards the Galactic bulge, more than twice the distance of the previous most distant planet found by Kepler. The sky-projected separation of the planet from its host is found to be $4.2\pm0.3~$au which, for circular orbits, deprojects to a host separation $a = 4.4^{+1.9}_{-0.4}~$au and orbital period $P = 13^{+9}_{-2}~$yr. This makes K2-2016-BLG-0005Lb a close Jupiter analogue orbiting a low-mass host star. According to current planet formation models, this system is very close to the host mass threshold below which Jupiters are not expected to form. Upcoming space-based exoplanet microlensing surveys by NASA's Nancy Grace Roman Space Telescope and, possibly, ESA's Euclid mission, will provide demanding tests of current planet formation models., Comment: 18 pages. Accepted for publication in MNRAS Published by Oxford University Press on behalf of the Royal Astronomical Society

The SuperCDMS Collaboration is currently building SuperCDMS SNOLAB, a dark matter search focused on nucleon-coupled dark matter in the 1-5 GeV/c$^2$ mass range. Looking to the future, the Collaboration has developed a set of experience-based upgrade scenarios, as well as novel directions, to extend the search for dark matter using the SuperCDMS technology in the SNOLAB facility. The experienced-based scenarios are forecasted to probe many square decades of unexplored dark matter parameter space below 5 GeV/c$^2$, covering over 6 decades in mass: 1-100 eV/c$^2$ for dark photons and axion-like particles, 1-100 MeV/c$^2$ for dark-photon-coupled light dark matter, and 0.05-5 GeV/c$^2$ for nucleon-coupled dark matter. They will reach the neutrino fog in the 0.5-5 GeV/c$^2$ mass range and test a variety of benchmark models and sharp targets. The novel directions involve greater departures from current SuperCDMS technology but promise even greater reach in the long run, and their development must begin now for them to be available in a timely fashion. The experienced-based upgrade scenarios rely mainly on dramatic improvements in detector performance based on demonstrated scaling laws and reasonable extrapolations of current performance. Importantly, these improvements in detector performance obviate significant reductions in background levels beyond current expectations for the SuperCDMS SNOLAB experiment. Given that the dominant limiting backgrounds for SuperCDMS SNOLAB are cosmogenically created radioisotopes in the detectors, likely amenable only to isotopic purification and an underground detector life-cycle from before crystal growth to detector testing, the potential cost and time savings are enormous and the necessary improvements much easier to prototype., Comment: contribution to Snowmass 2021; v2 updated (assorted corrections and improvements to forecasts) October 2022; v3 updated (corrected SuperCDMS SNOLAB sensitivity curves in upgrade forecast plots in body of text) April 2023

Context. Brown dwarfs are poorly understood transition objects between stars and planets, with several competing mechanisms having been proposed for their formation. Mass measurements are generally difficult for isolated objects but also for brown dwarfs orbiting low-mass stars, which are often too faint for spectroscopic follow-up. Aims. Microlensing provides an alternative tool for the discovery and investigation of such faint systems. Here we present the analysis of the microlensing event OGLE-2019-BLG-0033/MOA-2019-BLG-035, which is due to a binary system composed of a brown dwarf orbiting a red dwarf. Methods. Thanks to extensive ground observations and the availability of space observations from Spitzer, it has been possible to obtain accurate estimates of all microlensing parameters, including parallax, source radius and orbital motion of the binary lens. Results. After accurate modeling, we find that the lens is composed of a red dwarf with mass $M_1 = 0.149 \pm 0.010M_\odot$ and a brown dwarf with mass $M_2 = 0.0463 \pm 0.0031M_\odot$, at a projected separation of $a_\perp = 0.585$ au. The system has a peculiar velocity that is typical of old metal-poor populations in the thick disk. Percent precision in the mass measurement of brown dwarfs has been achieved only in a few microlensing events up to now, but will likely become common with the Roman space telescope., Comment: 13 pages, 6 figures

In this paper, we present a re-analysis of SuperCDMS data using a profile likelihood approach to search for sub-GeV dark matter particles (DM) through two inelastic scattering channels: bremsstrahlung radiation and the Migdal effect. By considering possible inelastic scattering channels, experimental sensitivity can be extended to DM masses that would otherwise be undetectable through the DM-nucleon elastic scattering channel, given the energy threshold of current experiments. We exclude DM masses down to $220~\textrm{MeV}/c^2$ at $2.7 \times 10^{-30}~\textrm{cm}^2$ via the bremsstrahlung channel. The Migdal channel search excludes DM masses down to $30~\textrm{MeV}/c^2$ at $5.0 \times 10^{-30}~\textrm{cm}^2$., Comment: This paper is being withdrawn due to an error in data selection during the analysis. Although incorrect, the limits are roughly representative of the sensitivity. The new corrected version of the result will be uploaded once ready

Two photo-neutron sources, $^{88}$Y$^{9}$Be and $^{124}$Sb$^{9}$Be, have been used to investigate the ionization yield of nuclear recoils in the CDMSlite germanium detectors by the SuperCDMS collaboration. This work evaluates the yield for nuclear recoil energies between 1 keV and 7 keV at a temperature of $\sim$ 50 mK. We use a Geant4 simulation to model the neutron spectrum assuming a charge yield model that is a generalization of the standard Lindhard model and consists of two energy dependent parameters. We perform a likelihood analysis using the simulated neutron spectrum, modeled background, and experimental data to obtain the best fit values of the yield model. The ionization yield between recoil energies of 1 keV and 7 keV is shown to be significantly lower than predicted by the standard Lindhard model for germanium. There is a general lack of agreement among different experiments using a variety of techniques studying the low-energy range of the nuclear recoil yield, which is most critical for interpretation of direct dark matter searches. This suggests complexity in the physical process that many direct detection experiments use to model their primary signal detection mechanism and highlights the need for further studies to clarify underlying systematic effects that have not been well understood up to this point.

Many low-threshold experiments observe sharply rising event rates of yet unknown origins below a few hundred eV, and larger than expected from known backgrounds. Due to the significant impact of this excess on the dark matter or neutrino sensitivity of these experiments, a collective effort has been started to share the knowledge about the individual observations. For this, the EXCESS Workshop was initiated. In its first iteration in June 2021, ten rare event search collaborations contributed to this initiative via talks and discussions. The contributing collaborations were CONNIE, CRESST, DAMIC, EDELWEISS, MINER, NEWS-G, NUCLEUS, RICOCHET, SENSEI and SuperCDMS. They presented data about their observed energy spectra and known backgrounds together with details about the respective measurements. In this paper, we summarize the presented information and give a comprehensive overview of the similarities and differences between the distinct measurements. The provided data is furthermore publicly available on the workshop's data repository together with a plotting tool for visualization., Comment: 44 pages, 20 figures; Editors: A. Fuss, M. Kaznacheeva, F. Reindl, F. Wagner; updated copyright statements and funding information

We consider the approximation of the inverse square root of regularly accretive operators in Hilbert spaces. The approximation is of rational type and comes from the use of the Gauss-Legendre rule applied to a special integral formulation of the problem. We derive sharp error estimates, based on the use of the numerical range, and provide some numerical experiments. For practical purposes, the finite dimensional case is also considered. In this setting, the convergence is shown to be of exponential type.

We report the first unambiguous detection and mass measurement of an isolated stellar-mass black hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source star of the long-duration (t_E~270 days), high-magnification microlensing event MOA-2011-BLG-191/OGLE-2011-BLG-0462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of the Galactic bulge. HST imaging, conducted at eight epochs over an interval of six years, reveals a clear relativistic astrometric deflection of the background star's apparent position. Ground-based photometry of MOA-11-191/OGLE-11-462 shows a parallactic signature of the effect of the Earth's motion on the microlensing light curve. Combining the HST astrometry with the ground-based light curve and the derived parallax, we obtain a lens mass of 7.1 +/- 1.3 Msun and a distance of 1.58 +/- 0.18 kpc. We show that the lens emits no detectable light, which, along with having a mass higher than is possible for a white dwarf or neutron star, confirms its BH nature. Our analysis also provides an absolute proper motion for the BH. The proper motion is offset from the mean motion of Galactic-disk stars at similar distances by an amount corresponding to a transverse space velocity of ~45 km/s, suggesting that the BH received a 'natal kick' from its supernova explosion. Previous mass determinations for stellar-mass BHs have come from radial-velocity measurements of Galactic X-ray binaries, and from gravitational radiation emitted by merging BHs in binary systems in external galaxies. Our mass measurement is the first for an isolated stellar-mass BH using any technique., Comment: 37 pages, Published in ApJ

OGLE-2016-BLG-1093 is a planetary microlensing event that is part of the statistical $Spitzer$ microlens parallax sample. The precise measurement of the microlens parallax effect for this event, combined with the measurement of finite source effects, leads to a direct measurement of the lens masses and system distance: $M_{\rm host} = 0.38$--$0.57\, M_{\odot}$, $m_p = 0.59$--$0.87\, M_{\rm Jup}$, and the system is located at the Galactic bulge ($D_L \sim 8.1$ kpc). Because this was a high-magnification event, we are also able to empirically show that the "cheap-space parallax" concept Gould & Yee (2012) produces well-constrained (and consistent) results for $|\pi_{\rm E}|$. This demonstrates that this concept can be extended to many two-body lenses. Finally, we briefly explore systematics in the $Spitzer$ light curve in this event and show that their potential impact is strongly mitigated by the color-constraint., Comment: 9 Figures, 3 Tables, submitted to the AAS journal

We present the photometric analysis of Gaia19bld, a high-magnification ($A\approx60$) microlensing event located in the southern Galactic plane, which exhibited finite source and microlensing parallax effects. Due to a prompt detection by the Gaia satellite and the very high brightness of $I = 9.05~$mag at the peak, it was possible to collect a complete and unique set of multi-channel follow-up observations, which allowed us to determine all parameters vital for the characterisation of the lens and the source in the microlensing event. Gaia19bld was discovered by the Gaia satellite and was subsequently intensively followed up with a network of ground-based observatories and the Spitzer Space Telescope. We collected multiple high-resolution spectra with Very Large Telescope (VLT)/X-Shooter to characterise the source star. The event was also observed with VLT Interferometer (VLTI)/PIONIER during the peak. Here we focus on the photometric observations and model the light curve composed of data from Gaia, Spitzer, and multiple optical, ground-based observatories. We find the best-fitting solution with parallax and finite source effects. We derived the limit on the luminosity of the lens based on the blended light model and spectroscopic distance. We compute the mass of the lens to be $1.13 \pm 0.03~M_{\odot}$ and derive its distance to be $5.52^{+0.35}_{-0.64}~\mathrm{kpc}$. The lens is likely a main sequence star, however its true nature has yet to be verified by future high-resolution observations. Our results are consistent with interferometric measurements of the angular Einstein radius, emphasising that interferometry can be a new channel for determining the masses of objects that would otherwise remain undetectable, including stellar-mass black holes., Comment: accepted to Astronomy&Astrophysics

In this work we develop the Gaussian quadrature rule for weight functions involving fractional powers, exponentials and Bessel functions of the first kind. Besides the computation based on the use of the standard and the modified Chebyshev algorithm, here we present a very stable algorithm based on the preconditioning of the moment matrix. Numerical experiments are provided and a geophysical application is considered.

The past few decades have seen major developments in the design and operation of cryogenic particle detectors. This technology offers an extremely good energy resolution - comparable to semiconductor detectors - and a wide choice of target materials, making low temperature calorimetric detectors ideal for a variety of particle physics applications. Rare event searches have continued to require ever greater exposures, which has driven them to ever larger cryogenic detectors, with the CUORE experiment being the first to reach a tonne-scale, mK-cooled, experimental mass. CUORE, designed to search for neutrinoless double beta decay, has been operational since 2017 at a temperature of about 10 mK. This result has been attained by the use of an unprecedentedly large cryogenic infrastructure called the CUORE cryostat: conceived, designed and commissioned for this purpose. In this article the main characteristics and features of the cryogenic facility developed for the CUORE experiment are highlighted. A brief introduction of the evolution of the field and of the past cryogenic facilities are given. The motivation behind the design and development of the CUORE cryogenic facility is detailed as are the steps taken toward realization, commissioning, and operation of the CUORE cryostat. The major challenges overcome by the collaboration and the solutions implemented throughout the building of the cryogenic facility will be discussed along with the potential improvements for future facilities. The success of CUORE has opened the door to a new generation of large-scale cryogenic facilities in numerous fields of science. Broader implications of the incredible feat achieved by the CUORE collaboration on the future cryogenic facilities in various fields ranging from neutrino and dark matter experiments to quantum computing will be examined., Comment: 45 pages, 14 figures

In this paper we are interested in the approximation of fractional powers of self-adjoint positive operators. Starting from the integral representation of the operators, we apply the trapezoidal rule combined with a single-exponential and a double-exponential transform of the integrand function. For the first approach our aim is only to review some theoretical aspects in order to refine the choice of the parameters that allow a faster convergence. As for the double exponential transform, in this work we show how to improve the existing error estimates for the scalar case and also extend the analysis to operators. We report some numerical experiments to show the reliability of the estimates obtained.

Swimming organisms can escape their predators by creating and harnessing unsteady flow fields through their body motions. Stochastic optimization and flow simulations have identified escape patterns that are consistent with those observed in natural larval swimmers. However, these patterns have been limited by the specification of a particular cost function and depend on a prescribed functional form of the body motion. Here, we deploy reinforcement learning to discover swimmer escape patterns for larval fish under energy constraints. The identified patterns include the C-start mechanism, in addition to more energetically efficient escapes. We find that maximizing distance with limited energy requires swimming via short bursts of accelerating motion interlinked with phases of gliding. The present, data efficient, reinforcement learning algorithm results in an array of patterns that reveal practical flow optimization principles for efficient swimming and the methodology can be transferred to the control of aquatic robotic devices operating under energy constraints.