-
Accelerating Time-to-Science by Streaming Detector Data Directly into Perlmutter Compute Nodes
Authors:
Samuel S. Welborn,
Bjoern Enders,
Chris Harris,
Peter Ercius,
Deborah J. Bard
Abstract:
Recent advancements in detector technology have significantly increased the size and complexity of experimental data, and high-performance computing (HPC) provides a path towards more efficient and timely data processing. However, movement of large data sets from acquisition systems to HPC centers introduces bottlenecks owing to storage I/O at both ends. This manuscript introduces a streaming work…
▽ More
Recent advancements in detector technology have significantly increased the size and complexity of experimental data, and high-performance computing (HPC) provides a path towards more efficient and timely data processing. However, movement of large data sets from acquisition systems to HPC centers introduces bottlenecks owing to storage I/O at both ends. This manuscript introduces a streaming workflow designed for an high data rate electron detector that streams data directly to compute node memory at the National Energy Research Scientific Computing Center (NERSC), thereby avoiding storage I/O. The new workflow deploys ZeroMQ-based services for data production, aggregation, and distribution for on-the-fly processing, all coordinated through a distributed key-value store. The system is integrated with the detector's science gateway and utilizes the NERSC Superfacility API to initiate streaming jobs through a web-based frontend. Our approach achieves up to a 14-fold increase in data throughput and enhances predictability and reliability compared to a I/O-heavy file-based transfer workflow. Our work highlights the transformative potential of streaming workflows to expedite data analysis for time-sensitive experiments.
△ Less
Submitted 2 July, 2024; v1 submitted 21 March, 2024;
originally announced March 2024.
-
Workflows Community Summit 2022: A Roadmap Revolution
Authors:
Rafael Ferreira da Silva,
Rosa M. Badia,
Venkat Bala,
Debbie Bard,
Peer-Timo Bremer,
Ian Buckley,
Silvina Caino-Lores,
Kyle Chard,
Carole Goble,
Shantenu Jha,
Daniel S. Katz,
Daniel Laney,
Manish Parashar,
Frederic Suter,
Nick Tyler,
Thomas Uram,
Ilkay Altintas,
Stefan Andersson,
William Arndt,
Juan Aznar,
Jonathan Bader,
Bartosz Balis,
Chris Blanton,
Kelly Rosa Braghetto,
Aharon Brodutch
, et al. (80 additional authors not shown)
Abstract:
Scientific workflows have become integral tools in broad scientific computing use cases. Science discovery is increasingly dependent on workflows to orchestrate large and complex scientific experiments that range from execution of a cloud-based data preprocessing pipeline to multi-facility instrument-to-edge-to-HPC computational workflows. Given the changing landscape of scientific computing and t…
▽ More
Scientific workflows have become integral tools in broad scientific computing use cases. Science discovery is increasingly dependent on workflows to orchestrate large and complex scientific experiments that range from execution of a cloud-based data preprocessing pipeline to multi-facility instrument-to-edge-to-HPC computational workflows. Given the changing landscape of scientific computing and the evolving needs of emerging scientific applications, it is paramount that the development of novel scientific workflows and system functionalities seek to increase the efficiency, resilience, and pervasiveness of existing systems and applications. Specifically, the proliferation of machine learning/artificial intelligence (ML/AI) workflows, need for processing large scale datasets produced by instruments at the edge, intensification of near real-time data processing, support for long-term experiment campaigns, and emergence of quantum computing as an adjunct to HPC, have significantly changed the functional and operational requirements of workflow systems. Workflow systems now need to, for example, support data streams from the edge-to-cloud-to-HPC enable the management of many small-sized files, allow data reduction while ensuring high accuracy, orchestrate distributed services (workflows, instruments, data movement, provenance, publication, etc.) across computing and user facilities, among others. Further, to accelerate science, it is also necessary that these systems implement specifications/standards and APIs for seamless (horizontal and vertical) integration between systems and applications, as well as enabling the publication of workflows and their associated products according to the FAIR principles. This document reports on discussions and findings from the 2022 international edition of the Workflows Community Summit that took place on November 29 and 30, 2022.
△ Less
Submitted 31 March, 2023;
originally announced April 2023.
-
The LBNL Superfacility Project Report
Authors:
Deborah Bard,
Cory Snavely,
Lisa Gerhardt,
Jason Lee,
Becci Totzke,
Katie Antypas,
William Arndt,
Johannes Blaschke,
Suren Byna,
Ravi Cheema,
Shreyas Cholia,
Mark Day,
Bjoern Enders,
Aditi Gaur,
Annette Greiner,
Taylor Groves,
Mariam Kiran,
Quincey Koziol,
Tom Lehman,
Kelly Rowland,
Chris Samuel,
Ashwin Selvarajan,
Alex Sim,
David Skinner,
Laurie Stephey
, et al. (2 additional authors not shown)
Abstract:
The Superfacility model is designed to leverage HPC for experimental science. It is more than simply a model of connected experiment, network, and HPC facilities; it encompasses the full ecosystem of infrastructure, software, tools, and expertise needed to make connected facilities easy to use. The three-year Lawrence Berkeley National Laboratory (LBNL) Superfacility project was initiated in 2019…
▽ More
The Superfacility model is designed to leverage HPC for experimental science. It is more than simply a model of connected experiment, network, and HPC facilities; it encompasses the full ecosystem of infrastructure, software, tools, and expertise needed to make connected facilities easy to use. The three-year Lawrence Berkeley National Laboratory (LBNL) Superfacility project was initiated in 2019 to coordinate work being performed at LBNL to support this model, and to provide a coherent and comprehensive set of science requirements to drive existing and new work.
A key component of the project was the in-depth engagements with eight science teams that represent challenging use cases across the DOE Office of Science. By the close of the project, we met our project goal by enabling our science application engagements to demonstrate automated pipelines that analyze data from remote facilities at large scale, without routine human intervention. In several cases, we have gone beyond demonstrations and now provide production-level services. To achieve this goal, the Superfacility team developed tools, infrastructure, and policies for near-real-time computing support, dynamic high-performance networking, data management and movement tools, API-driven automation, HPC-scale notebooks via Jupyter, authentication using Federated Identity and container-based edge services supported.
The lessons we learned during this project provide a valuable model for future large, complex, cross-disciplinary collaborations. There is a pressing need for a coherent computing infrastructure across national facilities, and LBNL's Superfacility project is a unique model for success in tackling the challenges that will be faced in hardware, software, policies, and services across multiple science domains.
△ Less
Submitted 27 June, 2022; v1 submitted 23 June, 2022;
originally announced June 2022.
-
Accelerating X-Ray Tracing for Exascale Systems using Kokkos
Authors:
Felix Wittwer,
Nicholas K. Sauter,
Derek Mendez,
Billy K. Poon,
Aaron S. Brewster,
James M. Holton,
Michael E. Wall,
William E. Hart,
Deborah J. Bard,
Johannes P. Blaschke
Abstract:
The upcoming exascale computing systems Frontier and Aurora will draw much of their computing power from GPU accelerators. The hardware for these systems will be provided by AMD and Intel, respectively, each supporting their own GPU programming model. The challenge for applications that harness one of these exascale systems will be to avoid lock-in and to preserve performance portability.
We rep…
▽ More
The upcoming exascale computing systems Frontier and Aurora will draw much of their computing power from GPU accelerators. The hardware for these systems will be provided by AMD and Intel, respectively, each supporting their own GPU programming model. The challenge for applications that harness one of these exascale systems will be to avoid lock-in and to preserve performance portability.
We report here on our results of using Kokkos to accelerate a real-world application on NERSC's Perlmutter Phase 1 (using NVIDIA A100 accelerators) and the testbed system for OLCF's Frontier (using AMD MI250X). By porting to Kokkos, we were able to successfully run the same X-ray tracing code on both systems and achieved speed-ups between 13% and 66% compared to the original CUDA code. These results are a highly encouraging demonstration of using Kokkos to accelerate production science code.
△ Less
Submitted 16 May, 2022;
originally announced May 2022.
-
Climate of the Field: Snowmass 2021
Authors:
Erin V. Hansen,
Erica Smith,
Deborah Bard,
Matthew Bellis,
Jessica Esquivel,
Tiffany R. Lewis,
Cameron Geddes,
Cindy Joe,
Alex G. Kim,
Asmita Patel,
Vitaly Pronskikh
Abstract:
How are formal policies put in place to create an inclusive, equitable, safe environment? How do these differ between different communities of practice (institutions, labs, collaborations, working groups)? What policies towards a more equitable community are working? For those that aren't working, what external support is needed in order to make them more effective? We present a discussion of the…
▽ More
How are formal policies put in place to create an inclusive, equitable, safe environment? How do these differ between different communities of practice (institutions, labs, collaborations, working groups)? What policies towards a more equitable community are working? For those that aren't working, what external support is needed in order to make them more effective? We present a discussion of the current climate of the field in high energy particle physics and astrophysics (HEPA), as well as current efforts toward making the community a more diverse, inclusive, and equitable environment. We also present issues facing both institutions and HEPA collaborations, with a set of interviews with a selection of HEPA collaboration DEI leaders. We encourage the HEPA community and the institutions & agencies that support it to think critically about the prioritization of people in HEPA over the coming decade, and what resources and policies need to be in place in order to protect and elevate minoritized populations within the HEPA community.
△ Less
Submitted 29 September, 2022; v1 submitted 7 April, 2022;
originally announced April 2022.
-
Data Preservation for Cosmology
Authors:
Marcelo Alvarez,
Stephen Bailey,
Deborah Bard,
Lisa Gerhardt,
Julien Guy,
Stéphanie Juneau,
Anthony Kremin,
Brian Nord,
David Schlegel,
Laurie Stephey,
Rollin Thomas,
Benjamin Weaver
Abstract:
We describe the needs and opportunities for preserving cosmology datasets and simulations, and facilitating their joint analysis beyond the lifetime of individual projects. We recommend that DOE fund a new cosmology data archive center to coordinate this work across the multiple DOE computing facilities.
We describe the needs and opportunities for preserving cosmology datasets and simulations, and facilitating their joint analysis beyond the lifetime of individual projects. We recommend that DOE fund a new cosmology data archive center to coordinate this work across the multiple DOE computing facilities.
△ Less
Submitted 15 March, 2022;
originally announced March 2022.
-
Snowmass 2021 CMB-S4 White Paper
Authors:
Kevork Abazajian,
Arwa Abdulghafour,
Graeme E. Addison,
Peter Adshead,
Zeeshan Ahmed,
Marco Ajello,
Daniel Akerib,
Steven W. Allen,
David Alonso,
Marcelo Alvarez,
Mustafa A. Amin,
Mandana Amiri,
Adam Anderson,
Behzad Ansarinejad,
Melanie Archipley,
Kam S. Arnold,
Matt Ashby,
Han Aung,
Carlo Baccigalupi,
Carina Baker,
Abhishek Bakshi,
Debbie Bard,
Denis Barkats,
Darcy Barron,
Peter S. Barry
, et al. (331 additional authors not shown)
Abstract:
This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan.
This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan.
△ Less
Submitted 15 March, 2022;
originally announced March 2022.
-
Snowmass2021 Cosmic Frontier: Modeling, statistics, simulations, and computing needs for direct dark matter detection
Authors:
Yonatan Kahn,
Maria Elena Monzani,
Kimberly J. Palladino,
Tyler Anderson,
Deborah Bard,
Daniel Baxter,
Micah Buuck,
Concetta Cartaro,
Juan I. Collar,
Miriam Diamond,
Alden Fan,
Simon Knapen,
Scott Kravitz,
Rafael F. Lang,
Benjamin Nachman,
Ibles Olcina Samblas,
Igor Ostrovskiy,
Aditya Parikh,
Quentin Riffard,
Amy Roberts,
Kelly Stifter,
Matthew Szydagis,
Christopher Tunnell,
Belina von Krosigk,
Dennis Wright
, et al. (12 additional authors not shown)
Abstract:
This paper summarizes the modeling, statistics, simulation, and computing needs of direct dark matter detection experiments in the next decade.
This paper summarizes the modeling, statistics, simulation, and computing needs of direct dark matter detection experiments in the next decade.
△ Less
Submitted 27 December, 2022; v1 submitted 15 March, 2022;
originally announced March 2022.
-
Software and Computing for Small HEP Experiments
Authors:
Dave Casper,
Maria Elena Monzani,
Benjamin Nachman,
Costas Andreopoulos,
Stephen Bailey,
Deborah Bard,
Wahid Bhimji,
Giuseppe Cerati,
Grigorios Chachamis,
Jacob Daughhetee,
Miriam Diamond,
V. Daniel Elvira,
Alden Fan,
Krzysztof Genser,
Paolo Girotti,
Scott Kravitz,
Robert Kutschke,
Vincent R. Pascuzzi,
Gabriel N. Perdue,
Erica Snider,
Elizabeth Sexton-Kennedy,
Graeme Andrew Stewart,
Matthew Szydagis,
Eric Torrence,
Christopher Tunnell
Abstract:
This white paper briefly summarized key conclusions of the recent US Community Study on the Future of Particle Physics (Snowmass 2021) workshop on Software and Computing for Small High Energy Physics Experiments.
This white paper briefly summarized key conclusions of the recent US Community Study on the Future of Particle Physics (Snowmass 2021) workshop on Software and Computing for Small High Energy Physics Experiments.
△ Less
Submitted 27 December, 2022; v1 submitted 15 March, 2022;
originally announced March 2022.
-
Real-Time XFEL Data Analysis at SLAC and NERSC: a Trial Run of Nascent Exascale Experimental Data Analysis
Authors:
Johannes P. Blaschke,
Aaron S. Brewster,
Daniel W. Paley,
Derek Mendez,
Asmit Bhowmick,
Nicholas K. Sauter,
Wilko Kröger,
Murali Shankar,
Bjoern Enders,
Deborah Bard
Abstract:
X-ray scattering experiments using Free Electron Lasers (XFELs) are a powerful tool to determine the molecular structure and function of unknown samples (such as COVID-19 viral proteins). XFEL experiments are a challenge to computing in two ways: i) due to the high cost of running XFELs, a fast turnaround time from data acquisition to data analysis is essential to make informed decisions on experi…
▽ More
X-ray scattering experiments using Free Electron Lasers (XFELs) are a powerful tool to determine the molecular structure and function of unknown samples (such as COVID-19 viral proteins). XFEL experiments are a challenge to computing in two ways: i) due to the high cost of running XFELs, a fast turnaround time from data acquisition to data analysis is essential to make informed decisions on experimental protocols; ii) data collection rates are growing exponentially, requiring new scalable algorithms. Here we report our experiences analyzing data from two experiments at the Linac Coherent Light Source (LCLS) during September 2020. Raw data were analyzed on NERSC's Cori XC40 system, using the Superfacility paradigm: our workflow automatically moves raw data between LCLS and NERSC, where it is analyzed using the software package CCTBX. We achieved real time data analysis with a turnaround time from data acquisition to full molecular reconstruction in as little as 10 min -- sufficient time for the experiment's operators to make informed decisions. By hosting the data analysis on Cori, and by automating LCLS-NERSC interoperability, we achieved a data analysis rate which matches the data acquisition rate. Completing data analysis with 10 mins is a first for XFEL experiments and an important milestone if we are to keep up with data collection trends.
△ Less
Submitted 31 December, 2023; v1 submitted 21 June, 2021;
originally announced June 2021.
-
Quantum Advantage on Proof of Work
Authors:
Dan A. Bard,
Joseph J. Kearney,
Carlos A. Perez-Delgado
Abstract:
Proof-of-Work (PoW) is a fundamental underlying technology behind most major blockchain cryptocurrencies. It has been previously pointed out that quantum devices provide a computational advantage in performing PoW in the context of Bitcoin. Here we make the case that this quantum advantage extends not only to all existing PoW mechanisms, but to any possible PoW as well. This has strong consequence…
▽ More
Proof-of-Work (PoW) is a fundamental underlying technology behind most major blockchain cryptocurrencies. It has been previously pointed out that quantum devices provide a computational advantage in performing PoW in the context of Bitcoin. Here we make the case that this quantum advantage extends not only to all existing PoW mechanisms, but to any possible PoW as well. This has strong consequences regarding both quantum-based attacks on the integrity of the entirety of the blockchain, as well as more legitimate uses of quantum computation for the purpose of mining Bitcoin and other cryptocurrencies. For the first case, we estimate when these quantum attacks will become feasible, for various cryptocurrencies, and discuss the impact of such attacks. For the latter, we derive a precise formula to calculate the economic incentive for switching to quantum-based cryptocurrency miners. Using this formula, we analyze several test scenarios, and conclude that investing in quantum hardware for cryptocurrency mining has the potential to pay off immensely.
△ Less
Submitted 4 May, 2021;
originally announced May 2021.
-
CMB-S4: Forecasting Constraints on Primordial Gravitational Waves
Authors:
CMB-S4 Collaboration,
:,
Kevork Abazajian,
Graeme E. Addison,
Peter Adshead,
Zeeshan Ahmed,
Daniel Akerib,
Aamir Ali,
Steven W. Allen,
David Alonso,
Marcelo Alvarez,
Mustafa A. Amin,
Adam Anderson,
Kam S. Arnold,
Peter Ashton,
Carlo Baccigalupi,
Debbie Bard,
Denis Barkats,
Darcy Barron,
Peter S. Barry,
James G. Bartlett,
Ritoban Basu Thakur,
Nicholas Battaglia,
Rachel Bean,
Chris Bebek
, et al. (212 additional authors not shown)
Abstract:
CMB-S4---the next-generation ground-based cosmic microwave background (CMB) experiment---is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. Among the science cases pursued with CMB-S4, the quest for detecting p…
▽ More
CMB-S4---the next-generation ground-based cosmic microwave background (CMB) experiment---is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. Among the science cases pursued with CMB-S4, the quest for detecting primordial gravitational waves is a central driver of the experimental design. This work details the development of a forecasting framework that includes a power-spectrum-based semi-analytic projection tool, targeted explicitly towards optimizing constraints on the tensor-to-scalar ratio, $r$, in the presence of Galactic foregrounds and gravitational lensing of the CMB. This framework is unique in its direct use of information from the achieved performance of current Stage 2--3 CMB experiments to robustly forecast the science reach of upcoming CMB-polarization endeavors. The methodology allows for rapid iteration over experimental configurations and offers a flexible way to optimize the design of future experiments given a desired scientific goal. To form a closed-loop process, we couple this semi-analytic tool with map-based validation studies, which allow for the injection of additional complexity and verification of our forecasts with several independent analysis methods. We document multiple rounds of forecasts for CMB-S4 using this process and the resulting establishment of the current reference design of the primordial gravitational-wave component of the Stage-4 experiment, optimized to achieve our science goals of detecting primordial gravitational waves for $r > 0.003$ at greater than $5σ$, or, in the absence of a detection, of reaching an upper limit of $r < 0.001$ at $95\%$ CL.
△ Less
Submitted 27 August, 2020;
originally announced August 2020.
-
CosmoFlow: Using Deep Learning to Learn the Universe at Scale
Authors:
Amrita Mathuriya,
Deborah Bard,
Peter Mendygral,
Lawrence Meadows,
James Arnemann,
Lei Shao,
Siyu He,
Tuomas Karna,
Daina Moise,
Simon J. Pennycook,
Kristyn Maschoff,
Jason Sewall,
Nalini Kumar,
Shirley Ho,
Mike Ringenburg,
Prabhat,
Victor Lee
Abstract:
Deep learning is a promising tool to determine the physical model that describes our universe. To handle the considerable computational cost of this problem, we present CosmoFlow: a highly scalable deep learning application built on top of the TensorFlow framework. CosmoFlow uses efficient implementations of 3D convolution and pooling primitives, together with improvements in threading for many el…
▽ More
Deep learning is a promising tool to determine the physical model that describes our universe. To handle the considerable computational cost of this problem, we present CosmoFlow: a highly scalable deep learning application built on top of the TensorFlow framework. CosmoFlow uses efficient implementations of 3D convolution and pooling primitives, together with improvements in threading for many element-wise operations, to improve training performance on Intel(C) Xeon Phi(TM) processors. We also utilize the Cray PE Machine Learning Plugin for efficient scaling to multiple nodes. We demonstrate fully synchronous data-parallel training on 8192 nodes of Cori with 77% parallel efficiency, achieving 3.5 Pflop/s sustained performance. To our knowledge, this is the first large-scale science application of the TensorFlow framework at supercomputer scale with fully-synchronous training. These enhancements enable us to process large 3D dark matter distribution and predict the cosmological parameters $Ω_M$, $σ_8$ and n$_s$ with unprecedented accuracy.
△ Less
Submitted 9 November, 2018; v1 submitted 14 August, 2018;
originally announced August 2018.
-
Scaling GRPC Tensorflow on 512 nodes of Cori Supercomputer
Authors:
Amrita Mathuriya,
Thorsten Kurth,
Vivek Rane,
Mustafa Mustafa,
Lei Shao,
Debbie Bard,
Prabhat,
Victor W Lee
Abstract:
We explore scaling of the standard distributed Tensorflow with GRPC primitives on up to 512 Intel Xeon Phi (KNL) nodes of Cori supercomputer with synchronous stochastic gradient descent (SGD), and identify causes of scaling inefficiency at higher node counts. To our knowledge, this is the first exploration of distributed GRPC Tensorflow scalability on a HPC supercomputer at such large scale with s…
▽ More
We explore scaling of the standard distributed Tensorflow with GRPC primitives on up to 512 Intel Xeon Phi (KNL) nodes of Cori supercomputer with synchronous stochastic gradient descent (SGD), and identify causes of scaling inefficiency at higher node counts. To our knowledge, this is the first exploration of distributed GRPC Tensorflow scalability on a HPC supercomputer at such large scale with synchronous SGD. We studied scaling of two convolution neural networks - ResNet-50, a state-of-the-art deep network for classification with roughly 25.5 million parameters, and HEP-CNN, a shallow topology with less than 1 million parameters for common scientific usages. For ResNet-50, we achieve >80% scaling efficiency on up to 128 workers, using 32 parameter servers (PS tasks) with a steep decline down to 23% for 512 workers using 64 PS tasks. Our analysis of the efficiency drop points to low network bandwidth utilization due to combined effect of three factors. (a) Heterogeneous distributed parallelization algorithm which uses PS tasks as centralized servers for gradient averaging is suboptimal for utilizing interconnect bandwidth. (b) Load imbalance among PS tasks hinders their efficient scaling. (c) Underlying communication primitive GRPC is currently inefficient on Cori high-speed interconnect. The HEP-CNN demands less interconnect bandwidth, and shows >80% weak scaling efficiency for up to 256 nodes with only 1 PS task. Our findings are applicable to other deep learning networks. Big networks with millions of parameters stumble upon the issues discussed here. Shallower networks like HEP-CNN with relatively lower number of parameters can efficiently enjoy weak scaling even with a single parameter server.
△ Less
Submitted 26 December, 2017;
originally announced December 2017.
-
Galactos: Computing the Anisotropic 3-Point Correlation Function for 2 Billion Galaxies
Authors:
Brian Friesen,
Md. Mostofa Ali Patwary,
Brian Austin,
Nadathur Satish,
Zachary Slepian,
Narayanan Sundaram,
Deborah Bard,
Daniel J Eisenstein,
Jack Deslippe,
Pradeep Dubey,
Prabhat
Abstract:
The nature of dark energy and the complete theory of gravity are two central questions currently facing cosmology. A vital tool for addressing them is the 3-point correlation function (3PCF), which probes deviations from a spatially random distribution of galaxies. However, the 3PCF's formidable computational expense has prevented its application to astronomical surveys comprising millions to bill…
▽ More
The nature of dark energy and the complete theory of gravity are two central questions currently facing cosmology. A vital tool for addressing them is the 3-point correlation function (3PCF), which probes deviations from a spatially random distribution of galaxies. However, the 3PCF's formidable computational expense has prevented its application to astronomical surveys comprising millions to billions of galaxies. We present Galactos, a high-performance implementation of a novel, O(N^2) algorithm that uses a load-balanced k-d tree and spherical harmonic expansions to compute the anisotropic 3PCF. Our implementation is optimized for the Intel Xeon Phi architecture, exploiting SIMD parallelism, instruction and thread concurrency, and significant L1 and L2 cache reuse, reaching 39% of peak performance on a single node. Galactos scales to the full Cori system, achieving 9.8PF (peak) and 5.06PF (sustained) across 9636 nodes, making the 3PCF easily computable for all galaxies in the observable universe.
△ Less
Submitted 31 August, 2017;
originally announced September 2017.
-
CosmoGAN: creating high-fidelity weak lensing convergence maps using Generative Adversarial Networks
Authors:
Mustafa Mustafa,
Deborah Bard,
Wahid Bhimji,
Zarija Lukić,
Rami Al-Rfou,
Jan M. Kratochvil
Abstract:
Inferring model parameters from experimental data is a grand challenge in many sciences, including cosmology. This often relies critically on high fidelity numerical simulations, which are prohibitively computationally expensive. The application of deep learning techniques to generative modeling is renewing interest in using high dimensional density estimators as computationally inexpensive emulat…
▽ More
Inferring model parameters from experimental data is a grand challenge in many sciences, including cosmology. This often relies critically on high fidelity numerical simulations, which are prohibitively computationally expensive. The application of deep learning techniques to generative modeling is renewing interest in using high dimensional density estimators as computationally inexpensive emulators of fully-fledged simulations. These generative models have the potential to make a dramatic shift in the field of scientific simulations, but for that shift to happen we need to study the performance of such generators in the precision regime needed for science applications. To this end, in this work we apply Generative Adversarial Networks to the problem of generating weak lensing convergence maps. We show that our generator network produces maps that are described by, with high statistical confidence, the same summary statistics as the fully simulated maps.
△ Less
Submitted 22 May, 2019; v1 submitted 7 June, 2017;
originally announced June 2017.
-
Probabilistic Cosmological Mass Mapping from Weak Lensing Shear
Authors:
Michael D. Schneider,
Karen Y. Ng,
William A. Dawson,
Philip J. Marshall,
Joshua Meyers,
Deborah J. Bard
Abstract:
We infer gravitational lensing shear and convergence fields from galaxy ellipticity catalogs under a spatial process prior for the lensing potential. We demonstrate the performance of our algorithm with simulated Gaussian-distributed cosmological lensing shear maps and a reconstruction of the mass distribution of the merging galaxy cluster Abell 781 using galaxy ellipticities measured with the Dee…
▽ More
We infer gravitational lensing shear and convergence fields from galaxy ellipticity catalogs under a spatial process prior for the lensing potential. We demonstrate the performance of our algorithm with simulated Gaussian-distributed cosmological lensing shear maps and a reconstruction of the mass distribution of the merging galaxy cluster Abell 781 using galaxy ellipticities measured with the Deep Lens Survey. Given interim posterior samples of lensing shear or convergence fields on the sky, we describe an algorithm to infer cosmological parameters via lens field marginalization. In the most general formulation of our algorithm we make no assumptions about weak shear or Gaussian distributed shape noise or shears. Because we require solutions and matrix determinants of a linear system of dimension that scales with the number of galaxies, we expect our algorithm to require parallel high-performance computing resources for application to ongoing wide field lensing surveys.
△ Less
Submitted 21 October, 2016;
originally announced October 2016.
-
ASCR/HEP Exascale Requirements Review Report
Authors:
Salman Habib,
Robert Roser,
Richard Gerber,
Katie Antypas,
Katherine Riley,
Tim Williams,
Jack Wells,
Tjerk Straatsma,
A. Almgren,
J. Amundson,
S. Bailey,
D. Bard,
K. Bloom,
B. Bockelman,
A. Borgland,
J. Borrill,
R. Boughezal,
R. Brower,
B. Cowan,
H. Finkel,
N. Frontiere,
S. Fuess,
L. Ge,
N. Gnedin,
S. Gottlieb
, et al. (29 additional authors not shown)
Abstract:
This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 ti…
▽ More
This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems.
△ Less
Submitted 31 March, 2016; v1 submitted 30 March, 2016;
originally announced March 2016.
-
First Observation of CP Violation in B0->D(*)CP h0 Decays by a Combined Time-Dependent Analysis of BaBar and Belle Data
Authors:
The BaBar,
Belle Collaborations,
:,
A. Abdesselam,
I. Adachi,
A. Adametz,
T. Adye,
H. Ahmed,
H. Aihara,
S. Akar,
M. S. Alam,
J. Albert,
S. Al Said,
R. Andreassen,
C. Angelini,
F. Anulli,
K. Arinstein,
N. Arnaud,
D. M. Asner,
D. Aston,
V. Aulchenko,
T. Aushev,
R. Ayad,
V. Babu,
I. Badhrees
, et al. (450 additional authors not shown)
Abstract:
We report a measurement of the time-dependent CP asymmetry of B0->D(*)CP h0 decays, where the light neutral hadron h0 is a pi0, eta or omega meson, and the neutral D meson is reconstructed in the CP eigenstates K+ K-, K0S pi0 or K0S omega. The measurement is performed combining the final data samples collected at the Y(4S) resonance by the BaBar and Belle experiments at the asymmetric-energy B fac…
▽ More
We report a measurement of the time-dependent CP asymmetry of B0->D(*)CP h0 decays, where the light neutral hadron h0 is a pi0, eta or omega meson, and the neutral D meson is reconstructed in the CP eigenstates K+ K-, K0S pi0 or K0S omega. The measurement is performed combining the final data samples collected at the Y(4S) resonance by the BaBar and Belle experiments at the asymmetric-energy B factories PEP-II at SLAC and KEKB at KEK, respectively. The data samples contain ( 471 +/- 3 ) x 10^6 BB pairs recorded by the BaBar detector and ( 772 +/- 11 ) x 10^6, BB pairs recorded by the Belle detector. We measure the CP asymmetry parameters -eta_f S = +0.66 +/- 0.10 (stat.) +/- 0.06 (syst.) and C = -0.02 +/- 0.07 (stat.) +/- 0.03 (syst.). These results correspond to the first observation of CP violation in B0->D(*)CP h0 decays. The hypothesis of no mixing-induced CP violation is excluded in these decays at the level of 5.4 standard deviations.
△ Less
Submitted 20 July, 2015; v1 submitted 15 May, 2015;
originally announced May 2015.
-
GREAT3 results I: systematic errors in shear estimation and the impact of real galaxy morphology
Authors:
Rachel Mandelbaum,
Barnaby Rowe,
Robert Armstrong,
Deborah Bard,
Emmanuel Bertin,
James Bosch,
Dominique Boutigny,
Frederic Courbin,
William A. Dawson,
Annamaria Donnarumma,
Ian Fenech Conti,
Raphael Gavazzi,
Marc Gentile,
Mandeep S. S. Gill,
David W. Hogg,
Eric M. Huff,
M. James Jee,
Tomasz Kacprzak,
Martin Kilbinger,
Thibault Kuntzer,
Dustin Lang,
Wentao Luo,
Marisa C. March,
Philip J. Marshall,
Joshua E. Meyers
, et al. (18 additional authors not shown)
Abstract:
We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-v…
▽ More
We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the Sérsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.
△ Less
Submitted 30 December, 2022; v1 submitted 4 December, 2014;
originally announced December 2014.
-
Hierarchical probabilistic inference of cosmic shear
Authors:
Michael D. Schneider,
David W. Hogg,
Philip J. Marshall,
William A. Dawson,
Joshua Meyers,
Deborah J. Bard,
Dustin Lang
Abstract:
Point estimators for the shearing of galaxy images induced by gravitational lensing involve a complex inverse problem in the presence of noise, pixelization, and model uncertainties. We present a probabilistic forward modeling approach to gravitational lensing inference that has the potential to mitigate the biased inferences in most common point estimators and is practical for upcoming lensing su…
▽ More
Point estimators for the shearing of galaxy images induced by gravitational lensing involve a complex inverse problem in the presence of noise, pixelization, and model uncertainties. We present a probabilistic forward modeling approach to gravitational lensing inference that has the potential to mitigate the biased inferences in most common point estimators and is practical for upcoming lensing surveys. The first part of our statistical framework requires specification of a likelihood function for the pixel data in an imaging survey given parameterized models for the galaxies in the images. We derive the lensing shear posterior by marginalizing over all intrinsic galaxy properties that contribute to the pixel data (i.e., not limited to galaxy ellipticities) and learn the distributions for the intrinsic galaxy properties via hierarchical inference with a suitably flexible conditional probabilitiy distribution specification. We use importance sampling to separate the modeling of small imaging areas from the global shear inference, thereby rendering our algorithm computationally tractable for large surveys. With simple numerical examples we demonstrate the improvements in accuracy from our importance sampling approach, as well as the significance of the conditional distribution specification for the intrinsic galaxy properties when the data are generated from an unknown number of distinct galaxy populations with different morphological characteristics.
△ Less
Submitted 10 November, 2014;
originally announced November 2014.
-
Masked areas in shear peak statistics: a forward modeling approach
Authors:
Deborah Bard,
Jan M. Kratochvil,
William Dawson
Abstract:
The statistics of shear peaks have been shown to provide valuable cosmological information beyond the power spectrum, and will be an important constraint of models of cosmology with the large survey areas provided by forthcoming astronomical surveys. Surveys include masked areas due to bright stars, bad pixels etc, which must be accounted for in producing constraints on cosmology from shear maps.…
▽ More
The statistics of shear peaks have been shown to provide valuable cosmological information beyond the power spectrum, and will be an important constraint of models of cosmology with the large survey areas provided by forthcoming astronomical surveys. Surveys include masked areas due to bright stars, bad pixels etc, which must be accounted for in producing constraints on cosmology from shear maps. We advocate a forward-modeling approach, where the impact of masking (and other survey artifacts) are accounted for in the theoretical prediction of cosmological parameters, rather than removed from survey data. We use masks based on the Deep Lens Survey, and explore the impact of up to 37% of the survey area being masked on LSST and DES-scale surveys. By reconstructing maps of aperture mass, the masking effect is smoothed out, resulting in up to 14% smaller statistical uncertainties compared to simply reducing the survey area by the masked area. We show that, even in the presence of large survey masks, the bias in cosmological parameter estimation produced in the forward-modeling process is ~1%, dominated by bias caused by limited simulation volume. We also explore how this potential bias scales with survey area and find that small survey areas are more significantly impacted by the differences in cosmological structure in the data and simulated volumes, due to cosmic variance.
△ Less
Submitted 20 October, 2014;
originally announced October 2014.
-
The Physics of the B Factories
Authors:
A. J. Bevan,
B. Golob,
Th. Mannel,
S. Prell,
B. D. Yabsley,
K. Abe,
H. Aihara,
F. Anulli,
N. Arnaud,
T. Aushev,
M. Beneke,
J. Beringer,
F. Bianchi,
I. I. Bigi,
M. Bona,
N. Brambilla,
J. B rodzicka,
P. Chang,
M. J. Charles,
C. H. Cheng,
H. -Y. Cheng,
R. Chistov,
P. Colangelo,
J. P. Coleman,
A. Drutskoy
, et al. (2009 additional authors not shown)
Abstract:
This work is on the Physics of the B Factories. Part A of this book contains a brief description of the SLAC and KEK B Factories as well as their detectors, BaBar and Belle, and data taking related issues. Part B discusses tools and methods used by the experiments in order to obtain results. The results themselves can be found in Part C.
Please note that version 3 on the archive is the auxiliary…
▽ More
This work is on the Physics of the B Factories. Part A of this book contains a brief description of the SLAC and KEK B Factories as well as their detectors, BaBar and Belle, and data taking related issues. Part B discusses tools and methods used by the experiments in order to obtain results. The results themselves can be found in Part C.
Please note that version 3 on the archive is the auxiliary version of the Physics of the B Factories book. This uses the notation alpha, beta, gamma for the angles of the Unitarity Triangle. The nominal version uses the notation phi_1, phi_2 and phi_3. Please cite this work as Eur. Phys. J. C74 (2014) 3026.
△ Less
Submitted 31 October, 2015; v1 submitted 24 June, 2014;
originally announced June 2014.
-
Evidence for the decay B0 --> omega omega and search for B0 --> omega phi
Authors:
The BABAR Collaboration,
J. P. Lees,
V. Poireau,
V. Tisserand,
E. Grauges,
A. Palano,
G. Eigen,
B. Stugu,
D. N. Brown,
L. T. Kerth,
Yu. G. Kolomensky,
M. J. Lee,
G. Lynch,
H. Koch,
T. Schroeder,
C. Hearty,
T. S. Mattison,
J. A. McKenna,
R. Y. So,
A. Khan,
V. E. Blinov,
A. R. Buzykaev,
V. P. Druzhinin,
V. B. Golubev,
E. A. Kravchenko
, et al. (312 additional authors not shown)
Abstract:
We describe searches for B meson decays to the charmless vector-vector final states omega omega and omega phi with 471 x 10^6 B Bbar pairs produced in e+ e- annihilation at sqrt(s) = 10.58 GeV using the BABAR detector at the PEP-II collider at the SLAC National Accelerator Laboratory. We measure the branching fraction B(B0 --> omega omega) = (1.2 +- 0.3 +0.3-0.2) x 10^-6, where the first uncertain…
▽ More
We describe searches for B meson decays to the charmless vector-vector final states omega omega and omega phi with 471 x 10^6 B Bbar pairs produced in e+ e- annihilation at sqrt(s) = 10.58 GeV using the BABAR detector at the PEP-II collider at the SLAC National Accelerator Laboratory. We measure the branching fraction B(B0 --> omega omega) = (1.2 +- 0.3 +0.3-0.2) x 10^-6, where the first uncertainty is statistical and the second is systematic, corresponding to a significance of 4.4 standard deviations. We also determine the upper limit B(B0 --> omega phi) < 0.7 x 10^-6 at 90% confidence level. These measurements provide the first evidence for the decay B0 --> omega omega, and an improvement of the upper limit for the decay B0 --> omega phi.
△ Less
Submitted 6 March, 2014; v1 submitted 29 November, 2013;
originally announced December 2013.
-
Facilities for Dark Energy Investigations
Authors:
David Weinberg,
Deborah Bard,
Kyle Dawson,
Olivier Dore,
Joshua Frieman,
Karl Gebhardt,
Michael Levi,
Jason Rhodes
Abstract:
The discovery of cosmic acceleration has inspired ambitious experimental and observational efforts to understand its origin. Many of these take the form of large astronomical surveys, sometimes using new, special-purpose instrumentation, and in some cases entirely new facilities. This Report summarizes some of the major ongoing and planned dark energy experiments, focusing on those in which the U.…
▽ More
The discovery of cosmic acceleration has inspired ambitious experimental and observational efforts to understand its origin. Many of these take the form of large astronomical surveys, sometimes using new, special-purpose instrumentation, and in some cases entirely new facilities. This Report summarizes some of the major ongoing and planned dark energy experiments, focusing on those in which the U.S. community has a leading or significant supporting role. It provides background information for other Reports from the "Dark Energy and CMB" working group of the APS Division of Particles and Fields long-term planning exercise. We provide 1-2 page summaries of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) and its SDSS-IV successor eBOSS, the Dark Energy Survey (DES), the Subaru Hyper-Suprime Camera (HSC) and Prime Focus Spectrograph (PFS), the Hobby-Eberley Telescope Dark Energy Experiment (HEDTEX), the Dark Energy Spectroscopic Instrument (DESI), the Large Synpotic Survey Telescope (LSST), and the Euclid and AFTA-WFIRST (Wide Field Infrared Survey Telescope) space missions. Over the next one to two decades, these experiments seek to improve the precision of measurements of the cosmic expansion history and cosmic structure growth by factors of 10-100, while simultaneously increasing their redshift range and tightening control and cross-checks of sytematic effects. These measurements will dramatically improve the empirical constraints on theoretical explanations of cosmic acceleration and could lead to breakthroughs in our understanding of the nature of gravity and spacetime or the forms of energy present in the universe.
△ Less
Submitted 20 September, 2013;
originally announced September 2013.
-
Effect of Measurement Errors on Predicted Cosmological Constraints from Shear Peak Statistics with LSST
Authors:
D. Bard,
J. M. Kratochvil,
C. Chang,
M. May,
S. M. Kahn,
Y. AlSayyad,
Z. Ahmad,
J. Bankert,
A. Connolly,
R. R. Gibson,
K. Gilmore,
E. Grace,
Z. Haiman,
M. Hannel,
K. M. Huffenberger,
J. G. Jernigan,
L. Jones,
S. Krughoff,
S. Lorenz,
S. Marshall,
A. Meert,
S. Nagarajan,
E. Peng,
J. Peterson,
A. P. Rasmussen
, et al. (4 additional authors not shown)
Abstract:
The statistics of peak counts in reconstructed shear maps contain information beyond the power spectrum, and can improve cosmological constraints from measurements of the power spectrum alone if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic Survey Tel…
▽ More
The statistics of peak counts in reconstructed shear maps contain information beyond the power spectrum, and can improve cosmological constraints from measurements of the power spectrum alone if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic Survey Telescope (LSST). We use the LSST image simulator in combination with cosmological N-body simulations to model realistic shear maps for different cosmological models. We include both galaxy shape noise and, for the first time, measurement errors on galaxy shapes. We find that the measurement errors considered have relatively little impact on the constraining power of shear peak counts for LSST.
△ Less
Submitted 4 January, 2013;
originally announced January 2013.
-
Search for direct CP-violation in singly-Cabibbo suppressed D+- --> K+ K- pi+- decays
Authors:
BaBar Collaboration,
J. P. Lees,
V. Poireau,
V. Tisserand,
J. Garra Tico,
E. Grauges,
A. Palano,
G. Eigen,
B. Stugu,
D. N. Brown,
L. T. Kerth,
Yu. G. Kolomensky,
G. Lynch,
H. Koch,
T. Schroeder,
D. J. Asgeirsson,
C. Hearty,
T. S. Mattison,
J. A. McKenna,
R. Y. So,
A. Khan,
V. E. Blinov,
A. R. Buzykaev,
V. P. Druzhinin,
V. B. Golubev
, et al. (338 additional authors not shown)
Abstract:
We report on a search for direct CP asymmetry in the singly Cabibbo-suppressed decay D+- --> K+ K- pi+- using a data sample of 476 fb-1 accumulated with the BaBar detector running at and just below the Y(4S) resonance. The CP-violating decay rate asymmetry A_CP is determined to be (0.35 +- 0.30 +- 0.15)%. Model-dependent and model-independent Dalitz plot analysis techniques are used to search for…
▽ More
We report on a search for direct CP asymmetry in the singly Cabibbo-suppressed decay D+- --> K+ K- pi+- using a data sample of 476 fb-1 accumulated with the BaBar detector running at and just below the Y(4S) resonance. The CP-violating decay rate asymmetry A_CP is determined to be (0.35 +- 0.30 +- 0.15)%. Model-dependent and model-independent Dalitz plot analysis techniques are used to search for CP-violating asymmetries in the various intermediate states.
△ Less
Submitted 21 February, 2013; v1 submitted 8 December, 2012;
originally announced December 2012.
-
Cosmological Calculations on the GPU
Authors:
Deborah Bard,
Matthew Bellis,
Mark T. Allen,
Hasmik Yepremyan,
Jan M. Kratochvil
Abstract:
Cosmological measurements require the calculation of nontrivial quantities over large datasets. The next generation of survey telescopes (such as DES, PanSTARRS, and LSST) will yield measurements of billions of galaxies. The scale of these datasets, and the nature of the calculations involved, make cosmological calculations ideal models for implementation on graphics processing units (GPUs). We co…
▽ More
Cosmological measurements require the calculation of nontrivial quantities over large datasets. The next generation of survey telescopes (such as DES, PanSTARRS, and LSST) will yield measurements of billions of galaxies. The scale of these datasets, and the nature of the calculations involved, make cosmological calculations ideal models for implementation on graphics processing units (GPUs). We consider two cosmological calculations, the two-point angular correlation function and the aperture mass statistic, and aim to improve the calculation time by constructing code for calculating them on the GPU. Using CUDA, we implement the two algorithms on the GPU and compare the calculation speeds to comparable code run on the CPU. We obtain a code speed-up of between 10 - 180x faster, compared to performing the same calculation on the CPU. The code has been made publicly available.
△ Less
Submitted 6 December, 2012; v1 submitted 17 August, 2012;
originally announced August 2012.
-
Branching fraction and form-factor shape measurements of exclusive charmless semileptonic B decays, and determination of |V_{ub}|
Authors:
BABAR Collaboration,
J. P. Lees,
V. Poireau,
V. Tisserand,
J. Garra Tico,
E. Grauges,
A. Palano,
G. Eigen,
B. Stugu,
D. N. Brown,
L. T. Kerth,
Yu. G. Kolomensky,
G. Lynch,
H. Koch,
T. Schroeder,
D. J. Asgeirsson,
C. Hearty,
T. S. Mattison,
J. A. McKenna,
R. Y. So,
A. Khan,
V. E. Blinov,
A. R. Buzykaev,
V. P. Druzhinin,
V. B. Golubev
, et al. (340 additional authors not shown)
Abstract:
We report the results of a study of the exclusive charmless semileptonic decays, B^0 --> pi^- l^+ nu, B^+ --> pi^0 l^+ nu, B^+ --> omega l^+ nu, B^+ --> eta l^+ nu and B^+ --> eta^' l^+ nu, (l = e or mu) undertaken with approximately 462x10^6 B\bar{B} pairs collected at the Upsilon(4S) resonance with the BABAR detector. The analysis uses events in which the signal B decays are reconstructed with a…
▽ More
We report the results of a study of the exclusive charmless semileptonic decays, B^0 --> pi^- l^+ nu, B^+ --> pi^0 l^+ nu, B^+ --> omega l^+ nu, B^+ --> eta l^+ nu and B^+ --> eta^' l^+ nu, (l = e or mu) undertaken with approximately 462x10^6 B\bar{B} pairs collected at the Upsilon(4S) resonance with the BABAR detector. The analysis uses events in which the signal B decays are reconstructed with a loose neutrino reconstruction technique. We obtain partial branching fractions in several bins of q^2, the square of the momentum transferred to the lepton-neutrino pair, for B^0 --> pi^- l^+ nu, B^+ --> pi^0 l^+ nu, B^+ --> omega l^+ nu and B^+ --> eta l^+ nu. From these distributions, we extract the form-factor shapes f_+(q^2) and the total branching fractions BF(B^0 --> pi^- l^+ nu) = (1.45 +/- 0.04_{stat} +/- 0.06_{syst})x10^-4 (combined pi^- and pi^0 decay channels assuming isospin symmetry), BF(B^+ --> omega l^+ nu) = (1.19 +/- 0.16_{stat} +/- 0.09_{syst})x10^-4 and BF(B^+ --> eta l^+ nu) = (0.38 +/- 0.05_{stat} +/- 0.05_{syst})x10^-4. We also measure BF(B^+ --> eta^' l^+ nu) = (0.24 +/- 0.08_{stat} +/- 0.03_{syst})x10^-4. We obtain values for the magnitude of the CKM matrix element V_{ub} by direct comparison with three different QCD calculations in restricted q^2 ranges of B --> pi l^+ nu decays. From a simultaneous fit to the experimental data over the full q^2 range and the FNAL/MILC lattice QCD predictions, we obtain |V_{ub}| = (3.25 +/- 0.31)x10^-3, where the error is the combined experimental and theoretical uncertainty.
△ Less
Submitted 6 November, 2012; v1 submitted 6 August, 2012;
originally announced August 2012.
-
Observation of Time Reversal Violation in the B0 Meson System
Authors:
The BABAR Collaboration,
J. P. Lees,
V. Poireau,
V. Tisserand,
J. Garra Tico,
E. Grauges,
A. Palanoab,
G. Eigen,
B. Stugu,
D. N. Brown,
L. T. Kerth,
Yu. G. Kolomensky,
G. Lynch,
H. Koch,
T. Schroeder,
D. J. Asgeirsson,
C. Hearty,
T. S. Mattison,
J. A. McKenna,
R. Y. So,
A. Khan,
V. E. Blinov,
A. R. Buzykaev,
V. P. Druzhinin,
V. B. Golubev
, et al. (343 additional authors not shown)
Abstract:
Although CP violation in the B meson system has been well established by the B factories, there has been no direct observation of time reversal violation. The decays of entangled neutral B mesons into definite flavor states ($B^0$ or $\bar{B}^0$), and $J/ψK_S^0$ or $c\bar{c} K_S^0$ final states (referred to as $B_+$ or $B_-$), allow comparisons between the probabilities of four pairs of T-conjugat…
▽ More
Although CP violation in the B meson system has been well established by the B factories, there has been no direct observation of time reversal violation. The decays of entangled neutral B mesons into definite flavor states ($B^0$ or $\bar{B}^0$), and $J/ψK_S^0$ or $c\bar{c} K_S^0$ final states (referred to as $B_+$ or $B_-$), allow comparisons between the probabilities of four pairs of T-conjugated transitions, for example, $\bar{B}^0 \rightarrow B_-$ and $B_- \rightarrow \bar{B}^0$, as a function of the time difference between the two B decays. Using 468 million $B\bar{B}$ pairs produced in $Υ(4S)$ decays collected by the BABAR detector at SLAC, we measure T-violating parameters in the time evolution of neutral B mesons, yielding $ΔS_T^+ = -1.37 \pm 0.14 (stat.) \pm 0.06 (syst.)$ and $ΔS_T^- = 1.17 \pm 0.18 (stat.) \pm 0.11 (syst.)$. These nonzero results represent the first direct observation of T violation through the exchange of initial and final states in transitions that can only be connected by a T-symmetry transformation.
△ Less
Submitted 2 January, 2013; v1 submitted 24 July, 2012;
originally announced July 2012.
-
Atmospheric PSF Interpolation for Weak Lensing in Short Exposure Imaging Data
Authors:
C. Chang,
P. J. Marshall,
J. G. Jernigan,
J. R. Peterson,
S. M. Kahn,
S. F. Gull,
Y. AlSayyad,
Z. Ahmad,
J. Bankert,
D. Bard,
A. Connolly,
R. R. Gibson,
K. Gilmore,
E. Grace,
M. Hannel,
M. A. Hodge,
L. Jones,
S. Krughoff,
S. Lorenz,
S. Marshall,
A. Meert,
S. Nagarajan,
E. Peng,
A. P. Rasmussen,
M. Shmakova
, et al. (3 additional authors not shown)
Abstract:
A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause…
▽ More
A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause the PSF to vary on angular scales similar or larger than a single CCD sensor, the atmosphere generates stochastic structures on a wide range of angular scales. It thus becomes a challenge to infer the multi-scale, complex atmospheric PSF patterns by interpolating the sparsely sampled stars in the field. In this paper we present a new method, PSFent, for interpolating the PSF shape parameters, based on reconstructing underlying shape parameter maps with a multi-scale maximum entropy algorithm. We demonstrate, using images from the LSST Photon Simulator, the performance of our approach relative to a 5th-order polynomial fit (representing the current standard) and a simple boxcar smoothing technique. Quantitatively, PSFent predicts more accurate PSF models in all scenarios and the residual PSF errors are spatially less correlated. This improvement in PSF interpolation leads to a factor of 3.5 lower systematic errors in the shear power spectrum on scales smaller than $\sim13'$, compared to polynomial fitting. We estimate that with PSFent and for stellar densities greater than $\simeq1/{\rm arcmin}^{2}$, the spurious shear correlation from PSF interpolation, after combining a complete 10-year dataset from LSST, is lower than the corresponding statistical uncertainties on the cosmic shear power spectrum, even under a conservative scenario.
△ Less
Submitted 12 November, 2012; v1 submitted 6 June, 2012;
originally announced June 2012.
-
Spurious Shear in Weak Lensing with LSST
Authors:
C. Chang,
S. M. Kahn,
J. G. Jernigan,
J. R. Peterson,
Y. AlSayyad,
Z. Ahmad,
J. Bankert,
D. Bard,
A. Connolly,
R. R. Gibson,
K. Gilmore,
E. Grace,
M. Hannel,
M. A. Hodge,
M. J. Jee,
L. Jones,
S. Krughoff,
S. Lorenz,
P. J. Marshall,
S. Marshall,
A. Meert,
S. Nagarajan,
E. Peng,
A. P. Rasmussen,
M. Shmakova
, et al. (3 additional authors not shown)
Abstract:
The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image $\sim$ 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to $r\sim27.5$, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controll…
▽ More
The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image $\sim$ 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to $r\sim27.5$, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controlled to a level similar to the statistical errors.
This work is the first attempt to quantitatively estimate the absolute level and statistical properties of the systematic errors on weak lensing shear measurements due to the most important physical effects in the LSST system via high fidelity ray-tracing simulations. We identify and isolate the different sources of algorithm-independent, \textit{additive} systematic errors on shear measurements for LSST and predict their impact on the final cosmic shear measurements using conventional weak lensing analysis techniques. We find that the main source of the errors comes from an inability to adequately characterise the atmospheric point spread function (PSF) due to its high frequency spatial variation on angular scales smaller than $\sim10'$ in the single short exposures, which propagates into a spurious shear correlation function at the $10^{-4}$--$10^{-3}$ level on these scales. With the large multi-epoch dataset that will be acquired by LSST, the stochastic errors average out, bringing the final spurious shear correlation function to a level very close to the statistical errors. Our results imply that the cosmological constraints from LSST will not be severely limited by these algorithm-independent, additive systematic effects.
△ Less
Submitted 16 October, 2012; v1 submitted 6 June, 2012;
originally announced June 2012.
-
Initial-State Radiation Measurement of the e+e- -> pi+pi-pi+pi- Cross Section
Authors:
J. P. Lees,
V. Poireau,
V. Tisserand,
J. Garra Tico,
E. Grauges,
M. Martinelli,
D. A. Milanes,
A. Palano,
M. Pappagallo,
G. Eigen,
B. Stugu,
D. N. Brown,
L. T. Kerth,
Yu. G. Kolomensky,
G. Lynch,
H. Koch,
T. Schroeder,
D. J. Asgeirsson,
C. Hearty,
T. S. Mattison,
J. A. McKenna,
A. Khan,
V. E. Blinov,
A. R. Buzykaev,
V. P. Druzhinin
, et al. (362 additional authors not shown)
Abstract:
We study the process e+e- -> pi+pi-pi+pi-gamma, with a photon emitted from the initial-state electron or positron, using 454.3 fb^-1 of data collected with the BABAR detector at SLAC, corresponding to approximately 260,000 signal events. We use these data to extract the non-radiative sigma(e+e- ->pi+pi-pi+pi-) cross section in the energy range from 0.6 to 4.5 Gev. The total uncertainty of the cros…
▽ More
We study the process e+e- -> pi+pi-pi+pi-gamma, with a photon emitted from the initial-state electron or positron, using 454.3 fb^-1 of data collected with the BABAR detector at SLAC, corresponding to approximately 260,000 signal events. We use these data to extract the non-radiative sigma(e+e- ->pi+pi-pi+pi-) cross section in the energy range from 0.6 to 4.5 Gev. The total uncertainty of the cross section measurement in the peak region is less than 3%, higher in precision than the corresponding results obtained from energy scan data.
△ Less
Submitted 27 August, 2012; v1 submitted 26 January, 2012;
originally announced January 2012.
-
A Measurement of the Semileptonic Branching Fraction of the B_s Meson
Authors:
The BaBar Collaboration,
J. P. Lees,
V. Poireau,
V. Tisserand,
J. Garra Tico,
E. Grauges,
M. Martinelli,
D. A. Milanes,
A. Palano,
M. Pappagallo,
G. Eigen,
B. Stugu,
D. N. Brown,
L. T. Kerth,
Yu. G. Kolomensky,
G. Lynch,
H. Koch,
T. Schroeder,
D. J. Asgeirsson,
C. Hearty,
T. S. Mattison,
J. A. McKenna,
A. Khan,
V. E. Blinov,
A. R. Buzykaev
, et al. (362 additional authors not shown)
Abstract:
We report a measurement of the inclusive semileptonic branching fraction of the B_s meson using data collected with the BaBar detector in the center-of-mass (CM) energy region above the Upsilon(4S) resonance. We use the inclusive yield of phi mesons and the phi yield in association with a high-momentum lepton to perform a simultaneous measurement of the semileptonic branching fraction and the pr…
▽ More
We report a measurement of the inclusive semileptonic branching fraction of the B_s meson using data collected with the BaBar detector in the center-of-mass (CM) energy region above the Upsilon(4S) resonance. We use the inclusive yield of phi mesons and the phi yield in association with a high-momentum lepton to perform a simultaneous measurement of the semileptonic branching fraction and the production rate of B_s mesons relative to all B mesons as a function of CM energy. The inclusive semileptonic branching fraction of the B_s meson is determined to be B(B_s to l nu X)=9.5 (+2.5/-2.0)(stat)(+1.1/-1.9)(syst)%, where l indicates the average of e and mu.
△ Less
Submitted 10 January, 2012; v1 submitted 25 October, 2011;
originally announced October 2011.
-
Search for hadronic decays of a light Higgs boson in the radiative decay Upsilon --> gamma A0
Authors:
The BABAR Collaboration,
J. P. Lees,
V. Poireau,
V. Tisserand,
J. Garra Tico,
E. Grauges,
M. Martinelli,
D. A. Milanes,
A. Palano,
M. Pappagallo,
G. Eigen,
B. Stugu,
D. N. Brown,
L. T. Kerth,
Yu. G. Kolomensky,
G. Lynch,
H. Koch,
T. Schroeder,
D. J. Asgeirsson,
C. Hearty,
T. S. Mattison,
J. A. McKenna,
R. Y. So,
A. Khan,
V. E. Blinov
, et al. (364 additional authors not shown)
Abstract:
We search for hadronic decays of a light Higgs boson (A0) produced in radiative decays of an Upsilon(2S) or Upsilon(3S) meson, Upsilon --> gamma A0. The data have been recorded by the BABAR experiment at the Upsilon(3S) and Upsilon(2S) center of mass energies, and include (121.3 \pm 1.2) x 10^6 Upsilon(3S) and (98.3 \pm 0.9) x 10^6 Upsilon(2S) mesons. No significant signal is observed. We set 90%…
▽ More
We search for hadronic decays of a light Higgs boson (A0) produced in radiative decays of an Upsilon(2S) or Upsilon(3S) meson, Upsilon --> gamma A0. The data have been recorded by the BABAR experiment at the Upsilon(3S) and Upsilon(2S) center of mass energies, and include (121.3 \pm 1.2) x 10^6 Upsilon(3S) and (98.3 \pm 0.9) x 10^6 Upsilon(2S) mesons. No significant signal is observed. We set 90% confidence level upper limits on the product branching fractions B(Upsilon(nS)-->gamma A0) x B(A0-->hadrons) (n=2 or 3) that range from 1 x 10^{-6} for an A0 mass of 0.3 GeV/c^2 to 8 x 10^{-5} at 7 GeV/c^2.
△ Less
Submitted 17 August, 2011;
originally announced August 2011.
-
Observation of the baryonic B decay B0bar --> Lambda_c^+ anti-Lambda K-
Authors:
The BABAR Collaboration,
J. P. Lees,
V. Poireau,
V. Tisserand,
J. Garra Tico,
E. Grauges,
M. Martinelli,
D. A. Milanes,
A. Palano,
M. Pappagallo,
G. Eigen,
B. Stugu,
L. Sun,
D. N. Brown,
L. T. Kerth,
Yu. G. Kolomensky,
G. Lynch,
H. Koch,
T. Schroeder,
D. J. Asgeirsson,
C. Hearty,
T. S. Mattison,
J. A. McKenna,
A. Khan,
V. E. Blinov
, et al. (363 additional authors not shown)
Abstract:
We report the observation of the baryonic B decay B0bar --> Lambda_c^+ anti-Lambda K- with a significance larger than 7 standard deviations based on 471x10^6$ BBbar pairs collected with the BABAR detector at the PEP-II storage ring at SLAC. We measure the branching fraction for the decay B0bar --> Lambda_c^+ anti-Lambda K- to be (3.8 \pm 0.8_{stat} \pm 0.2_{sys} \pm 1.0_{Lambda_c^+})x10^{-5}. The…
▽ More
We report the observation of the baryonic B decay B0bar --> Lambda_c^+ anti-Lambda K- with a significance larger than 7 standard deviations based on 471x10^6$ BBbar pairs collected with the BABAR detector at the PEP-II storage ring at SLAC. We measure the branching fraction for the decay B0bar --> Lambda_c^+ anti-Lambda K- to be (3.8 \pm 0.8_{stat} \pm 0.2_{sys} \pm 1.0_{Lambda_c^+})x10^{-5}. The uncertainties are statistical, systematic, and due to the uncertainty in the Lambda_c^+ branching fraction. We find that the Lambda_c^+ K^- invariant mass distribution shows an enhancement above 3.5 GeV/c^2.
△ Less
Submitted 16 August, 2011;
originally announced August 2011.
-
Study of radiative bottomonium transitions using converted photons
Authors:
The BABAR Collaboration,
J. P. Lees,
V. Poireau,
E. Prencipe,
V. Tisserand,
J. GarraTico,
E. Grauges,
M. Martinelli,
D. A. Milanes,
A. Palano,
M. Pappagallo,
G. Eigen,
B. Stugu,
L. Sun,
D. N. Brown,
L. T. Kerth,
Yu. G. Kolomensky,
G. Lynch,
H. Koch,
T. Schroeder,
D. J. Asgeirsson,
C. Hearty,
T. S. Mattison,
J. A. McKenna,
A. Khan
, et al. (367 additional authors not shown)
Abstract:
We use 111+/-1 million Upsilon(3S) and 89+/-1 million Upsilon(2S) events recorded by the BaBar detector at the PEP-II B-factory at SLAC to perform a study of radiative transitions between bottomonium states using photons that have been converted to e+e- pairs by the detector material. We observe Upsilon(3S) -> gamma chi_b0,2(1P) decay, make precise measurements of the branching fractions for chi_b…
▽ More
We use 111+/-1 million Upsilon(3S) and 89+/-1 million Upsilon(2S) events recorded by the BaBar detector at the PEP-II B-factory at SLAC to perform a study of radiative transitions between bottomonium states using photons that have been converted to e+e- pairs by the detector material. We observe Upsilon(3S) -> gamma chi_b0,2(1P) decay, make precise measurements of the branching fractions for chi_b1,2(1P,2P) -> gamma Upsilon(1S) and chi_b1,2(2P) -> gamma Upsilon(2S) decays, and search for radiative decay to the eta_b(1S) and eta_b(2S) states.
△ Less
Submitted 21 December, 2011; v1 submitted 27 April, 2011;
originally announced April 2011.
-
Evidence for the h_b(1P) meson in the decay Upsilon(3S) --> pi0 h_b(1P)
Authors:
The BABAR Collaboration,
J. P. Lees,
V. Poireau,
E. Prencipe,
V. Tisserand,
J. Garra Tico,
E. Grauges,
M. Martinelli,
D. A. Milanes,
A. Palano,
M. Pappagallo,
G. Eigen,
B. Stugu,
L. Sun,
D. N. Brown,
L. T. Kerth,
Yu. G. Kolomensky,
G. Lynch,
I. L. Osipenkov,
H. Koch,
T. Schroeder,
D. J. Asgeirsson,
C. Hearty,
T. S. Mattison,
J. A. McKenna
, et al. (386 additional authors not shown)
Abstract:
Using a sample of 122 million Upsilon(3S) events recorded with the BaBar detector at the PEP-II asymmetric-energy e+e- collider at SLAC, we search for the $h_b(1P)$ spin-singlet partner of the P-wave chi_{bJ}(1P) states in the sequential decay Upsilon(3S) --> pi0 h_b(1P), h_b(1P) --> gamma eta_b(1S). We observe an excess of events above background in the distribution of the recoil mass against the…
▽ More
Using a sample of 122 million Upsilon(3S) events recorded with the BaBar detector at the PEP-II asymmetric-energy e+e- collider at SLAC, we search for the $h_b(1P)$ spin-singlet partner of the P-wave chi_{bJ}(1P) states in the sequential decay Upsilon(3S) --> pi0 h_b(1P), h_b(1P) --> gamma eta_b(1S). We observe an excess of events above background in the distribution of the recoil mass against the pi0 at mass 9902 +/- 4(stat.) +/- 2(syst.) MeV/c^2. The width of the observed signal is consistent with experimental resolution, and its significance is 3.1sigma, including systematic uncertainties. We obtain the value (4.3 +/- 1.1(stat.) +/- 0.9(syst.)) x 10^{-4} for the product branching fraction BF(Upsilon(3S)-->pi0 h_b) x BF(h_b-->gamma eta_b).
△ Less
Submitted 17 October, 2011; v1 submitted 22 February, 2011;
originally announced February 2011.
-
Radiative B Decays
Authors:
D. Bard
Abstract:
Recent results in radiative B decays from the Belle and BaBar collaborations are discussed. New measurements of the decay rate and CP asymmetries in b -> s gamma and b -> d gamma decays are presented, together with measurements of the photon spectrum in b -> s gamma.
Recent results in radiative B decays from the Belle and BaBar collaborations are discussed. New measurements of the decay rate and CP asymmetries in b -> s gamma and b -> d gamma decays are presented, together with measurements of the photon spectrum in b -> s gamma.
△ Less
Submitted 8 April, 2009; v1 submitted 31 March, 2009;
originally announced March 2009.
-
LSST: from Science Drivers to Reference Design and Anticipated Data Products
Authors:
Željko Ivezić,
Steven M. Kahn,
J. Anthony Tyson,
Bob Abel,
Emily Acosta,
Robyn Allsman,
David Alonso,
Yusra AlSayyad,
Scott F. Anderson,
John Andrew,
James Roger P. Angel,
George Z. Angeli,
Reza Ansari,
Pierre Antilogus,
Constanza Araujo,
Robert Armstrong,
Kirk T. Arndt,
Pierre Astier,
Éric Aubourg,
Nicole Auza,
Tim S. Axelrod,
Deborah J. Bard,
Jeff D. Barr,
Aurelian Barrau,
James G. Bartlett
, et al. (288 additional authors not shown)
Abstract:
(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the…
▽ More
(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pachón in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$σ$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $δ<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.
△ Less
Submitted 23 May, 2018; v1 submitted 15 May, 2008;
originally announced May 2008.