Search

[Abridged] The Nancy Grace Roman Space Telescope will be one of several flagship survey facilities operating over the next decade starting $\sim$2025. The deep near-IR imaging that Roman will deliver will be highly complementary to the capabilities of other survey telescopes that will operate contemporaneously, particularly those that can provide data at different wavelengths and messengers, or different time intervals. Combining data from multiple facilities can provide important astrophysical insights, provided the data acquisition is carefully scheduled, and careful plans are made for appropriate joint data analyses. In this White Paper, we discuss the broad range of science that would be enabled by coordinating Roman observations of the Galactic Bulge with those of the Vera C. Rubin Observatory. Specifically, we discuss how Roman's characterization of lensing events caused by exoplanets, stellar systems and stellar remnants can be enhanced by data from Rubin. The same data will also be highly advantageous for the determination of stellar properties, and for distinguishing exoplanetary transits. It will enable more accurate period-color-luminosity relationships to be measured for RR~Lyrae throughout the Milky Way Bulge and Bar, probing galactic structure and dynamics. But we stress that this is only a sample of the full potential and advocate for a more complete study to be made as a joint effort between these major projects. We note that we do not suggest any changes beyond the established Science Requirements for the RGBTDS, in terms of survey footprint or filter selection., Comment: Submitted to the 2023 call for White Papers regarding the Roman Mission Core Community Surveys

Galactic science encompasses a wide range of subjects in the study of the Milky Way and Magellanic Clouds, from Young Stellar Objects to X-ray Binaries. Mapping these populations, and exploring transient phenomena within them, are among the primary science goals of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). While early versions of the survey strategy dedicated relatively few visits to the Galactic Plane region, more recent strategies under consideration envision higher cadence within selected regions of high scientific interest. The range of galactic science presents a challenge in evaluating which strategies deliver the highest scientific returns. Here we present metrics designed to evaluate Rubin survey strategy simulations based on the cadence of observations they deliver within regions of interest to different topics in galactic science, using variability categories defined by timescale. We also compare the fractions of exposures obtained in each filter with those recommended for the different science goals. We find that the baseline_v2.x simulations deliver observations of the high-priority regions at sufficiently high cadence to reliably detect variability on timescales >10 d or more. Follow-up observations may be necessary to properly characterize variability, especially transients, on shorter timescales. Combining the regions of interest for all the science cases considered, we identify those areas of the Galactic Plane and Magellanic Clouds of highest priority. We recommend that these refined survey footprints be used in future simulations to explore rolling cadence scenarios, and to optimize the sequence of observations in different bandpasses., Comment: Accepted for publication in ApJS. 26 pages and 36 figures

In this work, we expand on the XENON1T nuclear recoil searches to study the individual signals of dark matter interactions from operators up to dimension-eight in a Chiral Effective Field Theory (ChEFT) and a model of inelastic dark matter (iDM). We analyze data from two science runs of the XENON1T detector totaling 1\,tonne$\times$year exposure. For these analyses, we extended the region of interest from [4.9, 40.9]$\,$keV$_{\text{NR}}$ to [4.9, 54.4]$\,$keV$_{\text{NR}}$ to enhance our sensitivity for signals that peak at nonzero energies. We show that the data is consistent with the background-only hypothesis, with a small background over-fluctuation observed peaking between 20 and 50$\,$keV$_{\text{NR}}$, resulting in a maximum local discovery significance of 1.7\,$\sigma$ for the Vector$\otimes$Vector$_{\text{strange}}$ ($VV_s$) ChEFT channel for a dark matter particle of 70$\,$GeV/c$^2$, and $1.8\,\sigma$ for an iDM particle of 50$\,$GeV/c$^2$ with a mass splitting of 100$\,$keV/c$^2$. For each model, we report 90\,\% confidence level (CL) upper limits. We also report upper limits on three benchmark models of dark matter interaction using ChEFT where we investigate the effect of isospin-breaking interactions. We observe rate-driven cancellations in regions of the isospin-breaking couplings, leading to up to 6 orders of magnitude weaker upper limits with respect to the isospin-conserving case.

The XENON collaboration has published stringent limits on specific dark matter -nucleon recoil spectra from dark matter recoiling on the liquid xenon detector target. In this paper, we present an approximate likelihood for the XENON1T 1 tonne-year nuclear recoil search applicable to any nuclear recoil spectrum. Alongside this paper, we publish data and code to compute upper limits using the method we present. The approximate likelihood is constructed in bins of reconstructed energy, profiled along the signal expectation in each bin. This approach can be used to compute an approximate likelihood and therefore most statistical results for any nuclear recoil spectrum. Computing approximate results with this method is approximately three orders of magnitude faster than the likelihood used in the original publications of XENON1T, where limits were set for specific families of recoil spectra. Using this same method, we include toy Monte Carlo simulation-derived binwise likelihoods for the upcoming XENONnT experiment that can similarly be used to assess the sensitivity to arbitrary nuclear recoil signatures in its eventual 20 tonne-year exposure., Comment: Accepted by European Physical Journal C

We report on a blinded analysis of low-energy electronic-recoil data from the first science run of the XENONnT dark matter experiment. Novel subsystems and the increased 5.9 tonne liquid xenon target reduced the background in the (1, 30) keV search region to $(15.8 \pm 1.3)$ events/(tonne$\times$year$\times$keV), the lowest ever achieved in a dark matter detector and $\sim$5 times lower than in XENON1T. With an exposure of 1.16 tonne-years, we observe no excess above background and set stringent new limits on solar axions, an enhanced neutrino magnetic moment, and bosonic dark matter., Comment: 10 pages, 6 figures. Data points are available in https://zenodo.org/record/7311940

