Showing 1–35 of 35 results for author: Armstrong, E

Constraining electron number density in the Sun via Earth-based neutrino flavor data

Authors: Caroline Laber-Smith, Eve Armstrong, A. Baha Balantekin, Elizabeth K. Jones, Lily Newkirk, Amol V. Patwardhan, Sarah Ranginwala, M. Margarette Sanchez, Hansen Torres

Abstract: Neutrino flavor transformation offers a window into the physics of various astrophysical environments, including our Sun and the more exotic environs of core-collapse supernovae and binary neutron-star mergers. Here, we apply an inference framework - specifically: statistical data assimilation (SDA) - to neutrino flavor evolution in the Sun. We take a model for solar neutrino flavor evolution, tog… ▽ More Neutrino flavor transformation offers a window into the physics of various astrophysical environments, including our Sun and the more exotic environs of core-collapse supernovae and binary neutron-star mergers. Here, we apply an inference framework - specifically: statistical data assimilation (SDA) - to neutrino flavor evolution in the Sun. We take a model for solar neutrino flavor evolution, together with Earth-based neutrino measurements, to infer solar properties. Specifically, we ask what signature of the radially-varying solar electron number density $n_e(r)$ is contained within these Earth-based measurements. Currently, the best estimates of $n_e(r)$ come from the standard solar model. We seek to ascertain, through novel application of the SDA method, whether estimates of the same from neutrino data can serve as independent constraints. △ Less

Submitted 9 April, 2024; originally announced April 2024.

New constraints on ultraheavy dark matter from the LZ experiment

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, J. W. Bargemann, A. Baxter, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger, B. Boxer, C. A. J. Brew , et al. (174 additional authors not shown)

Abstract: Searches for dark matter with liquid xenon time projection chamber experiments have traditionally focused on the region of the parameter space that is characteristic of weakly interacting massive particles, ranging from a few GeV/$c^2$ to a few TeV/$c^2$. Models of dark matter with a mass much heavier than this are well motivated by early production mechanisms different from the standard thermal f… ▽ More Searches for dark matter with liquid xenon time projection chamber experiments have traditionally focused on the region of the parameter space that is characteristic of weakly interacting massive particles, ranging from a few GeV/$c^2$ to a few TeV/$c^2$. Models of dark matter with a mass much heavier than this are well motivated by early production mechanisms different from the standard thermal freeze-out, but they have generally been less explored experimentally. In this work, we present a re-analysis of the first science run (SR1) of the LZ experiment, with an exposure of $0.9$ tonne$\times$year, to search for ultraheavy particle dark matter. The signal topology consists of multiple energy deposits in the active region of the detector forming a straight line, from which the velocity of the incoming particle can be reconstructed on an event-by-event basis. Zero events with this topology were observed after applying the data selection calibrated on a simulated sample of signal-like events. New experimental constraints are derived, which rule out previously unexplored regions of the dark matter parameter space of spin-independent interactions beyond a mass of 10$^{17}$ GeV/$c^2$. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 9 pages, 7 figures

Journal ref: Phys. Rev. D 109, 112010 (2024)

First Constraints on WIMP-Nucleon Effective Field Theory Couplings in an Extended Energy Region From LUX-ZEPLIN

Authors: LZ Collaboration, J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, J. W. Bargemann, A. Baxter, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger , et al. (175 additional authors not shown)

Abstract: Following the first science results of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time projection chamber operating from the Sanford Underground Research Facility in Lead, South Dakota, USA, we report the initial limits on a model-independent non-relativistic effective field theory describing the complete set of possible interactions of a weakly interacting massive particle (WIMP) with a n… ▽ More Following the first science results of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time projection chamber operating from the Sanford Underground Research Facility in Lead, South Dakota, USA, we report the initial limits on a model-independent non-relativistic effective field theory describing the complete set of possible interactions of a weakly interacting massive particle (WIMP) with a nucleon. These results utilize the same 5.5 t fiducial mass and 60 live days of exposure collected for the LZ spin-independent and spin-dependent analyses while extending the upper limit of the energy region of interest by a factor of 7.5 to 270 keVnr. No significant excess in this high energy region is observed. Using a profile-likelihood ratio analysis, we report 90% confidence level exclusion limits on the coupling of each individual non-relativistic WIMP-nucleon operator for both elastic and inelastic interactions in the isoscalar and isovector bases. △ Less

Submitted 26 February, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 17 pages 11 figures

Journal ref: Phys. Rev. D 109, 092003 (2024)

Markov Chain Monte Carlo for Bayesian Parametric Galaxy Modeling in LSST

Authors: James J. Buchanan, Michael D. Schneider, Kerianne Pruett, Robert E. Armstrong

Abstract: We apply Markov Chain Monte Carlo (MCMC) to the problem of parametric galaxy modeling, estimating posterior distributions of galaxy properties such as ellipticity and brightness for more than 100,000 images of galaxies taken from DC2, a simulated telescope survey resembling the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We use a physically informed prior and apply selection… ▽ More We apply Markov Chain Monte Carlo (MCMC) to the problem of parametric galaxy modeling, estimating posterior distributions of galaxy properties such as ellipticity and brightness for more than 100,000 images of galaxies taken from DC2, a simulated telescope survey resembling the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We use a physically informed prior and apply selection corrections to the likelihood. The resulting posterior samples enable rigorous probabilistic inference of galaxy model parameters and their uncertainties. These posteriors are one key ingredient in a fully probabilistic description of galaxy catalogs, which can ultimately enable a refined Bayesian estimate of cosmological parameters. We systematically examine the reliability of the posterior mean as a point estimator of galaxy parameters, and of the posterior width as a measure of uncertainty, under some common modeling approximations. We implement the probabilistic modeling and MCMC inference using the JIF (Joint Image Framework) tool, which we make freely available online. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 34 pages, 16 figures

A search for new physics in low-energy electron recoils from the first LZ exposure

Authors: The LZ Collaboration, J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, J. W. Bargemann, A. Baxter, K. Beattie, P. Beltrame, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, G. M. Blockinger , et al. (178 additional authors not shown)

Abstract: The LUX-ZEPLIN (LZ) experiment is a dark matter detector centered on a dual-phase xenon time projection chamber. We report searches for new physics appearing through few-keV-scale electron recoils, using the experiment's first exposure of 60 live days and a fiducial mass of 5.5t. The data are found to be consistent with a background-only hypothesis, and limits are set on models for new physics inc… ▽ More The LUX-ZEPLIN (LZ) experiment is a dark matter detector centered on a dual-phase xenon time projection chamber. We report searches for new physics appearing through few-keV-scale electron recoils, using the experiment's first exposure of 60 live days and a fiducial mass of 5.5t. The data are found to be consistent with a background-only hypothesis, and limits are set on models for new physics including solar axion electron coupling, solar neutrino magnetic moment and millicharge, and electron couplings to galactic axion-like particles and hidden photons. Similar limits are set on weakly interacting massive particle (WIMP) dark matter producing signals through ionized atomic states from the Migdal effect. △ Less

