Showing 1–50 of 99 results for author: Kent, C

Towards improved property prediction of two-dimensional (2D) materials using many-body Quantum Monte Carlo methods

Authors: Daniel Wines, Jeonghwan Ahn, Anouar Benali, Paul R. C. Kent, Jaron T. Krogel, Yongkyung Kwon, Lubos Mitas, Fernando A. Reboredo, Brenda Rubenstein, Kayahan Saritas, Hyeondeok Shin, Ivan Štich, Can Ataca

Abstract: The field of two-dimensional (2D) materials has grown dramatically in the last two decades. 2D materials can be utilized for a variety of next-generation optoelectronic, spintronic, clean energy, and quantum computation applications. These 2D structures, which are often exfoliated from layered van der Waals (vdW) materials, possess highly inhomogeneous electron densities and can possess short- and… ▽ More The field of two-dimensional (2D) materials has grown dramatically in the last two decades. 2D materials can be utilized for a variety of next-generation optoelectronic, spintronic, clean energy, and quantum computation applications. These 2D structures, which are often exfoliated from layered van der Waals (vdW) materials, possess highly inhomogeneous electron densities and can possess short- and long-range electron correlations. The complexities of 2D materials make them challenging to study with standard mean-field electronic structure methods such as density functional theory (DFT), which relies on approximations for the unknown exchange-correlation functional. In order to overcome the limitations of DFT, highly accurate many-body electronic structure approaches such as Diffusion Monte Carlo (DMC) can be utilized. In the past decade, DMC has been used to calculate accurate magnetic, electronic, excitonic, and topological properties in addition to accurately capturing interlayer interactions and cohesion and adsorption energetics of 2D materials. This approach has been applied to 2D systems of wide interest including graphene, phosphorene, MoS$_2$, CrI$_3$, VSe$_2$, GaSe, GeSe, borophene, and several others. In this review article, we highlight some successful recent applications of DMC to 2D systems for improved property predictions beyond standard DFT. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Analysis of the BraTS 2023 Intracranial Meningioma Segmentation Challenge

Authors: Dominic LaBella, Ujjwal Baid, Omaditya Khanna, Shan McBurney-Lin, Ryan McLean, Pierre Nedelec, Arif Rashid, Nourel Hoda Tahon, Talissa Altes, Radhika Bhalerao, Yaseen Dhemesh, Devon Godfrey, Fathi Hilal, Scott Floyd, Anastasia Janas, Anahita Fathi Kazerooni, John Kirkpatrick, Collin Kent, Florian Kofler, Kevin Leu, Nazanin Maleki, Bjoern Menze, Maxence Pajot, Zachary J. Reitman, Jeffrey D. Rudie , et al. (96 additional authors not shown)

Abstract: We describe the design and results from the BraTS 2023 Intracranial Meningioma Segmentation Challenge. The BraTS Meningioma Challenge differed from prior BraTS Glioma challenges in that it focused on meningiomas, which are typically benign extra-axial tumors with diverse radiologic and anatomical presentation and a propensity for multiplicity. Nine participating teams each developed deep-learning… ▽ More We describe the design and results from the BraTS 2023 Intracranial Meningioma Segmentation Challenge. The BraTS Meningioma Challenge differed from prior BraTS Glioma challenges in that it focused on meningiomas, which are typically benign extra-axial tumors with diverse radiologic and anatomical presentation and a propensity for multiplicity. Nine participating teams each developed deep-learning automated segmentation models using image data from the largest multi-institutional systematically expert annotated multilabel multi-sequence meningioma MRI dataset to date, which included 1000 training set cases, 141 validation set cases, and 283 hidden test set cases. Each case included T2, T2/FLAIR, T1, and T1Gd brain MRI sequences with associated tumor compartment labels delineating enhancing tumor, non-enhancing tumor, and surrounding non-enhancing T2/FLAIR hyperintensity. Participant automated segmentation models were evaluated and ranked based on a scoring system evaluating lesion-wise metrics including dice similarity coefficient (DSC) and 95% Hausdorff Distance. The top ranked team had a lesion-wise median dice similarity coefficient (DSC) of 0.976, 0.976, and 0.964 for enhancing tumor, tumor core, and whole tumor, respectively and a corresponding average DSC of 0.899, 0.904, and 0.871, respectively. These results serve as state-of-the-art benchmarks for future pre-operative meningioma automated segmentation algorithms. Additionally, we found that 1286 of 1424 cases (90.3%) had at least 1 compartment voxel abutting the edge of the skull-stripped image edge, which requires further investigation into optimal pre-processing face anonymization steps. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 16 pages, 11 tables, 10 figures, MICCAI

Force-free identification of minimum-energy pathways and transition states for stochastic electronic structure theories

Authors: Gopal R. Iyer, Noah Whelpley, Juha Tiihonen, Paul R. C. Kent, Jaron T. Krogel, Brenda M. Rubenstein

Abstract: Stochastic electronic structure theories, e.g., Quantum Monte Carlo methods, enable highly accurate total energy calculations which in principle can be used to construct highly accurate potential energy surfaces. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and trans… ▽ More Stochastic electronic structure theories, e.g., Quantum Monte Carlo methods, enable highly accurate total energy calculations which in principle can be used to construct highly accurate potential energy surfaces. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged-elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method - previously developed for QMC structural optimization - to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and on saddle points, we show that it is possible to identify MEPs and TSs using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia molecule and an SN2 reaction. We validate our results using Density Functional Theory- and coupled cluster-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities - free energy differences, rate constants, and equilibrium constants - that incorporates stochastically-optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations, at significantly reduced computational expense. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Neural Population Decoding and Imbalanced Multi-Omic Datasets For Cancer Subtype Diagnosis

Authors: Charles Theodore Kent, Leila Bagheriye, Johan Kwisthout

Abstract: Recent strides in the field of neural computation has seen the adoption of Winner Take All (WTA) circuits to facilitate the unification of hierarchical Bayesian inference and spiking neural networks as a neurobiologically plausible model of information processing. Current research commonly validates the performance of these networks via classification tasks, particularly of the MNIST dataset. Howe… ▽ More Recent strides in the field of neural computation has seen the adoption of Winner Take All (WTA) circuits to facilitate the unification of hierarchical Bayesian inference and spiking neural networks as a neurobiologically plausible model of information processing. Current research commonly validates the performance of these networks via classification tasks, particularly of the MNIST dataset. However, researchers have not yet reached consensus about how best to translate the stochastic responses from these networks into discrete decisions, a process known as population decoding. Despite being an often underexamined part of SNNs, in this work we show that population decoding has a significanct impact on the classification performance of WTA networks. For this purpose, we apply a WTA network to the problem of cancer subtype diagnosis from multi omic data, using datasets from The Cancer Genome Atlas (TCGA). In doing so we utilise a novel implementation of gene similarity networks, a feature encoding technique based on Kohoens self organising map algorithm. We further show that the impact of selecting certain population decoding methods is amplified when facing imbalanced datasets. △ Less

Submitted 6 January, 2024; originally announced January 2024.

Comments: This paper has been accepted in BIOINFORMATICS 2024 (BIOSTEC 2024)

Path-homotopy is equivalent to $\mathbb{R}$-tree reduction

Authors: Jeremy Brazas, Gregory R. Conner, Paul Fabel, Curtis Kent

Abstract: Suppose a path $α$ factors through an $\mathbb{R}$-tree $T$ as $α=q \circ p $. Let $r$ parameterize the unique geodesic in $T$ joining the endpoints of $p$. Then we say that the path $β=q \circ r$ is obtained from $α$ by "geodesic $\mathbb{R}$-tree reduction." Essentially, $β$ is obtained from $α$ by deleting one-dimensional back-tracking. In this paper, we show that any two homotopic paths are ge… ▽ More Suppose a path $α$ factors through an $\mathbb{R}$-tree $T$ as $α=q \circ p $. Let $r$ parameterize the unique geodesic in $T$ joining the endpoints of $p$. Then we say that the path $β=q \circ r$ is obtained from $α$ by "geodesic $\mathbb{R}$-tree reduction." Essentially, $β$ is obtained from $α$ by deleting one-dimensional back-tracking. In this paper, we show that any two homotopic paths are geodesic $\mathbb{R}$-tree reductions of some single common path. Hence, the equivalence relation on paths generated by geodesic $\mathbb{R}$-tree reduction is precisely path-homotopy. The common path is explicitly constructed and is necessarily space-filling in the image of a given path-homotopy. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 16 pages, 6 page appendix, 6 figures

MSC Class: 54F15; 55P10; 54F50; 55R65; 57M10; 20E08

A Metacognitive Approach to Out-of-Distribution Detection for Segmentation

Authors: Meghna Gummadi, Cassandra Kent, Karl Schmeckpeper, Eric Eaton

Abstract: Despite outstanding semantic scene segmentation in closed-worlds, deep neural networks segment novel instances poorly, which is required for autonomous agents acting in an open world. To improve out-of-distribution (OOD) detection for segmentation, we introduce a metacognitive approach in the form of a lightweight module that leverages entropy measures, segmentation predictions, and spatial contex… ▽ More Despite outstanding semantic scene segmentation in closed-worlds, deep neural networks segment novel instances poorly, which is required for autonomous agents acting in an open world. To improve out-of-distribution (OOD) detection for segmentation, we introduce a metacognitive approach in the form of a lightweight module that leverages entropy measures, segmentation predictions, and spatial context to characterize the segmentation model's uncertainty and detect pixel-wise OOD data in real-time. Additionally, our approach incorporates a novel method of generating synthetic OOD data in context with in-distribution data, which we use to fine-tune existing segmentation models with maximum entropy training. This further improves the metacognitive module's performance without requiring access to OOD data while enabling compatibility with established pre-trained models. Our resulting approach can reliably detect OOD instances in a scene, as shown by state-of-the-art performance on OOD detection for semantic segmentation benchmarks. △ Less

Submitted 4 October, 2023; originally announced November 2023.

Software engineering to sustain a high-performance computing scientific application: QMCPACK

Authors: William F. Godoy, Steven E. Hahn, Michael M. Walsh, Philip W. Fackler, Jaron T. Krogel, Peter W. Doak, Paul R. C. Kent, Alfredo A. Correa, Ye Luo, Mark Dewing