We present results on the search for double-electron capture ($2\nu\text{ECEC}$) of $^{124}$Xe and neutrinoless double-$\beta$ decay ($0\nu\beta\beta$) of $^{136}$Xe in XENON1T. We consider captures from the K- up to the N-shell in the $2\nu\text{ECEC}$ signal model and measure a total half-life of $T_{1/2}^{2\nu\text{ECEC}}=(1.1\pm0.2_\text{stat}\pm0.1_\text{sys})\times 10^{22}\;\text{yr}$ with a $0.87\;\text{kg}\times\text{yr}$ isotope exposure. The statistical significance of the signal is $7.0\,\sigma$. We use XENON1T data with $36.16\;\text{kg}\times\text{yr}$ of $^{136}$Xe exposure to search for $0\nu\beta\beta$. We find no evidence of a signal and set a lower limit on the half-life of $T_{1/2}^{0\nu\beta\beta} > 1.2 \times 10^{24}\;\text{yr}\; \text{at}\; 90\,\%\;\text{CL}$. This is the best result from a dark matter detector without an enriched target to date. We also report projections on the sensitivity of XENONnT to $0\nu\beta\beta$. Assuming a $275\;\text{kg}\times\text{yr}$ $^{136}$Xe exposure, the expected sensitivity is $T_{1/2}^{0\nu\beta\beta} > 2.1 \times 10^{25}\;\text{yr}\; \text{at}\; 90\,\%\;\text{CL}$, corresponding to an effective Majorana mass range of $\langle m_{\beta\beta} \rangle < (0.19 - 0.59)\;\text{eV/c}^2$., Comment: 23 pages, 14 figures, 9 tables, version accepted for publication in Phys. Rev. C

A novel online distillation technique was developed for the XENON1T dark matter experiment to reduce intrinsic background components more volatile than xenon, such as krypton or argon, while the detector was operating. The method is based on a continuous purification of the gaseous volume of the detector system using the XENON1T cryogenic distillation column. A krypton-in-xenon concentration of $(360 \pm 60)$ ppq was achieved. It is the lowest concentration measured in the fiducial volume of an operating dark matter detector to date. A model was developed and fit to the data to describe the krypton evolution in the liquid and gas volumes of the detector system for several operation modes over the time span of 550 days, including the commissioning and science runs of XENON1T. The online distillation was also successfully applied to remove Ar-37 after its injection for a low energy calibration in XENON1T. This makes the usage of Ar-37 as a regular calibration source possible in the future. The online distillation can be applied to next-generation experiments to remove krypton prior to, or during, any science run. The model developed here allows further optimization of the distillation strategy for future large scale detectors.

Delayed single- and few-electron emissions plague dual-phase time projection chambers, limiting their potential to search for light-mass dark matter. This paper examines the origins of these events in the XENON1T experiment. Characterization of the intensity of delayed electron backgrounds shows that the resulting emissions are correlated, in time and position, with high-energy events and can effectively be vetoed. In this work we extend previous S2-only analyses down to a single electron. From this analysis, after removing the correlated backgrounds, we observe rates < 30 events/(electron*kg*day) in the region of interest spanning 1 to 5 electrons. We derive 90% confidence upper limits for dark matter-electron scattering, first direct limits on the electric dipole, magnetic dipole, and anapole interactions, and bosonic dark matter models, where we exclude new parameter space for dark photons and solar dark photons., Comment: 20 pages, 17 figures, Updated to correct published Solar Dark Photon limit

The selection of low-radioactive construction materials is of the utmost importance for rare-event searches and thus critical to the XENONnT experiment. Results of an extensive radioassay program are reported, in which material samples have been screened with gamma-ray spectroscopy, mass spectrometry, and $^{222}$Rn emanation measurements. Furthermore, the cleanliness procedures applied to remove or mitigate surface contamination of detector materials are described. Screening results, used as inputs for a XENONnT Monte Carlo simulation, predict a reduction of materials background ($\sim$17%) with respect to its predecessor XENON1T. Through radon emanation measurements, the expected $^{222}$Rn activity concentration in XENONnT is determined to be 4.2$\,(^{+0.5}_{-0.7})\,\mu$Bq/kg, a factor three lower with respect to XENON1T. This radon concentration will be further suppressed by means of the novel radon distillation system.

We suggest a new experiment sensitive to a possible difference between the amount of CP violation as measured on the surface of the Earth and in a lower gravity environment. Our proposed experiment is model independent and could yield a $5\sigma$ measurement within tens of days, indicating a dependence of the level of CP violation in the neutral kaon system on the local gravitational potential., Comment: Contribution to the 2021 Gravitation session of the 55th Rencontres de Moriond

The Muon g-2 Experiment at Fermi National Accelerator Laboratory (FNAL) has measured the muon anomalous precession frequency $\omega_a$ to an uncertainty of 434 parts per billion (ppb), statistical, and 56 ppb, systematic, with data collected in four storage ring configurations during its first physics run in 2018. When combined with a precision measurement of the magnetic field of the experiment's muon storage ring, the precession frequency measurement determines a muon magnetic anomaly of $a_{\mu}({\rm FNAL}) = 116\,592\,040(54) \times 10^{-11}$ (0.46 ppm). This article describes the multiple techniques employed in the reconstruction, analysis and fitting of the data to measure the precession frequency. It also presents the averaging of the results from the eleven separate determinations of \omega_a, and the systematic uncertainties on the result., Comment: 29 pages, 19 figures. Published in Physical Review D

The Fermi National Accelerator Laboratory has measured the anomalous precession frequency $a^{}_\mu = (g^{}_\mu-2)/2$ of the muon to a combined precision of 0.46 parts per million with data collected during its first physics run in 2018. This paper documents the measurement of the magnetic field in the muon storage ring. The magnetic field is monitored by nuclear magnetic resonance systems and calibrated in terms of the equivalent proton spin precession frequency in a spherical water sample at 34.7$^\circ$C. The measured field is weighted by the muon distribution resulting in $\tilde{\omega}'^{}_p$, the denominator in the ratio $\omega^{}_a$/$\tilde{\omega}'^{}_p$ that together with known fundamental constants yields $a^{}_\mu$. The reported uncertainty on $\tilde{\omega}'^{}_p$ for the Run-1 data set is 114 ppb consisting of uncertainty contributions from frequency extraction, calibration, mapping, tracking, and averaging of 56 ppb, and contributions from fast transient fields of 99 ppb., Comment: Added one citation and corrected missing normalization in Eqs (35) and (36)