Submitted 9 September, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 13 pages, 10 figures. See https://tinyurl.com/LZDataReleaseRun1ER for a data release related to this paper

Journal ref: Phys. Rev. D 108, 072006 (2023)

I Murdered Conan O'Brien and Nobody Will Ever Know -- an exercise in inference sabotage

Authors: Eve Armstrong

Abstract: I employ an optimization-based inference methodology together with an Ising model, in an intentionally ineffectual manner, to get away with murdering an obstreperous scientific collaborator. The antics of this collaborator, hereafter "Conan O'Brien," were impeding the publication of an important manuscript. With my tenure date looming, I found myself desperate. Luckily, I study inference, a comput… ▽ More I employ an optimization-based inference methodology together with an Ising model, in an intentionally ineffectual manner, to get away with murdering an obstreperous scientific collaborator. The antics of this collaborator, hereafter "Conan O'Brien," were impeding the publication of an important manuscript. With my tenure date looming, I found myself desperate. Luckily, I study inference, a computational means to find a solution to a physical problem, based on available measurements (say, a dead body) and a dynamical model assumed to give rise to those measurements (a murderer). If the measurements are insufficient and/or the model is incomplete, one obtains multiple "degenerate" solutions to the problem. Degenerate solutions are all equally valid given the information available, and thus render meaningless the notion of one "correct" solution. Typically in scientific research, degeneracy is undesirable. Here I describe the opposite situation: a quest to create degenerate solutions in which to cloak myself. Or even better: to render measurements incompatible with a solution in which I am the murderer. Moreover, I show how one may sabotage an inference procedure to commit an untraceable crime. I sit here now, typing victoriously, a free woman. Because you won't believe me anyway. And even if you do, you'll never prove a thing. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: 10 pages, 3 figures, 2 tables

Inference finds consistency between a neutrino flavor evolution model and Earth-based solar neutrino measurements

Authors: Caroline Laber-Smith, A. A. Ahmetaj, Eve Armstrong, A. Baha Balantekin, Amol V. Patwardhan, M. Margarette Sanchez, Sherry Wong

Abstract: We continue examining statistical data assimilation (SDA), an inference methodology, to infer solutions to neutrino flavor evolution, for the first time using real - rather than simulated - data. The model represents neutrinos streaming from the Sun's center and undergoing a Mikheyev-Smirnov-Wolfenstein (MSW) resonance in flavor space, due to the radially-varying electron number density. The model… ▽ More We continue examining statistical data assimilation (SDA), an inference methodology, to infer solutions to neutrino flavor evolution, for the first time using real - rather than simulated - data. The model represents neutrinos streaming from the Sun's center and undergoing a Mikheyev-Smirnov-Wolfenstein (MSW) resonance in flavor space, due to the radially-varying electron number density. The model neutrino energies are chosen to correspond to experimental bins in the Sudbury Neutrino Observatory (SNO) and Borexino experiments, which measure electron-flavor survival probability at Earth. The procedure successfully finds consistency between the observed fluxes and the model, if the MSW resonance - that is, flavor evolution due to solar electrons - is included in the dynamical equations representing the model. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Light curve completion and forecasting using fast and scalable Gaussian processes (MuyGPs)

Authors: Imène R. Goumiri, Alec M. Dunton, Amanda L. Muyskens, Benjamin W. Priest, Robert E. Armstrong

Abstract: Temporal variations of apparent magnitude, called light curves, are observational statistics of interest captured by telescopes over long periods of time. Light curves afford the exploration of Space Domain Awareness (SDA) objectives such as object identification or pose estimation as latent variable inference problems. Ground-based observations from commercial off the shelf (COTS) cameras remain… ▽ More Temporal variations of apparent magnitude, called light curves, are observational statistics of interest captured by telescopes over long periods of time. Light curves afford the exploration of Space Domain Awareness (SDA) objectives such as object identification or pose estimation as latent variable inference problems. Ground-based observations from commercial off the shelf (COTS) cameras remain inexpensive compared to higher precision instruments, however, limited sensor availability combined with noisier observations can produce gappy time-series data that can be difficult to model. These external factors confound the automated exploitation of light curves, which makes light curve prediction and extrapolation a crucial problem for applications. Traditionally, image or time-series completion problems have been approached with diffusion-based or exemplar-based methods. More recently, Deep Neural Networks (DNNs) have become the tool of choice due to their empirical success at learning complex nonlinear embeddings. However, DNNs often require large training data that are not necessarily available when looking at unique features of a light curve of a single satellite. In this paper, we present a novel approach to predicting missing and future data points of light curves using Gaussian Processes (GPs). GPs are non-linear probabilistic models that infer posterior distributions over functions and naturally quantify uncertainty. However, the cubic scaling of GP inference and training is a major barrier to their adoption in applications. In particular, a single light curve can feature hundreds of thousands of observations, which is well beyond the practical realization limits of a conventional GP on a single machine. Consequently, we employ MuyGPs, a scalable framework for hyperparameter estimation of GP models that uses nearest neighbors sparsification and local cross-validation. MuyGPs... △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: 14 pages, 7 figures, accepted to AMOS 2022 conference

First Dark Matter Search Results from the LUX-ZEPLIN (LZ) Experiment

Authors: J. Aalbers, D. S. Akerib, C. W. Akerlof, A. K. Al Musalhi, F. Alder, A. Alqahtani, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, S. Azadi, A. J. Bailey, A. Baker, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, M. J. Barry, J. Barthel, D. Bauer, A. Baxter , et al. (322 additional authors not shown)

Abstract: The LUX-ZEPLIN experiment is a dark matter detector centered on a dual-phase xenon time projection chamber operating at the Sanford Underground Research Facility in Lead, South Dakota, USA. This Letter reports results from LUX-ZEPLIN's first search for weakly interacting massive particles (WIMPs) with an exposure of 60~live days using a fiducial mass of 5.5 t. A profile-likelihood ratio analysis s… ▽ More The LUX-ZEPLIN experiment is a dark matter detector centered on a dual-phase xenon time projection chamber operating at the Sanford Underground Research Facility in Lead, South Dakota, USA. This Letter reports results from LUX-ZEPLIN's first search for weakly interacting massive particles (WIMPs) with an exposure of 60~live days using a fiducial mass of 5.5 t. A profile-likelihood ratio analysis shows the data to be consistent with a background-only hypothesis, setting new limits on spin-independent WIMP-nucleon, spin-dependent WIMP-neutron, and spin-dependent WIMP-proton cross sections for WIMP masses above 9 GeV/c$^2$. The most stringent limit is set for spin-independent scattering at 36 GeV/c$^2$, rejecting cross sections above 9.2$\times 10^{-48}$ cm$^2$ at the 90% confidence level. △ Less

Submitted 2 August, 2023; v1 submitted 8 July, 2022; originally announced July 2022.

Comments: 9 pages, 8 figures. See https://doi.org/10.1103/PhysRevLett.131.041002 for a data release related to this paper