Abstract: We provide an overview of the software engineering efforts and their impact in QMCPACK, a production-level ab-initio Quantum Monte Carlo open-source code targeting high-performance computing (HPC) systems. Aspects included are: (i) strategic expansion of continuous integration (CI) targeting CPUs, using GitHub Actions runners, and NVIDIA and AMD GPUs in pre-exascale systems, using self-hosted hard… ▽ More We provide an overview of the software engineering efforts and their impact in QMCPACK, a production-level ab-initio Quantum Monte Carlo open-source code targeting high-performance computing (HPC) systems. Aspects included are: (i) strategic expansion of continuous integration (CI) targeting CPUs, using GitHub Actions runners, and NVIDIA and AMD GPUs in pre-exascale systems, using self-hosted hardware; (ii) incremental reduction of memory leaks using sanitizers, (iii) incorporation of Docker containers for CI and reproducibility, and (iv) refactoring efforts to improve maintainability, testing coverage, and memory lifetime management. We quantify the value of these improvements by providing metrics to illustrate the shift towards a predictive, rather than reactive, sustainable maintenance approach. Our goal, in documenting the impact of these efforts on QMCPACK, is to contribute to the body of knowledge on the importance of research software engineering (RSE) for the sustainability of community HPC codes and scientific discovery at scale. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: Accepted at the first US-RSE Conference, USRSE2023, https://us-rse.org/usrse23/, 8 pages, 3 figures, 4 tables

Robotic Manipulation Datasets for Offline Compositional Reinforcement Learning

Authors: Marcel Hussing, Jorge A. Mendez, Anisha Singrodia, Cassandra Kent, Eric Eaton

Abstract: Offline reinforcement learning (RL) is a promising direction that allows RL agents to pre-train on large datasets, avoiding the recurrence of expensive data collection. To advance the field, it is crucial to generate large-scale datasets. Compositional RL is particularly appealing for generating such large datasets, since 1)~it permits creating many tasks from few components, 2)~the task structure… ▽ More Offline reinforcement learning (RL) is a promising direction that allows RL agents to pre-train on large datasets, avoiding the recurrence of expensive data collection. To advance the field, it is crucial to generate large-scale datasets. Compositional RL is particularly appealing for generating such large datasets, since 1)~it permits creating many tasks from few components, 2)~the task structure may enable trained agents to solve new tasks by combining relevant learned components, and 3)~the compositional dimensions provide a notion of task relatedness. This paper provides four offline RL datasets for simulated robotic manipulation created using the $256$ tasks from CompoSuite [Mendez at al., 2022a]. Each dataset is collected from an agent with a different degree of performance, and consists of $256$ million transitions. We provide training and evaluation settings for assessing an agent's ability to learn compositional task policies. Our benchmarking experiments show that current offline RL methods can learn the training tasks to some extent and that compositional methods outperform non-compositional methods. Yet current methods are unable to extract the compositional structure to generalize to unseen tasks, highlighting a need for future research in offline compositional RL. △ Less

Submitted 15 July, 2024; v1 submitted 13 July, 2023; originally announced July 2023.

Comments: Published as a conference paper at the First Reinforcement Learning Conference (RLC)

JARVIS-Leaderboard: A Large Scale Benchmark of Materials Design Methods

Authors: Kamal Choudhary, Daniel Wines, Kangming Li, Kevin F. Garrity, Vishu Gupta, Aldo H. Romero, Jaron T. Krogel, Kayahan Saritas, Addis Fuhr, Panchapakesan Ganesh, Paul R. C. Kent, Keqiang Yan, Yuchao Lin, Shuiwang Ji, Ben Blaiszik, Patrick Reiser, Pascal Friederich, Ankit Agrawal, Pratyush Tiwary, Eric Beyerle, Peter Minch, Trevor David Rhone, Ichiro Takeuchi, Robert B. Wexler, Arun Mannodi-Kanakkithodi , et al. (13 additional authors not shown)

Abstract: Lack of rigorous reproducibility and validation are major hurdles for scientific development across many fields. Materials science in particular encompasses a variety of experimental and theoretical approaches that require careful benchmarking. Leaderboard efforts have been developed previously to mitigate these issues. However, a comprehensive comparison and benchmarking on an integrated platform… ▽ More Lack of rigorous reproducibility and validation are major hurdles for scientific development across many fields. Materials science in particular encompasses a variety of experimental and theoretical approaches that require careful benchmarking. Leaderboard efforts have been developed previously to mitigate these issues. However, a comprehensive comparison and benchmarking on an integrated platform with multiple data modalities with both perfect and defect materials data is still lacking. This work introduces JARVIS-Leaderboard, an open-source and community-driven platform that facilitates benchmarking and enhances reproducibility. The platform allows users to set up benchmarks with custom tasks and enables contributions in the form of dataset, code, and meta-data submissions. We cover the following materials design categories: Artificial Intelligence (AI), Electronic Structure (ES), Force-fields (FF), Quantum Computation (QC) and Experiments (EXP). For AI, we cover several types of input data, including atomic structures, atomistic images, spectra, and text. For ES, we consider multiple ES approaches, software packages, pseudopotentials, materials, and properties, comparing results to experiment. For FF, we compare multiple approaches for material property predictions. For QC, we benchmark Hamiltonian simulations using various quantum algorithms and circuits. Finally, for experiments, we use the inter-laboratory approach to establish benchmarks. There are 1281 contributions to 274 benchmarks using 152 methods with more than 8 million data-points, and the leaderboard is continuously expanding. The JARVIS-Leaderboard is available at the website: https://pages.nist.gov/jarvis_leaderboard △ Less

Submitted 26 March, 2024; v1 submitted 20 June, 2023; originally announced June 2023.

The ASNR-MICCAI Brain Tumor Segmentation (BraTS) Challenge 2023: Intracranial Meningioma

Authors: Dominic LaBella, Maruf Adewole, Michelle Alonso-Basanta, Talissa Altes, Syed Muhammad Anwar, Ujjwal Baid, Timothy Bergquist, Radhika Bhalerao, Sully Chen, Verena Chung, Gian-Marco Conte, Farouk Dako, James Eddy, Ivan Ezhov, Devon Godfrey, Fathi Hilal, Ariana Familiar, Keyvan Farahani, Juan Eugenio Iglesias, Zhifan Jiang, Elaine Johanson, Anahita Fathi Kazerooni, Collin Kent, John Kirkpatrick, Florian Kofler , et al. (35 additional authors not shown)

Abstract: Meningiomas are the most common primary intracranial tumor in adults and can be associated with significant morbidity and mortality. Radiologists, neurosurgeons, neuro-oncologists, and radiation oncologists rely on multiparametric MRI (mpMRI) for diagnosis, treatment planning, and longitudinal treatment monitoring; yet automated, objective, and quantitative tools for non-invasive assessment of men… ▽ More Meningiomas are the most common primary intracranial tumor in adults and can be associated with significant morbidity and mortality. Radiologists, neurosurgeons, neuro-oncologists, and radiation oncologists rely on multiparametric MRI (mpMRI) for diagnosis, treatment planning, and longitudinal treatment monitoring; yet automated, objective, and quantitative tools for non-invasive assessment of meningiomas on mpMRI are lacking. The BraTS meningioma 2023 challenge will provide a community standard and benchmark for state-of-the-art automated intracranial meningioma segmentation models based on the largest expert annotated multilabel meningioma mpMRI dataset to date. Challenge competitors will develop automated segmentation models to predict three distinct meningioma sub-regions on MRI including enhancing tumor, non-enhancing tumor core, and surrounding nonenhancing T2/FLAIR hyperintensity. Models will be evaluated on separate validation and held-out test datasets using standardized metrics utilized across the BraTS 2023 series of challenges including the Dice similarity coefficient and Hausdorff distance. The models developed during the course of this challenge will aid in incorporation of automated meningioma MRI segmentation into clinical practice, which will ultimately improve care of patients with meningioma. △ Less

Submitted 12 May, 2023; originally announced May 2023.

A Domain-Agnostic Approach for Characterization of Lifelong Learning Systems

Authors: Megan M. Baker, Alexander New, Mario Aguilar-Simon, Ziad Al-Halah, Sébastien M. R. Arnold, Ese Ben-Iwhiwhu, Andrew P. Brna, Ethan Brooks, Ryan C. Brown, Zachary Daniels, Anurag Daram, Fabien Delattre, Ryan Dellana, Eric Eaton, Haotian Fu, Kristen Grauman, Jesse Hostetler, Shariq Iqbal, Cassandra Kent, Nicholas Ketz, Soheil Kolouri, George Konidaris, Dhireesha Kudithipudi, Erik Learned-Miller, Seungwon Lee , et al. (22 additional authors not shown)

Abstract: Despite the advancement of machine learning techniques in recent years, state-of-the-art systems lack robustness to "real world" events, where the input distributions and tasks encountered by the deployed systems will not be limited to the original training context, and systems will instead need to adapt to novel distributions and tasks while deployed. This critical gap may be addressed through th… ▽ More Despite the advancement of machine learning techniques in recent years, state-of-the-art systems lack robustness to "real world" events, where the input distributions and tasks encountered by the deployed systems will not be limited to the original training context, and systems will instead need to adapt to novel distributions and tasks while deployed. This critical gap may be addressed through the development of "Lifelong Learning" systems that are capable of 1) Continuous Learning, 2) Transfer and Adaptation, and 3) Scalability. Unfortunately, efforts to improve these capabilities are typically treated as distinct areas of research that are assessed independently, without regard to the impact of each separate capability on other aspects of the system. We instead propose a holistic approach, using a suite of metrics and an evaluation framework to assess Lifelong Learning in a principled way that is agnostic to specific domains or system techniques. Through five case studies, we show that this suite of metrics can inform the development of varied and complex Lifelong Learning systems. We highlight how the proposed suite of metrics quantifies performance trade-offs present during Lifelong Learning system development - both the widely discussed Stability-Plasticity dilemma and the newly proposed relationship between Sample Efficient and Robust Learning. Further, we make recommendations for the formulation and use of metrics to guide the continuing development of Lifelong Learning systems and assess their progress in the future. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: To appear in Neural Networks

DFT+U and Quantum Monte Carlo study of electronic and optical properties of AgNiO$_2$ and AgNi$_{1-x}$Co$_{x}$O$_2$ delafossite

Authors: Hyeondeok Shin, Panchapakesan Ganesh, Paul R. C. Kent, Anouar Benali, Anand Bhattacharya, Ho Nyung Lee, Olle Heinonen, Jaron T. Krogel