We present the first results of the Fermilab Muon g-2 Experiment for the positive muon magnetic anomaly $a_\mu \equiv (g_\mu-2)/2$. The anomaly is determined from the precision measurements of two angular frequencies. Intensity variation of high-energy positrons from muon decays directly encodes the difference frequency $\omega_a$ between the spin-precession and cyclotron frequencies for polarized muons in a magnetic storage ring. The storage ring magnetic field is measured using nuclear magnetic resonance probes calibrated in terms of the equivalent proton spin precession frequency ${\tilde{\omega}'^{}_p}$ in a spherical water sample at 34.7$^{\circ}$C. The ratio $\omega_a / {\tilde{\omega}'^{}_p}$, together with known fundamental constants, determines $a_\mu({\rm FNAL}) = 116\,592\,040(54)\times 10^{-11}$ (0.46\,ppm). The result is 3.3 standard deviations greater than the standard model prediction and is in excellent agreement with the previous Brookhaven National Laboratory (BNL) E821 measurement. After combination with previous measurements of both $\mu^+$ and $\mu^-$, the new experimental average of $a_\mu({\rm Exp}) = 116\,592\,061(41)\times 10^{-11}$ (0.35\,ppm) increases the tension between experiment and theory to 4.2 standard deviations, Comment: 10 pages; 4 figures

This paper presents the beam dynamics systematic corrections and their uncertainties for the Run-1 data set of the Fermilab Muon g-2 Experiment. Two corrections to the measured muon precession frequency $\omega_a^m$ are associated with well-known effects owing to the use of electrostatic quadrupole (ESQ) vertical focusing in the storage ring. An average vertically oriented motional magnetic field is felt by relativistic muons passing transversely through the radial electric field components created by the ESQ system. The correction depends on the stored momentum distribution and the tunes of the ring, which has relatively weak vertical focusing. Vertical betatron motions imply that the muons do not orbit the ring in a plane exactly orthogonal to the vertical magnetic field direction. A correction is necessary to account for an average pitch angle associated with their trajectories. A third small correction is necessary because muons that escape the ring during the storage time are slightly biased in initial spin phase compared to the parent distribution. Finally, because two high-voltage resistors in the ESQ network had longer than designed RC time constants, the vertical and horizontal centroids and envelopes of the stored muon beam drifted slightly, but coherently, during each storage ring fill. This led to the discovery of an important phase-acceptance relationship that requires a correction. The sum of the corrections to $\omega_a^m$ is 0.50 $\pm$ 0.09 ppm; the uncertainty is small compared to the 0.43 ppm statistical precision of $\omega_a^m$., Comment: 35 pages, 29 figures. Accepted by Phys. Rev. Accel. Beams

We report on a search for nuclear recoil signals from solar $^8$B neutrinos elastically scattering off xenon nuclei in XENON1T data, lowering the energy threshold from 2.6 keV to 1.6 keV. We develop a variety of novel techniques to limit the resulting increase in backgrounds near the threshold. No significant $^8$B neutrino-like excess is found in an exposure of 0.6 t $\times$ y. For the first time, we use the non-detection of solar neutrinos to constrain the light yield from 1-2 keV nuclear recoils in liquid xenon, as well as non-standard neutrino-quark interactions. Finally, we improve upon world-leading constraints on dark matter-nucleus interactions for dark matter masses between 3 GeV/c$^2$ and 11 GeV/c$^2$ by as much as an order of magnitude.

We report the results of a search for the inelastic scattering of weakly interacting massive particles (WIMPs) in the XENON1T dark matter experiment. Scattering off $^{129}$Xe is the most sensitive probe of inelastic WIMP interactions, with a signature of a 39.6 keV de-excitation photon detected simultaneously with the nuclear recoil. Using an exposure of 0.89 tonne-years, we find no evidence of inelastic WIMP scattering with a significance of more than 2$\sigma$. A profile-likelihood ratio analysis is used to set upper limits on the cross-section of WIMP-nucleus interactions. We exclude new parameter space for WIMPs heavier than 100 GeV/c${}^2$, with the strongest upper limit of $3.3 \times 10^{-39}$ cm${}^2$ for 130 GeV/c${}^2$ WIMPs at 90\% confidence level., Comment: 9 pages, 4 figures

The selection of low-radioactive construction materials is of utmost importance for the success of low-energy rare event search experiments. Besides radioactive contaminants in the bulk, the emanation of radioactive radon atoms from material surfaces attains increasing relevance in the effort to further reduce the background of such experiments. In this work, we present the $^{222}$Rn emanation measurements performed for the XENON1T dark matter experiment. Together with the bulk impurity screening campaign, the results enabled us to select the radio-purest construction materials, targeting a $^{222}$Rn activity concentration of 10 $\mu$Bq/kg in 3.2 t of xenon. The knowledge of the distribution of the $^{222}$Rn sources allowed us to selectively eliminate critical components in the course of the experiment. The predictions from the emanation measurements were compared to data of the $^{222}$Rn activity concentration in XENON1T. The final $^{222}$Rn activity concentration of (4.5 $\pm$ 0.1) $\mu$Bq/kg in the target of XENON1T is the lowest ever achieved in a xenon dark matter experiment., Comment: 14 pages, 3 figures

Do external galaxies host planetary systems? Many lines of reasoning suggest that the answer must be 'yes'. In the foreseeable future, however, the question cannot be answered by the methods most successful in our own Galaxy. We report on a different approach which focuses on bright X-ray sources (XRSs). M51-ULS-1b is the first planet candidate to be found because it produces a full, short-lived eclipse of a bright XRS. M51-ULS-1b has a most probable radius slightly smaller than Saturn. It orbits one of the brightest XRSs in the external galaxy M51, the Whirlpool Galaxy, located 8.6 Megaparsecs from Earth. It is the first candidate for a planet in an external galaxy. The binary it orbits, M51-ULS-1, is young and massive. One of the binary components is a stellar remnant, either a neutron star (NS) or black hole (BH), and the other is a massive star. X-ray transits can now be used to discover more planets in external galaxies and also planets orbiting XRSs inside the Milky Way., Comment: 19 pages, 13 figures, 2 tables. Submitted