Journal ref: Phys. Rev. Lett. 131, 041002 (2023)

The Three Little Pigs and the Big Bad Wolf: Case Studies of Peer Review

Authors: Eve Armstrong

Abstract: I present for your appraisal three independent cases of the manuscript referee process conducted by a venerable peer-reviewed scientific journal. Each case involves a little pig, who submitted for consideration a theoretical plan for a house to be constructed presently, in a faraway land. An anonymous big bad wolf was assigned by the journal to assess the merit of these manuscripts. The pigs propo… ▽ More I present for your appraisal three independent cases of the manuscript referee process conducted by a venerable peer-reviewed scientific journal. Each case involves a little pig, who submitted for consideration a theoretical plan for a house to be constructed presently, in a faraway land. An anonymous big bad wolf was assigned by the journal to assess the merit of these manuscripts. The pigs proposed three distinct construction frameworks, which varied in physical and mathematical sophistication. The first little pig submitted a model of straw, based on the numerical method of toe-counting. His design included odd features, such as spilled millet and cloven-hoofprints on the window sill -- possibly a ploy to distract the wolf from the manuscript's facile mathematical foundation. The second little pig used a more advanced approach, employing Newton's classical laws of motion, to propose a house of sticks. This pig included in her manuscript copious citations by a specific wolf, possibly aiming to ensure acceptance by flattering the wolf whom she anticipated would be the referee. The third little pig described an ostentatious house of bricks based on an elaborate dynamical systems and stability analysis, possibly scheming to dazzle and impress. The big bad wolf did not appear moved by any of the pigs' tactics. His recommendations were, for straw: the minor revision of water-proofing; for sticks: the major revision of fire-proofing, given concerns surrounding climate change; for bricks: unequivocal rejection, accompanied by multiple derogatory comments regarding "high-and-mighty theorists." I describe each case in detail, and suggest that the wolf's reports might be driven as much by self interest as the manuscripts themselves -- namely, that at the time the wolf wrote the reviews, he was rather hungry. Finally, I examine morals learned, if any. △ Less

Submitted 31 March, 2022; originally announced March 2022.

Comments: 7 pages, 6 figures. http://www.amazon.com/author/evearmstrong

A Next-Generation Liquid Xenon Observatory for Dark Matter and Neutrino Physics

Authors: J. Aalbers, K. Abe, V. Aerne, F. Agostini, S. Ahmed Maouloud, D. S. Akerib, D. Yu. Akimov, J. Akshat, A. K. Al Musalhi, F. Alder, S. K. Alsum, L. Althueser, C. S. Amarasinghe, F. D. Amaro, A. Ames, T. J. Anderson, B. Andrieu, N. Angelides, E. Angelino, J. Angevaare, V. C. Antochi, D. Antón Martin, B. Antunovic, E. Aprile, H. M. Araújo , et al. (572 additional authors not shown)

Abstract: The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neut… ▽ More The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neutrinos through neutrinoless double-beta decay and through a variety of astrophysical sources. A next-generation xenon-based detector will therefore be a true multi-purpose observatory to significantly advance particle physics, nuclear physics, astrophysics, solar physics, and cosmology. This review article presents the science cases for such a detector. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 77 pages, 40 figures, 1262 references

Report number: INT-PUB-22-003

Journal ref: J. Phys. G: Nucl. Part. Phys. 50 (2023) 013001

Inference of bipolar neutrino flavor oscillations near a core-collapse supernova, based on multiple measurements at Earth

Authors: Eve Armstrong, Amol V. Patwardhan, A. A. Ahmetaj, M. Margarette Sanchez, Sophia Miskiewicz, Marcus Ibrahim, Ishaan Singh

Abstract: Neutrinos in compact-object environments, such as core-collapse supernovae, can experience various kinds of collective effects in flavor space, engendered by neutrino-neutrino interactions. These include "bipolar" collective oscillations, which are exhibited by neutrino ensembles where different flavors dominate at different energies. Considering the importance of neutrinos in the dynamics and nuc… ▽ More Neutrinos in compact-object environments, such as core-collapse supernovae, can experience various kinds of collective effects in flavor space, engendered by neutrino-neutrino interactions. These include "bipolar" collective oscillations, which are exhibited by neutrino ensembles where different flavors dominate at different energies. Considering the importance of neutrinos in the dynamics and nucleosynthesis in these environments, it is desirable to ascertain whether an Earth-based detection could contain signatures of bipolar oscillations that occurred within a supernova envelope. To that end, we continue examining a cost-function formulation of statistical data assimilation (SDA) to infer solutions to a small-scale model of neutrino flavor transformation. SDA is an inference paradigm designed to optimize a model with sparse data. Our model consists of two mono-energetic neutrino beams with different energies emanating from a source and coherently interacting with each other and with a matter background, with time-varying interaction strengths. We attempt to infer flavor transformation histories of these beams using simulated measurements of the flavor content at locations in vacuum (that is, far from the source), which could in principle correspond to earth-based detectors. Within the scope of this small-scale model, we found that: (i) based on such measurements, the SDA procedure is able to infer \textit{whether} bipolar oscillations had occurred within the protoneutron star envelope, and (ii) if the measurements are able to sample the full amplitude of the neutrino oscillations in vacuum, then the amplitude of the prior bipolar oscillations is also well predicted. This result intimates that the inference paradigm can well complement numerical integration codes, via its ability to infer flavor evolution at physically inaccessible locations. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: 12 pages, 7 figures, 1 table

Cosmogenic production of $^{37}$Ar in the context of the LUX-ZEPLIN experiment

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, A. Baker, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, D. Bauer, A. Baxter, K. Beattie, E. P. Bernard, A. Bhatti, A. Biekert, T. P. Biesiadzinski , et al. (183 additional authors not shown)