Abstract: As the only semimetallic $d^{10}$-based delafossite, AgNiO$_2$ has received a great deal of attention due to both its unique semimetallicity and its antiferromagnetism in the NiO$_2$ layer that is coupled with a lattice distortion. In contrast, other delafossites such as AgCoO$_2$ are insulating. Here we study how the electronic structure of AgNi$_{1-x}$Co$_{x}$O$_2$ alloys vary with Ni/Co concent… ▽ More As the only semimetallic $d^{10}$-based delafossite, AgNiO$_2$ has received a great deal of attention due to both its unique semimetallicity and its antiferromagnetism in the NiO$_2$ layer that is coupled with a lattice distortion. In contrast, other delafossites such as AgCoO$_2$ are insulating. Here we study how the electronic structure of AgNi$_{1-x}$Co$_{x}$O$_2$ alloys vary with Ni/Co concentration, in order to investigate the electronic properties and phase stability of the intermetallics. While the electronic and magnetic structure of delafossites have been studied using Density Functional Theory (DFT), earlier studies have not included corrections for strong on-site Coulomb interactions. In order to treat these interactions accurately, in this study we use Quantum Monte Carlo (QMC) simulations to obtain accurate estimates for the electronic and magnetic properties of AgNiO$_2$. By comparison to DFT results we show that these electron correlations are critical to account for. We show that Co doping on the magnetic Ni sites results in a metal-insulator transition near $x\sim 0.33$, and reentrant behavior near $x\sim 0.66$ △ Less

Submitted 21 July, 2023; v1 submitted 9 April, 2022; originally announced April 2022.

Comments: 25 Pages, 12 figures

Interfacial charge transfer and interaction in the MXene/TiO2 heterostructures

Authors: Lihua Xu, Tao Wu, Paul R. C. Kent, De-en Jiang

Abstract: Hybrid materials of MXenes (2D carbides and nitrides) and transition-metal oxides (TMOs) have shown great promise in electrical energy storage and 2D heterostructures have been proposed as the next-generation electrode materials to expand the limits of current technology. Here we use first principles density functional theory to investigate the interfacial structure, energetics, and electronic pro… ▽ More Hybrid materials of MXenes (2D carbides and nitrides) and transition-metal oxides (TMOs) have shown great promise in electrical energy storage and 2D heterostructures have been proposed as the next-generation electrode materials to expand the limits of current technology. Here we use first principles density functional theory to investigate the interfacial structure, energetics, and electronic properties of the heterostructures of MXenes (Tin+1CnT2; T=terminal groups) and anatase TiO2. We find that the greatest work-function differences are between OH-terminated-MXene (1.6 eV) and anatase TiO2(101) (6.4 eV), resulting in the largest interfacial electron transfer (~0.9 e/nm2 across the interface) from MXene to the TiO2 layer. This interface also has the strongest adhesion and further strengthened by hydrogen bond formation. For O-, F-, or mixed O-/F- terminated Tin+1Cn MXenes, electron transfer is minimal and interfacial adhesion is weak for their heterostructures with TiO2. The strong dependence of the interfacial properties of the MXene/TiO2 heterostructures on the surface chemistry of the MXenes will be useful to tune the heterostructures for electric-energy-storage applications. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: 13 pages, 5 figures

Journal ref: Phys. Rev. Mater., 2021, 5, 054007

Inverse limit slender groups

Authors: Gregory Conner, Wolfgang Herfort, Curtis Kent, Peter Pavesic

Abstract: Classically, an abelian group $G$ is said to be slender if every homomorphism from the countable product $\mathbb Z^{\mathbb N}$ to $G$ factors through the projection to some finite product $\mathbb Z^n$. Various authors have proposed generalizations to non-commutative groups, resulting in a plethora of similar but not completely equivalent concepts. In the first part of this work we present a uni… ▽ More Classically, an abelian group $G$ is said to be slender if every homomorphism from the countable product $\mathbb Z^{\mathbb N}$ to $G$ factors through the projection to some finite product $\mathbb Z^n$. Various authors have proposed generalizations to non-commutative groups, resulting in a plethora of similar but not completely equivalent concepts. In the first part of this work we present a unified treatment of these concepts and examine how are they related. In the second part of the paper we study slender groups in the context of co-small objects in certain categories, and give several new applications including the proof that certain homology groups of Barratt-Milnor spaces are cotorsion groups and a universal coefficients theorem for Čech cohomology with coefficients in a slender group. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Frank-Wolfe Methods in Probability Space

Authors: Carson Kent, Jose Blanchet, Peter Glynn

Abstract: We introduce a new class of Frank-Wolfe algorithms for minimizing differentiable functionals over probability measures. This framework can be shown to encompass a diverse range of tasks in areas such as artificial intelligence, reinforcement learning, and optimization. Concrete computational complexities for these algorithms are established and demonstrate that these methods enjoy convergence in r… ▽ More We introduce a new class of Frank-Wolfe algorithms for minimizing differentiable functionals over probability measures. This framework can be shown to encompass a diverse range of tasks in areas such as artificial intelligence, reinforcement learning, and optimization. Concrete computational complexities for these algorithms are established and demonstrate that these methods enjoy convergence in regimes that go beyond convexity and require minimal regularity of the underlying functional. Novel techniques used to obtain these results also lead to the development of new complexity bounds and duality theorems for a family of distributionally robust optimization problems. The performance of our method is demonstrated on several nonparametric estimation problems. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Origin of Metal-Insulator Transitions in Correlated Perovskite Metals

Authors: M. Chandler Bennett, Guoxiang Hu, Guangming Wang, Olle Heinonen, Paul R. C. Kent, Jaron T. Krogel, Panchapakesan Ganesh

Abstract: The mechanisms that drive metal-to-insulator transitions (MIT) in correlated solids are not fully understood. For example, the perovskite (PV) SrCoO3 is a FM metal while the oxygen-deficient (n-doped) brownmillerite (BM) SrCoO2.5 is an anti-ferromagnetic (AFM) insulator. Given the magnetic and structural transitions that accompany the MIT, the driver for such a MIT transition is unclear. We also o… ▽ More The mechanisms that drive metal-to-insulator transitions (MIT) in correlated solids are not fully understood. For example, the perovskite (PV) SrCoO3 is a FM metal while the oxygen-deficient (n-doped) brownmillerite (BM) SrCoO2.5 is an anti-ferromagnetic (AFM) insulator. Given the magnetic and structural transitions that accompany the MIT, the driver for such a MIT transition is unclear. We also observe that the perovskite metals LaNiO3, SrFeO3, and SrCoO3 also undergo MIT when n-doped via high-to-low valence compositional changes. Also, pressurizing the insulating BM SrCoO2.5 phase, drives a gap closing. Using DFT and correlated diffusion Monte Carlo approaches we demonstrate that the ABO3 perovskites most prone to MIT are self hole-doped materials, reminiscent of a negative charge-transfer system. Upon n-doping away from the negative-charge transfer metallic phase, an underlying charge-lattice (or e-phonon) coupling drives the system to a bond-disproportionated gapped state, thereby achieving ligand hole passivation at certain sites only, leading to charge-disproportionated states. The size of the gap opened is correlated with the size of the hole-filling at these ligand sites. This suggests that the interactions driving the gap opening to realize a MIT even in correlated metals is the charge-transfer energy, but it couples with the underlying phonons to enable the transition to the insulating phase. Other orderings (magnetic, charge, etc.) driven by weaker interactions are secondary and may assist gap openings at small dopings, but its the charge-transfer energy that predominantly determines the bandgap, with a negative energy preferring the metallic phase. This n-doping can be achieved by modulations in stoichiometry or composition or pressure. Hence, controlling the amount of the ligand-hole is key in controlling MIT. We compare our predictions to experiments where possible. △ Less

Submitted 25 January, 2022; v1 submitted 17 March, 2021; originally announced March 2021.

Uncountable groups and the geometry of inverse limits of coverings

Authors: Gregory Conner, Wolfgang Herfort, Curtis Kent, Petar Pavesic

Abstract: In this paper we develop a new approach to the study of uncountable fundamental groups by using Hurewicz fibrations with the unique path-lifting property (lifting spaces for short) as a replacement for covering spaces. In particular, we consider the inverse limit of a sequence of covering spaces of $X$. It is known that the path-connectivity of the inverse limit can be expressed by means of the de… ▽ More In this paper we develop a new approach to the study of uncountable fundamental groups by using Hurewicz fibrations with the unique path-lifting property (lifting spaces for short) as a replacement for covering spaces. In particular, we consider the inverse limit of a sequence of covering spaces of $X$. It is known that the path-connectivity of the inverse limit can be expressed by means of the derived inverse limit functor $\varprojlim^1$, which is, however, notoriously difficult to compute when the $π_1(X)$ is uncountable.To circumvent this difficulty, we express the set of path-components of the inverse limit, $\widehat X$, in terms of the functors $\varprojlim$ and $\varprojlim^1$ applied to sequences of countable groups arising from polyhedral approximations of $X$. A consequence of our computation is that path-connectedness of lifting space implies that $π_1(\tilde X)$ supplements $π_1(X)$ in $\checkπ_1(X)$ where $\checkπ_1(X)$ is the inverse limit of fundamental groups of polyhedral approximations of $X$. As an application we show that $\mathcal G\cdot \ker_{\mathbb Z}(\widehat F)= \widehat F\ne\mathcal G\cdot \ker_{B(1,n)}(\widehat F)$, where $\widehat F$ is the canonical inverse limit of finite rank free groups, $\mathcal G$ is the fundamental group of the Hawaiian Earring, and $\ker_A(\widehat F)$ is the intersection of kernels of homomorphisms from $\widehat{F}$ to $A$. △ Less

Submitted 27 January, 2021; originally announced January 2021.

Towards a Systematic Improvement of the Fixed-Node Approximation in Diffusion Monte Carlo for Solids -- A Case Study In Diamond

Authors: Anouar Benali, Kevin Gasperich, Kenneth D. Jordan, Thomas Applencourt, Ye Luo, M. Chandler Bennett, Jaron T. Krogel, Luke Shulenburger, Paul R. C. Kent, Pierre-François Loos, Anthony Scemama, Michel Caffarel