XENONnT is a dark matter direct detection experiment, utilizing 5.9 t of instrumented liquid xenon, located at the INFN Laboratori Nazionali del Gran Sasso. In this work, we predict the experimental background and project the sensitivity of XENONnT to the detection of weakly interacting massive particles (WIMPs). The expected average differential background rate in the energy region of interest, corresponding to (1, 13) keV and (4, 50) keV for electronic and nuclear recoils, amounts to $12.3 \pm 0.6$ (keV t y)$^{-1}$ and $(2.2\pm 0.5)\times 10^{-3}$ (keV t y)$^{-1}$, respectively, in a 4 t fiducial mass. We compute unified confidence intervals using the profile construction method, in order to ensure proper coverage. With the exposure goal of 20 t$\,$y, the expected sensitivity to spin-independent WIMP-nucleon interactions reaches a cross-section of $1.4\times10^{-48}$ cm$^2$ for a 50 GeV/c$^2$ mass WIMP at 90% confidence level, more than one order of magnitude beyond the current best limit, set by XENON1T. In addition, we show that for a 50 GeV/c$^2$ WIMP with cross-sections above $2.6\times10^{-48}$ cm$^2$ ($5.0\times10^{-48}$ cm$^2$) the median XENONnT discovery significance exceeds 3$\sigma$ (5$\sigma$). The expected sensitivity to the spin-dependent WIMP coupling to neutrons (protons) reaches $2.2\times10^{-43}$ cm$^2$ ($6.0\times10^{-42}$ cm$^2$).

We report results from searches for new physics with low-energy electronic recoil data recorded with the XENON1T detector. With an exposure of 0.65 t-y and an unprecedentedly low background rate of $76\pm2$ events/(t y keV) between 1 and 30 keV, the data enables sensitive searches for solar axions, an enhanced neutrino magnetic moment, and bosonic dark matter. An excess over known backgrounds is observed at low energies and most prominent between 2 and 3 keV. The solar axion model has a 3.4$\sigma$ significance, and a 3D 90% confidence surface is reported for axion couplings to electrons, photons, and nucleons. This surface is inscribed in the cuboid defined by $g_{ae}<3.8 \times 10^{-12}$, $g_{ae}g_{an}^{eff}<4.8\times 10^{-18}$, and $g_{ae}g_{a\gamma}<7.7\times10^{-22} GeV^{-1}$, and excludes either $g_{ae}=0$ or $g_{ae}g_{a\gamma}=g_{ae}g_{an}^{eff}=0$. The neutrino magnetic moment signal is similarly favored over background at 3.2$\sigma$ and a confidence interval of $\mu_{\nu} \in (1.4,2.9)\times10^{-11}\mu_B$ (90% C.L.) is reported. Both results are in strong tension with stellar constraints. The excess can also be explained by $\beta$ decays of tritium at 3.2$\sigma$ with a trace amount that can neither be confirmed nor excluded with current knowledge of its production and reduction mechanisms. The significances of the solar axion and neutrino magnetic moment hypotheses are reduced to 2.0$\sigma$ and 0.9$\sigma$, respectively, if an unconstrained tritium component is included in the fitting. With respect to bosonic dark matter, the excess favors a monoenergetic peak at ($2.3\pm0.2$) keV (68% C.L.) with a 3.0$\sigma$ global (4.0$\sigma$ local) significance. We also consider the possibility that $^{37}$Ar may be present in the detector and yield a 2.82 keV peak. Contrary to tritium, the $^{37}$Ar concentration can be tightly constrained and is found to be negligible., Comment: 26 pages, 15 figures. v2 added Ar37 background discussion and best-fit mass of bosonic dark matter, v3 updated Ar37 discussion, tritium estimation, and solar axion energy spectrum. Data in Fig. 2, 4, and 15, including unbinned energy points in Fig. 4, are available in 10.5281/zenodo.4088778

The Muon $g-2$ experiment, E989, is currently taking data at Fermilab with the aim of reducing the experimental error on the muon anomaly by a factor of four and possibly clarifying the current discrepancy with the theoretical prediction. A central component of this four-fold improvement in precision is the laser calibration system of the calorimeters, which has to monitor the gain variations of the photo-sensors with a 0.04\% precision on the short-term ($\sim 1\,$ms). This is about one order of magnitude better than what has ever been achieved for the calibration of a particle physics calorimeter. The system is designed to monitor also long-term gain variations, mostly due to temperature effects, with a precision below the per mille level. This article reviews the design, the implementation and the performance of the Muon $g-2$ laser calibration system, showing how the experimental requirements have been met., Comment: 33 pages,24 figures. Matches the published version

We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at https://doi.org/10.5281/zenodo.2612896. We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models.

NASA's WFIRST mission will perform a wide-field, NIR survey of the Galactic Bulge to search for exoplanets via the microlensing techniques. As the mission is due to launch in the mid-2020s, around half-way through the LSST Main Survey, we have a unique opportunity to explore synergistic science from two landmark programs. LSST can survey the entire footprint of the WFIRST microlensing survey in a single Deep Drilling Field. Here we explore the great scientific potential of this proposal and recommend the most effective observing strategies., Comment: White Paper submitted to the 2018 call for LSST cadence proposals. arXiv admin note: text overlap with arXiv:1812.03137

The overwhelming majority of objects visible to LSST lie within the Galactic Plane. Though many previous surveys have avoided this region for fear of stellar crowding, LSST's spatial resolution combined with its state-of-the-art Difference Image Analysis mean that it can conduct a high cadence survey of most of the Galaxy for the first time. Here we outline the many areas of science that would greatly benefit from an LSST survey that included the Galactic Plane, Magellanic Clouds and Bulge at a cadence of 2-3 d. Particular highlights include measuring the mass spectrum of black holes, and mapping the population of exoplanets in the Galaxy in relation to variations in star forming environments. But the same survey data will provide a goldmine for a wide range of science, and we explore possible survey strategies which maximize the scientific return for a number of fields including young stellar objects, cataclysmic variables and Neptune Trojans., Comment: White Paper submitted in response to the 2018 call for LSST cadence proposals

