Normalized to: Chan, C.
[1]
oai:arXiv.org:2003.02845 [pdf] - 2122644
A Novel Multi-Dimensional Boltzmann Neutrino Transport Scheme for
Core-Collapse Supernovae
Submitted: 2020-03-05, last modified: 2020-06-26
We introduce a new discrete-ordinate scheme for solving the general
relativistic Boltzmann transport equation in the context of core-collapse
supernovae. Our algorithm avoids the need to spell out the complicated
advection terms in energy and angle that arise when the transport equation is
formulated in spherical polar coordinates, in the comoving frame, or in a
general relativistic spacetime. We instead approach the problem by calculating
the advection of neutrinos across momentum space using an intuitive
particle-like approach that has excellent conservation properties and fully
accounts for Lorentz boosts, general relativistic effects, and grid geometry
terms. In order to avoid the need for a global implicit solution, time
integration is performed using a locally implicit Lax-Wendroff scheme that
correctly reproduces the diffusion limit. This will facilitate the use of our
method on massively parallel distributed-memory architectures. We have verified
the accuracy and stability of our scheme with a suite of test problems in
spherical symmetry and axisymmetry. To demonstrate that the new algorithm works
stably in core-collapse supernova simulations, we have coupled it to the
general relativistic hydrodynamics code CoCoNuT and present a first
demonstration run of a $20M_\odot$ progenitor with a reduced set of neutrino
opacities.
[3]
oai:arXiv.org:2005.14251 [pdf] - 2104086
A Plasmoid Model for the Sgr A* Flares Observed with GRAVITY and Chandra
Submitted: 2020-05-28
The Galactic Center black hole Sgr A* shows significant variability and
flares in the submillimeter, infrared, and X-ray wavelengths. Owing to its
exquisite resolution in the IR bands, the GRAVITY experiment for the first time
spatially resolved the locations of three flares and showed that a bright
region moves in ellipse-like trajectories close to but offset from the black
hole over the course of each event. We present a model for plasmoids that form
during reconnection events and orbit in the coronal region around a black hole
to explain these observations. We utilize general-relativistic radiative
transfer calculations that include effects from finite light travel time,
plasmoid motion, particle acceleration, and synchrotron cooling and obtain a
rich structure in the flare lightcurves. This model can naturally account for
the observed motion of the bright regions observed by the GRAVITY experiment
and the offset between the center of the centroid motion and the position of
the black hole. It also explains why some flares may be double-peaked while
others have only a single peak and uncovers a correlation between the structure
in the lightcurve and the location of the flare. Finally, we make predictions
for future observations of flares from the inner accretion flow of Sgr A* that
will provide a test of this model.
[4]
oai:arXiv.org:2003.04320 [pdf] - 2100423
The impact of fallback on the compact remnants and chemical yields of
core-collapse supernovae
Submitted: 2020-03-09, last modified: 2020-05-24
Fallback in core-collapse supernovae plays a crucial role in determining the
properties of the compact remnants and of the ejecta composition. We perform
three-dimensional simulations of mixing and fallback for selected non-rotating
supernova models to study how explosion energy and asymmetries correlate with
the remnant mass, remnant kick, and remnant spin. We find that the strongest
kick and spin is imparted by partial fallback in an asymmetric explosion. Black
hole (BH) kicks of several hundred $\mathrm{km}\,\mathrm{s}^{-1}$ and spin
parameters of $\mathord{\sim}0.25$ can be obtained in this scenario. If the
initial explosion energy barely exceeds the envelope binding energy, stronger
fallback results, and the remnant kick and spin remain small. If the explosion
energy is high with respect to the envelope binding energy, there is little
fallback with a small effect on the remnant kick, but the spin-up by fallback
can be substantial. For a non-rotating $12\,\mathrm{M}_\odot$ progenitor, we
find that the neutron star (NS) is spun up to millisecond periods. The high
specific angular momentum of the fallback material can also lead to disk
formation around black holes. Fallback may thus be a pathway towards
millisecond-magnetar or collapsar-type engines for hypernovae and gamma-ray
bursts that does not require rapid progenitor rotation. Within our small set of
simulations, none reproduced the peculiar layered fallback necessary to explain
the metal-rich iron-poor composition of many carbon-enhanced metal-poor (CEMP)
stars. Models with different explosion energy and different realisations of
asymmetries may, however, be compatible with CEMP abundance patterns.
[5]
oai:arXiv.org:2005.09632 [pdf] - 2097551
Markov Chains for Horizons (MARCH). I. Identifying Biases in Fitting
Theoretical Models to Event Horizon Telescope Observations
Psaltis, Dimitrios;
Ozel, Feryal;
Medeiros, Lia;
Christian, Pierre;
Kim, Junhan;
Chan, Chi-kwan;
Conway, Landen J.;
Raithel, Carolyn A.;
Marrone, Dan;
Lauer, Tod R.
Submitted: 2020-05-19
We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel
tempering for fitting theoretical models of horizon-scale images of black holes
to the interferometric data from the Event Horizon Telescope (EHT). The
algorithm implements forms of the noise distribution in the data that are
accurate for all signal-to-noise ratios. In addition to being trivially
parallelizable, the algorithm is optimized for high performance, achieving 1
million MCMC chain steps in under 20 seconds on a single processor. We use
synthetic data for the 2017 EHT coverage of M87 that are generated based on
analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model
images to explore several potential sources of biases in fitting models to
sparse interferometric data. We demonstrate that a very small number of data
points that lie near salient features of the interferometric data exert
disproportionate influence on the inferred model parameters. We also show that
the preferred orientations of the EHT baselines introduce significant biases in
the inference of the orientation of the model images. Finally, we discuss
strategies that help identify the presence and severity of such biases in
realistic applications.
[6]
oai:arXiv.org:2004.13029 [pdf] - 2085279
The Origin of hotspots around Sgr A*: Orbital or pattern motion
Submitted: 2020-04-27
The Gravity Collaboration detected a near-infrared hotspot moving around Sgr
A* during the 2018 July 22 flare. They fitted the partial loop the hotspot made
on the sky with a circular Keplerian orbit of radius $\simeq7.5\,r_{\rm g}$
around the supermassive black hole (BH), where $r_{\rm g}$ is the gravitational
radius. However, because the hotspot traversed the loop in a short time, models
in which the hotspot tracks the motion of some fluid element tend to produce a
best-fit trajectory smaller than the observed loop. This is true for a circular
Keplerian orbit, even when BH spin is accounted for, and for motion along a
RIAF streamline. A marginally bound geodesic suffers from the same problem; in
addition, it is not clear what the origin of an object following the geodesic
would be. The observed hotspot motion is more likely a pattern motion. Circular
motion with $r\simeq12.5\,r_{\rm g}$ and a super-Keplerian speed $\simeq0.8\,c$
is a good fit. Such motion must be pattern motion because it cannot be
explained by physical forces. The pattern speed is compatible with
magnetohydrodynamic perturbations, provided that the magnetic field is
sufficiently strong. Circular pattern motion of radius $\sim20\, r_{\rm g}$ on
a plane above the BH is an equally good alternative; in this case, the hotspot
may be caused by a precessing outflow interacting with a surrounding disk. As
all our fits have relatively large radii, we cannot constrain the BH spin using
these observations.
[8]
oai:arXiv.org:2004.06234 [pdf] - 2077492
Light curves of tidal disruption events in active galactic nuclei
Submitted: 2020-04-13
The black hole of an active galactic nucleus is encircled by an accretion
disk. The surface density of the disk is always too low to affect the tidal
disruption of a star, but it can be high enough that a vigorous interaction
results when the debris stream returns to pericenter and punches through the
disk. Shocks excited in the disk dissipate the kinetic energy of the disk
interior to the impact point and expedite inflow toward the black hole.
Radiatively efficient disks with luminosity $\gtrsim10^{-3}$ Eddington have
high enough surface density that the initial stream-disk interaction leads to
energy dissipation at a super-Eddington rate. Because of the rapid inflow, only
part of this dissipated energy emerges as radiation, while the rest is advected
into the black hole. Dissipation, inflow, and cooling balance to keep the
bolometric luminosity at an Eddington-level plateau that lasts tens of days.
After the plateau, the luminosity decreases in proportion to the disk surface
density, with a power-law index between $-3$ and $-2$ at earlier times, and
possibly a steeper index at later times.
[9]
oai:arXiv.org:2004.06210 [pdf] - 2077487
Discretization and Filtering Effects on Black Hole Images Obtained with
the Event Horizon Telescope
Submitted: 2020-04-13
Interferometers, such as the Event Horizon Telescope (EHT), do not directly
observe the images of sources but rather measure their Fourier components at
discrete spatial frequencies up to a maximum value set by the longest baseline
in the array. Construction of images from the Fourier components or analysis of
them with high-resolution models requires careful treatment of fine source
structure nominally beyond the array resolution. The primary EHT targets, Sgr
A* and M87, are expected to have black-hole shadows with sharp edges and
strongly filamentary emission from the surrounding plasma on scales much
smaller than those probed by the currently largest baselines. We show that for
aliasing not to affect images reconstructed with regularized maximum likelihood
methods and model images that are directly compared to the data, the sampling
of these images (i.e., their pixel spacing) needs to be significantly finer
than the scale probed by the largest baseline in the array. Using GRMHD
simulations of black-hole images, we estimate the maximum allowable pixel
spacing to be approximately equal to (1/8)GM/c^2; for both of the primary EHT
targets, this corresponds to an angular pixel size of <0.5 microarcseconds.
With aliasing under control, we then advocate use of the second-order
Butterworth filter with a cut-off scale equal to the maximum array baseline as
optimal for visualizing the reconstructed images. In contrast to the
traditional Gaussian filters, this Butterworth filter retains most of the power
at the scales probed by the array while suppressing the fine image details for
which no data exist.
[10]
oai:arXiv.org:2004.01161 [pdf] - 2074562
SYMBA: An end-to-end VLBI synthetic data generation pipeline
Roelofs, F.;
Janssen, M.;
Natarajan, I.;
Deane, R.;
Davelaar, J.;
Olivares, H.;
Porth, O.;
Paine, S. N.;
Bouman, K. L.;
Tilanus, R. P. J.;
van Bemmel, I. M.;
Falcke, H.;
Akiyama, K.;
Alberdi, A.;
Alef, W.;
Asada, K.;
Azulay, R.;
Baczko, A.;
Ball, D.;
Baloković, M.;
Barrett, J.;
Bintley, D.;
Blackburn, L.;
Boland, W.;
Bower, G. C.;
Bremer, M.;
Brinkerink, C. D.;
Brissenden, R.;
Britzen, S.;
Broderick, A. E.;
Broguiere, D.;
Bronzwaer, T.;
Byun, D.;
Carlstrom, J. E.;
Chael, A.;
Chan, C.;
Chatterjee, S.;
Chatterjee, K.;
Chen, M.;
Chen, Y.;
Cho, I.;
Christian, P.;
Conway, J. E.;
Cordes, J. M.;
Crew, G. B.;
Cui, Y.;
De Laurentis, M.;
Dempsey, J.;
Desvignes, G.;
Dexter, J.;
Doeleman, S. S.;
Eatough, R. P.;
Fish, V. L.;
Fomalont, E.;
Fraga-Encinas, R.;
Friberg, P.;
Fromm, C. M.;
Gómez, J. L.;
Galison, P.;
Gammie, C. F.;
García, R.;
Gentaz, O.;
Georgiev, B.;
Goddi, C.;
Gold, R.;
Gu, M.;
Gurwell, M.;
Hada, K.;
Hecht, M. H.;
Hesper, R.;
Ho, L. C.;
Ho, P.;
Honma, M.;
Huang, C. L.;
Huang, L.;
Hughes, D. H.;
Ikeda, S.;
Inoue, M.;
Issaoun, S.;
James, D. J.;
Jannuzi, B. T.;
Jeter, B.;
Jiang, W.;
Johnson, M. D.;
Jorstad, S.;
Jung, T.;
Karami, M.;
Karuppusamy, R.;
Kawashima, T.;
Keating, G. K.;
Kettenis, M.;
Kim, J.;
Kim, J.;
Kim, J.;
Kino, M.;
Koay, J. Yi;
Koch, P. M.;
Koyama, S.;
Kramer, M.;
Kramer, C.;
Krichbaum, T. P.;
Kuo, C.;
Lauer, T. R.;
Lee, S.;
Li, Y.;
Li, Z.;
Lindqvist, M.;
Lico, R.;
Liu, K.;
Liuzzo, E.;
Lo, W.;
Lobanov, A. P.;
Loinard, L.;
Lonsdale, C.;
Lu, R.;
MacDonald, N. R.;
Mao, J.;
Markoff, S.;
Marrone, D. P.;
Marscher, A. P.;
Martí-Vidal, I.;
Matsushita, S.;
Matthews, L. D.;
Medeiros, L.;
Menten, K. M.;
Mizuno, Y.;
Mizuno, I.;
Moran, J. M.;
Moriyama, K.;
Moscibrodzka, M.;
Müller, C.;
Nagai, H.;
Nagar, N. M.;
Nakamura, M.;
Narayan, R.;
Narayanan, G.;
Neri, R.;
Ni, C.;
Noutsos, A.;
Okino, H.;
Ortiz-León, G. N.;
Oyama, T.;
Özel, F.;
Palumbo, D. C. M.;
Patel, N.;
Pen, U.;
Pesce, D. W.;
Piétu, V.;
Plambeck, R.;
PopStefanija, A.;
Prather, B.;
Preciado-López, J. A.;
Psaltis, D.;
Pu, H.;
Ramakrishnan, V.;
Rao, R.;
Rawlings, M. G.;
Raymond, A. W.;
Rezzolla, L.;
Ripperda, B.;
Rogers, A.;
Ros, E.;
Rose, M.;
Roshanineshat, A.;
Rottmann, H.;
Roy, A. L.;
Ruszczyk, C.;
Ryan, B. R.;
Rygl, K. L. J.;
Sánchez, S.;
Sánchez-Arguelles, D.;
Sasada, M.;
Savolainen, T.;
Schloerb, F. Peter;
Schuster, K.;
Shao, L.;
Shen, Z.;
Small, D.;
Sohn, B. Won;
SooHoo, J.;
Tazaki, F.;
Tiede, P.;
Titus, M.;
Toma, K.;
Torne, P.;
Trent, T.;
Trippe, S.;
Tsuda, S.;
van Langevelde, H. J.;
van Rossum, D. R.;
Wagner, J.;
Wardle, J.;
Ward-Thompson, D.;
Weintroub, J.;
Wex, N.;
Wharton, R.;
Wielgus, M.;
Wong, G. N.;
Wu, Q.;
Young, A.;
Young, K.;
Younsi, Z.;
Yuan, F.;
Yuan, Y.;
Zensus, J. A.;
Zhao, G.;
Zhao, S.;
Zhu, Z.
Submitted: 2020-04-02
Realistic synthetic observations of theoretical source models are essential
for our understanding of real observational data. In using synthetic data, one
can verify the extent to which source parameters can be recovered and evaluate
how various data corruption effects can be calibrated. These studies are
important when proposing observations of new sources, in the characterization
of the capabilities of new or upgraded instruments, and when verifying
model-based theoretical predictions in a comparison with observational data. We
present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a
novel synthetic data generation pipeline for Very Long Baseline Interferometry
(VLBI) observations. SYMBA takes into account several realistic atmospheric,
instrumental, and calibration effects. We used SYMBA to create synthetic
observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has
recently captured the first image of a black hole shadow. After testing SYMBA
with simple source and corruption models, we study the importance of including
all corruption and calibration effects. Based on two example general
relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed
case studies to assess the attainable image quality with the current and future
EHT array for different weather conditions. The results show that the effects
of atmospheric and instrumental corruptions on the measured visibilities are
significant. Despite these effects, we demonstrate how the overall structure of
the input models can be recovered robustly after performing calibration steps.
With the planned addition of new stations to the EHT array, images could be
reconstructed with higher angular resolution and dynamic range. In our case
study, these improvements allowed for a distinction between a thermal and a
non-thermal GRMHD model based on salient features in reconstructed images.
[11]
oai:arXiv.org:2003.11141 [pdf] - 2070172
An Effective Field Theory Analysis of the First LUX Dark Matter Search
Akerib, D. S.;
Alsum, S.;
Araújo, H. M.;
Bai, X.;
Balajthy, J.;
Baxter, A.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Boxer, B.;
Brás, P.;
Burdin, S.;
Byram, D.;
Carmona-Benitez, M. C.;
Chan, C.;
Cutter, J. E.;
de Viveiros, L.;
Druszkiewicz, E.;
Fan, A.;
Fiorucci, S.;
Gaitskell, R. J.;
Ghag, C.;
Gilchriese, M. G. D.;
Gwilliam, C.;
Hall, C. R.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Jacobsen, R. G.;
Jahangir, O.;
Ji, W.;
Kamdin, K.;
Kazkaz, K.;
Khaitan, D.;
Korolkova, E. V.;
Kravitz, S.;
Kudryavtsev, V. A.;
Larsen, N. A.;
Leason, E.;
Lenardo, B. G.;
Lesko, K. T.;
Liao, J.;
Lin, J.;
Lindote, A.;
Lopes, M. I.;
Manalaysay, A.;
Mannino, R. L.;
Marangou, N.;
McKinsey, D. N.;
Mei, D. -M.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Naylor, A.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
Nilima, A.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pease, E. K.;
Riffard, Q.;
Rischbieter, G. R. C.;
Rhyne, C.;
Rossiter, P.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, R.;
Taylor, W. C.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tvrznikova, L.;
Utku, U.;
Uvarov, S.;
Vacheret, A.;
Velan, V.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Woodward, D.;
Xu, J.;
Zhang, C.
Submitted: 2020-03-24
The Large Underground Xenon (LUX) dark matter search was a 250-kg active mass
dual-phase time projection chamber that operated by detecting light and
ionization signals from particles incident on a xenon target. In December 2015,
LUX reported a minimum 90% upper C.L. of 6e-46 cm^2 on the spin-independent
WIMP-nucleon elastic scattering cross section based on a 1.4e4 kg*day exposure
in its first science run. Tension between experiments and the absence of a
definitive positive detection suggest it would be prudent to search for WIMPs
outside the standard spin-independent/spin-dependent paradigm. Recent
theoretical work has identified a complete basis of 14 independent effective
field theory (EFT) operators to describe WIMP-nucleon interactions. In addition
to spin-independent and spin-dependent nuclear responses, these operators can
produce novel responses such as angular-momentum-dependent and spin-orbit
couplings. Here we report on a search for all 14 of these EFT couplings with
data from LUX's first science run. Limits are placed on each coupling as a
function of WIMP mass.
[13]
oai:arXiv.org:1904.12261 [pdf] - 2036238
Tidal disruption events in active galactic nuclei
Submitted: 2019-04-28, last modified: 2020-01-23
A fraction of tidal disruption events (TDEs) occur in active galactic nuclei
(AGNs) whose black holes possess accretion disks; these TDEs can be confused
with common AGN flares. The disruption itself is unaffected by the disk, but
the evolution of the bound debris stream is modified by its collision with the
disk when it returns to pericenter. The outcome of the collision is largely
determined by the ratio of the stream mass current to the azimuthal mass
current of the disk rotating underneath the stream footprint, which in turn
depends on the mass and luminosity of the AGN. To characterize TDEs in AGNs, we
simulated a suite of stream--disk collisions with various mass current ratios.
The collision excites shocks in the disk, leading to inflow and energy
dissipation orders of magnitude above Eddington; however, much of the radiation
is trapped in the inflow and advected into the black hole, so the actual
bolometric luminosity may be closer to Eddington. The emergent spectrum may not
be thermal, TDE-like, or AGN-like. The rapid inflow causes the disk interior to
the impact point to be depleted within a fraction of the mass return time. If
the stream is heavy enough to penetrate the disk, part of the outgoing material
eventually hits the disk again, dissipating its kinetic energy in the second
collision; another part becomes unbound, emitting synchrotron radiation as it
shocks with surrounding gas.
[14]
oai:arXiv.org:1907.06272 [pdf] - 2050176
Extending light WIMP searches to single scintillation photons in LUX
Akerib, D. S.;
Alsum, S.;
Araújo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Baxter, A.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Boxer, B.;
Brás, P.;
Burdin, S.;
Byram, D.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Chiller, A. A.;
Chiller, C.;
Currie, A.;
Cutter, J. E.;
de Viveiros, L.;
Dobi, A.;
Dobson, J. E. Y.;
Druszkiewicz, E.;
Edwards, B. N.;
Faham, C. H.;
Fallon, S. R.;
Fan, A.;
Fiorucci, S.;
Gaitskell, R. J.;
Gehman, V. M.;
Genovesi, J.;
Ghag, C.;
Gibson, K. R.;
Gilchriese, M. G. D.;
Grace, E.;
Gwilliam, C.;
Hall, C. R.;
Hanhardt, M.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Jacobsen, R. G.;
Jahangir, O.;
Ji, W.;
Kamdin, K.;
Kazka, K.;
Khaitan, D.;
Knoche, R.;
Korolkova, E. V.;
Kravitz, S.;
Kudryavtsev, V. A.;
Larsen, N. A.;
Leason, E.;
Lee, C.;
Lenardo, B. G.;
Lesko, K. T.;
Levy, C.;
Liao, J.;
Lin, J.;
Lindote, A.;
Lopes, M. I.;
López-Paredes, B.;
Manalaysay, A.;
Mannino, R. L.;
Marangou, N.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. M.;
Mock, J.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Naylor, A.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
Nilima, A.;
O'Sullivan, K.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pease, E. K.;
Reichhart, L.;
Riffard, Q.;
Rischbieter, G. R. C.;
Rossiter, P.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Stephenson, S.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, R.;
Taylor, W. C.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Utku, U.;
Uvarov, S.;
Vacheret, A.;
Velan, V.;
Verbus, J. R.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Woodward, D.;
Xu, J.;
Yazdani, K.;
Young, S. K.;
Zhang, C.
Submitted: 2019-07-14, last modified: 2019-12-27
We present a novel analysis technique for liquid xenon time projection
chambers that allows for a lower threshold by relying on events with a prompt
scintillation signal consisting of single detected photons. The energy
threshold of the LUX dark matter experiment is primarily determined by the
smallest scintillation response detectable, which previously required a 2-fold
coincidence signal in its photomultiplier arrays, enforced in data analysis.
The technique presented here exploits the double photoelectron emission effect
observed in some photomultiplier models at vacuum ultraviolet wavelengths. We
demonstrate this analysis using an electron recoil calibration dataset and
place new constraints on the spin-independent scattering cross section of
weakly interacting massive particles (WIMPs) down to 2.5 GeV/c$^2$ WIMP mass
using the 2013 LUX dataset. This new technique is promising to enhance light
WIMP and astrophysical neutrino searches in next-generation liquid xenon
experiments.
[16]
oai:arXiv.org:1802.06039 [pdf] - 2061504
Projected WIMP sensitivity of the LUX-ZEPLIN (LZ) dark matter experiment
Akerib, D. S.;
Akerlof, C. W.;
Alsum, S. K.;
Araújo, H. M.;
Arthurs, M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Balashov, S.;
Bauer, D.;
Belle, J.;
Beltrame, P.;
Benson, T.;
Bernard, E. P.;
Biesiadzinski, T. P.;
Boast, K. E.;
Boxer, B.;
Brás, P.;
Buckley, J. H.;
Bugaev, V. V.;
Burdin, S.;
Busenitz, J. K.;
Carels, C.;
Carlsmith, D. L.;
Carlson, B.;
Carmona-Benitez, M. C.;
Chan, C.;
Cherwinka, J. J.;
Cole, A.;
Cottle, A.;
Craddock, W. W.;
Currie, A.;
Cutter, J. E.;
Dahl, C. E.;
de Viveiros, L.;
Dobi, A.;
Dobson, J. E. Y.;
Druszkiewicz, E.;
Edberg, T. K.;
Edwards, W. R.;
Fan, A.;
Fayer, S.;
Fiorucci, S.;
Fruth, T.;
Gaitskell, R. J.;
Genovesi, J.;
Ghag, C.;
Gilchriese, M. G. D.;
van der Grinten, M. G. D.;
Hall, C. R.;
Hans, S.;
Hanzel, K.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hillbrand, S.;
Hjemfelt, C.;
Hoff, M. D.;
Hor, J. Y-K.;
Huang, D. Q.;
Ignarra, C. M.;
Ji, W.;
Kaboth, A. C.;
Kamdin, K.;
Keefner, J.;
Khaitan, D.;
Khazov, A.;
Kim, Y. D.;
Kocher, C. D.;
Korolkova, E. V.;
Kraus, H.;
Krebs, H. J.;
Kreczko, L.;
Krikler, B.;
Kudryavtsev, V. A.;
Kyre, S.;
Lee, J.;
Lenardo, B. G.;
Leonard, D. S.;
Lesko, K. T.;
Levy, C.;
Li, J.;
Liao, J.;
Liao, F. -T.;
Lin, J.;
Lindote, A.;
Linehan, R.;
Lippincott, W. H.;
Liu, X.;
Lopes, M. I.;
Paredes, B. López;
Lorenzon, W.;
Luitz, S.;
Lyle, J. M.;
Majewski, P.;
Manalaysay, A.;
Mannino, R. L.;
Maupin, C.;
McKinsey, D. N.;
Meng, Y.;
Miller, E. H.;
Mock, J.;
Monzani, M. E.;
Morad, J. A.;
Morrison, E.;
Mount, B. J.;
Murphy, A. St. J.;
Nelson, H. N.;
Neves, F.;
Nikoleyczik, J.;
O'Sullivan, K.;
Olcina, I.;
Olevitch, M. A.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Patton, S. J.;
Pease, E. K.;
Penning, B.;
Piepke, A.;
Powell, S.;
Preece, R. M.;
Pushkin, K.;
Ratcliff, B. N.;
Reichenbacher, J.;
Rhyne, C. A.;
Richards, A.;
Rodrigues, J. P.;
Rosero, R.;
Rossiter, P.;
Saba, J. S.;
Sarychev, M.;
Schnee, R. W.;
Schubnell, M.;
Scovell, P. R.;
Shaw, S.;
Shutt, T. A.;
Silk, J. J.;
Silva, C.;
Skarpaas, K.;
Skulski, W.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Stancu, I.;
Stark, M. R.;
Stiegler, T. M.;
Stifter, K.;
Szydagis, M.;
Taylor, W. C.;
Taylor, R.;
Taylor, D. J.;
Temples, D.;
Terman, P. A.;
Thomas, K. J.;
Timalsina, M.;
To, W. H.;
Tomás, A.;
Tope, T. E.;
Tripathi, M.;
Tull, C. E.;
Tvrznikova, L.;
Utku, U.;
Va'vra, J.;
Vacheret, A.;
Verbus, J. R.;
Voirin, E.;
Waldron, W. L.;
Watson, J. R.;
Webb, R. C.;
White, D. T.;
Whitis, T. J.;
Wisniewski, W. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Woodward, D.;
Worm, S. D.;
Yeh, M.;
Yin, J.;
Young, I.
Submitted: 2018-02-16, last modified: 2019-12-02
LUX-ZEPLIN (LZ) is a next generation dark matter direct detection experiment
that will operate 4850 feet underground at the Sanford Underground Research
Facility (SURF) in Lead, South Dakota, USA. Using a two-phase xenon detector
with an active mass of 7 tonnes, LZ will search primarily for low-energy
interactions with Weakly Interacting Massive Particles (WIMPs), which are
hypothesized to make up the dark matter in our galactic halo. In this paper,
the projected WIMP sensitivity of LZ is presented based on the latest
background estimates and simulations of the detector. For a 1000 live day run
using a 5.6 tonne fiducial mass, LZ is projected to exclude at 90% confidence
level spin-independent WIMP-nucleon cross sections above $1.6 \times 10^{-48}$
cm$^{2}$ for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP. Additionally, a $5\sigma$
discovery potential is projected reaching cross sections below the existing and
projected exclusion limits of similar experiments that are currently operating.
For spin-dependent WIMP-neutron(-proton) scattering, a sensitivity of $2.7
\times 10^{-43}$ cm$^{2}$ ($8.1 \times 10^{-42}$ cm$^{2}$) for a 40
$\mathrm{GeV}/c^{2}$ mass WIMP is expected. With underground installation well
underway, LZ is on track for commissioning at SURF in 2020.
[18]
oai:arXiv.org:1910.10193 [pdf] - 1985370
First M87 Event Horizon Telescope Results and the Role of ALMA
Goddi, Ciriaco;
Crew, Geoff;
Impellizzeri, Violette;
Marti-Vidal, Ivan;
Matthews, Lynn D.;
Messias, Hugo;
Rottmann, Helge;
Alef, Walter;
Blackburn, Lindy;
Bronzwaer, Thomas;
Chan, Chi-Kwan;
Davelaar, Jordy;
Deane, Roger;
Dexter, Jason;
Doeleman, Shep;
Falcke, Heino;
Fish, Vincent L.;
Fraga-Encinas, Raquel;
Fromm, Christian M.;
Herrero-Illana, Ruben;
Issaoun, Sara;
James, David;
Janssen, Michael;
Kramer, Michael;
Krichbaum, Thomas P.;
De Laurentis, Mariafelicia;
Liuzzo, Elisabetta;
Mizuno, Yosuke;
Moscibrodzka, Monika;
Natarajan, Iniyan;
Porth, Oliver;
Rezzolla, Luciano;
Rygl, Kazi;
Roelofs, Freek;
Ros, Eduardo;
Roy, Alan L.;
Shao, Lijing;
van Langevelde, Huib Jan;
van Bemmel, Ilse;
Tilanus, Remo;
Torne, Pablo;
Wielgus, Maciek;
Younsi, Ziri;
Zensus, J. Anton
Submitted: 2019-10-22
In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the
first image of the candidate super-massive black hole (SMBH) at the centre of
the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image
shows a ring of glowing plasma with a dark patch at the centre, which is
interpreted as the shadow of the black hole. This breakthrough result, which
represents a powerful confirmation of Einstein's theory of gravity, or general
relativity, was made possible by assembling a global network of radio
telescopes operating at millimetre wavelengths that for the first time included
the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA
as an anchor station has enabled a giant leap forward by increasing the
sensitivity limits of the EHT by an order of magnitude, effectively turning it
into an imaging array. The published image demonstrates that it is now possible
to directly study the event horizon shadows of SMBHs via electromagnetic
radiation, thereby transforming this elusive frontier from a mathematical
concept into an astrophysical reality. The expansion of the array over the next
few years will include new stations on different continents - and eventually
satellites in space. This will provide progressively sharper and
higher-fidelity images of SMBH candidates, and potentially even movies of the
hot plasma orbiting around SMBHs. These improvements will shed light on the
processes of black hole accretion and jet formation on event-horizon scales,
thereby enabling more precise tests of general relativity in the truly strong
field regime.
[19]
oai:arXiv.org:1910.08376 [pdf] - 1982431
The Growing Importance of a Tech Savvy Astronomy and Astrophysics
Workforce
Norman, Dara;
Cruz, Kelle;
Desai, Vandana;
Lundgren, Britt;
Bellm, Eric;
Economou, Frossie;
Smith, Arfon;
Bauer, Amanda;
Nord, Brian;
Schafer, Chad;
Narayan, Gautham;
Li, Ting;
Tollerud, Erik;
Sipocz, Brigitta;
Stevance, Heloise;
Pickering, Timothy;
Sinha, Manodeep;
Harrington, Joseph;
Kartaltepe, Jeyhan;
Vohl, Dany;
Price-Whelan, Adrian;
Cherinka, Brian;
Chan, Chi-kwan;
Weiner, Benjamin;
Modjaz, Maryam;
Bianco, Federica;
Kerzendorf, Wolfgang;
Laginja, Iva;
Dong, Chuanfei
Submitted: 2019-10-17
Fundamental coding and software development skills are increasingly necessary
for success in nearly every aspect of astronomical and astrophysical research
as large surveys and high resolution simulations become the norm. However,
professional training in these skills is inaccessible or impractical for many
members of our community. Students and professionals alike have been expected
to acquire these skills on their own, apart from formal classroom curriculum or
on-the-job training. Despite the recognized importance of these skills, there
is little opportunity to develop them - even for interested researchers. To
ensure a workforce capable of taking advantage of the computational resources
and the large volumes of data coming in the next decade, we must identify and
support ways to make software development training widely accessible to
community members, regardless of affiliation or career level. To develop and
sustain a technology capable astronomical and astrophysical workforce, we
recommend that agencies make funding and other resources available in order to
encourage, support and, in some cases, require progress on necessary training,
infrastructure and policies. In this white paper, we focus on recommendations
for how funding agencies can lead in the promotion of activities to support the
astronomy and astrophysical workforce in the 2020s.
[20]
oai:arXiv.org:1811.11241 [pdf] - 1981692
Results of a Search for Sub-GeV Dark Matter Using 2013 LUX Data
Akerib, D. S.;
Alsum, S.;
Araújo, H. M.;
Bai, X.;
Balajthy, J.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Boxer, B.;
Brás, P.;
Burdin, S.;
Byram, D.;
Carmona-Benitez, M. C.;
Chan, C.;
Cutter, J. E.;
Davison, T. J. R.;
Druszkiewicz, E.;
Fallon, S. R.;
Fan, A.;
Fiorucci, S.;
Gaitskell, R. J.;
Genovesi, J.;
Ghag, C.;
Gilchriese, M. G. D.;
Gwilliam, C.;
Hall, C. R.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Jacobsen, R. G.;
Jahangir, O.;
Ji, W.;
Kamdin, K.;
Kazkaz, K.;
Khaitan, D.;
Knoche, R.;
Korolkova, E. V.;
Kravitz, S.;
Kudryavtsev, V. A.;
Lenardo, B. G.;
Lesko, K. T.;
Liao, J.;
Lin, J.;
Lindote, A.;
Lopes, M. I.;
Manalaysay, A.;
Mannino, R. L.;
Marangou, N.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. -M.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Naylor, A.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pease, E. K.;
Riffard, Q.;
Rischbieter, G. R. C.;
Rhyne, C.;
Rossiter, P.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, W. C.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Utku, U.;
Uvarov, S.;
Velan, V.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Woodward, D.;
Xu, J.;
Yazdani, K.;
Zhang, C.
Submitted: 2018-11-27, last modified: 2019-10-16
The scattering of dark matter (DM) particles with sub-GeV masses off nuclei
is difficult to detect using liquid xenon-based DM search instruments because
the energy transfer during nuclear recoils is smaller than the typical detector
threshold. However, the tree-level DM-nucleus scattering diagram can be
accompanied by simultaneous emission of a Bremsstrahlung photon or a so-called
"Migdal" electron. These provide an electron recoil component to the
experimental signature at higher energies than the corresponding nuclear
recoil. The presence of this signature allows liquid xenon detectors to use
both the scintillation and the ionization signals in the analysis where the
nuclear recoil signal would not be otherwise visible. We report constraints on
spin-independent DM-nucleon scattering for DM particles with masses of 0.4-5
GeV/c$^2$ using 1.4$\times10^4$ kg$\cdot$day of search exposure from the 2013
data from the Large Underground Xenon (LUX) experiment for four different
classes of mediators. This analysis extends the reach of liquid xenon-based DM
search instruments to lower DM masses than has been achieved previously.
[21]
oai:arXiv.org:1903.08832 [pdf] - 1950552
EHT-HOPS pipeline for millimeter VLBI data reduction
Blackburn, Lindy;
Chan, Chi-kwan;
Crew, Geoffrey B.;
Fish, Vincent L.;
Issaoun, Sara;
Johnson, Michael D.;
Wielgus, Maciek;
Akiyama, Kazunori;
Barrett, John;
Bouman, Katherine L.;
Cappallo, Roger;
Chael, Andrew A.;
Janssen, Michael;
Lonsdale, Colin J.;
Doeleman, Sheperd S.
Submitted: 2019-03-21, last modified: 2019-08-28
We present the design and implementation of an automated data calibration and
reduction pipeline for very-long-baseline interferometric (VLBI) observations
taken at millimeter wavelengths. These short radio-wavelengths provide the best
imaging resolution available from ground-based VLBI networks such as the Event
Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA), but
require specialized processing due to the strong effects from atmospheric
opacity and turbulence as well as the heterogeneous nature of existing global
arrays. The pipeline builds upon a calibration suite (HOPS) originally designed
for precision geodetic VLBI. To support the reduction of data for astronomical
observations, we have developed an additional framework for global phase and
amplitude calibration which provides output in a standard data format for
astronomical imaging and analysis. We test the pipeline on observations taken
at 3.5 mm (86 GHz) by the GMVA joined by the phased Atacama Large
Millimeter/submillimeter Array in April 2017, and demonstrate the benefits from
the specialized processing of high frequency VLBI data with respect to
classical analysis techniques.
[22]
oai:arXiv.org:1908.03991 [pdf] - 1986239
Prospects for Wideband VLBI Correlation in the Cloud
Submitted: 2019-08-12
This paper proposes a cloud architecture for the correlation of wide
bandwidth VLBI data. Cloud correlation facilitates processing of entire
experiments in parallel using flexibly allocated and practically unlimited
compute resources. This approach offers a potential improvement over dedicated
correlation clusters, which are constrained by a fixed number of installed
processor nodes and playback units. Additionally, cloud storage offers an
alternative to maintaining a fleet of hard-disk drives that might be utilized
intermittently. We describe benchmarks of VLBI correlation using the DiFX-2.5.2
software on the Google Cloud Platform to assess cloud-based correlation
performance. The number of virtual CPUs per Virtual Machine was varied to
determine the optimum configuration of cloud resources. The number of stations
was varied to determine the scaling of correlation time with VLBI arrays of
different sizes. Data transfer rates from Google Cloud Storage to the Virtual
Machines performing the correlation were also measured. We also present an
example cloud correlation configuration. Current cloud service and equipment
pricing data is used to compile cost estimates allowing an approximate economic
comparison between cloud and cluster processing. The economic comparisons are
based on cost figures which are a moving target, and are highly dependent on
factors such as the utilization of cluster and media, which are a challenge to
estimate. Our model suggests that shifting to the cloud is an alternative path
for high data rate, low duty cycle wideband VLBI correlation that should
continue to be explored. In the production phase of VLBI correlation, the cloud
has the potential to significantly reduce data processing times and allow the
processing of more science experiments in a given year for the petabyte-scale
data sets increasingly common in both astronomy and geodesy VLBI applications.
[23]
oai:arXiv.org:1904.04923 [pdf] - 1929749
The Event Horizon General Relativistic Magnetohydrodynamic Code
Comparison Project
Porth, Oliver;
Chatterjee, Koushik;
Narayan, Ramesh;
Gammie, Charles F.;
Mizuno, Yosuke;
Anninos, Peter;
Baker, John G.;
Bugli, Matteo;
Chan, Chi-kwan;
Davelaar, Jordy;
Del Zanna, Luca;
Etienne, Zachariah B.;
Fragile, P. Chris;
Kelly, Bernard J.;
Liska, Matthew;
Markoff, Sera;
McKinney, Jonathan C.;
Mishra, Bhupendra;
Noble, Scott C.;
Olivares, Héctor;
Prather, Ben;
Rezzolla, Luciano;
Ryan, Benjamin R.;
Stone, James M.;
Tomei, Niccolò;
White, Christopher J.;
Younsi, Ziri;
Collaboration, The Event Horizon Telescope
Submitted: 2019-04-09, last modified: 2019-08-05
Recent developments in compact object astrophysics, especially the discovery
of merging neutron stars by LIGO, the imaging of the black hole in M87 by the
Event Horizon Telescope (EHT) and high precision astrometry of the Galactic
Center at close to the event horizon scale by the GRAVITY experiment motivate
the development of numerical source models that solve the equations of general
relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for
the evolution of a magnetized accretion flow where turbulence is promoted by
the magnetorotational instability from a set of nine GRMHD codes: Athena++,
BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL.
Agreement between the codes improves as resolution increases, as measured by a
consistently applied, specially developed set of code performance metrics. We
conclude that the community of GRMHD codes is mature, capable, and consistent
on these test problems.
[24]
oai:arXiv.org:1907.06981 [pdf] - 1917113
Astro2020 APC White Paper: Elevating the Role of Software as a Product
of the Research Enterprise
Smith, Arfon M.;
Norman, Dara;
Cruz, Kelle;
Desai, Vandana;
Bellm, Eric;
Lundgren, Britt;
Economou, Frossie;
Nord, Brian D.;
Schafer, Chad;
Narayan, Gautham;
Harrington, Joseph;
Tollerud, Erik;
Sipőcz, Brigitta;
Pickering, Timothy;
Peeples, Molly S.;
Berriman, Bruce;
Teuben, Peter;
Rodriguez, David;
Gradvohl, Andre;
Shamir, Lior;
Allen, Alice;
Brownstein, Joel R.;
Ginsburg, Adam;
Sinha, Manodeep;
Hummels, Cameron;
Smith, Britton;
Stevance, Heloise;
Price-Whelan, Adrian;
Cherinka, Brian;
Chan, Chi-kwan;
Kartaltepe, Jeyhan;
Turk, Matthew;
Weiner, Benjamin;
Modjaz, Maryam;
Nemiroff, Robert J.;
Kerzendorf, Wolfgang;
Laginja, Iva;
Dong, Chuanfei;
Merín, Bruno;
Sobeck, Jennifer;
Buzasi, Derek;
Faherty, Jacqueline K;
Momcheva, Ivelina;
Connolly, Andrew;
Golkhou, V. Zach
Submitted: 2019-07-14
Software is a critical part of modern research, and yet there are
insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and
measure the impact of research software. The majority of academic fields rely
on a one-dimensional credit model whereby academic articles (and their
associated citations) are the dominant factor in the success of a researcher's
career. In the petabyte era of astronomical science, citing software and
measuring its impact enables academia to retain and reward researchers that
make significant software contributions. These highly skilled researchers must
be retained to maximize the scientific return from petabyte-scale datasets.
Evolving beyond the one-dimensional credit model requires overcoming several
key challenges, including the current scholarly ecosystem and scientific
culture issues. This white paper will present these challenges and suggest
practical solutions for elevating the role of software as a product of the
research enterprise.
[25]
oai:arXiv.org:1907.04460 [pdf] - 1914294
Astro2020: Training the Future Generation of Computational Researchers
Besla, Gurtina;
Huppenkothen, Daniela;
Lloyd-Ronning, Nicole;
Schneider, Evan;
Behroozi, Peter;
Burkhart, Blakesley;
Chan, C. K.;
Jacobson, Seth A.;
Morrison, Sarah;
Nam, Hai Ah;
Naoz, Smadar;
Peter, Annika;
Ramirez-Ruiz, Enrico
Submitted: 2019-07-09
The current disparity in computational knowledge is a critical hindrance to
the diversity and success of the field. Recommendations are outlined for
policies and funding models to enable the growth and retention of a new
generation of computational researchers that reflect the demographics of the
undergraduate population in Astronomy and Physics.
[26]
oai:arXiv.org:1903.12372 [pdf] - 1916847
Improved Measurements of the \b{eta}-Decay Response of Liquid Xenon with
the LUX Detector
Akerib, D. S.;
Alsum, S.;
Araújo, H. M.;
Bai, X.;
Balajthy, J.;
Baxter, A.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Boxer, B.;
Brás, P.;
Burdin, S.;
Byram, D.;
Carmona-Benitez, M. C.;
Chan, C.;
Cutter, J. E.;
de Viveiros, L.;
Druszkiewicz, E.;
Fallon, S. R.;
Fan, A.;
Fiorucci, S.;
Gaitskell, R. J.;
Genovesi, J.;
Ghag, C.;
Gilchriese, M. G. D.;
Gwilliam, C.;
Hall, C. R.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Jacobsen, R. G.;
Jahangir, O.;
Ji, W.;
Kamdin, K.;
Kazkaz, K.;
Khaitan, D.;
Korolkova, E. V.;
Kravitz, S.;
Kudryavtsev, V. A.;
Leason, E.;
Lenardo, B. G.;
Lesko, K. T.;
Liao, J.;
Lin, J.;
Lindote, A.;
Lopes, M. I.;
Manalaysay, A.;
Mannino, R. L.;
Marangou, N.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. M.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Naylor, A.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
Nilima, A.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pease, E. K.;
Riffard, Q.;
Rischbieter, G. R. C.;
Rhyne, C.;
Rossiter, P.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, R.;
Taylor, W. C.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Utku, U.;
Uvarov, S.;
Vacheret, A.;
Velan, V.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Woodward, D.;
Xu, J.;
Zhang, C.
Submitted: 2019-03-29, last modified: 2019-06-07
We report results from an extensive set of measurements of the \b{eta}-decay
response in liquid xenon.These measurements are derived from high-statistics
calibration data from injected sources of both $^{3}$H and $^{14}$C in the LUX
detector. The mean light-to-charge ratio is reported for 13 electric field
values ranging from 43 to 491 V/cm, and for energies ranging from 1.5 to 145
keV.
[27]
oai:arXiv.org:1906.00153 [pdf] - 1938443
Fallback accretion powered supernova light curves based on a
neutrino-driven explosion simulation of a 40 Msun star
Submitted: 2019-06-01
We present synthetic light curves of fallback-powered supernovae based on a
neutrino-driven explosion of a 40 Msun zero-metallicity star with significant
fallback accretion onto a black hole that was previously simulated by Chan et
al. (2018) until shock breakout. Here, we investigate the light curve
properties of the explosion after shock breakout for various fallback models.
Without extra power from fallback accretion, the light curve is that of a Type
IIP supernova with a plateau magnitude of around -14 mag and a plateau duration
of 40 days. With extra power for the light curve from fallback accretion,
however, we find that the transient can be significantly more luminous. The
light-curve shape can be SN 1987A-like or Type IIP-like, depending on the
efficiency of the fallback engine. If the accretion disk forms soon after the
collapse and more than 1% of the accretion energy is released as the central
engine, fallback accretion powered supernovae become as luminous as
superluminous supernovae. We suggest that Type II superluminous supernovae with
broad hydrogen features could be related to such hydrogen-rich supernovae
powered by fallback accretion. In the future, such superluminous supernovae
powered by fallback accretion might be found among the supernovae from the
first stars in addition to pair-instability supernovae and pulstational
pair-instability supernovae.
[28]
oai:arXiv.org:1806.06623 [pdf] - 1814930
An FFT-based Solution Method for the Poisson Equation on 3D Spherical
Polar Grids
Submitted: 2018-06-18, last modified: 2018-12-04
The solution of the Poisson equation is a ubiquitous problem in computational
astrophysics. Most notably, the treatment of self-gravitating flows involves
the Poisson equation for the gravitational field. In hydrodynamics codes using
spherical polar grids, one often resorts to a truncated spherical harmonics
expansion for an approximate solution. Here we present a non-iterative method
that is similar in spirit, but uses the full set of eigenfunctions of the
discretized Laplacian to obtain an exact solution of the discretized Poisson
equation. This allows the solver to handle density distributions for which the
truncated multipole expansion fails, such as off-center point masses. In three
dimensions, the operation count of the new method is competitive with a naive
implementation of the truncated spherical harmonics expansion with $N_\ell
\approx 15$ multipoles. We also discuss the parallel implementation of the
algorithm. The serial code and a template for the parallel solver are made
publicly available.
[29]
oai:arXiv.org:1811.05483 [pdf] - 1822945
Three-Dimensional Simulations of Neutrino-Driven Core-Collapse
Supernovae from Low-Mass Single and Binary Star Progenitors
Submitted: 2018-11-13
We present a suite of seven 3D supernova simulations of non-rotating low-mass
progenitors using multi-group neutrino transport. Our simulations cover single
star progenitors with zero-age main sequence masses between $9.6 M_\odot$ and
$12.5 M_\odot$ and (ultra)stripped-envelope progenitors with initial helium
core masses between $2.8 M_\odot$ and $3.5 M_\odot$. We find explosion energies
between $0.1\,\mathrm{Bethe}$ and $0.4\,\mathrm{Bethe}$, which are still rising
by the end of the simulations. Although less energetic than typical events, our
models are compatible with observations of less energetic explosions of
low-mass progenitors. In six of our models, the mass outflow rate already
exceeds the accretion rate onto the proto-neutron star, and the mass and
angular momentum of the compact remnant have closely approached their final
value, barring the possibility of later fallback. While the proto-neutron star
is still accelerated by the gravitational tug of the asymmetric ejecta, the
acceleration can be extrapolated to obtain estimates for the final kick
velocity. We obtain gravitational neutron star masses between $1.22 M_\odot$
and $1.44 M_\odot$, kick velocities between $11\, \mathrm{km}\,
\mathrm{s}^{-1}$ and $695\, \mathrm{km}\, \mathrm{s}^{-1}$, and spin periods
from $20\, \mathrm{ms}$ to $2.7\,\mathrm{s}$, which suggests that typical
neutron star birth properties can be naturally obtained in the neutrino-driven
paradigm. We find a loose correlation between the explosion energy and the kick
velocity. There is no indication of spin-kick alignment, but a correlation
between the kick velocity and the neutron star angular momentum, which needs to
be investigated further as a potential point of tension between models and
observations.
[30]
oai:arXiv.org:1807.07113 [pdf] - 1758152
Search for annual and diurnal rate modulations in the LUX experiment
Akerib, D. S.;
Alsum, S.;
Araújo, H. M.;
Bai, X.;
Balajthy, J.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Boxer, B.;
Brás, P.;
Burdin, S.;
Byram, D.;
Carmona-Benitez, M. C.;
Chan, C.;
Cutter, J. E.;
Davison, T. J. R.;
Druszkiewicz, E.;
Fallon, S. R.;
Fan, A.;
Fiorucci, S.;
Gaitskell, R. J.;
Genovesi, J.;
Ghag, C.;
Gilchriese, M. G. D.;
Gwilliam, C.;
Hall, C. R.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Jacobsen, R. G.;
Ji, W.;
Kamdin, K.;
Kazkaz, K.;
Khaitan, D.;
Knoche, R.;
Korolkova, E. V.;
Kravitz, S.;
Kudryavtsev, V. A.;
Lenardo, B. G.;
Lesko, K. T.;
Liao, J.;
Lin, J.;
Lindote, A.;
Lopes, M. I.;
Manalaysay, A.;
Mannino, R. L.;
Marangou, N.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. -M.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pease, E. K.;
Rischbieter, G. R. C.;
Rhyne, C.;
Rossiter, P.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, W. C.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Utku, U.;
Uvarov, S.;
Velan, V.;
Verbus, J. R.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Woodward, D.;
Xu, J.;
Yazdani, K.;
Zhang, C.
Submitted: 2018-07-18, last modified: 2018-09-27
Various dark matter models predict annual and diurnal modulations of dark
matter interaction rates in Earth-based experiments as a result of the Earth's
motion in the halo. Observation of such features can provide generic evidence
for detection of dark matter interactions. This paper reports a search for both
annual and diurnal rate modulations in the LUX dark matter experiment using
over 20 calendar months of data acquired between 2013 and 2016. This search
focuses on electron recoil events at low energies, where leptophilic dark
matter interactions are expected to occur and where the DAMA experiment has
observed a strong rate modulation for over two decades. By using the innermost
volume of the LUX detector and developing robust cuts and corrections, we
obtained a stable event rate of 2.3$\pm$0.2~cpd/keV$_{\text{ee}}$/tonne, which
is among the lowest in all dark matter experiments. No statistically
significant annual modulation was observed in energy windows up to
26~keV$_{\text{ee}}$. Between 2 and 6~keV$_{\text{ee}}$, this analysis
demonstrates the most sensitive annual modulation search up to date, with
9.2$\sigma$ tension with the DAMA/LIBRA result. We also report no observation
of diurnal modulations above 0.2~cpd/keV$_{\text{ee}}$/tonne amplitude between
2 and 6~keV$_{\text{ee}}$.
[31]
oai:arXiv.org:1601.06799 [pdf] - 1732514
GRMHD simulations of visibility amplitude variability for Event Horizon
Telescope images of Sgr A*
Submitted: 2016-01-25, last modified: 2018-08-06
Synthesis imaging of the black hole in the center of the Milky Way, Sgr A*,
with the Event Horizon Telescope (EHT) rests on the assumption of a stationary
image. We explore the limitations of this assumption using high-cadence GRMHD
simulations of Sgr A*. We employ analytic models that capture the basic
characteristics of the images to understand the origin of the variability in
the simulated visibility amplitudes. We find that, in all simulations, the
visibility amplitudes for baselines oriented perpendicular to the spin axis of
the black hole typically decrease smoothly over baseline lengths that are
comparable to those of the EHT. On the other hand, the visibility amplitudes
for baselines oriented parallel to the spin axis show significant structure
with one or more minima. This suggests that fitting EHT observations with
geometric models will lead to reasonably accurate determination of the
orientation of the black-hole on the plane of the sky. However, in the
disk-dominated models, the locations and depths of the minima in the visibility
amplitudes depend primarily on the width and asymmetry of the crescent-like
images and are highly variable. In the jet-dominated models, the locations of
the minima are determined by the separation of the two image components but
their depths depend primarily on the relative brightness of the two components
and are also variable. This suggests that using time-independent models to
infer additional black-hole parameters, such as the shadow size or the spin
magnitude, will be severely affected by the variability of the accretion flow.
[32]
oai:arXiv.org:1702.03930 [pdf] - 1767316
Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code
for astrophysics
Price, Daniel J.;
Wurster, James;
Tricco, Terrence S.;
Nixon, Chris;
Toupin, Stéven;
Pettitt, Alex;
Chan, Conrad;
Mentiplay, Daniel;
Laibe, Guillaume;
Glover, Simon;
Dobbs, Clare;
Nealon, Rebecca;
Liptai, David;
Worpel, Hauke;
Bonnerot, Clément;
Dipierro, Giovanni;
Ballabio, Giulia;
Ragusa, Enrico;
Federrath, Christoph;
Iaconi, Roberto;
Reichardt, Thomas;
Forgan, Duncan;
Hutchison, Mark;
Constantino, Thomas;
Ayliffe, Ben;
Hirsh, Kieran;
Lodato, Giuseppe
Submitted: 2017-02-13, last modified: 2018-06-15
We present Phantom, a fast, parallel, modular and low-memory smoothed
particle hydrodynamics and magnetohydrodynamics code developed over the last
decade for astrophysical applications in three dimensions. The code has been
developed with a focus on stellar, galactic, planetary and high energy
astrophysics and has already been used widely for studies of accretion discs
and turbulence, from the birth of planets to how black holes accrete. Here we
describe and test the core algorithms as well as modules for
magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas
mixtures, physical viscosity, external forces including numerous galactic
potentials as well as implementations of Lense-Thirring precession,
Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby
made publicly available.
[33]
oai:arXiv.org:1712.01882 [pdf] - 2036220
Magnetorotational instability in eccentric disks
Submitted: 2017-12-05, last modified: 2018-03-08
Eccentric disks arise in such astrophysical contexts as tidal disruption
events, but it is unknown whether the magnetorotational instability (MRI),
which powers accretion in circular disks, operates in eccentric disks as well.
We examine the linear evolution of unstratified, incompressible MRI in an
eccentric disk orbiting a point mass. We consider vertical modes of wavenumber
$k$ on a background flow with uniform eccentricity $e$ and vertical Alfv\'en
speed $v_\mathrm A$ along an orbit with mean motion $n$. We find two mode
families, one with dominant magnetic components, the other with dominant
velocity components; the former is unstable at $(1-e)^3f^2\lesssim3$, where
$f\equiv kv_\mathrm A/n$, the latter at $e\gtrsim0.8$. For $f^2\lesssim3$, MRI
behaves much like in circular disks, but the growth per orbit declines slowly
with increasing $e$; for $f^2\gtrsim3$, modes grow by parametric amplification,
which is resonant for $0<e\ll1$. MRI growth and the attendant angular momentum
and energy transport happen chiefly near pericenter, where orbital shear
dominates magnetic tension.
[34]
oai:arXiv.org:1710.00838 [pdf] - 1614836
Black hole formation and fallback during the supernova explosion of a
$40 \,\mathrm{M}_\odot$ star
Submitted: 2017-10-02
Fallback in core-collapse supernovae is considered a major ingredient for
explaining abundance anomalies in metal-poor stars and the natal kicks and
spins of black holes (BHs). We present a first 3D simulation of BH formation
and fallback in an "aborted" neutrino-driven explosion of a $40$ solar mass
zero-metallicity progenitor from collapse to shock breakout. We follow the
phase up to BH formation using the relativistic CoCoNuT-FMT code. For the
subsequent evolution to shock breakout we apply the moving-mesh code Arepo to
core-collapse supernovae for the first time. Our simulation shows that despite
early BH formation, neutrino-heated bubbles can survive for tens of seconds
before being accreted, leaving them sufficient time to transfer part of their
energy to sustain the shock wave as is propagates through the envelope.
Although the initial net energy ($\sim 2$ Bethe) of the neutrino-heated ejecta
barely equals the binding energy of the envelope, $11\,\mathrm{M}_\odot$ of
hydrogen are still expelled with an energy of $0.23$ Bethe. We find no
significant mixing and only a modest BH kick and spin, but speculate that
stronger effects could occur for slightly more energetic explosions or
progenitors with less tightly bound envelopes.
[35]
oai:arXiv.org:1709.00800 [pdf] - 1611913
Ultra-Low Energy Calibration of LUX Detector using $^{127}$Xe Electron
Capture
LUX Collaboration;
Akerib, D. S.;
Alsum, S.;
Araújo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Brás, P.;
Byram, D.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Currie, A.;
Cutter, J. E.;
Davison, T. J. R.;
Dobi, A.;
Druszkiewicz, E.;
Edwards, B. N.;
Fallon, S. R.;
Fan, A.;
Fiorucci, S.;
Gaitskell, R. J.;
Genovesi, J.;
Ghag, C.;
Gilchriese, M. G. D.;
Hall, C. R.;
Hanhardt, M.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Jacobsen, R. G.;
Ji, W.;
Kamdin, K.;
Kazkaz, K.;
Khaitan, D.;
Knoche, R.;
Larsen, N. A.;
Lenardo, B. G.;
Lesko, K. T.;
Lindote, A.;
Lopes, M. I.;
Manalaysay, A.;
Mannino, R. L.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. M.;
Mock, J.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
O'Sullivan, K.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pease, E. K.;
Rhyne, C.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, W. C.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Uvarov, S.;
Velan, V.;
Verbus, J. R.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Xu, J.;
Yazdani, K.;
Young, S. K.;
Zhang, C.
Submitted: 2017-09-03
We report an absolute calibration of the ionization yields($\textit{Q$_y$})$
and fluctuations for electronic recoil events in liquid xenon at discrete
energies between 186 eV and 33.2 keV. The average electric field applied across
the liquid xenon target is 180 V/cm. The data are obtained using low energy
$^{127}$Xe electron capture decay events from the 95.0-day first run from LUX
(WS2013) in search of Weakly Interacting Massive Particles (WIMPs). The
sequence of gamma-ray and X-ray cascades associated with $^{127}$I
de-excitations produces clearly identified 2-vertex events in the LUX detector.
We observe the K- (binding energy, 33.2 keV), L- (5.2 keV), M- (1.1 keV), and
N- (186 eV) shell cascade events and verify that the relative ratio of observed
events for each shell agrees with calculations. The N-shell cascade analysis
includes single extracted electron (SE) events and represents the lowest-energy
electronic recoil $\textit{in situ}$ measurements that have been explored in
liquid xenon.
[36]
oai:arXiv.org:1706.00240 [pdf] - 1584139
Geometrically thick obscuration by radiation-driven outflow from
magnetized tori of active galactic nuclei
Submitted: 2017-06-01, last modified: 2017-08-18
Near-Eddington radiation from active galactic nuclei (AGNs) has significant
dynamical influence on the surrounding dusty gas, plausibly furnishing AGNs
with geometrically thick obscuration. We investigate this paradigm with
radiative magnetohydrodynamics simulations. The simulations solve the
magnetohydrodynamics equations simultaneously with the infrared (IR) and
ultraviolet (UV) radiative transfer (RT) equations; no approximate closure is
used for RT. We find that our torus, when given a suitable sub-Keplerian
angular momentum profile, spontaneously evolves toward a state in which its
opening angle, density distribution, and flow pattern change only slowly. This
"steady" state lasts for as long as there is gas resupply toward the inner
edge. The torus is best described as a mid-plane inflow and a high-latitude
outflow. The outflow is launched from the torus inner edge by UV radiation and
expands in solid angle as it ascends; IR radiation continues to drive the
wide-angle outflow outside the central hole. The dusty outflow obscures the
central source in soft X-rays, the IR, and the UV over three quarters of solid
angle, and each decade in column density covers roughly equal solid angle
around the central source; these obscuration properties are similar to what
observations imply.
[37]
oai:arXiv.org:1704.02297 [pdf] - 1582421
First Searches for Axions and Axion-Like Particles with the LUX
Experiment
Akerib, D. S.;
Alsum, S.;
Aquino, C.;
Araújo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Brás, P.;
Byram, D.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Chiller, A. A.;
Chiller, C.;
Currie, A.;
Cutter, J. E.;
Davison, T. J. R.;
Dobi, A.;
Dobson, J. E. Y.;
Druszkiewicz, E.;
Edwards, B. N.;
Faham, C. H.;
Fallon, S. R.;
Fiorucci, S.;
Gaitskell, R. J.;
Gehman, V. M.;
Ghag, C.;
Gibson, K. R.;
Gilchriese, M. G. D.;
Hall, C. R.;
Hanhardt, M.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Jacobsen, R. G.;
Ji, W.;
Kamdin, K.;
Kazkaz, K.;
Khaitan, D.;
Knoche, R.;
Larsen, N. A.;
Lee, C.;
Lenardo, B. G.;
Lesko, K. T.;
Lindote, A.;
Lopes, M. I.;
Manalaysay, A.;
Mannino, R. L.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. M.;
Mock, J.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
O'Sullivan, K.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pease, E. K.;
Reichhart, L.;
Rhyne, C.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Stephenson, S.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, W. C.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Uvarov, S.;
Velan, V.;
Verbus, J. R.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Xu, J.;
Yazdani, K.;
Young, S. K.;
Zhang, C.
Submitted: 2017-04-07, last modified: 2017-06-30
The first searches for axions and axion-like particles with the Large
Underground Xenon (LUX) experiment are presented. Under the assumption of an
axio-electric interaction in xenon, the coupling constant between axions and
electrons, gAe is tested, using data collected in 2013 with an exposure
totalling 95 live-days $\times$ 118 kg. A double-sided, profile likelihood
ratio statistic test excludes gAe larger than 3.5 $\times$ 10$^{-12}$ (90%
C.L.) for solar axions. Assuming the DFSZ theoretical description, the upper
limit in coupling corresponds to an upper limit on axion mass of 0.12
eV/c$^{2}$, while for the KSVZ description masses above 36.6 eV/c$^{2}$ are
excluded. For galactic axion-like particles, values of gAe larger than 4.2
$\times$ 10$^{-13}$ are excluded for particle masses in the range 1-16
keV/c$^{2}$. These are the most stringent constraints to date for these
interactions.
[38]
oai:arXiv.org:1705.03380 [pdf] - 1583256
Limits on spin-dependent WIMP-nucleon cross section obtained from the
complete LUX exposure
LUX Collaboration;
Akerib, D. S.;
Alsum, S.;
Araújo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Brás, P.;
Byram, D.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Chiller, A. A.;
Chiller, C.;
Currie, A.;
Cutter, J. E.;
Davison, T. J. R.;
Dobi, A.;
Dobson, J. E. Y.;
Druszkiewicz, E.;
Edwards, B. N.;
Faham, C. H.;
Fallon, S. R.;
Fiorucci, S.;
Gaitskell, R. J.;
Gehman, V. M.;
Ghag, C.;
Gilchriese, M. G. D.;
Hall, C. R.;
Hanhardt, M.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Jacobsen, R. G.;
Ji, W.;
Kamdin, K.;
Kazkaz, K.;
Khaitan, D.;
Knoche, R.;
Larsen, N. A.;
Lee, C.;
Lenardo, B. G.;
Lesko, K. T.;
Lindote, A.;
Lopes, M. I.;
Manalaysay, A.;
Mannino, R. L.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. M.;
Mock, J.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
O'Sullivan, K.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pease, E. K.;
Reichhart, L.;
Rhyne, C.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Stephenson, S.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, W. C.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Uvarov, S.;
Velan, V.;
Verbus, J. R.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Xu, J.;
Yazdani, K.;
Young, S. K.;
Zhang, C.
Submitted: 2017-05-09, last modified: 2017-06-23
We present experimental constraints on the spin-dependent WIMP-nucleon
elastic cross sections from the total 129.5 kg-year exposure acquired by the
Large Underground Xenon experiment (LUX), operating at the Sanford Underground
Research Facility in Lead, South Dakota (USA). A profile likelihood ratio
analysis allows 90% CL upper limits to be set on the WIMP-neutron (WIMP-proton)
cross section of $\sigma_n$ = 1.6$\times 10^{-41}$ cm$^{2}$ ($\sigma_p$ =
5$\times 10^{-40}$ cm$^{2}$) at 35 GeV$c^{-2}$, almost a sixfold improvement
over the previous LUX spin-dependent results. The spin-dependent WIMP-neutron
limit is the most sensitive constraint to date.
[39]
oai:arXiv.org:1706.07062 [pdf] - 1779469
GRay2: A General Purpose Geodesic Integrator for Kerr Spacetimes
Submitted: 2017-06-21
Fast and accurate integration of geodesics in Kerr spacetimes is an important
tool in modeling the orbits of stars and the transport of radiation in the
vicinities of black holes. Most existing integration algorithms employ
Boyer-Lindquist coordinates, which have coordinate singularities at the event
horizon and along the poles. Handling the singularities requires special
numerical treatment in these regions, often slows down the calculations, and
may lead to inaccurate geodesics. We present here a new general-purpose
geodesic integrator, GRay2, that overcomes these issues by employing the
Cartesian form of Kerr-Schild coordinates. By performing particular
mathematical manipulations of the geodesic equations and several optimizations,
we develop an implementation of the Cartesian Kerr-Schild coordinates that
outperforms calculations that use the seemingly simpler equations in
Boyer-Lindquist coordinates. We also employ the OpenCL framework, which allows
GRay2 to run on multi-core CPUs as well as on a wide range of GPU hardware
accelerators, making the algorithm more versatile. We report numerous
convergence tests and benchmark results for GRay2 for both time-like (particle)
and null (photon) geodesics.
[40]
oai:arXiv.org:1705.06293 [pdf] - 1637475
The Properties of Reconnection Current Sheets in GRMHD Simulations of
Radiatively Inefficient Accretion Flows
Submitted: 2017-05-17
Non-ideal MHD effects may play a significant role in determining the
dynamics, thermal properties, and observational signatures of radiatively
inefficient accretion flows onto black holes. In particular, particle
acceleration during magnetic reconnection events may influence black hole
spectra and flaring properties. We use representative GRMHD simulations of
black hole accretion flows to identify and explore the structures and
properties of current sheets as potential sites of magnetic reconnection. In
the case of standard and normal (SANE) disks, we find that, in the reconnection
sites, the plasma beta ranges from $0.1$ to $1000$, the magnetization ranges
from $10^{-4}$ to $1$, and the guide fields are weak compared to the
reconnecting fields. In magnetically arrested (MAD) disks, we find typical
values for plasma beta from $10^{-2}$ to $10^3$, magnetizations from $10^{-3}$
to $10$, and typically stronger guide fields, with strengths comparable to or
greater than the reconnecting fields. These are critical parameters that govern
the electron energy distribution resulting from magnetic reconnection and can
be used in the context of plasma simulations to provide microphysics inputs to
global simulations. We also find that ample magnetic energy is available in the
reconnection regions to power the fluence of bright X-ray flares observed from
the black hole in the center of the Milky Way.
[41]
oai:arXiv.org:1511.03428 [pdf] - 1563900
The Mass Distribution of Population III Stars
Submitted: 2015-11-11, last modified: 2017-04-03
Extremely metal-poor stars are uniquely informative on the nature of massive
Population III stars. Modulo a few elements that vary with stellar evolution,
the present-day photospheric abundances observed in extremely metal-poor stars
are representative of their natal gas cloud composition. For this reason, the
chemistry of extremely metal-poor stars closely reflects the nucleosynthetic
yields of supernovae from massive Population III stars. Here we collate
detailed abundances of 53 extremely metal-poor stars from the literature and
infer the masses of their Population III progenitors. We fit a simple initial
mass function to a subset of 29 of theinferred Population III star masses, and
find that the mass distribution is well-represented by a power law IMF with
exponent $\alpha = 2.35^{+0.29}_{-0.24}$. The inferred maximum progenitor mass
for supernovae from massive Population III stars is $M_{\rm{max}} =
87^{+13}_{-33}$ M$_\odot$, and we find no evidence in our sample for a
contribution from stars with masses above $\sim$120 M$_\odot$. The minimum mass
is strongly consistent with the theoretical lower mass limit for Population III
supernovae. We conclude that the IMF for massive Population III stars is
consistent with the initial mass function of present-day massive stars and
there may well have formed stars much below the supernova mass limit that could
have survived to the present day.
[43]
oai:arXiv.org:1608.07648 [pdf] - 1531370
Results from a search for dark matter in the complete LUX exposure
Akerib, D. S.;
Alsum, S.;
Araújo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Bramante, R.;
Brás, P.;
Byram, D.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Chiller, A. A.;
Chiller, C.;
Currie, A.;
Cutter, J. E.;
Davison, T. J. R.;
Dobi, A.;
Dobson, J. E. Y.;
Druszkiewicz, E.;
Edwards, B. N.;
Faham, C. H.;
Fiorucci, S.;
Gaitskell, R. J.;
Gehman, V. M.;
Ghag, C.;
Gibson, K. R.;
Gilchriese, M. G. D.;
Hall, C. R.;
Hanhardt, M.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Ihm, M.;
Jacobsen, R. G.;
Ji, W.;
Kamdin, K.;
Kazkaz, K.;
Khaitan, D.;
Knoche, R.;
Larsen, N. A.;
Lee, C.;
Lenardo, B. G.;
Lesko, K. T.;
Lindote, A.;
Lopes, M. I.;
Manalaysay, A.;
Mannino, R. L.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. M.;
Mock, J.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
O`Sullivan, K.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pease, E. K.;
Phelps, P.;
Reichhart, L.;
Rhyne, C.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Stephenson, S.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, W. C.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Uvarov, S.;
Verbus, J. R.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Xu, J.;
Yazdani, K.;
Young, S. K.;
Zhang, C.
Submitted: 2016-08-26, last modified: 2017-01-13
We report constraints on spin-independent weakly interacting massive particle
(WIMP)-nucleon scattering using a 3.35e4 kg-day exposure of the Large
Underground Xenon (LUX) experiment. A dual-phase xenon time projection chamber
with 250 kg of active mass is operated at the Sanford Underground Research
Facility under Lead, South Dakota (USA). With roughly fourfold improvement in
sensitivity for high WIMP masses relative to our previous results, this search
yields no evidence of WIMP nuclear recoils. At a WIMP mass of 50 GeV/c^2,
WIMP-nucleon spin-independent cross sections above 2.2e-46 cm^2 are excluded at
the 90% confidence level. When combined with the previously reported LUX
exposure, this exclusion strengthens to 1.1e-46 cm^2 at 50 GeV/c^2.
[44]
oai:arXiv.org:1602.00692 [pdf] - 1530453
Bayesian techniques for comparing time-dependent GRMHD simulations to
variable Event Horizon Telescope observations
Submitted: 2016-02-01, last modified: 2016-12-16
The Event Horizon Telescope (EHT) is a millimeter-wavelength,
very-long-baseline interferometry (VLBI) experiment that is capable of
observing black holes with horizon-scale resolution. Early observations have
revealed variable horizon-scale emission in the Galactic Center black hole,
Sagittarius A* (Sgr A*). Comparing such observations to time-dependent general
relativistic magnetohydrodynamic (GRMHD) simulations requires statistical tools
that explicitly consider the variability in both the data and the models. We
develop here a Bayesian method to compare time-resolved simulation images to
variable VLBI data, in order to infer model parameters and perform model
comparisons. We use mock EHT data based on GRMHD simulations to explore the
robustness of this Bayesian method and contrast it to approaches that do not
consider the effects of variability. We find that time-independent models lead
to offset values of the inferred parameters with artificially reduced
uncertainties. Moreover, neglecting the variability in the data and the models
often leads to erroneous model selections. We finally apply our method to the
early EHT data on Sgr A*.
[45]
oai:arXiv.org:1608.05381 [pdf] - 1503451
Low-energy (0.7-74 keV) nuclear recoil calibration of the LUX dark
matter experiment using D-D neutron scattering kinematics
LUX Collaboration;
Akerib, D. S.;
Alsum, S.;
Araújo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Bradley, A.;
Bramante, R.;
Brás, P.;
Byram, D.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Chapman, J. J.;
Chiller, A. A.;
Chiller, C.;
Currie, A.;
Cutter, J. E.;
Davison, T. J. R.;
de Viveiros, L.;
Dobi, A.;
Dobson, J. E. Y.;
Druszkiewicz, E.;
Edwards, B. N.;
Faham, C. H.;
Fiorucci, S.;
Gaitskell, R. J.;
Gehman, V. M.;
Ghag, C.;
Gibson, K. R.;
Gilchriese, M. G. D.;
Hall, C. R.;
Hanhardt, M.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Ihm, M.;
Jacobsen, R. G.;
Ji, W.;
Kamdin, K.;
Kazkaz, K.;
Khaitan, D.;
Knoche, R.;
Larsen, N. A.;
Lee, C.;
Lenardo, B. G.;
Lesko, K. T.;
Lindote, A.;
Lopes, M. I.;
Malling, D. C.;
Manalaysay, A.;
Mannino, R. L.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. M.;
Mock, J.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
O'Sullivan, K.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pangilinan, M.;
Pease, E. K.;
Phelps, P.;
Reichhart, L.;
Rhyne, C. A.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solmaz, M.;
Solovov, V. N.;
Sorensen, P.;
Stephenson, S.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, W. C.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Uvarov, S.;
Verbus, J. R.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Xu, J.;
Yazdani, K.;
Young, S. K.;
Zhang, C.
Submitted: 2016-08-18, last modified: 2016-10-26
The Large Underground Xenon (LUX) experiment is a dual-phase liquid xenon
time projection chamber (TPC) operating at the Sanford Underground Research
Facility in Lead, South Dakota. A calibration of nuclear recoils in liquid
xenon was performed $\textit{in situ}$ in the LUX detector using a collimated
beam of mono-energetic 2.45 MeV neutrons produced by a deuterium-deuterium
(D-D) fusion source. The nuclear recoil energy from the first neutron scatter
in the TPC was reconstructed using the measured scattering angle defined by
double-scatter neutron events within the active xenon volume. We measured the
absolute charge ($Q_{y}$) and light ($L_{y}$) yields at an average electric
field of 180 V/cm for nuclear recoil energies spanning 0.7 to 74 keV and 1.1 to
74 keV, respectively. This calibration of the nuclear recoil signal yields will
permit the further refinement of liquid xenon nuclear recoil signal models and,
importantly for dark matter searches, clearly demonstrates measured ionization
and scintillation signals in this medium at recoil energies down to
$\mathcal{O}$(1 keV).
[46]
oai:arXiv.org:1610.06339 [pdf] - 1553908
Combined Nucleosynthetic Yields of Multiple First Stars
Submitted: 2016-10-20
Modern numerical simulations of the formation of the first stars predict that
the first stars formed in multiples. In those cases, the chemical yields of
multiple supernova explosions may have contributed to the formation of a next
generation star. We match the chemical abundances of the oldest observed stars
in the universe to a database of theoretical supernova models, to show that it
is likely that the first stars formed from the ashes of two or more
progenitors.
[47]
oai:arXiv.org:1610.03505 [pdf] - 1580386
Variability in GRMHD simulations of Sgr A$^*$: Implications for EHT
closure phase observations
Submitted: 2016-10-11
The observable quantities that carry the most information regarding the
structures of the images of black holes in the interferometric observations
with the Event Horizon Telescope are the closure phases along different
baseline triangles. We use long time span, high cadence, GRMHD+radiative
transfer models of Sgr A$^*$ to investigate the expected variability of closure
phases in such observations. We find that, in general, closure phases along
small baseline triangles show little variability, except in the cases when one
of the triangle vertices crosses one of a small regions of low visibility
amplitude. The closure phase variability increases with the size of the
baseline triangle, as larger baselines probe the small-scale structures of the
images, which are highly variable. On average, the jet-dominated MAD models
show less closure phase variability than the disk-dominated SANE models, even
in the large baseline triangles, because the images from the latter are more
sensitive to the turbulence in the accretion flow. Our results suggest that
image reconstruction techniques need to explicitly take into account the
closure phase variability, especially if the quality and quantity of data allow
for a detailed characterization of the nature of variability. This also implies
that, if image reconstruction techniques that rely on the assumption of a
static image are utilized, regions of the $u-v$ space that show a high level of
variability will need to be identified and excised.
[48]
oai:arXiv.org:1609.02134 [pdf] - 1531508
Observational Constraints on First-Star Nucleosynthesis. II.
Spectroscopy of an Ultra Metal-Poor CEMP-no Star
Submitted: 2016-09-07
We report on the first high-resolution spectroscopic analysis of HE0020-1741,
a bright (V=12.9), ultra metal-poor ([Fe/H] = -4.1), carbon-enhanced ([C/Fe] =
+1.7) star selected from the Hamburg/ESO Survey. This star exhibits low
abundances of neutron-capture elements ([Ba/Fe] = -1.1), and an absolute carbon
abundance A(C) = 6.1; based on either criterion, HE0020-1741 is sub-classified
as a CEMP-no star. We show that the light-element abundance pattern of
HE0020-1741 is consistent with predicted yields from a massive (M = 21.5 Mo),
primordial composition, supernova (SN) progenitor. We also compare the
abundance patterns of other ultra metal-poor stars from the literature with
available measures of C, N, Na, Mg, and Fe abundances with an extensive grid of
SN models (covering the mass range 10 Mo - 100 Mo), in order to probe the
nature of their likely stellar progenitors. Our results suggest that at least
two classes of progenitors are required at [Fe/H] < -4.0, as the abundance
patterns for more than half of the sample studied in this work (7 out of 12
stars) cannot be easily reproduced by the predicted yields.
[49]
oai:arXiv.org:1602.05968 [pdf] - 1462494
Particle Acceleration and the Origin of X-ray Flares in GRMHD
simulations of Sgr A*
Submitted: 2016-02-18, last modified: 2016-08-23
Significant X-ray variability and flaring has been observed from Sgr A* but
is poorly understood from a theoretical standpoint. We perform GRMHD
simulations that take into account a population of non-thermal electrons with
energy distributions and injection rates that are motivated by PIC simulations
of magnetic reconnection. We explore the effects of including these non-thermal
electrons on the predicted broadband variability of Sgr A* and find that X-ray
variability is a generic result of localizing non-thermal electrons to highly
magnetized regions, where particles are likely to be accelerated via magnetic
reconnection. The proximity of these high-field regions to the event horizon
forms a natural connection between IR and X-ray variability and accounts for
the rapid timescales associated with the X-ray flares. The qualitative nature
of this variability is consistent with observations, producing X-ray flares
that are always coincident with IR flares, but not vice versa, i.e., there are
a number of IR flares without X-ray counterparts.
[50]
oai:arXiv.org:1511.00333 [pdf] - 1431612
Radiation-driven outflows from and radiative support in dusty tori of
active galactic nuclei
Submitted: 2015-11-01, last modified: 2016-07-01
Substantial evidence points to dusty, geometrically thick tori obscuring the
central engines of active galactic nuclei (AGNs), but so far no mechanism
satisfactorily explains why cool dust in the torus remains in a puffy geometry.
Near-Eddington infrared (IR) and ultraviolet (UV) luminosities coupled with
high dust opacities at these frequencies suggest that radiation pressure on
dust can play a significant role in shaping the torus. To explore the possible
effects of radiation pressure, we perform three-dimensional radiative
hydrodynamics simulations of an initially smooth torus. Our code solves the
hydrodynamics equations, the time-dependent multi-angle group IR radiative
transfer (RT) equation, and the time-independent UV RT equation. We find a
highly dynamic situation. IR radiation is anisotropic, leaving primarily
through the central hole. The torus inner surface exhibits a break in
axisymmetry under the influence of radiation and differential rotation;
clumping follows. In addition, UV radiation pressure on dust launches a strong
wind along the inner surface; when scaled to realistic AGN parameters, this
outflow travels at $\sim 5000 (M/10^7 M_\odot)^{1/4} [L_\mathrm{UV}/(0.1
L_\mathrm E)]^{1/4} \mathrm{km}\,\mathrm s^{-1}$ and carries $\sim 0.1 (M/10^7
M_\odot)^{3/4} [L_\mathrm{UV}/(0.1 L_\mathrm E)]^{3/4}
M_\odot\,\mathrm{yr}^{-1}$, where $M$, $L_\mathrm{UV}$, and $L_\mathrm E$ are
the mass, UV luminosity, and Eddington luminosity of the central object
respectively.
[51]
oai:arXiv.org:1512.03506 [pdf] - 1406496
Improved Limits on Scattering of Weakly Interacting Massive Particles
from Reanalysis of 2013 LUX data
LUX Collaboration;
Akerib, D. S.;
Araújo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Bradley, A.;
Bramante, R.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Chapman, J. J.;
Chiller, A. A.;
Chiller, C.;
Currie, A.;
Cutter, J. E.;
Davison, T. J. R.;
de Viveiros, L.;
Dobi, A.;
Dobson, J. E. Y.;
Druszkiewicz, E.;
Edwards, B. N.;
Faham, C. H.;
Fiorucci, S.;
Gaitskell, R. J.;
Gehman, V. M.;
Ghag, C.;
Gibson, K. R.;
Gilchriese, M. G. D.;
Hall, C. R.;
Hanhardt, M.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Ihm, M.;
Jacobsen, R. G.;
Ji, W.;
Kazkaz, K.;
Khaitan, D.;
Knoche, R.;
Larsen, N. A.;
Lee, C.;
Lenardo, B. G.;
Lesko, K. T.;
Lindote, A.;
Lopes, M. I.;
Malling, D. C.;
Manalaysay, A.;
Mannino, R. L.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. M.;
Mock, J.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
O`Sullivan, K.;
Oliver-Mallory, K. C.;
Ott, R. A.;
Palladino, K. J.;
Pangilinan, M.;
Pease, E. K.;
Phelps, P.;
Reichhart, L.;
Rhyne, C.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solovov, V. N.;
Sorensen, P.;
Stephenson, S.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, W.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Uvarov, S.;
Verbus, J. R.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Yazdani, K.;
Young, S. K.;
Zhang, C.
Submitted: 2015-12-10, last modified: 2016-05-16
We present constraints on weakly interacting massive particles (WIMP)-nucleus
scattering from the 2013 data of the Large Underground Xenon dark matter
experiment, including $1.4\times10^{4}\;\mathrm{kg\; day}$ of search exposure.
This new analysis incorporates several advances: single-photon calibration at
the scintillation wavelength, improved event-reconstruction algorithms, a
revised background model including events originating on the detector walls in
an enlarged fiducial volume, and new calibrations from decays of an injected
tritium $\beta$ source and from kinematically constrained nuclear recoils down
to 1.1 keV. Sensitivity, especially to low-mass WIMPs, is enhanced compared to
our previous results which modeled the signal only above a 3 keV minimum
energy. Under standard dark matter halo assumptions and in the mass range above
4 $\mathrm{GeV}\,c^{-2}$, these new results give the most stringent direct
limits on the spin-independent WIMP-nucleon cross section. The 90% C.L. upper
limit has a minimum of 0.6 zb at 33 $\mathrm{GeV}\,c^{-2}$ WIMP mass.
[52]
oai:arXiv.org:1512.03133 [pdf] - 1402511
Tritium calibration of the LUX dark matter experiment
LUX Collaboration;
Akerib, D. S.;
Araújo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boulton, E. M.;
Bradley, A.;
Bramante, R.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Chapman, J. J.;
Chiller, A. A.;
Chiller, C.;
Currie, A.;
Cutter, J. E.;
Davison, T. J. R.;
de Viveiros, L.;
Dobi, A.;
Dobson, J. E. Y.;
Druszkiewicz, E.;
Edwards, B. N.;
Faham, C. H.;
Fiorucci, S.;
Gaitskell, R. J.;
Gehman, V. M.;
Ghag, C.;
Gibson, K. R.;
Gilchriese, M. G. D.;
Hall, C. R.;
Hanhardt, M.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hogan, D. P.;
Horn, M.;
Huang, D. Q.;
Ignarra, C. M.;
Ihm, M.;
Jacobsen, R. G.;
Ji, W.;
Kazkaz, K.;
Khaitan, D.;
Knoche, R.;
Larsen, N. A.;
Lee, C.;
Lenardo, B. G.;
Lesko, K. T.;
Lindote, A.;
Lopes, M. I.;
Malling, D. C.;
Manalaysay, A. G.;
Mannino, R. L.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. M.;
Mock, J.;
Moongweluwan, M.;
Morad, J. A.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H. N.;
Neves, F.;
O`Sullivan, K.;
Oliver-Mallory, K. C.;
Ott, R. A.;
Palladino, K. J.;
Pangilinan, M.;
Pease, E. K.;
Phelps, P.;
Reichhart, L.;
Rhyne, C.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Solovov, V. N.;
Sorensen, P.;
Stephenson, S.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D. J.;
Taylor, W.;
Tennyson, B. P.;
Terman, P. A.;
Tiedt, D. R.;
To, W. H.;
Tripathi, M.;
Tvrznikova, L.;
Uvarov, S.;
Verbus, J. R.;
Webb, R. C.;
White, J. T.;
Whitis, T. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Young, S. K.;
Zhang, C.
Submitted: 2015-12-09, last modified: 2016-05-05
We present measurements of the electron-recoil (ER) response of the LUX dark
matter detector based upon 170,000 highly pure and spatially-uniform tritium
decays. We reconstruct the tritium energy spectrum using the combined energy
model and find good agreement with expectations. We report the average charge
and light yields of ER events in liquid xenon at 180 V/cm and 105 V/cm and
compare the results to the NEST model. We also measure the mean charge
recombination fraction and its fluctuations, and we investigate the location
and width of the LUX ER band. These results provide input to a re-analysis of
the LUX Run3 WIMP search.
[53]
oai:arXiv.org:1602.05527 [pdf] - 1385498
Persistent Asymmetric Structure of Sagittarius A* on Event Horizon
Scales
Fish, Vincent L.;
Johnson, Michael D.;
Doeleman, Sheperd S.;
Broderick, Avery E.;
Psaltis, Dimitrios;
Lu, Ru-Sen;
Akiyama, Kazunori;
Alef, Walter;
Algaba, Juan Carlos;
Asada, Keiichi;
Beaudoin, Christopher;
Bertarini, Alessandra;
Blackburn, Lindy;
Blundell, Ray;
Bower, Geoffrey C.;
Brinkerink, Christiaan;
Cappallo, Roger;
Chael, Andrew A.;
Chamberlin, Richard;
Chan, Chi-Kwan;
Crew, Geoffrey B.;
Dexter, Jason;
Dexter, Matt;
Dzib, Sergio A.;
Falcke, Heino;
Freund, Robert;
Friberg, Per;
Greer, Christopher H.;
Gurwell, Mark A.;
Ho, Paul T. P.;
Honma, Mareki;
Inoue, Makoto;
Johannsen, Tim;
Kim, Junhan;
Krichbaum, Thomas P.;
Lamb, James;
León-Tavares, Jonathan;
Loeb, Abraham;
Loinard, Laurent;
MacMahon, David;
Marrone, Daniel P.;
Moran, James M.;
Mościbrodzka, Monika;
Ortiz-León, Gisela N.;
Oyama, Tomoaki;
Özel, Feryal;
Plambeck, Richard L.;
Pradel, Nicolas;
Primiani, Rurik A.;
Rogers, Alan E. E.;
Rosenfeld, Katherine;
Rottmann, Helge;
Roy, Alan L.;
Ruszczyk, Chester;
Smythe, Daniel L.;
SooHoo, Jason;
Spilker, Justin;
Stone, Jordan;
Strittmatter, Peter;
Tilanus, Remo P. J.;
Titus, Michael;
Vertatschitsch, Laura;
Wagner, Jan;
Wardle, John F. C.;
Weintroub, Jonathan;
Woody, David;
Wright, Melvyn;
Yamaguchi, Paul;
Young, André;
Young, Ken H.;
Zensus, J. Anton;
Ziurys, Lucy M.
Submitted: 2016-02-17
The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing
target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm
emission from this source on angular scales comparable to that of the general
relativistic shadow. Previous EHT observations have used visibility amplitudes
to infer the morphology of the millimeter-wavelength emission. Potentially much
richer source information is contained in the phases. We report on 1.3 mm phase
information on Sgr A* obtained with the EHT on a total of 13 observing nights
over 4 years. Closure phases, the sum of visibility phases along a closed
triangle of interferometer baselines, are used because they are robust against
phase corruptions introduced by instrumentation and the rapidly variable
atmosphere. The median closure phase on a triangle including telescopes in
California, Hawaii, and Arizona is nonzero. This result conclusively
demonstrates that the millimeter emission is asymmetric on scales of a few
Schwarzschild radii and can be used to break 180-degree rotational ambiguities
inherent from amplitude data alone. The stability of the sign of the closure
phase over most observing nights indicates persistent asymmetry in the image of
Sgr A* that is not obscured by refraction due to interstellar electrons along
the line of sight.
[54]
oai:arXiv.org:1411.1454 [pdf] - 1319513
A General Relativistic Null Hypothesis Test with Event Horizon Telescope
Observations of the black-hole shadow in Sgr A*
Submitted: 2014-11-05, last modified: 2015-11-10
The half opening angle of a Kerr black-hole shadow is always equal to
(5+-0.2)GM/Dc^2, where M is the mass of the black hole and D is its distance
from the Earth. Therefore, measuring the size of a shadow and verifying whether
it is within this 4% range constitutes a null hypothesis test of General
Relativity. We show that the black hole in the center of the Milky Way, Sgr A*,
is the optimal target for performing this test with upcoming observations using
the Event Horizon Telescope. We use the results of optical/IR monitoring of
stellar orbits to show that the mass-to-distance ratio for Sgr A* is already
known to an accuracy of +-4%. We investigate our prior knowledge of the
properties of the scattering screen between Sgr A and the Earth, the effects of
which will need to be corrected for in order for the black-hole shadow to
appear sharp against the background emission. Finally, we explore an edge
detection scheme for interferometric data and a pattern matching algorithm
based on the Hough/Radon transform and demonstrate that the shadow of the black
hole at 1.3 mm can be localized, in principle, to within ~9%. All these results
suggest that our prior knowledge of the properties of the black hole, of
scattering broadening, and of the accretion flow can only limit this General
Relativistic null hypothesis test with Event Horizon Telescope observations of
Sgr A* to 10%.
[55]
oai:arXiv.org:1509.02910 [pdf] - 1286987
LUX-ZEPLIN (LZ) Conceptual Design Report
The LZ Collaboration;
Akerib, D. S.;
Akerlof, C. W.;
Akimov, D. Yu.;
Alsum, S. K.;
Araújo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Balashov, S.;
Barry, M. J.;
Bauer, P.;
Beltrame, P.;
Bernard, E. P.;
Bernstein, A.;
Biesiadzinski, T. P.;
Boast, K. E.;
Bolozdynya, A. I.;
Boulton, E. M.;
Bramante, R.;
Buckley, J. H.;
Bugaev, V. V.;
Bunker, R.;
Burdin, S.;
Busenitz, J. K.;
Carels, C.;
Carlsmith, D. L.;
Carlson, B.;
Carmona-Benitez, M. C.;
Cascella, M.;
Chan, C.;
Cherwinka, J. J.;
Chiller, A. A.;
Chiller, C.;
Craddock, W. W.;
Currie, A.;
Cutter, J. E.;
da Cunha, J. P.;
Dahl, C. E.;
Dasu, S.;
Davison, T. J. R.;
de Viveiros, L.;
Dobi, A.;
Dobson, J. E. Y.;
Druszkiewicz, E.;
Edberg, T. K.;
Edwards, B. N.;
Edwards, W. R.;
Elnimr, M. M.;
Emmet, W. T.;
Faham, C. H.;
Fiorucci, S.;
Ford, P.;
Francis, V. B.;
Fu, C.;
Gaitskell, R. J.;
Gantos, N. J.;
Gehman, V. M.;
Gerhard, R. M.;
Ghag, C.;
Gilchriese, M. G. D.;
Gomber, B.;
Hall, C. R.;
Harris, A.;
Haselschwardt, S. J.;
Hertel, S. A.;
Hoff, M. D.;
Holbrook, B.;
Holtom, E.;
Huang, D. Q.;
Hurteau, T. W.;
Ignarra, C. M.;
Jacobsen, R. G.;
Ji, W.;
Ji, X.;
Johnson, M.;
Ju, Y.;
Kamdin, K.;
Kazkaz, K.;
Khaitan, D.;
Khazov, A.;
Khromov, A. V.;
Konovalov, A. M.;
Korolkova, E. V.;
Kraus, H.;
Krebs, H. J.;
Kudryavtsev, V. A.;
Kumpan, A. V.;
Kyre, S.;
Larsen, N. A.;
Lee, C.;
Lenardo, B. G.;
Lesko, K. T.;
Liao, F. -T.;
Lin, J.;
Lindote, A.;
Lippincott, W. H.;
Liu, J.;
Liu, X.;
Lopes, M. I.;
Lorenzon, W.;
Luitz, S.;
Majewski, P.;
Malling, D. C.;
Manalaysay, A. G.;
Manenti, L.;
Mannino, R. L.;
Markley, D. J.;
Martin, T. J.;
Marzioni, M. F.;
McKinsey, D. N.;
Mei, D. -M.;
Meng, Y.;
Miller, E. H.;
Mock, J.;
Monzani, M. E.;
Morad, J. A.;
Murphy, A. St. J.;
Nelson, H. N.;
Neves, F.;
Nikkel, J. A.;
O'Neill, F. G.;
O'Dell, J.;
O'Sullivan, K.;
Olevitch, M. A.;
Oliver-Mallory, K. C.;
Palladino, K. J.;
Pangilinan, M.;
Patton, S. J.;
Pease, E. K.;
Piepke, A.;
Powell, S.;
Preece, R. M.;
Pushkin, K.;
Ratcliff, B. N.;
Reichenbacher, J.;
Reichhart, L.;
Rhyne, C.;
Rodrigues, J. P.;
Rose, H. J.;
Rosero, R.;
Saba, J. S.;
Sarychev, M.;
Schnee, R. W.;
Schubnell, M. S. G.;
Scovell, P. R.;
Shaw, S.;
Shutt, T. A.;
Silva, C.;
Skarpaas, K.;
Skulski, W.;
Solovov, V. N.;
Sorensen, P.;
Sosnovtsev, V. V.;
Stancu, I.;
Stark, M. R.;
Stephenson, S.;
Stiegler, T. M.;
Sumner, T. J.;
Sundarnath, K.;
Szydagis, M.;
Taylor, D. J.;
Taylor, W.;
Tennyson, B. P.;
Terman, P. A.;
Thomas, K. J.;
Thomson, J. A.;
Tiedt, D. R.;
To, W. H.;
Tomás, A.;
Tripathi, M.;
Tull, C. E.;
Tvrznikova, L.;
Uvarov, S.;
Va'vra, J.;
van der Grinten, M. G. D.;
Verbus, J. R.;
Vuosalo, C. O.;
Waldron, W. L.;
Wang, L.;
Webb, R. C.;
Wei, W. -Z.;
While, M.;
White, D. T.;
Whitis, T. J.;
Wisniewski, W. J.;
Witherell, M. S.;
Wolfs, F. L. H.;
Woods, E.;
Woodward, D.;
Worm, S. D.;
Yeh, M.;
Yin, J.;
Young, S. K.;
Zhang, C.
Submitted: 2015-09-09, last modified: 2015-09-23
The design and performance of the LUX-ZEPLIN (LZ) detector is described as of
March 2015 in this Conceptual Design Report. LZ is a second-generation
dark-matter detector with the potential for unprecedented sensitivity to weakly
interacting massive particles (WIMPs) of masses from a few GeV/c2 to hundreds
of TeV/c2. With total liquid xenon mass of about 10 tonnes, LZ will be the most
sensitive experiment for WIMPs in this mass region by the end of the decade.
This report describes in detail the design of the LZ technical systems.
Expected backgrounds are quantified and the performance of the experiment is
presented. The LZ detector will be located at the Sanford Underground Research
Facility in South Dakota. The organization of the LZ Project and a summary of
the expected cost and current schedule are given.
[56]
oai:arXiv.org:1507.03656 [pdf] - 1265764
Metal-Poor Stars Observed with the Magellan Telescope. III. New
Extremely and Ultra Metal-Poor Stars from SDSS/SEGUE and Insights on the
Formation of Ultra Metal-Poor Stars
Submitted: 2015-07-13
We report the discovery of one extremely metal-poor (EMP; [Fe/H]<-3) and one
ultra metal-poor (UMP; [Fe/H]<-4) star selected from the SDSS/SEGUE survey.
These stars were identified as EMP candidates based on their medium-resolution
(R~2,000) spectra, and were followed-up with high-resolution (R~35,000)
spectroscopy with the Magellan-Clay Telescope. Their derived chemical
abundances exhibit good agreement with those of stars with similar
metallicities. We also provide new insights on the formation of the UMP stars,
based on comparison with a new set of theoretical models of supernovae
nucleosynthesis. The models were matched with 20 UMP stars found in the
literature, together with one of the program stars (SDSS J1204+1201), with
[Fe/H]=-4.34. From fitting their abundances, we find that the supernovae
progenitors, for stars where carbon and nitrogen are measured, had masses
ranging from 20.5 M_sun to 28 M_sun and explosion energies from 0.3 to
0.9x10^51 erg. These results are highly sensitive to the carbon and nitrogen
abundance determinations, which is one of the main drivers for future
high-resolution follow-up of UMP candidates. In addition, we are able to
reproduce the different CNO abundance patterns found in UMP stars with a single
progenitor type, by varying its mass and explosion energy.
[57]
oai:arXiv.org:1505.03756 [pdf] - 1183013
Nucleosynthesis in a Primordial Supernova: Carbon and Oxygen Abundances
in SMSS J031300.36-670839.31
Bessell, Michael;
Collett, Remo;
Keller, Stefan;
Frebel, Anna;
Heger, Alexander;
Casey, Andrew;
Masseron, Thomas;
Asplund, Martin;
Jacobson, Heather;
Lind, Karin;
Marino, Anna;
Norris, John;
Yong, David;
Da Costa, Gary;
Chan, Conrad;
Magic, Zazralt;
Schmidt, Brian;
Tisserand, Patrick
Submitted: 2015-05-14, last modified: 2015-05-15
SMSS J031300.36-670839.3 (hereafter SM0313-6708) is a sub-giant halo star,
with no detectable Fe lines and large overabundances of C and Mg relative to
Ca. We obtained VLT-UVES spectra extending to 3060 Angstroms showing strong OH
A-X band lines enabling an oxygen abundance to be derived. The OH A-X band
lines in SM0313-6708 are much stronger than the CH C-X band lines. Spectrum
synthesis fits indicate an [O/C] ratio of 0.02 +- 0.175. Our high S/N UVES data
also enabled us to lower the Fe abundance limit to [Fe/H]{3D},NLTE < -7.52 (3
sigma). These data support our previous suggestion that the star formed from
the iron-poor ejecta of a single massive star Population III supernova.
[58]
oai:arXiv.org:1505.01500 [pdf] - 1296148
Fast Variability and mm/IR flares in GRMHD Models of Sgr A* from
Strong-Field Gravitational Lensing
Submitted: 2015-05-06
We explore the variability properties of long, high cadence GRMHD simulations
across the electromagnetic spectrum using an efficient, GPU-based radiative
transfer algorithm. We focus on both disk- and jet-dominated simulations with
parameters that successfully reproduce the time-averaged spectral properties of
Sgr A* and the size of its image at 1.3mm. We find that the disk-dominated
models produce short timescale variability with amplitudes and power spectra
that closely resemble those inferred observationally. In contrast,
jet-dominated models generate only slow variability, at lower flux levels.
Neither set of models show any X-ray flares, which most likely indicate that
additional physics, such as particle acceleration mechanisms, need to be
incorporated into the GRMHD simulations to account for them. The disk-dominated
models show strong, short-lived mm/IR flares, with short (<~ 1hr) time lags
between the mm and IR wavelengths, that arise from strong-field gravitational
lensing of magnetic flux tubes near the horizon. Such events provide a natural
explanation for the observed IR flares with no X-ray counterparts.
[59]
oai:arXiv.org:1410.3492 [pdf] - 1222552
The Power of Imaging: Constraining the Plasma Properties of GRMHD
Simulations using EHT Observations of Sgr A*
Submitted: 2014-10-13
Recent advances in general relativistic magnetohydrodynamic simulations have
expanded and improved our understanding of the dynamics of black-hole accretion
disks. However, current simulations do not capture the thermodynamics of
electrons in the low density accreting plasma. This poses a significant
challenge in predicting accretion flow images and spectra from first
principles. Because of this, simplified emission models have often been used,
with widely different configurations (e.g., disk- versus jet-dominated
emission), and were able to account for the observed spectral properties of
accreting black-holes. Exploring the large parameter space introduced by such
models, however, requires significant computational power that exceeds
conventional computational facilities. In this paper, we use GRay, a fast
GPU-based ray-tracing algorithm, on the GPU cluster El Gato, to compute images
and spectra for a set of six general relativistic magnetohydrodynamic
simulations with different magnetic field configurations and black-hole spins.
We also employ two different parametric models for the plasma thermodynamics in
each of the simulations. We show that, if only the spectral properties of Sgr
A* are used, all twelve models tested here can fit the spectra equally well.
However, when combined with the measurement of the image size of the emission
using the Event Horizon Telescope, current observations rule out all models
with strong funnel emission, because the funnels are typically very extended.
Our study shows that images of accretion flows with horizon-scale resolution
offer a powerful tool in understanding accretion flows around black-holes and
their thermodynamic properties.
[60]
oai:arXiv.org:1403.1299 [pdf] - 862558
Radiogenic and Muon-Induced Backgrounds in the LUX Dark Matter Detector
Akerib, D. S.;
Araujo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Bernard, E.;
Bernstein, A.;
Bradley, A.;
Byram, D.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Chapman, J. J.;
Chiller, A. A.;
Chiller, C.;
Coffey, T.;
Currie, A.;
de Viveiros, L.;
Dobi, A.;
Dobson, J.;
Druszkiewicz, E.;
Edwards, B.;
Faham, C. H.;
Fiorucci, S.;
Flores, C.;
Gaitskell, R. J.;
Gehman, V. M.;
Ghag, C.;
Gibson, K. R.;
Gilchriese, M. G. D.;
Hall, C.;
Hertel, S. A.;
Horn, M.;
Huang, D. Q.;
Ihm, M.;
Jacobsen, R. G.;
Kazkaz, K.;
Knoche, R.;
Larsen, N. A.;
Lee, C.;
Lindote, A.;
Lopes, M. I.;
Malling, D. C.;
Mannino, R.;
McKinsey, D. N.;
Mei, D. -M.;
Mock, J.;
Moongweluwan, M.;
Morad, J.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H.;
Neves, F.;
Ott, R. A.;
Pangilinan, M.;
Parker, P. D.;
Pease, E. K.;
Pech, K.;
Phelps, P.;
Reichhart, L.;
Shutt, T.;
Silva, C.;
Solovov, V. N.;
Sorensen, P.;
O'Sullivan, K.;
Sumner, T. J.;
Szydagis, M.;
Taylor, D.;
Tennyson, B.;
Tiedt, D. R.;
Tripathi, M.;
Uvarov, S.;
Verbus, J. R.;
Walsh, N.;
Webb, R.;
White, J. T.;
Witherell, M. S.;
Wolfs, F. L. H.;
Woods, M.;
Zhang, C.
Submitted: 2014-03-05
The Large Underground Xenon (LUX) dark matter experiment aims to detect rare
low-energy interactions from Weakly Interacting Massive Particles (WIMPs). The
radiogenic backgrounds in the LUX detector have been measured and compared with
Monte Carlo simulation. Measurements of LUX high-energy data have provided
direct constraints on all background sources contributing to the background
model. The expected background rate from the background model for the 85.3 day
WIMP search run is
$(2.6\pm0.2_{\textrm{stat}}\pm0.4_{\textrm{sys}})\times10^{-3}$~events~keV$_{ee}^{-1}$~kg$^{-1}$~day$^{-1}$
in a 118~kg fiducial volume. The observed background rate is
$(3.6\pm0.4_{\textrm{stat}})\times10^{-3}$~events~keV$_{ee}^{-1}$~kg$^{-1}$~day$^{-1}$,
consistent with model projections. The expectation for the radiogenic
background in a subsequent one-year run is presented.
[61]
oai:arXiv.org:1402.3731 [pdf] - 788195
A Detailed Look at the First Results from the Large Underground Xenon
(LUX) Dark Matter Experiment
Szydagis, M.;
Akerib, D. S.;
Araujo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Bernard, E.;
Bernstein, A.;
Bradley, A.;
Byram, D.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Chapman, J. J.;
Chiller, A. A.;
Chiller, C.;
Coffey, T.;
Currie, A.;
de Viveiros, L.;
Dobi, A.;
Dobson, J.;
Druszkiewicz, E.;
Edwards, B.;
Faham, C. H.;
Fiorucci, S.;
Flores, C.;
Gaitskell, R. J.;
Gehman, V. M.;
Ghag, C.;
Gibson, K. R.;
Gilchriese, M. G. D.;
Hall, C.;
Hertel, S. A.;
Horn, M.;
Huang, D. Q.;
Ihm, M.;
Jacobsen, R. G.;
Kazkaz, K.;
Knoche, R.;
Larsen, N. A.;
Lee, C.;
Lindote, A.;
Lopes, M. I.;
Malling, D. C.;
Mannino, R.;
McKinsey, D. N.;
Mei, D. -M.;
Mock, J.;
Moongweluwan, M.;
Morad, J.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H.;
Neves, F.;
Ott, R. A.;
Pangilinan, M.;
Parker, P. D.;
Pease, E. K.;
Pech, K.;
Phelps, P.;
Reichhart, L.;
Shutt, T.;
Silva, C.;
Solovov, V. N.;
Sorensen, P.;
O'Sullivan, K.;
Sumner, T.;
Taylor, D.;
Tennyson, B.;
Tiedt, D. R.;
Tripathi, M.;
Uvarov, S.;
Verbus, J. R.;
Walsh, N.;
Webb, R.;
White, J. T.;
Witherell, M. S.;
Wolfs, F. L. H.;
Woods, M.;
Zhang, C.
Submitted: 2014-02-15, last modified: 2014-02-25
LUX, the world's largest dual-phase xenon time-projection chamber, with a
fiducial target mass of 118 kg and 10,091 kg-days of exposure thus far, is
currently the most sensitive direct dark matter search experiment. The initial
null-result limit on the spin-independent WIMP-nucleon scattering cross-section
was released in October 2013, with a primary scintillation threshold of 2 phe,
roughly 3 keVnr for LUX. The detector has been deployed at the Sanford
Underground Research Facility (SURF) in Lead, South Dakota, and is the first
experiment to achieve a limit on the WIMP cross-section lower than $10^{-45}$
cm$^{2}$. Here we present a more in-depth discussion of the novel energy scale
employed to better understand the nuclear recoil light and charge yields, and
of the calibration sources, including the new internal tritium source. We found
the LUX data to be in conflict with low-mass WIMP signal interpretations of
other results.
[62]
oai:arXiv.org:1310.8214 [pdf] - 795009
First results from the LUX dark matter experiment at the Sanford
Underground Research Facility
LUX Collaboration;
Akerib, D. S.;
Araujo, H. M.;
Bai, X.;
Bailey, A. J.;
Balajthy, J.;
Bedikian, S.;
Bernard, E.;
Bernstein, A.;
Bolozdynya, A.;
Bradley, A.;
Byram, D.;
Cahn, S. B.;
Carmona-Benitez, M. C.;
Chan, C.;
Chapman, J. J.;
Chiller, A. A.;
Chiller, C.;
Clark, K.;
Coffey, T.;
Currie, A.;
Curioni, A.;
Dazeley, S.;
de Viveiros, L.;
Dobi, A.;
Dobson, J.;
Dragowsky, E. M.;
Druszkiewicz, E.;
Edwards, B.;
Faham, C. H.;
Fiorucci, S.;
Flores, C.;
Gaitskell, R. J.;
Gehman, V. M.;
Ghag, C.;
Gibson, K. R.;
Gilchriese, M. G. D.;
Hall, C.;
Hanhardt, M.;
Hertel, S. A.;
Horn, M.;
Huang, D. Q.;
Ihm, M.;
Jacobsen, R. G.;
Kastens, L.;
Kazkaz, K.;
Knoche, R.;
Kyre, S.;
Lander, R.;
Larsen, N. A.;
Lee, C.;
Leonard, D. S.;
Lesko, K. T.;
Lindote, A.;
Lopes, M. I.;
Lyashenko, A.;
Malling, D. C.;
Mannino, R.;
McKinsey, D. N.;
Mei, D. -M.;
Mock, J.;
Moongweluwan, M.;
Morad, J.;
Morii, M.;
Murphy, A. St. J.;
Nehrkorn, C.;
Nelson, H.;
Neves, F.;
Nikkel, J. A.;
Ott, R. A.;
Pangilinan, M.;
Parker, P. D.;
Pease, E. K.;
Pech, K.;
Phelps, P.;
Reichhart, L.;
Shutt, T.;
Silva, C.;
Skulski, W.;
Sofka, C. J.;
Solovov, V. N.;
Sorensen, P.;
Stiegler, T.;
O`Sullivan, K.;
Sumner, T. J.;
Svoboda, R.;
Sweany, M.;
Szydagis, M.;
Taylor, D.;
Tennyson, B.;
Tiedt, D. R.;
Tripathi, M.;
Uvarov, S.;
Verbus, J. R.;
Walsh, N.;
Webb, R.;
White, J. T.;
White, D.;
Witherell, M. S.;
Wlasenko, M.;
Wolfs, F. L. H.;
Woods, M.;
Zhang, C.
Submitted: 2013-10-30, last modified: 2014-02-05
The Large Underground Xenon (LUX) experiment, a dual-phase xenon
time-projection chamber operating at the Sanford Underground Research Facility
(Lead, South Dakota), was cooled and filled in February 2013. We report results
of the first WIMP search dataset, taken during the period April to August 2013,
presenting the analysis of 85.3 live-days of data with a fiducial volume of 118
kg. A profile-likelihood analysis technique shows our data to be consistent
with the background-only hypothesis, allowing 90% confidence limits to be set
on spin-independent WIMP-nucleon elastic scattering with a minimum upper limit
on the cross section of $7.6 \times 10^{-46}$ cm$^{2}$ at a WIMP mass of 33
GeV/c$^2$. We find that the LUX data are in strong disagreement with low-mass
WIMP signal interpretations of the results from several recent direct detection
experiments.
[63]
oai:arXiv.org:1303.5057 [pdf] - 1165424
GRay: a Massively Parallel GPU-Based Code for Ray Tracing in
Relativistic Spacetimes
Submitted: 2013-03-20
We introduce GRay, a massively parallel integrator designed to trace the
trajectories of billions of photons in a curved spacetime. This GPU-based
integrator employs the stream processing paradigm, is implemented in CUDA
C/C++, and runs on nVidia graphics cards. The peak performance of GRay using
single precision floating-point arithmetic on a single GPU exceeds 300 GFLOP
(or 1 nanosecond per photon per time step). For a realistic problem, where the
peak performance cannot be reached, GRay is two orders of magnitude faster than
existing CPU-based ray tracing codes. This performance enhancement allows more
effective searches of large parameter spaces when comparing theoretical
predictions of images, spectra, and lightcurves from the vicinities of compact
objects to observations. GRay can also perform on-the-fly ray tracing within
general relativistic magnetohydrodynamic algorithms that simulate accretion
flows around compact objects. Making use of this algorithm, we calculate the
properties of the shadows of Kerr black holes and the photon rings that
surround them. We also provide accurate fitting formulae of their dependencies
on black hole spin and observer inclination, which can be used to interpret
upcoming observations of the black holes at the center of the Milky Way, as
well as M87, with the Event Horizon Telescope.
[64]
oai:arXiv.org:1211.6150 [pdf] - 595129
MHD Turbulence in Accretion Disk Boundary Layers
Submitted: 2012-11-26
The physical modeling of the accretion disk boundary layer, the region where
the disk meets the surface of the accreting star, usually relies on the
assumption that angular momentum transport is opposite to the radial angular
frequency gradient of the disk. The standard model for turbulent shear
viscosity, widely adopted in astrophysics, satisfies this assumption by
construction. However, this behavior is not supported by numerical simulations
of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular
momentum transport driven by the magnetorotational instability is inefficient
in this inner disk region. I will discuss the results of a recent study on the
generation of hydromagnetic stresses and energy density in the boundary layer
around a weakly magnetized star. Our findings suggest that although magnetic
energy density can be significantly amplified in this region, angular momentum
transport is rather inefficient. This seems consistent with the results
obtained in numerical simulations and suggests that the detailed structure of
turbulent MHD boundary layers could differ appreciably from those derived
within the standard framework of turbulent shear viscosity.
[65]
oai:arXiv.org:1111.4219 [pdf] - 516449
On Hydromagnetic Stresses in Accretion Disk Boundary Layers
Submitted: 2011-11-17, last modified: 2012-05-23
Detailed calculations of the physical structure of accretion disk boundary
layers, and thus their inferred observational properties, rely on the
assumption that angular momentum transport is opposite to the radial angular
frequency gradient of the disk. The standard model for turbulent shear
viscosity satisfies this assumption by construction. However, this behavior is
not supported by numerical simulations of turbulent magnetohydrodynamic (MHD)
accretion disks, which show that angular momentum transport driven by the
magnetorotational instability (MRI) is inefficient in disk regions where, as
expected in boundary layers, the angular frequency increases with radius. In
order to shed light into physically viable mechanisms for angular momentum
transport in this inner disk region, we examine the generation of hydromagnetic
stresses and energy density in differentially rotating backgrounds with angular
frequencies that increase outward in the shearing-sheet framework. We isolate
the modes that are unrelated to the standard MRI and provide analytic solutions
for the long-term evolution of the resulting shearing MHD waves. We show that,
although the energy density of these waves can be amplified significantly,
their associated stresses oscillate around zero, rendering them an inefficient
mechanism to transport significant angular momentum (inward). These findings
are consistent with the results obtained in numerical simulations of MHD
accretion disk boundary layers and challenge the standard assumption of
efficient angular momentum transport in the inner disk regions. This suggests
that the detailed structure of turbulent MHD accretion disk boundary layers
could differ appreciably from those derived within the standard framework of
turbulent shear viscosity.
[66]
oai:arXiv.org:1109.6937 [pdf] - 494909
Dynamics of Saturated Energy Condensation in Two-Dimensional Turbulence
Submitted: 2011-09-30, last modified: 2012-03-31
In two-dimensional forced Navier-Stokes turbulence, energy cascades to the
largest scales in the system to form a pair of coherent vortices known as the
Bose condensate. We show, both numerically and analytically, that the energy
condensation saturates and the system reaches a statistically stationary state.
The time scale of saturation is inversely proportional to the viscosity and the
saturation energy level is determined by both the viscosity and the force. We
further show that, without sufficient resolution to resolve the small-scale
enstrophy spectrum, numerical simulations can give a spurious result for the
saturation energy level. We also find that the movement of the condensate is
similar to the motion of an inertial particle with an effective drag force.
Furthermore, we show that the profile of the saturated coherent vortices can be
described by a Gaussian core with exponential wings.
[67]
oai:arXiv.org:0911.5351 [pdf] - 302001
A Class of Physically Motivated Closures for Radiation Hydrodynamics
Submitted: 2009-11-30, last modified: 2011-02-02
Radiative transfer and radiation hydrodynamics use the relativistic Boltzmann
equation to describe the kinetics of photons. It is difficult to solve the
six-dimensional time-dependent transfer equation unless the problem is highly
symmetric or in equilibrium. When the radiation field is smooth, it is natural
to take angular moments of the transfer equation to reduce the degrees of
freedom. However, low order moment equations contain terms that depend on
higher order moments. To close the system of moment equations, approximations
are made to truncate this hierarchy. Popular closures used in astrophysics
include flux limited diffusion and the M1 closure, which are rather ad hoc and
do not necessarily capture the correct physics. In this paper, we propose a new
class of closures for radiative transfer and radiation hydrodynamics. We start
from a different perspective and highlight the consistency of a fully
relativistic formalism. We present a generic framework to approximate radiative
transfer based on relativistic Grad's moment method. We then derive a 14-field
method that minimizes unphysical photon self-interaction.
[68]
oai:arXiv.org:0812.2031 [pdf] - 900519
Oscillations of the Inner Regions of Viscous Accretion Disks
Submitted: 2008-12-11
Although quasi-periodic oscillations (QPOs) have been discovered in different
X-ray sources, their origin is still a matter of debate. Analytical studies of
hydrodynamic accretion disks have shown three types of trapped global modes
with properties that appear to agree with the observations. However, these
studies take only linear effects into account and do not address the issues of
mode excitation and decay. Moreover, observations suggest that resonances
between modes play a crucial role. A systematic, numerical study of this
problem is therefore needed. In this paper, we use a pseudo-spectral algorithm
to perform a parameter study of the inner regions of hydrodynamic disks. By
assuming alpha-viscosity, we show that steady state solutions rarely exist. The
inner edges of the disks oscillate and excite axisymmetric waves. In addition,
the flows inside the inner edges are sometimes unstable to non-axisymmetric
perturbations. One-armed, or even two-armed, spirals are developed, which
provides a plausible explanation for the high-frequency QPOs observed from
accreting black holes. When the Reynolds numbers are above certain critical
values, the inner disks go through some transient turbulent states
characterized by strong trailing spirals; while large-scale leading spirals
developed in the outer disks. We compared our numerical results with standard
thin disk oscillation models. Although the non-axisymmetric features have their
analytical counterparts, more careful study is needed to explain the
axisymmetric oscillations.
[69]
oai:arXiv.org:astro-ph/0612404 [pdf] - 87742
The fundamental difference between shear alpha viscosity and turbulent
magnetorotational stresses
Submitted: 2006-12-15, last modified: 2008-02-18
Numerical simulations of turbulent, magnetized, differentially rotating flows
driven by the magnetorotational instability are often used to calculate the
effective values of alpha viscosity that is invoked in analytical models of
accretion discs. In this paper we use various dynamical models of turbulent
magnetohydrodynamic stresses, as well as numerical simulations of shearing
boxes, to show that angular momentum transport in MRI-driven accretion discs
cannot be described by the standard model for shear viscosity. In particular,
we demonstrate that turbulent magnetorotational stresses are not linearly
proportional to the local shear and vanish identically for angular velocity
profiles that increase outwards.
[70]
oai:arXiv.org:0801.4570 [pdf] - 9541
Viscous, Resistive Magnetorotational Modes
Submitted: 2008-01-29
We carry out a comprehensive analysis of the behavior of the
magnetorotational instability (MRI) in viscous, resistive plasmas. We find
exact, non-linear solutions of the non-ideal magnetohydrodynamic (MHD)
equations describing the local dynamics of an incompressible, differentially
rotating background threaded by a vertical magnetic field when disturbances
with wavenumbers perpendicular to the shear are considered. We provide a
geometrical description of these viscous, resistive MRI modes and show how
their physical structure is modified as a function of the Reynolds and magnetic
Reynolds numbers. We demonstrate that when finite dissipative effects are
considered, velocity and magnetic field disturbances are no longer orthogonal
(as it is the case in the ideal MHD limit) unless the magnetic Prandtl number
is unity. We generalize previous results found in the ideal limit and show that
a series of key properties of the mean Reynolds and Maxwell stresses also hold
for the viscous, resistive MRI. In particular, we show that the Reynolds stress
is always positive and the Maxwell stress is always negative. Therefore, even
in the presence of viscosity and resistivity, the total mean angular momentum
transport is always directed outwards. We also find that, for any combination
of the Reynolds and magnetic Reynolds numbers, magnetic disturbances dominate
both the energetics and the transport of angular momentum and that the total
mean energy density is an upper bound for the total mean stress responsible for
angular momentum transport. The ratios between the Maxwell and Reynolds
stresses and between magnetic and kinetic energy densities increase with
decreasing Reynolds numbers for any magnetic Reynolds number; the lowest limit
of both ratios is reached in the ideal MHD regime.
[71]
oai:arXiv.org:0705.0352 [pdf] - 942
Angular Momentum Transport in Accretion Disks: Scaling Laws in
MRI-driven Turbulence
Submitted: 2007-05-02, last modified: 2007-08-20
We present a scaling law that predicts the values of the stresses obtained in
numerical simulations of saturated MRI-driven turbulence in non-stratified
shearing boxes. It relates the turbulent stresses to the strength of the
vertical magnetic field, the sound speed, the vertical size of the box, and the
numerical resolution and predicts accurately the results of 35 numerical
simulations performed for a wide variety of physical conditions. We use our
result to show that the saturated stresses in simulations with zero net
magnetic flux depend linearly on the numerical resolution and would become
negligible if the resolution were set equal to the natural dissipation scale in
astrophysical disks. We conclude that, in order for MRI-driven turbulent
angular momentum transport to be able to account for the large value of the
effective alpha viscosity inferred observationally, the disk must be threaded
by a significant vertical magnetic field and the turbulent magnetic energy must
be in near equipartition with the thermal energy. This result has important
implications for the spectra of accretion disks and their stability.
[72]
oai:arXiv.org:astro-ph/0612742 [pdf] - 88080
Spectral Methods for Time-Dependent Studies of Accretion Flows. III.
Three-Dimensional, Self-gravitating, Magnetohydrodynamic Disks
Submitted: 2006-12-27
Accretion disks are three-dimensional, turbulent, often self-gravitating,
magnetohydrodynamic flows, which can be modeled in detail with numerical
simulations. In this paper, we present a new algorithm that is based on a
spectral decomposition method to simulate such flows. Because of the high order
of the method, we can solve the induction equation in terms of the magnetic
potential and, therefore, ensure trivially that the magnetic fields in the
numerical solution are divergence free. The spectral method also suffers
minimally from numerical dissipation and allows for an easy implementation of
models for sub-grid physics. Both properties make our method ideal for studying
MHD turbulent flows such as those found in accretion disks around compact
objects. We verify our algorithm with a series of standard tests and use it to
show the development of MHD turbulnce in a simulation of an accretion disk.
Finally, we study the evolution and saturation of the power spectrum of MHD
turbulence driven by the magnetorotational instability.
[73]
oai:arXiv.org:astro-ph/0611269 [pdf] - 316597
MHD Simulations of Accretion onto Sgr A*: Quiescent Fluctuations,
Outbursts, and Quasi-Periodicity
Submitted: 2006-11-08
High resolution observations of Sgr A* have revealed a wide variety of
phenomena, ranging from intense rapid flares to quasi-periodic oscillations,
making this object an ideal system to study the properties of low luminosity
accreting black holes. In this paper, we use a pseudo-spectral algorithm to
construct and evolve a three-dimensional magnetohydrodynamic model of the
accretion disk in Sgr A*. Assuming a hybrid thermal-nonthermal emission scheme,
we show that the MHD turbulence can by itself only produce factor of two
fluctuations in luminosity. These amplitudes in variation cannot explain the
magnitude of flares observed in this system. However, we also demonstrate that
density perturbations in the disk do produce outbursts qualitatively similar to
those observed by XMM-Newton in X-rays and ground-based facilities in the near
infrared. Quasi-periodic oscillations emerge naturally in the simulated
lightcurves. We attribute these to non-axisymmetric density perturbations that
emerge as the disk evolves back toward its quiescent state.
[74]
oai:arXiv.org:astro-ph/0610565 [pdf] - 85989
A Local Model for Angular Momentum Transport in Accretion Disks Driven
by the Magnetorotational Instability
Submitted: 2006-10-18
We develop a local model for the exponential growth and saturation of the
Reynolds and Maxwell stresses in turbulent flows driven by the
magnetorotational instability. We first derive equations that describe the
effects of the instability on the growth and pumping of the stresses. We
highlight the relevance of a new type of correlations that couples the
dynamical evolution of the Reynolds and Maxwell stresses and plays a key role
in developing and sustaining the magnetorotational turbulence. We then
supplement these equations with a phenomenological description of the triple
correlations that lead to a saturated turbulent state. We show that the
steady-state limit of the model describes successfully the correlations among
stresses found in numerical simulations of shearing boxes.
[75]
oai:arXiv.org:astro-ph/0603178 [pdf] - 80387
The signature of the magnetorotational instability in the Reynolds and
Maxwell stress tensors in accretion discs
Submitted: 2006-03-08, last modified: 2006-07-20
The magnetorotational instability is thought to be responsible for the
generation of magnetohydrodynamic turbulence that leads to enhanced outward
angular momentum transport in accretion discs. Here, we present the first
formal analytical proof showing that, during the exponential growth of the
instability, the mean (averaged over the disc scale-height) Reynolds stress is
always positive, the mean Maxwell stress is always negative, and hence the mean
total stress is positive and leads to a net outward flux of angular momentum.
More importantly, we show that the ratio of the Maxwell to the Reynolds
stresses during the late times of the exponential growth of the instability is
determined only by the local shear and does not depend on the initial spectrum
of perturbations or the strength of the seed magnetic. Even though we derived
these properties of the stress tensors for the exponential growth of the
instability in incompressible flows, numerical simulations of shearing boxes
show that this characteristic is qualitatively preserved under more general
conditions, even during the saturated turbulent state generated by the
instability.
[76]
oai:arXiv.org:astro-ph/0512448 [pdf] - 78648
Spectral Methods for Time-Dependent Studies of Accretion Flows. II.
Two-Dimensional Hydrodynamic Disks with Self-Gravity
Submitted: 2005-12-16
Spectral methods are well suited for solving hydrodynamic problems in which
the self-gravity of the flow needs to be considered. Because Poisson's equation
is linear, the numerical solution for the gravitational potential for each
individual mode of the density can be pre-computed, thus reducing substantially
the computational cost of the method. In this second paper, we describe two
different approaches to computing the gravitational field of a two-dimensional
flow with pseudo-spectral methods. For situations in which the density profile
is independent of the third coordinate (i.e., an infinite cylinder), we use a
standard Poisson solver in spectral space. On the other hand, for situations in
which the density profile is a delta function along the third coordinate (i.e.,
an infinitesimally thin disk), or any other function known a priori, we perform
a direct integration of Poisson's equation using a Green's functions approach.
We devise a number of test problems to verify the implementations of these two
methods. Finally, we use our method to study the stability of polytropic,
self-gravitating disks. We find that, when the polytropic index Gamma is <=
4/3, Toomre's criterion correctly describes the stability of the disk. However,
when Gamma > 4/3 and for large values of the polytropic constant K, the
numerical solutions are always stable, even when the linear criterion predicts
the contrary. We show that, in the latter case, the minimum wavelength of the
unstable modes is larger than the extent of the unstable region and hence the
local linear analysis is inapplicable.
[77]
oai:arXiv.org:astro-ph/0406073 [pdf] - 65230
Spectral Methods for Time-Dependent Studies of Accretion Flows. I.
Two-dimensional, Viscous, Hydrodynamic Disks
Submitted: 2004-06-03, last modified: 2005-03-26
We present a numerical method for studying the normal modes of accretion
flows around black holes. In this first paper, we focus on two-dimensional,
viscous, hydrodynamic disks, for which the linear modes have been calculated
analytically in previous investigations. We use pseudo-spectral methods and low
storage Runge-Kutta methods to solve the continuity equation, the Navier-Stokes
equation, and the energy equation. We devise a number of test problems to
verify the implementation. These tests demonstrate the ability of spectral
methods to handle accurately advection problems and to reproduce correctly the
stability criteria for differentially rotating hydrodynamic flows. They also
show that our implementation is able to handle sound wave correctly with
non-reflective boundary conditions, to recover the standard solution for a
viscous spreading ring, and produce correctly the Shakura-Sunyaev steady disk
solution. Finally, we have applied our algorithm to the problem of a
non-axisymmetric viscous spreading ring and verify that such configuration is
unstable to non-antisymmetric perturbations.
[78]
oai:arXiv.org:astro-ph/0309374 [pdf] - 59201
Recombination Induced Softening and Reheating of the Cosmic Plasma
Submitted: 2003-09-13, last modified: 2003-12-10
The atomic recombination process leads to a softening of the matter equation
of state as reflected by a reduced generalized adiabatic index, with
accompanying heat release. We study the effects of this recombination softening
and reheating of the cosmic plasma on the ionization history, visibility
function, Cold Dark Matter (CDM) transfer function, and the Cosmic Microwave
Background (CMB) spectra. The resulting modifications of the CMB spectrm is
1/10 of WMAP's current error and is comparable to PLANCK's error. Therefore,
this effect should be considered when data with higher accuracy are analysed.