Abstract: We estimate the amount of $^{37}$Ar produced in natural xenon via cosmic ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth's surface. We then calculate the resulting $^{37}$Ar concentration in a 10-tonne payload~(similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage and delivery to the… ▽ More We estimate the amount of $^{37}$Ar produced in natural xenon via cosmic ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth's surface. We then calculate the resulting $^{37}$Ar concentration in a 10-tonne payload~(similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage and delivery to the underground facility. Using the spallation model by Silberberg and Tsao, the sea level production rate of $^{37}$Ar in natural xenon is estimated to be 0.024~atoms/kg/day. Assuming the xenon is successively purified to remove radioactive contaminants in 1-tonne batches at a rate of 1~tonne/month, the average $^{37}$Ar activity after 10~tonnes are purified and transported underground is 0.058--0.090~$μ$Bq/kg, depending on the degree of argon removal during above-ground purification. Such cosmogenic $^{37}$Ar will appear as a noticeable background in the early science data, while decaying with a 35~day half-life. This newly-noticed production mechanism of $^{37}$Ar should be considered when planning for future liquid xenon-based experiments. △ Less

Submitted 22 March, 2022; v1 submitted 8 January, 2022; originally announced January 2022.

Inference solves a boundary-value collision problem, with relevance to neutrino flavor transformation

Authors: Eve Armstrong

Abstract: Understanding neutrino flavor transformation in dense environments such as core-collapse supernovae (CCSN) is critical for inferring nucleosynthesis and interpreting a detected neutrino signal. The role of direction-changing collisions in shaping the neutrino flavor field in these environments is important and poorly understood; it has not been treated self-consistently. There has been progress, v… ▽ More Understanding neutrino flavor transformation in dense environments such as core-collapse supernovae (CCSN) is critical for inferring nucleosynthesis and interpreting a detected neutrino signal. The role of direction-changing collisions in shaping the neutrino flavor field in these environments is important and poorly understood; it has not been treated self-consistently. There has been progress, via numerical integration, to include the effects of collisions in the dynamics of the neutrino flavor field. While this has led to important insights, integration is limited by its requirement that full initial conditions must be assumed known. On the contrary, feedback from collisions to the neutrino field is a boundary value problem. Numerical integration techniques are poorly equipped to handle that formulation. This paper demonstrates that an inference formulation of the problem can solve a simple collision-only model representing a CCSN core -- without full knowledge of initial conditions. Rather, the procedure solves a two-point boundary value problem with partial information at the bounds. The model is sufficiently simple that physical reasoning may be used as a confidence check on the inference-based solution, and the procedure recovers the expected model dynamics. This result demonstrates that inference can solve a problem that is artificially hidden from integration techniques -- a problem that is an important feature of flavor evolution in dense environments. Thus, it is worthwhile to explore means of augmenting the existing powerful integration tools with inference-based approaches. △ Less

Submitted 14 November, 2021; originally announced November 2021.

Comments: 9 pages, 3 figures

Gaussian Process Classification for Galaxy Blend Identification in LSST

Authors: James J. Buchanan, Michael D. Schneider, Robert E. Armstrong, Amanda L. Muyskens, Benjamin W. Priest, Ryan J. Dana

Abstract: A significant fraction of observed galaxies in the Rubin Observatory Legacy Survey of Space and Time (LSST) will overlap at least one other galaxy along the same line of sight, in a so-called "blend." The current standard method of assessing blend likelihood in LSST images relies on counting up the number of intensity peaks in the smoothed image of a blend candidate, but the reliability of this pr… ▽ More A significant fraction of observed galaxies in the Rubin Observatory Legacy Survey of Space and Time (LSST) will overlap at least one other galaxy along the same line of sight, in a so-called "blend." The current standard method of assessing blend likelihood in LSST images relies on counting up the number of intensity peaks in the smoothed image of a blend candidate, but the reliability of this procedure has not yet been comprehensively studied. Here we construct a realistic distribution of blended and unblended galaxies through high-fidelity simulations of LSST-like images, and from this we examine the blend classification accuracy of the standard peak-finding method. Furthermore, we develop a novel Gaussian process blend classifier model, and show that this classifier is competitive with both the peak-finding method as well as with a convolutional neural network model. Finally, whereas the peak-finding method does not naturally assign probabilities to its classification estimates, the Gaussian process model does, and we show that the Gaussian process classification probabilities are generally reliable. △ Less

Submitted 10 December, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 20 pages, 6 figures, version accepted by ApJ

Star-Galaxy Image Separation with Computationally Efficient Gaussian Process Classification

Authors: Amanda L. Muyskens, Imène R. Goumiri, Benjamin W. Priest, Michael D. Schneider, Robert E. Armstrong, Jason M. Bernstein, Ryan Dana

Abstract: We introduce a novel method for discerning optical telescope images of stars from those of galaxies using Gaussian processes (GPs). Although applications of GPs often struggle in high-dimensional data modalities such as optical image classification, we show that a low-dimensional embedding of images into a metric space defined by the principal components of the data suffices to produce high-qualit… ▽ More We introduce a novel method for discerning optical telescope images of stars from those of galaxies using Gaussian processes (GPs). Although applications of GPs often struggle in high-dimensional data modalities such as optical image classification, we show that a low-dimensional embedding of images into a metric space defined by the principal components of the data suffices to produce high-quality predictions from real large-scale survey data. We develop a novel method of GP classification hyperparameter training that scales approximately linearly in the number of image observations, which allows for application of GP models to large-size Hyper Suprime-Cam (HSC) Subaru Strategic Program data. In our experiments we evaluate the performance of a principal component analysis (PCA) embedded GP predictive model against other machine learning algorithms including a convolutional neural network and an image photometric morphology discriminator. Our analysis shows that our methods compare favorably with current methods in optical image classification while producing posterior distributions from the GP regression that can be used to quantify object classification uncertainty. We further describe how classification uncertainty can be used to efficiently parse large-scale survey imaging data to produce high-confidence object catalogs. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: 18 pages, 8 figures

Projected sensitivities of the LUX-ZEPLIN (LZ) experiment to new physics via low-energy electron recoils

Authors: The LZ Collaboration, D. S. Akerib, A. K. Al Musalhi, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, D. Bauer, A. Baxter, P. Beltrame, E. P. Bernard, A. Bernstein, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch , et al. (172 additional authors not shown)

Abstract: LUX-ZEPLIN (LZ) is a dark matter detector expected to obtain world-leading sensitivity to weakly interacting massive particles (WIMPs) interacting via nuclear recoils with a ~7-tonne xenon target mass. This manuscript presents sensitivity projections to several low-energy signals of the complementary electron recoil signal type: 1) an effective neutrino magnetic moment and 2) an effective neutrino… ▽ More LUX-ZEPLIN (LZ) is a dark matter detector expected to obtain world-leading sensitivity to weakly interacting massive particles (WIMPs) interacting via nuclear recoils with a ~7-tonne xenon target mass. This manuscript presents sensitivity projections to several low-energy signals of the complementary electron recoil signal type: 1) an effective neutrino magnetic moment and 2) an effective neutrino millicharge, both for pp-chain solar neutrinos, 3) an axion flux generated by the Sun, 4) axion-like particles forming the galactic dark matter, 5) hidden photons, 6) mirror dark matter, and 7) leptophilic dark matter. World-leading sensitivities are expected in each case, a result of the large 5.6t 1000d exposure and low expected rate of electron recoil backgrounds in the $<$100keV energy regime. A consistent signal generation, background model and profile-likelihood analysis framework is used throughout. △ Less

Submitted 18 May, 2021; v1 submitted 23 February, 2021; originally announced February 2021.

Comments: v2 updates exclusion sensitivities from single-sided to two-sided

Enhancing the sensitivity of the LUX-ZEPLIN (LZ) dark matter experiment to low energy signals

Authors: D. S. Akerib, A. K. Al Musalhi, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, D. Bauer, A. Baxter, P. Beltrame, E. P. Bernard, A. Bernstein, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, G. M. Blockinger , et al. (162 additional authors not shown)