Gravitational lensing provides a means to measure mass that does not rely on detecting and analysing light from the lens itself. Compact objects are ideal gravitational lenses, because they have relatively large masses and are dim. In this paper we describe the prospects for predicting lensing events generated by the local population of compact objects, consisting of 250 neutron stars, 5 black holes, and approximately 35,000 white dwarfs. By focusing on a population of nearby compact objects with measured proper motions and known distances from us, we can measure their masses by studying the characteristics of any lensing event they generate. Here we concentrate on shifts in the position of a background source due to lensing by a foreground compact object. With HST, JWST, and Gaia, measurable centroid shifts caused by lensing are relatively frequent occurrences. We find that 30-50 detectable events per decade are expected for white dwarfs. Because relatively few neutron stars and black holes have measured distances and proper motions, it is more difficult to compute realistic rates for them. However, we show that at least one isolated neutron star has likely produced detectable events during the past several decades. This work is particularly relevant to the upcoming data releases by the Gaia mission and also to data that will be collected by JWST. Monitoring predicted microlensing events will not only help to determine the masses of compact objects, but will also potentially discover dim companions to these stellar remnants, including orbiting exoplanets., Comment: 16 pages, 10 figures, 8 tables

We report the test of many of the key elements of the laser-based calibration system for muon g - 2 experiment E989 at Fermilab. The test was performed at the Laboratori Nazionali di Frascati's Beam Test Facility using a 450 MeV electron beam impinging on a small subset of the final g - 2 lead-fluoride crystal calorimeter system. The calibration system was configured as planned for the E989 experiment and uses the same type of laser and most of the final optical elements. We show results regarding the calorimeter's response calibration, the maximum equivalent electron energy which can be provided by the laser and the stability of the calibration system components., Comment: Submitted to Nuclear Instrument and Methods A

Globular clusters are ancient stellar populations with no star formation or core-collapse supernovae. Several lines of evidence suggest that globular clusters are rich in planets. If so, and if advanced civilizations can develop there, then the distances between these civilizations and other stars would be far smaller than typical distances between stars in the Galactic disk. The relative proximity would facilitate interstellar communication and travel. However, the very proximity that promotes interstellar travel also brings danger, since stellar interactions can destroy planetary systems. However, by modeling globular clusters and their stellar populations, we find that large regions of many globular clusters can be thought of as "sweet spots" where habitable-zone planetary orbits can be stable for long times. We also compute the ambient densities and fluxes in the regions within which habitable-zone planets can survive. Globular clusters are among the best targets for searches for extraterrestrial intelligence (SETI). We use the Drake equation to compare globular clusters to the Galactic disk, in terms of the likelihood of housing advanced communicating civilizations. We also consider free-floating planets, since wide-orbit planets can be ejected and travel freely through the cluster. A civilization spawned in a globular cluster may have opportunities to establish self-sustaining outposts, thereby reducing the probability that a single catastrophic event will destroy the civilization or its descendants. Although individual civilizations within a cluster may follow different evolutionary paths, or even be destroyed, the cluster may always host some advanced civilization, once a small number of them have managed to jump across interstellar space., Comment: 28 pages, 3 figures 1 table

The new muon (g-2) experiment E989 at Fermilab will be equipped with a laser calibration system for all the 1296 channels of the calorimeters. An integrating sphere and an alternative system based on an engineered diffuser have been considered as possible light distributors for the experiment. We present here a detailed comparison of the two based on temporal response, spatial uniformity, transmittance and time stability., Comment: accepted to Nucl.Instrum.Meth. A

We report on the discovery of a new X-ray pulsator, Swift J201424.9+152930 (Sw J2014). Owing to its X-ray modulation at 491 s, it was discovered in a systematic search for coherent signals in the archival data of the Swift X-ray Telescope. To investigate the nature of Sw J2014, we performed multi-wavelength follow-up observations with space-borne (Swift and XMM-Newton) and ground-based (the 1.5-m Loiano Telescope and the 3.6-m Telescopio Nazionale Galileo) instruments. The X-ray spectrum of Sw J2014 can be described by a hard and highly absorbed power law. The optical observations made it possible to single out the optical counterpart to this source, which displays several variable emission lines and total eclipses lasting ~20 min. Total eclipses of similar length were observed also in X-rays. The study of the eclipses, allowed us to infer a second periodicity of 3.44 h, which we interpret as the orbital period of a close binary system. We also found that the period has not significantly changed over a ~7 yr timespan. Based on the timing signatures of Sw J2014, and its optical and X-ray spectral properties, we suggest that it is a close binary hosting an accreting magnetic white dwarf. The system is therefore a cataclysmic variable of the intermediate polar type and one of the very few showing deep eclipses., Comment: 12 pages, 8 figures, 5 tables; minor changes to match the final MNRAS version

The Muon (g-2) Experiment, E989 at Fermilab, will measure the muon anomalous magnetic moment a factor-of-four more precisely than was done in E821 at the Brookhaven National Laboratory AGS. The E821 result appears to be greater than the Standard-Model prediction by more than three standard deviations. When combined with expected improvement in the Standard-Model hadronic contributions, E989 should be able to determine definitively whether or not the E821 result is evidence for physics beyond the Standard Model. After a review of the physics motivation and the basic technique, which will use the muon storage ring built at BNL and now relocated to Fermilab, the design of the new experiment is presented. This document was created in partial fulfillment of the requirements necessary to obtain DOE CD-2/3 approval., Comment: 666 pages