Abstract: While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for \textit{ab initio} electronic structure calculations, in practice the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros) must be used. This approximation introduces a variational error in the energy that potentially can be tested and sy… ▽ More While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for \textit{ab initio} electronic structure calculations, in practice the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros) must be used. This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. △ Less

Submitted 29 October, 2020; v1 submitted 22 July, 2020; originally announced July 2020.

Systematic Comparison and Cross-validation of Fixed-Node Diffusion Monte Carlo and Phaseless Auxiliary-Field Quantum Monte Carlo in Solids

Authors: Fionn D. Malone, Anouar Benali, Miguel A. Morales, Michel Caffarel, P. R. C. Kent, Luke Shulenburger

Abstract: Quantum Monte Carlo (QMC) methods are some of the most accurate methods for simulating correlated electronic systems. We investigate the compatibility, strengths and weaknesses of two such methods, namely, diffusion Monte Carlo (DMC) and auxiliary-field quantum Monte Carlo (AFQMC). The multi-determinant trial wave functions employed in both approaches are generated using the configuration interact… ▽ More Quantum Monte Carlo (QMC) methods are some of the most accurate methods for simulating correlated electronic systems. We investigate the compatibility, strengths and weaknesses of two such methods, namely, diffusion Monte Carlo (DMC) and auxiliary-field quantum Monte Carlo (AFQMC). The multi-determinant trial wave functions employed in both approaches are generated using the configuration interaction using a perturbative selection made iteratively (CIPSI) technique. Complete basis set full configuration interaction (CBS-FCI) energies estimated with CIPSI are used as a reference in this comparative study between DMC and AFQMC. By focusing on a set of canonical finite size solid state systems, we show that both QMC methods can be made to systematically converge towards the same energy once basis set effects and systematic biases have been removed. AFQMC shows a much smaller dependence on the trial wavefunction than DMC while simultaneously exhibiting a much larger basis set dependence. We outline some of the remaining challenges and opportunities for improving these approaches. △ Less

Submitted 10 July, 2020; originally announced July 2020.

Journal ref: Phys. Rev. B 102, 161104 (2020)

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion Quantum Monte Carlo

Authors: P. R. C. Kent, Abdulgani Annaberdiyev, Anouar Benali, M. Chandler Bennett, Edgar Josue Landinez Borda, Peter Doak, Kenneth D. Jordan, Jaron T. Krogel, Ilkka Kylanpaa, Joonho Lee, Ye Luo, Fionn D. Malone, Cody A. Melton, Lubos Mitas, Miguel A. Morales, Eric Neuscamman, Fernando A. Reboredo, Brenda Rubenstein, Kayahan Saritas, Shiv Upadhyay, Hongxia Hao, Guangming Wang, Shuai Zhang, Luning Zhao

Abstract: We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and me… ▽ More We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing accuracy. Advances in real space methods include techniques for accurate computation of band gaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods including GW and density functional based techniques. To provide an improved foundation for these calculations we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK. △ Less

Submitted 6 May, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Report number: SAND2020-2836 O

Journal ref: J. Chem. Phys. 152, 174105 (2020)

Doped NiO: the Mottness of a charge transfer insulator

Authors: Friederike Wrobel, Hyowon Park, Changhee Sohn, Haw-Wen Hsia, Jian-Min Zuo, Hyeondeok Shin, Ho Nyung Lee, P. Ganesh, Anouar Benali, Paul R. C. Kent, Olle Heinonen, Anand Bhattacharya

Abstract: The evolution of the electronic structures of strongly correlated insulators with doping has long been a central fundamental question in condensed matter physics; it is also of great practical relevance for applications. We have studied the evolution of NiO under hole {\em and} electron doping using high-quality thin film and a wide range of experimental and theoretical methods. The evolution is i… ▽ More The evolution of the electronic structures of strongly correlated insulators with doping has long been a central fundamental question in condensed matter physics; it is also of great practical relevance for applications. We have studied the evolution of NiO under hole {\em and} electron doping using high-quality thin film and a wide range of experimental and theoretical methods. The evolution is in both cases very smooth with dopant concentration. The band gap is asymmetric under electron and hole doping, consistent with a charge-transfer insulator picture, and is reduced faster under hole than electron doping. For both electron and hole doping, occupied states are introduced at the top of the valence band. The formation of deep donor levels under electron doping and the inability to pin otherwise empty states near the conduction band edge is indicative that local electron addition and removal energies are dominated by a Mott-like Hubbard $U$-interaction even though the global bandgap is predominantly a charge-transfer type gap. △ Less

Submitted 11 May, 2020; v1 submitted 30 January, 2020; originally announced January 2020.

Comments: Revised version

Journal ref: Phys. Rev. B 101, 195128 (2020)

An empirical, Bayesian approach to modelling the impact of weather on crop yield: maize in the US

Authors: Raphael Shirley, Edward Pope, Myles Bartlett, Seb Oliver, Novi Quadrianto, Peter Hurley, Steven Duivenvoorden, Phil Rooney, Adam B. Barrett, Chris Kent, James Bacon

Abstract: We apply an empirical, data-driven approach for describing crop yield as a function of monthly temperature and precipitation by employing generative probabilistic models with parameters determined through Bayesian inference. Our approach is applied to state-scale maize yield and meteorological data for the US Corn Belt from 1981 to 2014 as an exemplar, but would be readily transferable to other cr… ▽ More We apply an empirical, data-driven approach for describing crop yield as a function of monthly temperature and precipitation by employing generative probabilistic models with parameters determined through Bayesian inference. Our approach is applied to state-scale maize yield and meteorological data for the US Corn Belt from 1981 to 2014 as an exemplar, but would be readily transferable to other crops, locations and spatial scales. Experimentation with a number of models shows that maize growth rates can be characterised by a two-dimensional Gaussian function of temperature and precipitation with monthly contributions accumulated over the growing period. This approach accounts for non-linear growth responses to the individual meteorological variables, and allows for interactions between them. Our models correctly identify that temperature and precipitation have the largest impact on yield in the six months prior to the harvest, in agreement with the typical growing season for US maize (April to September). Maximal growth rates occur for monthly mean temperature 18-19$^\circ$C, corresponding to a daily maximum temperature of 24-25$^\circ$C (in broad agreement with previous work) and monthly total precipitation 115 mm. Our approach also provides a self-consistent way of investigating climate change impacts on current US maize varieties in the absence of adaptation measures. Keeping precipitation and growing area fixed, a temperature increase of $2^\circ$C, relative to 1981-2014, results in the mean yield decreasing by 8\%, while the yield variance increases by a factor of around 3. We thus provide a flexible, data-driven framework for exploring the impacts of natural climate variability and climate change on globally significant crops based on their observed behaviour. In concert with other approaches, this can help inform the development of adaptation strategies that will ensure food security under a changing climate. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Robots as Actors in a Film: No War, A Robot Story

Authors: Andreagiovanni Reina, Viktor Ioannou, Junjin Chen, Lu Lu, Charles Kent, James A. R. Marshall

Abstract: Will the Third World War be fought by robots? This short film is a light-hearted comedy that aims to trigger an interesting discussion and reflexion on the terrifying killer-robot stories that increasingly fill us with dread when we read the news headlines. The fictional scenario takes inspiration from current scientific research and describes a future where robots are asked by humans to join the… ▽ More Will the Third World War be fought by robots? This short film is a light-hearted comedy that aims to trigger an interesting discussion and reflexion on the terrifying killer-robot stories that increasingly fill us with dread when we read the news headlines. The fictional scenario takes inspiration from current scientific research and describes a future where robots are asked by humans to join the war. Robots are divided, sparking protests in robot society... will robots join the conflict or will they refuse to be employed in human warfare? Food for thought for engineers, roboticists and anyone imagining what the upcoming robot revolution could look like. We let robots pop on camera to tell a story, taking on the role of actors playing in the film, instructed through code on how to "act" for each scene. △ Less

Submitted 27 October, 2019; originally announced October 2019.

Interfacial and Electronic Properties of Heterostructures of MXene and Graphene

Authors: Rui Li, Weiwei Sun, Cheng Zhan, Paul R. C. Kent, De-en Jiang

Abstract: MXene-based heterostructures have received considerable interest owing to their unique properties. Herein, we examine various heterostructures of a prototypical MXene and graphene using density functional theory. We find that the adhesion energy, charge transfer, and band structure of these heterostructures are sensitive not only to the surface functional group, but also to the stacking order. Dif… ▽ More MXene-based heterostructures have received considerable interest owing to their unique properties. Herein, we examine various heterostructures of a prototypical MXene and graphene using density functional theory. We find that the adhesion energy, charge transfer, and band structure of these heterostructures are sensitive not only to the surface functional group, but also to the stacking order. Difference in work function dictates the direction and amount of electron transfer across the interface, which causes a shift in the Dirac point of the graphene bands in the heterostructures of monolayer graphene and monolayer MXene. In the heterostructures of bilayer graphene and monolayer MXene, the interface breaks the symmetry of the bilayer graphene; in the case of the AB-stacking bilayer, the electron transfer leads to an interfacial electric field that opens up a gap in the graphene bands at the K point. This internal polarization strengthens both the interfacial adhesions and the cohesion between the two graphene layers. The MXene-graphene-MXene and graphene-MXene-graphene sandwich structures behave as two mirror-symmetric MXene-graphene interfaces. Our first principles studies provide a comprehensive understanding for the interaction between a typical MXene and graphene. △ Less

Submitted 2 February, 2019; originally announced February 2019.

Comments: 17 pages, 7 figures

Journal ref: Phys. Rev. B 99, 085429 (2019)

Doping a Bad Metal: Origin of Suppression of Metal-Insulator Transition in Non-Stoichiometric VO$_2$

Authors: P. Ganesh, Frank Lechermann, Ilkka Kylanpaa, Jaron Krogel, Paul R. C. Kent, Olle Heinonen