Abstract: Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matt… ▽ More Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matter and astrophysical neutrinos, which will be applicable to other liquid xenon detectors. The energy threshold is determined by the number of detected S1 photons; typically, these must be recorded in three or more photomultiplier channels to avoid dark count coincidences that mimic real signals. To lower this threshold: a) we take advantage of the double photoelectron emission effect, whereby a single vacuum ultraviolet photon has a $\sim20\%$ probability of ejecting two photoelectrons from a photomultiplier tube photocathode; and b) we drop the requirement of an S1 signal altogether, and use only the ionization signal, which can be detected more efficiently. For both techniques we develop signal and background models for the nominal exposure, and explore accompanying systematic effects, including the dependence on the free electron lifetime in the liquid xenon. When incorporating double photoelectron signals, we predict a factor of $\sim 4$ sensitivity improvement to the dark matter-nucleon scattering cross-section at $2.5$ GeV/c$^2$, and a factor of $\sim1.6$ increase in the solar $^8$B neutrino detection rate. Dropping the S1 requirement may allow sensitivity gains of two orders of magnitude in both cases. Finally, we apply these techniques to even lower masses by taking into account the atomic Migdal effect; this could lower the dark matter particle mass threshold to $80$ MeV/c$^2$. △ Less

Submitted 21 January, 2021; originally announced January 2021.

Comments: 14 pages, 6 figures

Star-Galaxy Separation via Gaussian Processes with Model Reduction

Authors: Imène R. Goumiri, Amanda L. Muyskens, Michael D. Schneider, Benjamin W. Priest, Robert E. Armstrong

Abstract: Modern cosmological surveys such as the Hyper Suprime-Cam (HSC) survey produce a huge volume of low-resolution images of both distant galaxies and dim stars in our own galaxy. Being able to automatically classify these images is a long-standing problem in astronomy and critical to a number of different scientific analyses. Recently, the challenge of "star-galaxy" classification has been approached… ▽ More Modern cosmological surveys such as the Hyper Suprime-Cam (HSC) survey produce a huge volume of low-resolution images of both distant galaxies and dim stars in our own galaxy. Being able to automatically classify these images is a long-standing problem in astronomy and critical to a number of different scientific analyses. Recently, the challenge of "star-galaxy" classification has been approached with Deep Neural Networks (DNNs), which are good at learning complex nonlinear embeddings. However, DNNs are known to overconfidently extrapolate on unseen data and require a large volume of training images that accurately capture the data distribution to be considered reliable. Gaussian Processes (GPs), which infer posterior distributions over functions and naturally quantify uncertainty, haven't been a tool of choice for this task mainly because popular kernels exhibit limited expressivity on complex and high-dimensional data. In this paper, we present a novel approach to the star-galaxy separation problem that uses GPs and reap their benefits while solving many of the issues traditionally affecting them for classification of high-dimensional celestial image data. After an initial filtering of the raw data of star and galaxy image cutouts, we first reduce the dimensionality of the input images by using a Principal Components Analysis (PCA) before applying GPs using a simple Radial Basis Function (RBF) kernel on the reduced data. Using this method, we greatly improve the accuracy of the classification over a basic application of GPs while improving the computational efficiency and scalability of the method. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Inference of neutrino flavor evolution through data assimilation and neural differential equations

Authors: Ermal Rrapaj, Amol V. Patwardhan, Eve Armstrong, George Fuller

Abstract: The evolution of neutrino flavor in dense environments such as core-collapse supernovae and binary compact object mergers constitutes an important and unsolved problem. Its solution has potential implications for the dynamics and heavy-element nucleosynthesis in these environments. In this paper, we build upon recent work to explore inference-based techniques for estimation of model parameters and… ▽ More The evolution of neutrino flavor in dense environments such as core-collapse supernovae and binary compact object mergers constitutes an important and unsolved problem. Its solution has potential implications for the dynamics and heavy-element nucleosynthesis in these environments. In this paper, we build upon recent work to explore inference-based techniques for estimation of model parameters and neutrino flavor evolution histories. We combine data assimilation, ordinary differential equation solvers, and neural networks to craft an inference approach tailored for non-linear dynamical systems. Using this architecture, and a simple two-neutrino, two-flavor model, we test various optimization algorithms with the help of four experimental setups. We find that employing this new architecture, together with evolutionary optimization algorithms, accurately captures flavor histories in the four experiments. This work provides more options for extending inference techniques to large numbers of neutrinos. △ Less

Submitted 27 April, 2021; v1 submitted 1 October, 2020; originally announced October 2020.

Report number: N3AS-20-009

Journal ref: Phys. Rev. D 103, 043006 (2021)

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

Authors: Eve Armstrong, Amol V. Patwardhan, Ermal Rrapaj, Sina Fallah Ardizi, George M. Fuller

Abstract: The multi-messenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational… ▽ More The multi-messenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning (ML) formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two mono-energetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein (MSW) resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results intimate the promise of this "variational annealing" methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access. △ Less

Submitted 21 August, 2020; v1 submitted 13 June, 2020; originally announced June 2020.

Comments: 14 pages, 6 figures

Journal ref: Phys. Rev. D 102, 043013 (2020)

The LUX-ZEPLIN (LZ) Experiment

Authors: The LZ Collaboration, D. S. Akerib, C. W. Akerlof, D. Yu. Akimov, A. Alquahtani, S. K. Alsum, T. J. Anderson, N. Angelides, H. M. Araújo, A. Arbuckle, J. E. Armstrong, M. Arthurs, H. Auyeung, X. Bai, A. J. Bailey, J. Balajthy, S. Balashov, J. Bang, M. J. Barry, J. Barthel, D. Bauer, P. Bauer, A. Baxter, J. Belle, P. Beltrame , et al. (357 additional authors not shown)

Abstract: We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient n… ▽ More We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements. △ Less

Submitted 3 November, 2019; v1 submitted 20 October, 2019; originally announced October 2019.

Velocity-resolved reverberation mapping of five bright Seyfert 1 galaxies