The detections of both X-ray and radio emission from the cluster G1 in M31 have provided strong support for existing dynamical evidence for an intermediate mass black hole (IMBH) of mass 1.8 +/- 0.5 x 10^4 solar masses at the cluster center. However, given the relatively low significance and astrometric accuracy of the radio detection, and the non-simultaneity of the X-ray and radio measurements, this identification required further confirmation. Here we present deep, high angular resolution, strictly simultaneous X-ray and radio observations of G1. While the X-ray emission (L_X = 1.74^{+0.53}_{-0.44} x 10^{36} (d/750 kpc)^2 erg/s in the 0.5-10 keV band) remained fully consistent with previous observations, we detected no radio emission from the cluster center down to a 3-sigma upper limit of 4.7 microJy/beam. Our favored explanation for the previous radio detection is flaring activity from a black hole low mass X-ray binary (LMXB). We performed a new regression of the Fundamental Plane of black hole activity, valid for determining black hole mass from radio and X-ray observations of sub-Eddington black holes, finding log M_{BH} = (1.638 +/- 0.070)log L_R - (1.136 +/- 0.077)log L_X - (6.863 +/- 0.790), with an empirically-determined uncertainty of 0.44 dex. This constrains the mass of the X-ray source in G1, if a black hole, to be <9700 solar masses at 95% confidence, suggesting that it is a persistent LMXB. This annuls what was previously the most convincing evidence from radiation for an IMBH in the Local Group, though the evidence for an IMBH in G1 from velocity dispersion measurements remains unaffected by these results., Comment: 6 pages, 3 figures. Accepted for publication in ApJ Letters

We report the discovery of an X-ray/UV stellar flare from the source LMC 335, captured by XMM-Newton in the field of the Large Magellanic Cloud. The flare event was recorded continuously in X-ray for its first 10 hours from the precursor to the late decay phases. The observed fluxes increased by more than two orders of magnitude at its peak in X-ray and at least one in the UV as compared to quiescence. The peak 0.1-7.0 keV X-ray flux is derived from the two-temperature APEC model to be ~(8.4 +/- 0.6) x 10^-12 erg cm-2 s-1. Combining astrometric information from multiple X-ray observations in the quiescent and flare states, we identify the NIR counterpart of LMC 335 as the 2MASS source J05414534-6921512. The NIR color relations and spectroscopic parallax characterize the source as a Galactic K7-M4 dwarf at a foreground distance of (100 - 264) pc, implying a total energy output of the entire event of ~(0.4 - 2.9) x 10^35 erg. This report comprises detailed analyses of this late-K / early-M dwarf flare event that has the longest time coverage yet reported in the literature. The flare decay can be modeled with two exponential components with timescales of ~28 min and ~4 hours, with a single component decay firmly ruled out. The X-ray spectra during flare can be described by two components, a dominant high temperature component of ~40-60MK and a low temperature component of ~10MK, with a flare loop length of about 1.1-1.3 stellar radius., Comment: 35 pages, 6 figures, 5 tables, accepted for publication in ApJ

Type Ia supernovae (SNe Ia) play important roles in our study of the expansion and acceleration of the Universe, but because we do not know the exact nature or natures of the progenitors, there is a systematic uncertainty that must be resolved if SNe Ia are to become more precise cosmic probes. No progenitor system has ever been identified either in the pre- or post-explosion images of a Ia event. There have been recent claims for and against the detection of ex-companion stars in several SNe Ia remnants. These studies, however, usually ignore the angular momentum gain of the progenitor WD, which leads to a spin-up phase and a subsequent spin-down phase before explosion. For spin-down timescales greater than 10^5 years, the donor star could be too dim to detect by the time of explosion. Here we revisit the current limits on ex-companion stars to SNR 0509-67.5, a 400 year old remnant in the Large Magellanic Cloud (LMC). If the effects of possible angular momentum gain on the WD are included, a wide range of single-degenerate progenitor models are allowed for this remnant. We demonstrate that the current absence of evidence for ex-companion stars in this remnant, as well as other SNe Ia remnants, does not necessarily provide the evidence of absence for ex-companions. We discuss potential ways to identify such ex-companion stars through deep imaging observations., Comment: ApJ, accepted version

Recent observations indicate that many if not all galaxies host massive central black holes (BHs). In this paper we explore the influence of supermassive binary black holes (SMBBHs) on their actions as gravitational lenses. When lenses are modelled as singular isothermal ellipsoids, binary black holes change the critical curves and caustics differently as a function of distance. Each black hole can in principle create at least one additional image, which, if observed, provides evidence of black holes. By studying how SMBBHs affect the cumulative distribution of magnification for images created by black holes, we find that the cross section for at least one such additional image to have a magnification larger than $10^{-5}$ is comparable to the cross section for producing multiple-images in singular isothermal lenses. Such additional images may be detectable with high-resolution and large dynamic range maps of multiply-imaged systems from future facilities, such as the SKA. The probability of detecting at least one image (two images) with magnification above $10^{-3}$ is $\sim 0.2 \fBH$ ($\sim 0.05 \fBH$) in a multiply-imaged lens system, where $\fBH$ is the fraction of galaxies housing binary black holes. We also study the effects of SMBBHs on the core images when galaxies have shallower central density profiles (modelled as non-singular isothermal ellipsoids). We find that the cross section of the usually faint core images is further suppressed by SMBBHs. Thus their presence should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses., Comment: 9 pages, 5 figures, accepted for publication in MNRAS

In the single degenerate scenario for Type Ia supernova (SNeIa), a white dwarf (WD) must gain a significant amount of matter from a companion star. Because the accreted mass carries angular momentum, the WD is likely to achieve fast spin periods, which can increase the critical mass, $M_{crit}$, needed for explosion. When $M_{crit}$ is higher than the maximum mass achieved by the WD, the WD must spin down before it can explode. This introduces a delay between the time at which the WD has completed its epoch of mass gain and the time of the explosion. Matter ejected from the binary during mass transfer therefore has a chance to become diffuse, and the explosion occurs in a medium with a density similar to that of typical regions of the interstellar medium. Also, either by the end of the WD's mass increase or else by the time of explosion, the donor may exhaust its stellar envelope and become a WD. This alters, generally diminishing, explosion signatures related to the donor star. Nevertheless, the spin-up/spin-down model is highly predictive. Prior to explosion, progenitors can be super-$M_{Ch}$ WDs in either wide binaries with WD companions, or else in cataclysmic variables. These systems can be discovered and studied through wide-field surveys. Post explosion, the spin-up/spin-down model predicts a population of fast-moving WDs, low-mass stars, and even brown dwarfs. In addition, the spin-up/spin-down model provides a paradigm which may be able to explain both the similarities and the diversity observed among SNeIa., Comment: Submitted to ApJ Letters