Abstract: Rutile ($R$) phase VO$_2$ is a quintessential example of a strongly correlated bad-metal, which undergoes a metal-insulator transition (MIT) concomitant with a structural transition to a V-V dimerized monoclinic phase below T$_{MIT} \sim 340K$. It has been experimentally shown that one can control this transition by doping VO$_2$. In particular, doping with oxygen vacancies ($V_O$) has been shown… ▽ More Rutile ($R$) phase VO$_2$ is a quintessential example of a strongly correlated bad-metal, which undergoes a metal-insulator transition (MIT) concomitant with a structural transition to a V-V dimerized monoclinic phase below T$_{MIT} \sim 340K$. It has been experimentally shown that one can control this transition by doping VO$_2$. In particular, doping with oxygen vacancies ($V_O$) has been shown to completely suppress this MIT {\em without} any structural transition. We explain this suppression by elucidating the influence of oxygen-vacancies on the electronic-structure of the metallic $R$ phase VO$_2$, explicitly treating strong electron-electron correlations using dynamical mean-field theory (DMFT) as well as diffusion Monte Carlo (DMC) flavor of quantum Monte Carlo (QMC) techniques. We show that $V_O$'s tend to change the V-3$d$ filling away from its nominal half-filled value, with the $e_{g}^π$ orbitals competing with the otherwise dominant $a_{1g}$ orbital. Loss of this near orbital polarization of the $a_{1g}$ orbital is associated with a weakening of electron correlations, especially along the V-V dimerization direction. This removes a charge-density wave (CDW) instability along this direction above a critical doping concentration, which further suppresses the metal-insulator transition. Our study also suggests that the MIT is predominantly driven by a correlation-induced CDW instability along the V-V dimerization direction. △ Less

Submitted 2 November, 2018; originally announced November 2018.

Journal ref: Phys. Rev. B 101, 155129 (2020)

Towards Optimal Running Times for Optimal Transport

Authors: Jose Blanchet, Arun Jambulapati, Carson Kent, Aaron Sidford

Abstract: In this work, we provide faster algorithms for approximating the optimal transport distance, e.g. earth mover's distance, between two discrete probability distributions $μ, ν\in Δ^n$. Given a cost function $C : [n] \times [n] \to \mathbb{R}_{\geq 0}$ where $C(i,j) \leq 1$ quantifies the penalty of transporting a unit of mass from $i$ to $j$, we show how to compute a coupling $X$ between $r$ and… ▽ More In this work, we provide faster algorithms for approximating the optimal transport distance, e.g. earth mover's distance, between two discrete probability distributions $μ, ν\in Δ^n$. Given a cost function $C : [n] \times [n] \to \mathbb{R}_{\geq 0}$ where $C(i,j) \leq 1$ quantifies the penalty of transporting a unit of mass from $i$ to $j$, we show how to compute a coupling $X$ between $r$ and $c$ in time $\widetilde{O}\left(n^2 /ε\right)$ whose expected transportation cost is within an additive $ε$ of optimal. This improves upon the previously best known running time for this problem of $\widetilde{O}\left(\text{min}\left\{ n^{9/4}/ε, n^2/ε^2 \right\}\right)$. We achieve our results by providing reductions from optimal transport to canonical optimization problems for which recent algorithmic efforts have provided nearly-linear time algorithms. Leveraging nearly linear time algorithms for solving packing linear programs and for solving the matrix balancing problem, we obtain two separate proofs of our stated running time. Further, one of our algorithms is easily parallelized and can be implemented with depth $\widetilde{O}(1/ε)$. Moreover, we show that further algorithmic improvements to our result would be surprising in the sense that any improvement would yield an $o(n^{2.5})$ algorithm for \textit{maximum cardinality bipartite matching}, for which currently the only known algorithms for achieving such a result are based on fast-matrix multiplication. △ Less

Submitted 28 January, 2020; v1 submitted 17 October, 2018; originally announced October 2018.

Anomalous dielectric response at intermixed oxide heterointerfaces

Authors: Valentino R. Cooper, Houlong L. Zhuang, Lipeng Zhang, P. Ganesh, Haixuan Xu, Paul R. C. Kent

Abstract: Two-dimensional charge carrier accumulation at oxide heterointerfaces presents a paradigm shift for oxide electronics. Like a capacitor, interfacial charge buildup couples to an electric field across the dielectric medium. To prevent the so-called polar catastrophe, several charge screening mechanisms emerge, including polar distortions and interfacial intermixing which reduce the sharpness of the… ▽ More Two-dimensional charge carrier accumulation at oxide heterointerfaces presents a paradigm shift for oxide electronics. Like a capacitor, interfacial charge buildup couples to an electric field across the dielectric medium. To prevent the so-called polar catastrophe, several charge screening mechanisms emerge, including polar distortions and interfacial intermixing which reduce the sharpness of the interface. Here, we examine how atomic intermixing at oxide interfaces affect the balance between polar distortions and electric potential across the dielectric medium. We find that intermixing moves the peak charge distribution away from the oxide/oxide interface; thereby changing the direction of polar distortions away from this boundary with minimal effect on the electric field. This opposing electric field and polar distortions is equivalent to the transient phase transition tipping point observed in double well ferroelectrics; resulting in an anomalous dielectric response -- a possible signature of local negative differential capacitance, with implications for designing dissipationless oxide electronics. △ Less

Submitted 21 June, 2018; originally announced June 2018.

Oxygen vacancy formation energies in PbTiO$_3$/SrTiO$_3$ superlattice

Authors: Lipeng Zhang, Isaac Bredeson, Axiel Yaël Birenbaum, Paul R. C. Kent, Valentino R. Cooper, P. Ganesh, Haixuan Xu

Abstract: The defect stability in a prototypical perovskite oxide superlattice consisting of SrTiO$_3$ and PbTiO$_3$ (STO/PTO) is determined using first principles density functional theory calculations. Specifically, the oxygen vacancy formation energies E$_v$ in the paraelectric and ferroelectric phases of a superlattice with four atomic layers of STO and four layers of PTO (4STO/4PTO) are determined and… ▽ More The defect stability in a prototypical perovskite oxide superlattice consisting of SrTiO$_3$ and PbTiO$_3$ (STO/PTO) is determined using first principles density functional theory calculations. Specifically, the oxygen vacancy formation energies E$_v$ in the paraelectric and ferroelectric phases of a superlattice with four atomic layers of STO and four layers of PTO (4STO/4PTO) are determined and compared. The effects of charge state, octahedral rotation, polarization, and interfaces on the E$_v$ are examined. The formation energies vary layer-by-layer in the superlattices, with E$_v$ being higher in the ferroelectric phase than that in the paraelectric phase. The two interfaces constructed in these oxide superlattices, which are symmetrically equivalent in the paraelectric systems, exhibit very different formation energies in the ferroelectric superlattices and this can be seen to be driven by the coupling of ferroelectric and rotational modes. At equivalent lattice sites, E$_v$ of charged vacancies is generally lower than that of neutral vacancies. Octahedral rotations (a$^0$a$^0$c$^-$) in the FE superlattice have a significant effect on the E$_v$, increasing the formation energy of vacancies located near the interface but decreasing the formation energy of the oxygen vacancies located in the bulk-like regions of the STO and PTO constituent parts. The formation energy variations among different layers are found to be primarily caused by the difference in the local relaxation at each layer. These fundamental insights into the defect stability in perovskite superlattices can be used to tune defect properties via controlling the constituent materials of superlattices and interface engineering. △ Less

Submitted 29 May, 2018; v1 submitted 25 May, 2018; originally announced May 2018.

Comments: 20 pages, 5 figures, accepted at Physical Review Materials

Journal ref: Phys. Rev. Materials 2, 064409 (2018)

An efficient hybrid orbital representation for quantum Monte Carlo calculations

Authors: Ye Luo, Kenneth P. Esler, Paul R. C. Kent, Luke Shulenburger

Abstract: The scale and complexity of quantum system to which real-space quantum Monte Carlo (QMC) can be applied in part depends on the representation and memory usage of the trial wavefunction. B-splines, the computationally most efficient basis set, can have memory requirements exceeding the capacity of a single computational node. This situation has traditionally forced a difficult choice of either usin… ▽ More The scale and complexity of quantum system to which real-space quantum Monte Carlo (QMC) can be applied in part depends on the representation and memory usage of the trial wavefunction. B-splines, the computationally most efficient basis set, can have memory requirements exceeding the capacity of a single computational node. This situation has traditionally forced a difficult choice of either using slow internode communication or a potentially less accurate but smaller basis set such as Gaussians. Here, we introduce a hybrid representation of the single particle orbitals that combine a localized atomic basis set around atomic cores and B-splines in the interstitial regions to reduce the memory usage while retaining high speed of evaluation and either retaining or increasing overall accuracy. We present a benchmark calculation for NiO demonstrating a superior accuracy while using only one eighth the memory required for conventional B-splines. The hybrid orbital representation therefore expands the overall range of systems that can be practically studied with QMC. △ Less

Submitted 18 May, 2018; originally announced May 2018.

Journal ref: J. Chem. Phys. 149, 084107 (2018)

Asymptotic cones and boundaries of CAT(0) groups

Authors: Curtis Kent, Russell Ricks

Abstract: It is well known that the Tits boundary of a proper cocompact CAT(0) space embeds into every asymptotic cone of the space. We explore the relationships between the asymptotic cones of a CAT(0) space and its boundary under both the standard visual (i.e. cone) topology and the Tits metric. We show that the set of asymptotic cones of a proper cocompact CAT(0) space admits canonical connecting maps un… ▽ More It is well known that the Tits boundary of a proper cocompact CAT(0) space embeds into every asymptotic cone of the space. We explore the relationships between the asymptotic cones of a CAT(0) space and its boundary under both the standard visual (i.e. cone) topology and the Tits metric. We show that the set of asymptotic cones of a proper cocompact CAT(0) space admits canonical connecting maps under which the direct limit is isometric to the Euclidean cone on the Tits boundary. The resulting projection from any asymptotic cone to the Tits boundary is determined by the visual topology; on the other hand, the visual topology can be recovered from the connecting maps between asymptotic cones. We also demonstrate how maps between asymptotic cones induce maps between Tits boundaries. △ Less

Submitted 11 September, 2018; v1 submitted 17 May, 2018; originally announced May 2018.

Comments: Added a brief section on metric direct and inverse limits

A Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1075 additional authors not shown)

Abstract: The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational w… ▽ More The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually-unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generically-polarized gravitational waves. We find no evidence for a background of any polarization, and place the first direct bounds on the contributions of vector and scalar polarizations to the stochastic background. Under log-uniform priors for the energy in each polarization, we limit the energy-densities of tensor, vector, and scalar modes at 95% credibility to $Ω^T_0 < 5.6 \times 10^{-8}$, $Ω^V_0 < 6.4\times 10^{-8}$, and $Ω^S_0 < 1.1\times 10^{-7}$ at a reference frequency $f_0 = 25$ Hz. △ Less