Authors: G. De Rosa, M. M. Fausnaugh, C. J. Grier, B. M. Peterson, K. D. Denney, Keith Horne, M. C. Bentz, S. Ciroi, E. Dalla Bonta`, M. D. Joner, S. Kaspi, C. S. Kochanek, R. W. Pogge, S. G. Sergeev, M. Vestergaard, S. M. Adams, J. Antognini, C. Araya Salvo, E. Armstrong, J. Bae, A. J. Barth, T. G. Beatty, A. Bhattacharjee, G. A. Borman, T. A. Boroson , et al. (77 additional authors not shown)

Abstract: We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line reg… ▽ More We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007. △ Less

Submitted 3 August, 2018; v1 submitted 12 July, 2018; originally announced July 2018.

Comments: Accepted for publication on ApJ; 32 pages, 16 figures, 10 tables

An optimization-based approach to calculating neutrino flavor evolution

Authors: Eve Armstrong, Amol V. Patwardhan, Lucas Johns, Chad T. Kishimoto, Henry D. I. Abarbanel, George M. Fuller

Abstract: We assess the utility of an optimization-based data assimilation (D.A.) technique for treating the problem of nonlinear neutrino flavor transformation in core collapse supernovae. D.A. uses measurements obtained from a physical system to estimate the state variable evolution and parameter values of the associated model. Formulated as an optimization procedure, D.A. can offer an integration-blind a… ▽ More We assess the utility of an optimization-based data assimilation (D.A.) technique for treating the problem of nonlinear neutrino flavor transformation in core collapse supernovae. D.A. uses measurements obtained from a physical system to estimate the state variable evolution and parameter values of the associated model. Formulated as an optimization procedure, D.A. can offer an integration-blind approach to predicting model evolution, which offers an advantage for models that thwart solution via traditional numerical integration techniques. Further, D.A. performs most optimally for models whose equations of motion are nonlinearly coupled. In this exploratory work, we consider a simple steady-state model with two mono-energetic neutrino beams coherently interacting with each other and a background medium. As this model can be solved via numerical integration, we have an independent consistency check for D.A. solutions. We find that the procedure can capture key features of flavor evolution over the entire trajectory, even given measurements of neutrino flavor only at the endpoint, and with an assumed known initial flavor distribution. Further, the procedure permits an examination of the sensitivity of flavor evolution to estimates of unknown model parameters, locates degeneracies in parameter space, and can identify the specific measurements required to break those degeneracies. △ Less

Submitted 17 August, 2017; v1 submitted 11 December, 2016; originally announced December 2016.

Comments: 26 pages, 4 figures

Journal ref: Phys. Rev. D 96, 083008 (2017)

Orbital, Superhump, and Superorbital Periods in the Cataclysmic Variables AQ Mensae and IM Eridani

Authors: E. Armstrong, J. Patterson, E. Michelsen, J. Thorstensen, H. Uthas, T. Vanmunster, F. -J. Hambsch, G. Roberts, S. Dvorak

Abstract: We report photometric detections of orbital and superorbital signals, and negative orbital sidebands, in the light curves of the nova-like cataclysmic variables AQ Mensae and IM Eridani. The frequencies of the orbital, superorbital, and sideband signals are 7.0686 (3), 0.263 (3), and 7.332 (3) cycles per day (c/d) in AQ Mensae, and 6.870 (1), 0.354 (7), and 7.226 (1) c/d in IM Eridani. We also fin… ▽ More We report photometric detections of orbital and superorbital signals, and negative orbital sidebands, in the light curves of the nova-like cataclysmic variables AQ Mensae and IM Eridani. The frequencies of the orbital, superorbital, and sideband signals are 7.0686 (3), 0.263 (3), and 7.332 (3) cycles per day (c/d) in AQ Mensae, and 6.870 (1), 0.354 (7), and 7.226 (1) c/d in IM Eridani. We also find a spectroscopic orbital frequency in IM Eridani of 6.86649 (2) c/d. These observations can be reproduced by invoking an accretion disc that is tilted with respect to the orbital plane. This model works well for X-ray binaries, in which irradiation by a primary neutron star can account for the disc's tilt. A likely tilt mechanism has yet to be identified in CVs, yet the growing collection of observational evidence indicates that the phenomenon of tilt is indeed at work in this class of object. The results presented in this paper bring the number of CVs known to display signals associated with retrograde disc precession to twelve. We also find AQ Mensae to be an eclipsing system. The eclipse depths are highly variable, which suggests that the eclipses are grazing. This finding raises the possibility of probing variations in disc tilt by studying systematic variations in the eclipse profile. △ Less

Submitted 7 August, 2013; v1 submitted 20 July, 2013; originally announced July 2013.

Comments: 18 pages, 9 figures, 5 tables

Non-detection of the Tooth Fairy at Optical Wavelengths

Authors: Eve Armstrong

Abstract: We report a non-detection, to a limiting magnitude of V = 18.4 (9), of the elusive entity commonly described as the Tooth Fairy. We review various physical models and conclude that follow-up observations must precede an interpretation of our result. We report a non-detection, to a limiting magnitude of V = 18.4 (9), of the elusive entity commonly described as the Tooth Fairy. We review various physical models and conclude that follow-up observations must precede an interpretation of our result. △ Less

Submitted 20 March, 2016; v1 submitted 2 April, 2012; originally announced April 2012.

Comments: 5 pages, 3 figures

Journal ref: Journal of Irreproducible Results. 2014, Volume 52, Number 3, pages 22-25

The photometric period in ES Ceti

Authors: C. M. Copperwheat, T. R. Marsh, V. S. Dhillon, S. P. Littlefair, P. A. Woudt, B. Warner, J. Patterson, D. Steeghs, J. Kemp, E. Armstrong, R. Rea

Abstract: We present ULTRACAM photometry of ES Cet, an ultracompact binary with a 620s orbital period. The mass transfer in systems such as this one is thought to be driven by gravitational radiation, which causes the binary to evolve to longer periods since the semi-degenerate donor star expands in size as it loses mass. We supplement these ULTRACAM+WHT data with observations made with smaller telescopes a… ▽ More We present ULTRACAM photometry of ES Cet, an ultracompact binary with a 620s orbital period. The mass transfer in systems such as this one is thought to be driven by gravitational radiation, which causes the binary to evolve to longer periods since the semi-degenerate donor star expands in size as it loses mass. We supplement these ULTRACAM+WHT data with observations made with smaller telescopes around the world over a nine year baseline. All of the observations show variation on the orbital period, and by timing this variation we track the period evolution of this system. We do not detect any significant departure from a linear ephemeris, implying a donor star that is of small mass and close to a fully degenerate state. This finding favours the double white dwarf formation channel for this AM CVn star. An alternative explanation is that the system is in the relatively short-lived phase in which the mass transfer rate climbs towards its long-term value. △ Less

Submitted 21 February, 2011; v1 submitted 18 January, 2011; originally announced January 2011.

Comments: Accepted for publication in MNRAS. 8 pages with 2 figures and 1 table. This revision corrects some typographical errors in Table 1

The New Eclipsing Cataclysmic Variable SDSS 154453+2553

Authors: Julie N. Skinner, John R. Thorstensen, Eve Armstrong, Steve Brady

Abstract: The cataclysmic variable SDSS154453+2553 was recently identified in the Sloan Digital Sky Survey. We obtained spectra and photometry at the MDM Observatory, which revealed an eclipse with a 6.03 hour period. The Hα emission line exhibits a strong rotational disturbance during eclipse, indicating that it arises in an accretion disk. A contribution from an M-type companion is also observed. Time-ser… ▽ More The cataclysmic variable SDSS154453+2553 was recently identified in the Sloan Digital Sky Survey. We obtained spectra and photometry at the MDM Observatory, which revealed an eclipse with a 6.03 hour period. The Hα emission line exhibits a strong rotational disturbance during eclipse, indicating that it arises in an accretion disk. A contribution from an M-type companion is also observed. Time-series photometry during eclipse gives an ephemeris of 2454878.0062(15) + 0.251282(2)E. We present spectroscopy through the orbit and eclipse photometry. Our analysis of the secondary star indicates a distance of 800 {\pm} 180 pc. △ Less

Submitted 7 January, 2011; originally announced January 2011.

Comments: 6 pages, 3 figures, Accepted for publication in PASP

Optical and X-Ray Observations of GRB 060526: A Complex Afterglow Consistent with An Achromatic Jet Break

Authors: X. Dai, J. P. Halpern, N. D. Morgan, E. Armstrong, N. Mirabal, J. B. Haislip, D. E. Reichart, K. Z. Stanek

Abstract: We obtained 98 R-band and 18 B, r', i' images of the optical afterglow of GRB 060526 (z=3.21) with the MDM 1.3m, 2.4m, and the PROMPT telescopes in Cerro Tololo over the 5 nights following the burst trigger. Combining these data with other optical observations reported in GCN and the Swift-XRT observations, we compare the optical and X-ray afterglow light curves of GRB 060526. Both the optical a… ▽ More We obtained 98 R-band and 18 B, r', i' images of the optical afterglow of GRB 060526 (z=3.21) with the MDM 1.3m, 2.4m, and the PROMPT telescopes in Cerro Tololo over the 5 nights following the burst trigger. Combining these data with other optical observations reported in GCN and the Swift-XRT observations, we compare the optical and X-ray afterglow light curves of GRB 060526. Both the optical and X-ray afterglow light curves show rich features, such as flares and breaks. The densely sampled optical observations provide very good coverage at T>1.e4 sec. We observed a break at 2.4e5 sec in the optical afterglow light curve. Compared with the X-ray afterglow light curve, the break is consistent with an achromatic break supporting the beaming models of GRBs. However, the pre-break and post-break temporal decay slopes are difficult to explain in simple afterglow models. We estimated a jet angle of θ_j ~ 7 degrees and a prompt emission size of R_{prompt} ~ 2e14 cm. In addition, we detected several optical flares with amplitudes of Δm ~ 0.2, 0.6, and 0.2 mag. The X-ray afterglows detected by Swift have shown complicated decay patterns. Recently, many well-sampled optical afterglows also show decays with flares and multiple breaks. GRB 060526 provides an additional case of such a complex, well observed optical afterglow. The accumulated well-sampled afterglows indicate that most of the optical afterglows are complex. △ Less

Submitted 17 January, 2007; v1 submitted 10 September, 2006; originally announced September 2006.

Comments: Accepted by ApJ, 21 pages, 4 figures

Journal ref: Astrophys.J.658:509-513,2007

Superhumps in Cataclysmic Binaries. XXV. q_crit, epsilon(q), and Mass-Radius

Authors: Joseph Patterson, Jonathan Kemp, David Harvey, Robert Fried, Robert Rea, Berto Monard, Lewis Cook, David Skillman, Tonny Vanmunster, Greg Bolt, Eve Armstrong, Jennie McCormick, Thomas Krajci, Lasse Jensen, Jerry Gunn, Neil Butterworth, Jerry Foote, Marc Bos, Gianluca Masi, Paul Warhurst

Abstract: We report on successes and failures in searching for positive superhumps in cataclysmic variables, and show the superhumping fraction as a function of orbital period. Basically, all short-period systems do, all long-period systems don't, and a 50% success rate is found at P_orb=3.1+-0.2 hr. We can use this to measure the critical mass ratio for the creation of superhumps. With a mass-radius rela… ▽ More We report on successes and failures in searching for positive superhumps in cataclysmic variables, and show the superhumping fraction as a function of orbital period. Basically, all short-period systems do, all long-period systems don't, and a 50% success rate is found at P_orb=3.1+-0.2 hr. We can use this to measure the critical mass ratio for the creation of superhumps. With a mass-radius relation appropriate for cataclysmic variables, and an assumed mean white-dwarf mass of 0.75 M_sol, we find a mass ratio q_crit=0.35+-0.02. We also report superhump studies of several stars of independently known mass ratio: OU Virginis, XZ Eridani, UU Aquarii, and KV UMa (= XTE J1118+480). The latter two are of special interest, because they represent the most extreme mass ratios for which accurate superhump measurements have been made. We use these to improve the epsilon(q) calibration, by which we can infer the elusive q from the easy-to-measure epsilon (the fractional period excess of P_superhump over P_orb). This relation allows mass and radius estimates for the secondary star in any CV showing superhumps. The consequent mass-radius law shows an apparent discontinuity in radius near 0.2 M_sol, as predicted by the disrupted magnetic braking model for the 2.1-2.7 hour period gap. This is effectively the "empirical main sequence" for CV secondaries. △ Less

Submitted 15 July, 2005; originally announced July 2005.

Comments: PDF, 45 pages, 9 tables, 12 figures; accepted, in press, to appear November 2005, PASP; more info at http://cba.phys.columbia.edu/

The Dwarf Nova PQ Andromedae

Authors: Joseph Patterson, John Thorstensen, Eve Armstrong, Arne Henden, Robert Hynes

Abstract: We report a photometric study of the WZ Sagittae-type dwarf nova PQ Andromedae. The light curve shows strong (0.05 mag full amplitude) signals with periods of 1263(1) and 634(1) s, and a likely double-humped signal with P=80.6(2) min. We interpret the first two as nonradial pulsation periods of the underlying white dwarf, and the last as the orbital period of the underlying binary. We estimate a… ▽ More We report a photometric study of the WZ Sagittae-type dwarf nova PQ Andromedae. The light curve shows strong (0.05 mag full amplitude) signals with periods of 1263(1) and 634(1) s, and a likely double-humped signal with P=80.6(2) min. We interpret the first two as nonradial pulsation periods of the underlying white dwarf, and the last as the orbital period of the underlying binary. We estimate a distance of 150(50) pc from proper motions and the two standard candles available: the white dwarf and the dwarf-nova outburst. At this distance, the K magnitude implies that the secondary is probably fainter than any star on the main sequence -- indicating a mass below the Kumar limit at 0.075 M_sol. PQ And may be another "period bouncer", where evolution now drives the binary out to longer period. △ Less

Submitted 7 June, 2005; originally announced June 2005.

Comments: PDF, 13 pages, 2 figures; accepted, in press, to appear September 2005, PASP; more info at http://cba.phys.columbia.edu/

Is FIRST J102347.6+003841 Really a Cataclysmic Binary?

Authors: J. R. Thorstensen, E. Armstrong

Abstract: The radio source FIRST J102347.6+003841 was presented as the first radio-selected cataclysmic. In the discovery paper, Bond et al. (2002) show a spectrum consistent with a magnetic AM Her-type system and a light curve with rapid, irregular flickering. In contrast, Woudt, Warner, and Pretorius (2004) found a smoothly-varying light curve with a period near 4.75 h and one minimum per orbit, indicat… ▽ More The radio source FIRST J102347.6+003841 was presented as the first radio-selected cataclysmic. In the discovery paper, Bond et al. (2002) show a spectrum consistent with a magnetic AM Her-type system and a light curve with rapid, irregular flickering. In contrast, Woudt, Warner, and Pretorius (2004) found a smoothly-varying light curve with a period near 4.75 h and one minimum per orbit, indicating a dramatic change. We present time-resolved spectra showing a superficially normal, mid-G type photosphere, with no detectable emission lines. The absorption-line radial velocity varies sinusoidally, with semi-amplitude 268 +- 4 km/s, on the orbital period, which is refined to 0.198094(2) d. At this orbital period the secondary's spectral type is atypically early, suggesting an unusual evolutionary history. We also obtained BVI photometry around the orbit. The light curve resembles that given by Woudt et al., and the color modulation is consistent with a heating effect. A simple illumination model matches the observations strikingly well with a Roche-lobe filling secondary near 5650 kelvin being illuminated by a primary giving out around 2 solar luminosities. The modest amplitude of the observed modulation constrains the orbital inclination to be less than about 55 deg, unless the gravity darkening is artificially reduced. The resulting primary star mass is above the Chandrasekhar limit (assuming conventional gravity darkening). We examine the possibility that the compact object in this system is not a white dwarf, in which case this is not actually a cataclysmic variable. On close examination, FIRST J102347.6+003841 defies easy classification. △ Less

Submitted 23 April, 2005; originally announced April 2005.

Comments: 24 pages, 5 postscript and two JPG figures; Astronomical Journal, accepted

Journal ref: Astron.J.130:759-766,2005

High Speed Photometry of SDSS J013701.06-091234.9

Authors: M. L. Pretorius, P. A. Woudt, B. Warner, G. Bolt, J. Patterson, E. Armstrong

Abstract: We present high speed photometry of the Sloan Digital Sky Survey cataclysmic variable SDSS J013701.06-091234.9 in quiescence and during its 2003 December superoutburst. The orbital modulation at 79.71\pm0.01 min is double humped; the superhump period is 81.702\pm0.007 min. Towards the end of the outburst late superhumps with a period of 81.29\pm0.01 min were observed. We argue that this is a sys… ▽ More We present high speed photometry of the Sloan Digital Sky Survey cataclysmic variable SDSS J013701.06-091234.9 in quiescence and during its 2003 December superoutburst. The orbital modulation at 79.71\pm0.01 min is double humped; the superhump period is 81.702\pm0.007 min. Towards the end of the outburst late superhumps with a period of 81.29\pm0.01 min were observed. We argue that this is a system of very low mass transfer rate, and that it probably has a long outburst interval. △ Less

Submitted 11 May, 2004; originally announced May 2004.

Comments: 5 pages, 8 figures. Accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 352 (2004) 1056

Rapid Oscillations in Cataclysmic Variables. XVI. DW Cancri

Authors: Joseph Patterson, John Thorstensen, Tonny Vanmunster, Robert Fried, Brian Martin, Tut Campbell, Jeff Robertson, Jonathan Kemp, David Messier, Eve Armstrong

Abstract: We report photometry and spectroscopy of the novalike variable DW Cancri. The spectra show the usual broad H and He emission lines, with an excitation and continuum slope characteristic of a moderately high accretion rate. A radial-velocity search yields strong detections at two periods, 86.1015(3) min and 38.58377(6) min. We interpret these as respectively the orbital period P_orb of the binary… ▽ More We report photometry and spectroscopy of the novalike variable DW Cancri. The spectra show the usual broad H and He emission lines, with an excitation and continuum slope characteristic of a moderately high accretion rate. A radial-velocity search yields strong detections at two periods, 86.1015(3) min and 38.58377(6) min. We interpret these as respectively the orbital period P_orb of the binary, and the spin period P_spin of a magnetic white dwarf. The light curve also shows the spin period, plus an additional strong signal at 69.9133(10) min, which coincides with the difference frequency 1/P_spin-1/P_orb. These periods are stable over the 1 year baseline of measurement. This triply-periodic structure mimics the behavior of several well-credentialed members of the "DQ Herculis" (intermediate polar) class of cataclysmic variables. DQ Her membership is also suggested by the mysteriously strong sideband signal (at nu_spin-nu_orb), attesting to a strong pulsed flux at X-ray/EUV/UV wavelengths. DW Cnc is a new member of this class, and would be an excellent target for extended observation at these wavelengths. △ Less

Submitted 5 April, 2004; originally announced April 2004.

Comments: PDF, 28 pages, 6 tables, 9 figures; accepted, in press, to appear June 2004, PASP; more info at http://cba.phys.columbia.edu/

GRB 021004: A Possible Shell Nebula around a Wolf-Rayet Star Gamma-Ray Burst Progenitor

Authors: N. Mirabal, J. P. Halpern, Ryan Chornock, Alexei V. Filippenko, D. M. Terndrup, E. Armstrong, J. Kemp, J. R. Thorstensen, M. Tavarez, C. Espaillat

Abstract: The rapid localization of GRB 021004 by the HETE-2 satellite allowed nearly continuous monitoring of its early optical afterglow decay, as well as high-quality optical spectra that determined a redshift of z3=2.328 for its host galaxy, an active starburst galaxy with strong Lyman-alpha emission and several absorption lines. Spectral observations show multiple absorbers at z3A=2.323, z3B= 2.317,… ▽ More The rapid localization of GRB 021004 by the HETE-2 satellite allowed nearly continuous monitoring of its early optical afterglow decay, as well as high-quality optical spectra that determined a redshift of z3=2.328 for its host galaxy, an active starburst galaxy with strong Lyman-alpha emission and several absorption lines. Spectral observations show multiple absorbers at z3A=2.323, z3B= 2.317, and z3C= 2.293 blueshifted by 450, 990, and 3,155 km/s respectively relative to the host galaxy Lyman-alpha emission. We argue that these correspond to a fragmented shell nebula that has been radiatively accelerated by the gamma-ray burst (GRB) afterglow at a distance greater than 0.3 pc from a Wolf-Rayet star progenitor. The chemical abundance ratios indicate that the nebula is overabundant in carbon and silicon. The high level of carbon and silicon is consistent with a swept-up shell nebula gradually enriched by a WCL progenitor wind over the lifetime of the nebula prior to the GRB onset. The detection of statistically significant fluctuations and color changes about the jet-like optical decay further supports this interpretation since fluctuations must be present at some level due to inhomogeneities in a clumpy stellar wind medium or if the progenitor has undergone massive ejection prior to the GRB onset. This evidence suggests that the mass-loss process in a Wolf-Rayet star might lead naturally to an iron-core collapse with sufficient angular momentum that could serve as a suitable GRB progenitor. △ Less

Submitted 11 June, 2003; v1 submitted 27 March, 2003; originally announced March 2003.

Comments: Replaced with version accepted by ApJ; 40 pages, 9 figures

Journal ref: Astrophys.J. 595 (2003) 935-949