We report the most accurate X-ray position of the X-ray source in the giant globular cluster G1 in M31 by using the Chandra X-ray Observatory, Hubble Space Telescope (HST), and Canada-France-Hawaii Telescope (CFHT). G1 is clearly detected with Chandra and by cross-registering with HST and CFHT images, we derive a 1-sigma error radius of 0.15 arcsec, significantly smaller than the previous measurement by XMM-Newton. We conclude that the X-ray emission of G1 is likely to come from within the core radius of the cluster. We have considered a number of possibilities for the origin of the X-ray emission but can rule all but two scenarios out: it could be due to either accretion onto a central intermediate-mass black hole (IMBH), or an ordinary low-mass X-ray binary (LMXB). Based on the X-ray luminosity and the Bondi accretion rate, an IMBH accreting from the cluster gas seems unlikely and we suggest that the X-rays are due to accretion from a companion. Alternatively, the probability that a 1.5 solar masses cluster LMXB lies within the 95 per cent X-ray error circle is about 0.7. Therefore we cannot rule out a single LMXB as the origin of the X-ray emission. While we cannot distinguish between different models with current observations, future high-resolution and high-sensitivity radio imaging observations will reveal whether there is an IMBH at the centre of G1., Comment: 5 pages, 1 figure, accepted for publication in MNRAS

This paper is devoted to exploring how we can discover and study nearby (< 1-2 kpc) planetary and binary systems by observing their action as gravitational lenses. Lensing can extend the realm of nearby binaries and planets that can be systematically studied to include dark and dim binaries, and face-on systems. As with more traditional studies, which use light from the system, orbital parameters (including the total mass, mass ratio, and orbital separation) can be extracted from lensing data, Also in common with these traditional methods, individual systems can be targeted for study. We discuss the specific observing strategies needed in order to optimize the discovery and study of nearby planetary and binary systems by observing their actions as lenses., Comment: 12 pages, 4 figures, submitted to ApJL 14 September 2007

Nearby masses can have a high probability of lensing stars in a distant background field. High-probability lensing, or mesolensing, can therefore be used to dramatically increase our knowledge of dark and dim objects in the solar neighborhood, where it can discover and study members of the local dark population (free-floating planets, low-mass dwarfs, white dwarfs, neutron stars, and stellar mass black holes). We can measure the mass and transverse velocity of those objects discovered (or already known), and determine whether or not they are in binaries with dim companions. We explore these and other applications of mesolensing, including the study of forms of matter that have been hypothesized but not discovered, such as intermediate-mass black holes, dark matter objects free-streaming through the Galactic disk, and planets in the outermost regions of the solar system. In each case we discuss the feasibility of deriving results based on present-day monitoring systems, and also consider the vistas that will open with the advent of all-sky monitoring in the era of the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), and the Large Synoptic Survey Telescope (LSST)., Comment: To appear in the Astrophysical Journal; 10 pages, no figures

We report the results of Chandra and XMM-Newton observations of a new outburst of an ultra-luminous X-ray source in M101. During a Chandra monitoring observation of M101, M101 ULX-1 was found to be in outburst in 2004 December, the second outburst in 2004. The peak bolometric luminosity is about 3e40 erg/s (7e39 erg/s in 0.3-7 keV). The outburst spectra are very soft and can generally be fitted with a blackbody model with temperatures of 40-80 eV, similar to supersoft X-ray sources in the Milky Way and in the Magellanic Clouds. In one Chandra observation, the source spectrum appears to be harder with a temperature of 150 eV. Such a spectral state is rarely seen in M101 ULX-1 and no X-ray source in the Milky Way shows this kind of spectrum. However, such an unusual spectral state very likely belongs to a new class of X-ray sources, quasisoft X-ray sources, recently discovered in nearby galaxies. M101 ULX-1 returned to supersoft state in a subsequent XMM-Newton observation. Based on the two outbursts in 2004, the extremely high luminosity (L_bol=1e40-1e41 erg/s), very soft X-ray spectra (kT=40-150 eV), transient behavior, and state transition provide strong evidence that M101 ULX-1 harbors an intermediate-mass black hole., Comment: 4 pages, 2 figures, accepted for publication in ApJL

We report the results of Chandra and XMM-Newton observations of an ultra-luminous supersoft X-ray source in M101. M101 ULX-1 underwent 2 outbursts in 2004 during which the peak bolometric luminosities reached 1e41 erg/s. The outburst spectra were very soft and can generally be fitted with a blackbody model with temperatures of 50-160 eV. In two of the observations, absorption edges at 0.33 keV, 0.56 keV, 0.66 keV, and 0.88 keV were found. A cool accretion disk was also found in the 2004 December outburst. During the low luminosity state, a power-law tail was seen up to 7 keV. It is clear the source changed from a low/hard state to a high/soft state. In addition, it showed at least 5 outbursts between 1996 and 2004. This is the first ultra-luminous X-ray source for which recurrent outbursts with state transitions similar to Galactic X-ray binaries have been observed. From the Hubble Space Telescope data, we found an optical counterpart to the source. During the 2004 outbursts, we also performed radio and ground-based optical observations. All the results strongly suggest that the accreting object is a > 2800 solar mass black hole., Comment: 6 pages, 4 figures; paper presented at the XXII Texas Symposium on Relativistic Astrophysics, Stanford University, Dec. 13-17, 2004