Submitted 2 October, 2019; v1 submitted 27 February, 2018; originally announced February 2018.

Comments: Minor updates to match published version

Journal ref: Phys. Rev. Lett. 120, 201102 (2018)

QMCPACK : An open source ab initio Quantum Monte Carlo package for the electronic structure of atoms, molecules, and solids

Authors: Jeongnim Kim, Andrew Baczewski, Todd D. Beaudet, Anouar Benali, M. Chandler Bennett, Mark A. Berrill, Nick S. Blunt, Edgar Josue Landinez Borda, Michele Casula, David M. Ceperley, Simone Chiesa, Bryan K. Clark, Raymond C. Clay III, Kris T. Delaney, Mark Dewing, Kenneth P. Esler, Hongxia Hao, Olle Heinonen, Paul R. C. Kent, Jaron T. Krogel, Ilkka Kylanpaa, Ying Wai Li, M. Graham Lopez, Ye Luo, Fionn D. Malone , et al. (23 additional authors not shown)

Abstract: QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a s… ▽ More QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit (CPU) and graphical processing unit (GPU) systems. We detail the program's capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://www.qmcpack.org . △ Less

Submitted 4 April, 2018; v1 submitted 19 February, 2018; originally announced February 2018.

Journal ref: J. Phys.: Condens. Matter 30 195901 (2018)

Full Band All-sky Search for Periodic Gravitational Waves in the O1 LIGO Data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, B. Allen, G. Allen, A. Allocca, P. A. Altin , et al. (1077 additional authors not shown)

Abstract: We report on a new all-sky search for periodic gravitational waves in the frequency band 475-2000 Hz and with a frequency time derivative in the range of [-1.0e-8, +1e-9] Hz/s. Potential signals could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGO's first observational run O1. No gravitational wave… ▽ More We report on a new all-sky search for periodic gravitational waves in the frequency band 475-2000 Hz and with a frequency time derivative in the range of [-1.0e-8, +1e-9] Hz/s. Potential signals could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGO's first observational run O1. No gravitational wave signals were observed, and upper limits were placed on their strengths. For completeness, results from the separately published low frequency search 20-475 Hz are included as well. Our lowest upper limit on worst-case (linearly polarized) strain amplitude h_0 is 4e-25 near 170 Hz, while at the high end of our frequency range we achieve a worst-case upper limit of 1.3e-24. For a circularly polarized source (most favorable orientation), the smallest upper limit obtained is ~1.5e-25. △ Less

Submitted 14 February, 2018; originally announced February 2018.

Journal ref: Phys. Rev. D 97, 102003 (2018)

Analysis of scalar and fermion quantum field theory on anti-de Sitter space-time

Authors: Victor E. Ambrus, Carl Kent, Elizabeth Winstanley

Abstract: We study vacuum and thermal expectation values of quantum scalar and Dirac fermion fields on anti-de Sitter space-time. Anti-de Sitter space-time is maximally symmetric and this enables expressions for the scalar and fermion vacuum Feynman Green's functions to be derived in closed form. We employ Hadamard renormalization to find the vacuum expectation values. The thermal Feynman Green's functions… ▽ More We study vacuum and thermal expectation values of quantum scalar and Dirac fermion fields on anti-de Sitter space-time. Anti-de Sitter space-time is maximally symmetric and this enables expressions for the scalar and fermion vacuum Feynman Green's functions to be derived in closed form. We employ Hadamard renormalization to find the vacuum expectation values. The thermal Feynman Green's functions are constructed from the vacuum Feynman Green's functions using the imaginary time periodicity/anti-periodicity property for scalars/fermions. Focussing on massless fields with either conformal or minimal coupling to the space-time curvature (these two cases being the same for fermions) we compute the differences between the thermal and vacuum expectation values. We compare the resulting energy densities, pressures and pressure deviators with the corresponding classical quantities calculated using relativistic kinetic theory. △ Less

Submitted 30 April, 2018; v1 submitted 31 January, 2018; originally announced January 2018.

Comments: 30 pages, 10 figures, minor changes, to appear in the Proceedings of the Fourth Amazonian Symposium on Physics, Belem, Brazil, 18-22 September 2017

Constraints on cosmic strings using data from the first Advanced LIGO observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1020 additional authors not shown)

Abstract: Cosmic strings are topological defects which can be formed in GUT-scale phase transitions in the early universe. They are also predicted to form in the context of string theory. The main mechanism for a network of Nambu-Goto cosmic strings to lose energy is through the production of loops and the subsequent emission of gravitational waves, thus offering an experimental signature for the existence… ▽ More Cosmic strings are topological defects which can be formed in GUT-scale phase transitions in the early universe. They are also predicted to form in the context of string theory. The main mechanism for a network of Nambu-Goto cosmic strings to lose energy is through the production of loops and the subsequent emission of gravitational waves, thus offering an experimental signature for the existence of cosmic strings. Here we report on the analysis conducted to specifically search for gravitational-wave bursts from cosmic string loops in the data of Advanced LIGO 2015-2016 observing run (O1). No evidence of such signals was found in the data, and as a result we set upper limits on the cosmic string parameters for three recent loop distribution models. In this paper, we initially derive constraints on the string tension $Gμ$ and the intercommutation probability, using not only the burst analysis performed on the O1 data set, but also results from the previously published LIGO stochastic O1 analysis, pulsar timing arrays, cosmic microwave background and Big-Bang nucleosynthesis experiments. We show that these data sets are complementary in that they probe gravitational waves produced by cosmic string loops during very different epochs. Finally, we show that the data sets exclude large parts of the parameter space of the three loop distribution models we consider. △ Less

Submitted 2 May, 2018; v1 submitted 4 December, 2017; originally announced December 2017.

Comments: Physical Review D, in-press

Journal ref: Phys. Rev. D 97, 102002 (2018)

GW170608: Observation of a 19-solar-mass Binary Black Hole Coalescence

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1079 additional authors not shown)

Abstract: On June 8, 2017 at 02:01:16.49 UTC, a gravitational-wave signal from the merger of two stellar-mass black holes was observed by the two Advanced LIGO detectors with a network signal-to-noise ratio of 13. This system is the lightest black hole binary so far observed, with component masses $12^{+7}_{-2}\,M_\odot$ and $7^{+2}_{-2}\,M_\odot$ (90% credible intervals). These lie in the range of measured… ▽ More On June 8, 2017 at 02:01:16.49 UTC, a gravitational-wave signal from the merger of two stellar-mass black holes was observed by the two Advanced LIGO detectors with a network signal-to-noise ratio of 13. This system is the lightest black hole binary so far observed, with component masses $12^{+7}_{-2}\,M_\odot$ and $7^{+2}_{-2}\,M_\odot$ (90% credible intervals). These lie in the range of measured black hole masses in low-mass X-ray binaries, thus allowing us to compare black holes detected through gravitational waves with electromagnetic observations. The source's luminosity distance is $340^{+140}_{-140}$ Mpc, corresponding to redshift $0.07^{+0.03}_{-0.03}$. We verify that the signal waveform is consistent with the predictions of general relativity. △ Less

Submitted 15 November, 2017; originally announced November 2017.

Report number: LIGO Document P170608-v8

Search for post-merger gravitational waves from the remnant of the binary neutron star merger GW170817

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1083 additional authors not shown)

Abstract: The first observation of a binary neutron star coalescence by the Advanced LIGO and Advanced Virgo gravitational-wave detectors offers an unprecedented opportunity to study matter under the most extreme conditions. After such a merger, a compact remnant is left over whose nature depends primarily on the masses of the inspiralling objects and on the equation of state of nuclear matter. This could b… ▽ More The first observation of a binary neutron star coalescence by the Advanced LIGO and Advanced Virgo gravitational-wave detectors offers an unprecedented opportunity to study matter under the most extreme conditions. After such a merger, a compact remnant is left over whose nature depends primarily on the masses of the inspiralling objects and on the equation of state of nuclear matter. This could be either a black hole or a neutron star (NS), with the latter being either long-lived or too massive for stability implying delayed collapse to a black hole. Here, we present a search for gravitational waves from the remnant of the binary neutron star merger GW170817 using data from Advanced LIGO and Advanced Virgo. We search for short ($\lesssim1$ s) and intermediate-duration ($\lesssim 500$ s) signals, which includes gravitational-wave emission from a hypermassive NS or supramassive NS, respectively. We find no signal from the post-merger remnant. Our derived strain upper limits are more than an order of magnitude larger than those predicted by most models. For short signals, our best upper limit on the root-sum-square of the gravitational-wave strain emitted from 1--4 kHz is $h_{\rm rss}^{50\%}=2.1\times 10^{-22}$ Hz$^{-1/2}$ at 50% detection efficiency. For intermediate-duration signals, our best upper limit at 50% detection efficiency is $h_{\rm rss}^{50\%}=8.4\times 10^{-22}$ Hz$^{-1/2}$ for a millisecond magnetar model, and $h_{\rm rss}^{50\%}=5.9\times 10^{-22}$ Hz$^{-1/2}$ for a bar-mode model. These results indicate that post-merger emission from a similar event may be detectable when advanced detectors reach design sensitivity or with next-generation detectors. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Report number: LIGO-P1700318

Journal ref: ApJL, 851:L16 (2017)

Search for High-energy Neutrinos from Binary Neutron Star Merger GW170817 with ANTARES, IceCube, and the Pierre Auger Observatory

Authors: A. Albert, M. Andre, M. Anghinolfi, M. Ardid, J. -J. Aubert, J. Aublin, T. Avgitas, B. Baret, J. Barrios-Marti, S. Basa, B. Belhorma, V. Bertin, S. Biagi, R. Bormuth, S. Bourret, M. C. Bouwhuis, H. Branzacs, R. Bruijn, J. Brunner, J. Busto, A. Capone, L. Caramete, J. Carr, S. Celli, R. Cherkaoui El Moursli , et al. (1916 additional authors not shown)

Abstract: The Advanced LIGO and Advanced Virgo observatories recently discovered gravitational waves from a binary neutron star inspiral. A short gamma-ray burst (GRB) that followed the merger of this binary was also recorded by the Fermi Gamma-ray Burst Monitor (Fermi-GBM), and the Anticoincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), indicating par… ▽ More The Advanced LIGO and Advanced Virgo observatories recently discovered gravitational waves from a binary neutron star inspiral. A short gamma-ray burst (GRB) that followed the merger of this binary was also recorded by the Fermi Gamma-ray Burst Monitor (Fermi-GBM), and the Anticoincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), indicating particle acceleration by the source. The precise location of the event was determined by optical detections of emission following the merger. We searched for high-energy neutrinos from the merger in the GeV--EeV energy range using the ANTARES, IceCube, and Pierre Auger Observatories. No neutrinos directionally coincident with the source were detected within $\pm500$ s around the merger time. Additionally, no MeV neutrino burst signal was detected coincident with the merger. We further carried out an extended search in the direction of the source for high-energy neutrinos within the 14-day period following the merger, but found no evidence of emission. We used these results to probe dissipation mechanisms in relativistic outflows driven by the binary neutron star merger. The non-detection is consistent with model predictions of short GRBs observed at a large off-axis angle. △ Less

Submitted 9 November, 2017; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 22 pages, 2 figures

Report number: LIGO-P1700344

On the Progenitor of Binary Neutron Star Merger GW170817

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1073 additional authors not shown)

Abstract: On 2017 August 17 the merger of two compact objects with masses consistent with two neutron stars was discovered through gravitational-wave (GW170817), gamma-ray (GRB 170817A), and optical (SSS17a/AT 2017gfo) observations. The optical source was associated with the early-type galaxy NGC 4993 at a distance of just $\sim$40 Mpc, consistent with the gravitational-wave measurement, and the merger was… ▽ More On 2017 August 17 the merger of two compact objects with masses consistent with two neutron stars was discovered through gravitational-wave (GW170817), gamma-ray (GRB 170817A), and optical (SSS17a/AT 2017gfo) observations. The optical source was associated with the early-type galaxy NGC 4993 at a distance of just $\sim$40 Mpc, consistent with the gravitational-wave measurement, and the merger was localized to be at a projected distance of $\sim$2 kpc away from the galaxy's center. We use this minimal set of facts and the mass posteriors of the two neutron stars to derive the first constraints on the progenitor of GW170817 at the time of the second supernova (SN). We generate simulated progenitor populations and follow the three-dimensional kinematic evolution from the binary neutron star (BNS) birth to the merger time, accounting for pre-SN galactic motion, for considerably different input distributions of the progenitor mass, pre-SN semimajor axis, and SN-kick velocity. Though not considerably tight, we find these constraints to be comparable to those for Galactic BNS progenitors. The derived constraints are very strongly influenced by the requirement of keeping the binary bound after the second SN and having the merger occur relatively close to the center of the galaxy. These constraints are insensitive to the galaxy's star formation history, provided the stellar populations are older than 1 Gyr. △ Less

Submitted 1 December, 2017; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 18 pages, 8 figures, 2 tables

Report number: LIGO-P1700264

Journal ref: ApJL, 850, L40, 2017

GW170817: Implications for the Stochastic Gravitational-Wave Background from Compact Binary Coalescences

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1077 additional authors not shown)

Abstract: The LIGO Scientific and Virgo Collaborations have announced the first detection of gravitational waves from the coalescence of two neutron stars. The merger rate of binary neutron stars estimated from this event suggests that distant, unresolvable binary neutron stars create a significant astrophysical stochastic gravitational-wave background. The binary neutron star background will add to the bac… ▽ More The LIGO Scientific and Virgo Collaborations have announced the first detection of gravitational waves from the coalescence of two neutron stars. The merger rate of binary neutron stars estimated from this event suggests that distant, unresolvable binary neutron stars create a significant astrophysical stochastic gravitational-wave background. The binary neutron star background will add to the background from binary black holes, increasing the amplitude of the total astrophysical background relative to previous expectations. In the Advanced LIGO-Virgo frequency band most sensitive to stochastic backgrounds (near 25 Hz), we predict a total astrophysical background with amplitude $Ω_{\rm GW} (f=25 \text{Hz}) = 1.8_{-1.3}^{+2.7} \times 10^{-9}$ with $90\%$ confidence, compared with $Ω_{\rm GW} (f=25 \text{Hz}) = 1.1_{-0.7}^{+1.2} \times 10^{-9}$ from binary black holes alone. Assuming the most probable rate for compact binary mergers, we find that the total background may be detectable with a signal-to-noise-ratio of 3 after 40 months of total observation time, based on the expected timeline for Advanced LIGO and Virgo to reach their design sensitivity. △ Less

Submitted 30 September, 2019; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 12 pages, 3 figures

Report number: LIGO P1700272

Journal ref: Phys. Rev. Lett. 120, 091101 (2018)

Estimating the Contribution of Dynamical Ejecta in the Kilonova Associated with GW170817

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1078 additional authors not shown)

Abstract: The source of the gravitational-wave signal GW170817, very likely a binary neutron star merger, was also observed electromagnetically, providing the first multi-messenger observations of this type. The two week long electromagnetic counterpart had a signature indicative of an r-process-induced optical transient known as a kilonova. This Letter examines how the mass of the dynamical ejecta can be e… ▽ More The source of the gravitational-wave signal GW170817, very likely a binary neutron star merger, was also observed electromagnetically, providing the first multi-messenger observations of this type. The two week long electromagnetic counterpart had a signature indicative of an r-process-induced optical transient known as a kilonova. This Letter examines how the mass of the dynamical ejecta can be estimated without a direct electromagnetic observation of the kilonova, using gravitational-wave measurements and a phenomenological model calibrated to numerical simulations of mergers with dynamical ejecta. Specifically, we apply the model to the binary masses inferred from the gravitational-wave measurements, and use the resulting mass of the dynamical ejecta to estimate its contribution (without the effects of wind ejecta) to the corresponding kilonova light curves from various models. The distributions of dynamical ejecta mass range between $M_{ej} = 10^{-3} - 10^{-2} M_{\odot}$ for various equations of state, assuming the neutron stars are rotating slowly. In addition, we use our estimates of the dynamical ejecta mass and the neutron star merger rates inferred from GW170817 to constrain the contribution of events like this to the r-process element abundance in the Galaxy when ejecta mass from post-merger winds is neglected. We find that if $\gtrsim10\%$ of the matter dynamically ejected from BNS mergers is converted to r-process elements, GW170817-like BNS mergers could fully account for the amount of r-process material observed in the Milky Way. △ Less

Submitted 8 January, 2018; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 18 pages, 5 figures, accepted for publication in ApJ Letters

Report number: LIGO-P1700309

A gravitational-wave standard siren measurement of the Hubble constant

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca, P. A. Altin, A. Amato , et al. (1289 additional authors not shown)

Abstract: The detection of GW170817 in both gravitational waves and electromagnetic waves heralds the age of gravitational-wave multi-messenger astronomy. On 17 August 2017 the Advanced LIGO and Virgo detectors observed GW170817, a strong signal from the merger of a binary neutron-star system. Less than 2 seconds after the merger, a gamma-ray burst (GRB 170817A) was detected within a region of the sky consi… ▽ More The detection of GW170817 in both gravitational waves and electromagnetic waves heralds the age of gravitational-wave multi-messenger astronomy. On 17 August 2017 the Advanced LIGO and Virgo detectors observed GW170817, a strong signal from the merger of a binary neutron-star system. Less than 2 seconds after the merger, a gamma-ray burst (GRB 170817A) was detected within a region of the sky consistent with the LIGO-Virgo-derived location of the gravitational-wave source. This sky region was subsequently observed by optical astronomy facilities, resulting in the identification of an optical transient signal within $\sim 10$ arcsec of the galaxy NGC 4993. These multi-messenger observations allow us to use GW170817 as a standard siren, the gravitational-wave analog of an astronomical standard candle, to measure the Hubble constant. This quantity, which represents the local expansion rate of the Universe, sets the overall scale of the Universe and is of fundamental importance to cosmology. Our measurement combines the distance to the source inferred purely from the gravitational-wave signal with the recession velocity inferred from measurements of the redshift using electromagnetic data. This approach does not require any form of cosmic "distance ladder;" the gravitational wave analysis can be used to estimate the luminosity distance out to cosmological scales directly, without the use of intermediate astronomical distance measurements. We determine the Hubble constant to be $70.0^{+12.0}_{-8.0} \, \mathrm{km} \, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$ (maximum a posteriori and 68% credible interval). This is consistent with existing measurements, while being completely independent of them. Additional standard-siren measurements from future gravitational-wave sources will provide precision constraints of this important cosmological parameter. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: 26 pages, 5 figures, Nature in press. For more information see https://dcc.ligo.org/LIGO-P1700296/public

Report number: LIGO P1700296

First narrow-band search for continuous gravitational waves from known pulsars in advanced detector data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, B. Allen, G. Allen, A. Allocca, P. A. Altin , et al. (1074 additional authors not shown)

Abstract: Spinning neutron stars asymmetric with respect to their rotation axis are potential sources of continuous gravitational waves for ground-based interferometric detectors. In the case of known pulsars a fully coherent search, based on matched filtering, which uses the position and rotational parameters obtained from electromagnetic observations, can be carried out. Matched filtering maximizes the si… ▽ More Spinning neutron stars asymmetric with respect to their rotation axis are potential sources of continuous gravitational waves for ground-based interferometric detectors. In the case of known pulsars a fully coherent search, based on matched filtering, which uses the position and rotational parameters obtained from electromagnetic observations, can be carried out. Matched filtering maximizes the signal-to-noise (SNR) ratio, but a large sensitivity loss is expected in case of even a very small mismatch between the assumed and the true signal parameters. For this reason, {\it narrow-band} analyses methods have been developed, allowing a fully coherent search for gravitational waves from known pulsars over a fraction of a hertz and several spin-down values. In this paper we describe a narrow-band search of eleven pulsars using data from Advanced LIGO's first observing run. Although we have found several initial outliers, further studies show no significant evidence for the presence of a gravitational wave signal. Finally, we have placed upper limits on the signal strain amplitude lower than the spin-down limit for 5 of the 11 targets over the bands searched: in the case of J1813-1749 the spin-down limit has been beaten for the first time. For an additional 3 targets, the median upper limit across the search bands is below the spin-down limit. This is the most sensitive narrow-band search for continuous gravitational waves carried out so far. △ Less

Submitted 5 December, 2017; v1 submitted 6 October, 2017; originally announced October 2017.

Comments: 9 Figures, 7 tables, submitted to PRD

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