We report the results of Chandra and XMM-Newton observations of a new outburst of an ultraluminous supersoft X-ray source in M101. CXOU J140332.3+542103 was observed in a low luminosity state (L_x~1e37 erg/s) between 2004 January and May. The low-state X-ray spectra were relatively hard; the combined low state spectrum can be fitted with a combination of a power-law with photon index of 1.4 and a blackbody of 63 eV. During 2004 July, the source underwent a strong outburst and the peak 0.3-7 keV luminosity reached 3e40 erg/s, with a bolometric luminosity of about 1e41 erg/s. The outburst spectra were very soft and can generally be fitted with a blackbody model with temperatures of 50-100 eV. In two of the observations, absorption edges at 0.33 keV, 0.56 keV, 0.66 keV, and 0.88 keV were found. An XMM-Newton observation was also performed during the decay of the outburst, and a power-law tail was seen in addition to the supersoft spectrum. We consider different accretion models; one involving an intermediate-mass black hole can explain the observations., Comment: 4 pages, 1 figure; accepted for publication in ApJL

Chandra observations of nearby galaxies have revealed a number of X-ray sources characterized by high luminosities ($L_X > \sim10^{36}$ erg s$^{-1}$) and $k T$ in the range $100 - 350$ eV. These ``quasisoft X-ray sources'' (QSSs) are harder than luminous supersoft X-ray sources (SSSs), whose characteristic temperatures are tens of eV, but, with little or no emission above 2 keV, they are significantly softer than most canonical X-ray sources. They are likely to include a range of physical systems; some may be common systems in unusual states: neutron stars or stellar-mass black holes (BHs), or SNRs. Others may be accreting BHs of intermediate mass. We have analyzed {\it Chandra} data from $19 $ near-by galaxies to identify QSSs with sufficiently high count rates to allow spectral fits. Six of these galaxies have been studied in great detail; in these we find 89 SSSs and 122 QSSs. In this paper we present spectra for QSSs those with more than 50 counts. We also use data from six of these galaxies to study the broadband spectral distributions of QSSs and to compare them with those of other X-ray sources. Since QSSs in the Milky Way are likely to be missed by our selection procedures, we discuss other signatures by which they can be found. If some of the QSSs for which we have spectra are accreting BHs, the lower bounds for the BH mass range from roughly $30 M_\odot$ to $1000 M_\odot$., Comment: 9 pages, 3 figures, ApJ, submitted

We report on possible associations between 6 globular clusters in the Virgo elliptical galaxy NGC 4472 (M49) and bright ($L_x > 10^{38}$ erg s$^{-1}$) very soft X-ray sources (VSSs). Two of the VSSs have broad-band spectral properties consistent with those of luminous supersoft X-ray sources (SSSs). The other VSSs are somewhat harder harder, possibly with values of kT between roughly 150 eV and 250 eV. These sources may be too hot to be explained by the white dwarf models so promising for SSSs; they are members of the newly-established class of quasisoft sources (QSSs). We examine white dwarf, neutron star, and black hole models for the VSSs. One of the SSSs is hot and bright enough to be a possible progenitor of a Type Ia supernova, while the most natural model for the other VSSs is one in which the accretors are intermediate-mass black holes. Whatever their physical natures, these sources are unlike any X-ray sources in Galactic globular clusters. No Galactic globular cluster houses a quasisoft source and only one, M3, contains a dim ($L_x \sim 10^{36}$ erg s$^{-1}$) transient SSS., Comment: 5 pages, 1 figure, submitted to ApJL

We introduce a procedure to identify very soft X-ray sources (VSSs) in external galaxies. Our immediate goal was to formulate a systematic procedure to identify luminous supersoft X-ray sources (SSSs), so as to allow comparisons among galaxies and to study environmental effects. The focus of this paper is on the design of the selection algorithm and on its application to simulated data. In the companion paper we test it by applying it to sources discovered through Chandra observations of 4 galaxies. We find that, in its application to both simulated and real data, our procedure also selects somewhat harder sources, which we call quasisoft. Whereas values of kT for SSSs are typically tens of eV, some quasisoft sources (QSSs) may have kT as high as ~250-300 eV. The dominant spectral component of other QSSs may be as soft as SSS spectra, but the spectra may also include a low-luminosity harder component. We sketch physical models for both supersoft and quasisoft sources. Some SSSs are likely to be accreting white dwarfs; some of these may be progenitors of Type Ia supernovae. Most QSSs may be too hot to be white dwarfs. They, together with a subset of SSSs, may be neutron stars or, perhaps most likely, accreting intermediate-mass black holes., Comment: 15 pages, 4 figures, ApJ, submitted

We apply a uniform procedure to select very soft sources from point sources observed by Chandra in 4 galaxies. This sample includes one elliptical galaxy (NGC 4967), 2 face-on spirals (M101 and M83), and an interacting galaxy (M51). We have found very soft X-ray sources (VSSs) in every galaxy. Some of these fit the criteria for canonical supersoft sources (SSSs), while others are somewhat harder. These latter have characteristic values of kT < 300 eV; we refer to them as quasisoft sources (QSSs). We found a combined total of 149 VSSs in the 4 galaxies we considered; 77 were SSSs and 72 were QSSs. (See the paper for the original long abstract), Comment: 20 pages, 6 figures. Accepted for publication in ApJ

We present new optical and {\it Chandra} observations of the field containing the ultraluminous X-ray source NGC 1313 X-2. On an ESO 3.6 m image, the {\it Chandra} error box embraces a $R=21.6$ point-like object and excludes a previously proposed optical counterpart. The resulting X-ray/optical flux ratio of NGC 1313 X-2 is $\sim 500$. The value of $f_X/f_{opt}$, the X-ray variability history and the spectral distribution derived from a re-analysis of the {\it ROSAT}, {\it ASCA} and {\it XMM} data indicate a luminous X-ray binary in NGC 1313 as a likely explanation for NGC 1313 X-2. If the X-ray soft component observed in the {\it XMM} EPIC spectrum originates from an accretion disk, the inferred mass of the compact remnant is $\approx 100 M_\odot$, making it an intermediate mass black hole. The derived optical luminosity ($L\approx 10^5 L_\odot$) is consistent with that of a $\approx 15-20 M_\odot$ companion. The properties of the environment of NGC 1313 X-2 are briefly discussed., Comment: 12 pages, 6 figures; submitted to The Astrophysical Journal