GW170814: A Three-Detector Observation of Gravitational Waves from a Binary Black Hole Coalescence

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1085 additional authors not shown)

Abstract: On August 14, 2017 at 10:30:43 UTC, the Advanced Virgo detector and the two Advanced LIGO detectors coherently observed a transient gravitational-wave signal produced by the coalescence of two stellar mass black holes, with a false-alarm-rate of $\lesssim$ 1 in 27000 years. The signal was observed with a three-detector network matched-filter signal-to-noise ratio of 18. The inferred masses of the… ▽ More On August 14, 2017 at 10:30:43 UTC, the Advanced Virgo detector and the two Advanced LIGO detectors coherently observed a transient gravitational-wave signal produced by the coalescence of two stellar mass black holes, with a false-alarm-rate of $\lesssim$ 1 in 27000 years. The signal was observed with a three-detector network matched-filter signal-to-noise ratio of 18. The inferred masses of the initial black holes are $30.5_{-3.0}^{+5.7}$ Msun and $25.3_{-4.2}^{+2.8}$ Msun (at the 90% credible level). The luminosity distance of the source is $540_{-210}^{+130}~\mathrm{Mpc}$, corresponding to a redshift of $z=0.11_{-0.04}^{+0.03}$. A network of three detectors improves the sky localization of the source, reducing the area of the 90% credible region from 1160 deg$^2$ using only the two LIGO detectors to 60 deg$^2$ using all three detectors. For the first time, we can test the nature of gravitational wave polarizations from the antenna response of the LIGO-Virgo network, thus enabling a new class of phenomenological tests of gravity. △ Less

Submitted 13 October, 2017; v1 submitted 27 September, 2017; originally announced September 2017.

Journal ref: Phys. Rev. Lett. 119, 141101 (2017)

Homotopy type of planar continua

Authors: Curtis Kent

Abstract: We prove that the homotopy type of a map from a Peano continuum into a planar or one-dimensional space is determined by the induced homomorphism of fundamental groups. This provides a new proof that planar sets are aspherical and is used to prove that two planar or one-dimensional Peano continua are homotopy equivalent if and only if they have isomorphic fundamental groups. We prove that the homotopy type of a map from a Peano continuum into a planar or one-dimensional space is determined by the induced homomorphism of fundamental groups. This provides a new proof that planar sets are aspherical and is used to prove that two planar or one-dimensional Peano continua are homotopy equivalent if and only if they have isomorphic fundamental groups. △ Less

Submitted 26 September, 2017; originally announced September 2017.

Comments: 18 pages. Comments welcome

First search for nontensorial gravitational waves from known pulsars

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, P. A. Altin , et al. (1028 additional authors not shown)

Abstract: We present results from the first directed search for nontensorial gravitational waves. While general relativity allows for tensorial (plus and cross) modes only, a generic metric theory may, in principle, predict waves with up to six different polarizations. This analysis is sensitive to continuous signals of scalar, vector or tensor polarizations, and does not rely on any specific theory of grav… ▽ More We present results from the first directed search for nontensorial gravitational waves. While general relativity allows for tensorial (plus and cross) modes only, a generic metric theory may, in principle, predict waves with up to six different polarizations. This analysis is sensitive to continuous signals of scalar, vector or tensor polarizations, and does not rely on any specific theory of gravity. After searching data from the first observation run of the advanced LIGO detectors for signals at twice the rotational frequency of 200 known pulsars, we find no evidence of gravitational waves of any polarization. We report the first upper limits for scalar and vector strains, finding values comparable in magnitude to previously-published limits for tensor strain. Our results may be translated into constraints on specific alternative theories of gravity. △ Less

Submitted 16 October, 2019; v1 submitted 26 September, 2017; originally announced September 2017.

Comments: journal version

Report number: LIGO-P1700009

Journal ref: Phys. Rev. Lett. 120, 031104 (2018)

Delayed Slater determinant update algorithms for high efficiency quantum Monte Carlo

Authors: T. McDaniel, E. F. D'Azevedo, Y. W. Li, K. Wong, P. R. C. Kent

Abstract: Within ab initio Quantum Monte Carlo simulations, the leading numerical cost for large systems is the computation of the values of the Slater determinants in the trial wavefunction. Each Monte Carlo step requires finding the determinant of a dense matrix. This is most commonly iteratively evaluated using a rank-1 Sherman-Morrison updating scheme to avoid repeated explicit calculation of the invers… ▽ More Within ab initio Quantum Monte Carlo simulations, the leading numerical cost for large systems is the computation of the values of the Slater determinants in the trial wavefunction. Each Monte Carlo step requires finding the determinant of a dense matrix. This is most commonly iteratively evaluated using a rank-1 Sherman-Morrison updating scheme to avoid repeated explicit calculation of the inverse. The overall computational cost is therefore formally cubic in the number of electrons or matrix size. To improve the numerical efficiency of this procedure, we propose a novel multiple rank delayed update scheme. This strategy enables probability evaluation with application of accepted moves to the matrices delayed until after a predetermined number of moves, K. The accepted events are then applied to the matrices en bloc with enhanced arithmetic intensity and computational efficiency via matrix-matrix operations instead of matrix-vector operations. This procedure does not change the underlying Monte Carlo sampling or its statistical efficiency. For calculations on large systems and algorithms such as diffusion Monte Carlo where the acceptance ratio is high, order of magnitude improvements in the update time can be obtained on both multi-core CPUs and GPUs. △ Less

Submitted 2 August, 2017; originally announced August 2017.

First low-frequency Einstein@Home all-sky search for continuous gravitational waves in Advanced LIGO data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, B. Allen, G. Allen, A. Allocca, P. A. Altin , et al. (1017 additional authors not shown)

Abstract: We report results of a deep all-sky search for periodic gravitational waves from isolated neutron stars in data from the first Advanced LIGO observing run. This search investigates the low frequency range of Advanced LIGO data, between 20 and 100 Hz, much of which was not explored in initial LIGO. The search was made possible by the computing power provided by the volunteers of the Einstein@Home p… ▽ More We report results of a deep all-sky search for periodic gravitational waves from isolated neutron stars in data from the first Advanced LIGO observing run. This search investigates the low frequency range of Advanced LIGO data, between 20 and 100 Hz, much of which was not explored in initial LIGO. The search was made possible by the computing power provided by the volunteers of the Einstein@Home project. We find no significant signal candidate and set the most stringent upper limits to date on the amplitude of gravitational wave signals from the target population, corresponding to a sensitivity depth of 48.7 [1/$\sqrt{\textrm{Hz}}$]. At the frequency of best strain sensitivity, near 100 Hz, we set 90% confidence upper limits of $1.8 \times 10^{-25}$. At the low end of our frequency range, 20 Hz, we achieve upper limits of $3.9 \times 10^{-24}$. At 55 Hz we can exclude sources with ellipticities greater than $10^{-5}$ within 100 pc of Earth with fiducial value of the principal moment of inertia of $10^{38} \textrm{kg m}^2$. △ Less

Submitted 14 July, 2017; v1 submitted 9 July, 2017; originally announced July 2017.

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

All-sky Search for Periodic Gravitational Waves in the O1 LIGO Data

Authors: LIGO Scientific Collaboration, Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca, P. A. Altin , et al. (1020 additional authors not shown)

Abstract: We report on an all-sky search for periodic gravitational waves in the frequency band 20-475 Hz and with a frequency time derivative in the range of [-1.0, +0.1]e-8 Hz/s. Such a signal could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGO's first observational run, O1. No periodic gravitational wave si… ▽ More We report on an all-sky search for periodic gravitational waves in the frequency band 20-475 Hz and with a frequency time derivative in the range of [-1.0, +0.1]e-8 Hz/s. Such a signal could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGO's first observational run, O1. No periodic gravitational wave signals were observed, and upper limits were placed on their strengths. The lowest upper limits on worst-case (linearly polarized) strain amplitude h0 are 4e-25 near 170 Hz. For a circularly polarized source (most favorable orientation), the smallest upper limits obtained are 1.5e-25. These upper limits refer to all sky locations and the entire range of frequency derivative values. For a population-averaged ensemble of sky locations and stellar orientations, the lowest upper limits obtained for the strain amplitude are 2.5e-25. △ Less

Submitted 15 July, 2017; v1 submitted 9 July, 2017; originally announced July 2017.

Comments: Updated reference to just arXiv'ed Einstein@Home paper, fix e-mail in arXiv metadata

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

Upper Limits on Gravitational Waves from Scorpius X-1 from a Model-Based Cross-Correlation Search in Advanced LIGO Data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1024 additional authors not shown)

Abstract: We present the results of a semicoherent search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1, using data from the first Advanced LIGO observing run. The search method uses details of the modelled, parametrized continuous signal to combine coherently data separated by less than a specified coherence time, which can be adjusted to trade off sensitivity against compu… ▽ More We present the results of a semicoherent search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1, using data from the first Advanced LIGO observing run. The search method uses details of the modelled, parametrized continuous signal to combine coherently data separated by less than a specified coherence time, which can be adjusted to trade off sensitivity against computational cost. A search was conducted over the frequency range from 25 Hz to 2000 Hz, spanning the current observationally-constrained range of the binary orbital parameters. No significant detection candidates were found, and frequency-dependent upper limits were set using a combination of sensitivity estimates and simulated signal injections. The most stringent upper limit was set at 175 Hz, with comparable limits set across the most sensitive frequency range from 100 Hz to 200 Hz. At this frequency, the 95 pct upper limit on signal amplitude h0 is 2.3e-25 marginalized over the unknown inclination angle of the neutron star's spin, and 8.03e-26 assuming the best orientation (which results in circularly polarized gravitational waves). These limits are a factor of 3-4 stronger than those set by other analyses of the same data, and a factor of about 7 stronger than the best upper limits set using initial LIGO data. In the vicinity of 100 Hz, the limits are a factor of between 1.2 and 3.5 above the predictions of the torque balance model, depending on inclination angle, if the most likely inclination angle of 44 degrees is assumed, they are within a factor of 1.7. △ Less

Submitted 16 November, 2019; v1 submitted 9 June, 2017; originally announced June 2017.

Comments: 19 pages, 8 figures. Formatted with AASTeX 6.1. Published in The Astrophysical Journal

Report number: LIGO-P1600297

Journal ref: The Astrophysical Journal, 847:47 (14pp), 2017 September 20