-
Smoothing the Edges: Smooth Optimization for Sparse Regularization using Hadamard Overparametrization
Authors:
Chris Kolb,
Christian L. Müller,
Bernd Bischl,
David Rügamer
Abstract:
We present a framework for smooth optimization of explicitly regularized objectives for (structured) sparsity. These non-smooth and possibly non-convex problems typically rely on solvers tailored to specific models and regularizers. In contrast, our method enables fully differentiable and approximation-free optimization and is thus compatible with the ubiquitous gradient descent paradigm in deep l…
▽ More
We present a framework for smooth optimization of explicitly regularized objectives for (structured) sparsity. These non-smooth and possibly non-convex problems typically rely on solvers tailored to specific models and regularizers. In contrast, our method enables fully differentiable and approximation-free optimization and is thus compatible with the ubiquitous gradient descent paradigm in deep learning. The proposed optimization transfer comprises an overparameterization of selected parameters and a change of penalties. In the overparametrized problem, smooth surrogate regularization induces non-smooth, sparse regularization in the base parametrization. We prove that the surrogate objective is equivalent in the sense that it not only has identical global minima but also matching local minima, thereby avoiding the introduction of spurious solutions. Additionally, our theory establishes results of independent interest regarding matching local minima for arbitrary, potentially unregularized, objectives. We comprehensively review sparsity-inducing parametrizations across different fields that are covered by our general theory, extend their scope, and propose improvements in several aspects. Numerical experiments further demonstrate the correctness and effectiveness of our approach on several sparse learning problems ranging from high-dimensional regression to sparse neural network training.
△ Less
Submitted 26 April, 2024; v1 submitted 7 July, 2023;
originally announced July 2023.
-
Factorized Structured Regression for Large-Scale Varying Coefficient Models
Authors:
David Rügamer,
Andreas Bender,
Simon Wiegrebe,
Daniel Racek,
Bernd Bischl,
Christian L. Müller,
Clemens Stachl
Abstract:
Recommender Systems (RS) pervade many aspects of our everyday digital life. Proposed to work at scale, state-of-the-art RS allow the modeling of thousands of interactions and facilitate highly individualized recommendations. Conceptually, many RS can be viewed as instances of statistical regression models that incorporate complex feature effects and potentially non-Gaussian outcomes. Such structur…
▽ More
Recommender Systems (RS) pervade many aspects of our everyday digital life. Proposed to work at scale, state-of-the-art RS allow the modeling of thousands of interactions and facilitate highly individualized recommendations. Conceptually, many RS can be viewed as instances of statistical regression models that incorporate complex feature effects and potentially non-Gaussian outcomes. Such structured regression models, including time-aware varying coefficients models, are, however, limited in their applicability to categorical effects and inclusion of a large number of interactions. Here, we propose Factorized Structured Regression (FaStR) for scalable varying coefficient models. FaStR overcomes limitations of general regression models for large-scale data by combining structured additive regression and factorization approaches in a neural network-based model implementation. This fusion provides a scalable framework for the estimation of statistical models in previously infeasible data settings. Empirical results confirm that the estimation of varying coefficients of our approach is on par with state-of-the-art regression techniques, while scaling notably better and also being competitive with other time-aware RS in terms of prediction performance. We illustrate FaStR's performance and interpretability on a large-scale behavioral study with smartphone user data.
△ Less
Submitted 25 May, 2022;
originally announced May 2022.
-
Objective hearing threshold identification from auditory brainstem response measurements using supervised and self-supervised approaches
Authors:
Dominik Thalmeier,
Gregor Miller,
Elida Schneltzer,
Anja Hurt,
Martin Hrabě de Angelis,
Lore Becker,
Christian L. Müller,
Holger Maier
Abstract:
Hearing loss is a major health problem and psychological burden in humans. Mouse models offer a possibility to elucidate genes involved in the underlying developmental and pathophysiological mechanisms of hearing impairment. To this end, large-scale mouse phenotyping programs include auditory phenotyping of single-gene knockout mouse lines. Using the auditory brainstem response (ABR) procedure, th…
▽ More
Hearing loss is a major health problem and psychological burden in humans. Mouse models offer a possibility to elucidate genes involved in the underlying developmental and pathophysiological mechanisms of hearing impairment. To this end, large-scale mouse phenotyping programs include auditory phenotyping of single-gene knockout mouse lines. Using the auditory brainstem response (ABR) procedure, the German Mouse Clinic and similar facilities worldwide have produced large, uniform data sets of averaged ABR raw data of mutant and wildtype mice. In the course of standard ABR analysis, hearing thresholds are assessed visually by trained staff from series of signal curves of increasing sound pressure level. This is time-consuming and prone to be biased by the reader as well as the graphical display quality and scale. In an attempt to reduce workload and improve quality and reproducibility, we developed and compared two methods for automated hearing threshold identification from averaged ABR raw data: a supervised approach involving two combined neural networks trained on human-generated labels and a self-supervised approach, which exploits the signal power spectrum and combines random forest sound level estimation with a piece-wise curve fitting algorithm for threshold finding. We show that both models work well, outperform human threshold detection, and are suitable for fast, reliable, and unbiased hearing threshold detection and quality control. In a high-throughput mouse phenotyping environment, both methods perform well as part of an automated end-to-end screening pipeline to detect candidate genes for hearing involvement. Code for both models as well as data used for this work are freely available.
△ Less
Submitted 16 December, 2021;
originally announced December 2021.
-
Ocean Mover's Distance: Using Optimal Transport for Analyzing Oceanographic Data
Authors:
Sangwon Hyun,
Aditya Mishra,
Christopher L. Follett,
Bror Jonsson,
Gemma Kulk,
Gael Forget,
Marie-Fanny Racault,
Thomas Jackson,
Stephanie Dutkiewicz,
Christian L. Müller,
Jacob Bien
Abstract:
Remote sensing observations from satellites and global biogeochemical models have combined to revolutionize the study of ocean biogeochemical cycling, but comparing the two data streams to each other and across time remains challenging due to the strong spatial-temporal structuring of the ocean. Here, we show that the Wasserstein distance provides a powerful metric for harnessing these structured…
▽ More
Remote sensing observations from satellites and global biogeochemical models have combined to revolutionize the study of ocean biogeochemical cycling, but comparing the two data streams to each other and across time remains challenging due to the strong spatial-temporal structuring of the ocean. Here, we show that the Wasserstein distance provides a powerful metric for harnessing these structured datasets for better marine ecosystem and climate predictions. Wasserstein distance complements commonly used point-wise difference methods such as the root mean squared error, by quantifying differences in terms of spatial displacement in addition to magnitude. As a test case we consider Chlorophyll (a key indicator of phytoplankton biomass) in the North-East Pacific Ocean, obtained from model simulations, in situ measurements, and satellite observations. We focus on two main applications: 1) Comparing model predictions with satellite observations, and 2) temporal evolution of Chlorophyll both seasonally and over longer time frames. Wasserstein distance successfully isolates temporal and depth variability and quantifies shifts in biogeochemical province boundaries. It also exposes relevant temporal trends in satellite Chlorophyll consistent with climate change predictions. Our study shows that optimal transport vectors underlying Wasserstein distance provide a novel visualization tool for testing models and better understanding temporal dynamics in the ocean.
△ Less
Submitted 4 November, 2022; v1 submitted 16 November, 2021;
originally announced November 2021.
-
GGLasso -- a Python package for General Graphical Lasso computation
Authors:
Fabian Schaipp,
Christian L. Müller,
Oleg Vlasovets
Abstract:
We introduce GGLasso, a Python package for solving General Graphical Lasso problems. The Graphical Lasso scheme, introduced by (Friedman 2007) (see also (Yuan 2007; Banerjee 2008)), estimates a sparse inverse covariance matrix $Θ$ from multivariate Gaussian data $\mathcal{X} \sim \mathcal{N}(μ, Σ) \in \mathbb{R}^p$. Originally proposed by (Dempster 1972) under the name Covariance Selection, this e…
▽ More
We introduce GGLasso, a Python package for solving General Graphical Lasso problems. The Graphical Lasso scheme, introduced by (Friedman 2007) (see also (Yuan 2007; Banerjee 2008)), estimates a sparse inverse covariance matrix $Θ$ from multivariate Gaussian data $\mathcal{X} \sim \mathcal{N}(μ, Σ) \in \mathbb{R}^p$. Originally proposed by (Dempster 1972) under the name Covariance Selection, this estimation framework has been extended to include latent variables in (Chandrasekaran 2012). Recent extensions also include the joint estimation of multiple inverse covariance matrices, see, e.g., in (Danaher 2013; Tomasi 2018). The GGLasso package contains methods for solving a general problem formulation, including important special cases, such as, the single (latent variable) Graphical Lasso, the Group, and the Fused Graphical Lasso.
△ Less
Submitted 20 October, 2021;
originally announced October 2021.
-
latentcor: An R Package for estimating latent correlations from mixed data types
Authors:
Mingze Huang,
Christian L. Müller,
Irina Gaynanova
Abstract:
We present `latentcor`, an R package for correlation estimation from data with mixed variable types. Mixed variables types, including continuous, binary, ordinal, zero-inflated, or truncated data are routinely collected in many areas of science. Accurate estimation of correlations among such variables is often the first critical step in statistical analysis workflows. Pearson correlation as the de…
▽ More
We present `latentcor`, an R package for correlation estimation from data with mixed variable types. Mixed variables types, including continuous, binary, ordinal, zero-inflated, or truncated data are routinely collected in many areas of science. Accurate estimation of correlations among such variables is often the first critical step in statistical analysis workflows. Pearson correlation as the default choice is not well suited for mixed data types as the underlying normality assumption is violated. The concept of semi-parametric latent Gaussian copula models, on the other hand, provides a unifying way to estimate correlations between mixed data types. The R package `latentcor` comprises a comprehensive list of these models, enabling the estimation of correlations between any of continuous/binary/ternary/zero-inflated (truncated) variable types. The underlying implementation takes advantage of a fast multi-linear interpolation scheme with an efficient choice of interpolation grid points, thus giving the package a small memory footprint without compromising estimation accuracy. This makes latent correlation estimation readily available for modern high-throughput data analysis.
△ Less
Submitted 20 August, 2021;
originally announced August 2021.
-
Instrumental Variable Estimation for Compositional Treatments
Authors:
Elisabeth Ailer,
Christian L. Müller,
Niki Kilbertus
Abstract:
Many scientific datasets are compositional in nature. Important biological examples include species abundances in ecology, cell-type compositions derived from single-cell sequencing data, and amplicon abundance data in microbiome research. Here, we provide a causal view on compositional data in an instrumental variable setting where the composition acts as the cause. First, we crisply articulate p…
▽ More
Many scientific datasets are compositional in nature. Important biological examples include species abundances in ecology, cell-type compositions derived from single-cell sequencing data, and amplicon abundance data in microbiome research. Here, we provide a causal view on compositional data in an instrumental variable setting where the composition acts as the cause. First, we crisply articulate potential pitfalls for practitioners regarding the interpretation of compositional causes from the viewpoint of interventions and warn against attributing causal meaning to common summary statistics such as diversity indices in microbiome data analysis. We then advocate for and develop multivariate methods using statistical data transformations and regression techniques that take the special structure of the compositional sample space into account while still yielding scientifically interpretable results. In a comparative analysis on synthetic and real microbiome data we show the advantages and limitations of our proposal. We posit that our analysis provides a useful framework and guidance for valid and informative cause-effect estimation in the context of compositional data.
△ Less
Submitted 28 May, 2024; v1 submitted 21 June, 2021;
originally announced June 2021.
-
deepregression: a Flexible Neural Network Framework for Semi-Structured Deep Distributional Regression
Authors:
David Rügamer,
Chris Kolb,
Cornelius Fritz,
Florian Pfisterer,
Philipp Kopper,
Bernd Bischl,
Ruolin Shen,
Christina Bukas,
Lisa Barros de Andrade e Sousa,
Dominik Thalmeier,
Philipp Baumann,
Lucas Kook,
Nadja Klein,
Christian L. Müller
Abstract:
In this paper we describe the implementation of semi-structured deep distributional regression, a flexible framework to learn conditional distributions based on the combination of additive regression models and deep networks. Our implementation encompasses (1) a modular neural network building system based on the deep learning library \pkg{TensorFlow} for the fusion of various statistical and deep…
▽ More
In this paper we describe the implementation of semi-structured deep distributional regression, a flexible framework to learn conditional distributions based on the combination of additive regression models and deep networks. Our implementation encompasses (1) a modular neural network building system based on the deep learning library \pkg{TensorFlow} for the fusion of various statistical and deep learning approaches, (2) an orthogonalization cell to allow for an interpretable combination of different subnetworks, as well as (3) pre-processing steps necessary to set up such models. The software package allows to define models in a user-friendly manner via a formula interface that is inspired by classical statistical model frameworks such as \pkg{mgcv}. The packages' modular design and functionality provides a unique resource for both scalable estimation of complex statistical models and the combination of approaches from deep learning and statistics. This allows for state-of-the-art predictive performance while simultaneously retaining the indispensable interpretability of classical statistical models.
△ Less
Submitted 10 March, 2022; v1 submitted 6 April, 2021;
originally announced April 2021.
-
STENCIL-NET: Data-driven solution-adaptive discretization of partial differential equations
Authors:
Suryanarayana Maddu,
Dominik Sturm,
Bevan L. Cheeseman,
Christian L. Müller,
Ivo F. Sbalzarini
Abstract:
Numerical methods for approximately solving partial differential equations (PDE) are at the core of scientific computing. Often, this requires high-resolution or adaptive discretization grids to capture relevant spatio-temporal features in the PDE solution, e.g., in applications like turbulence, combustion, and shock propagation. Numerical approximation also requires knowing the PDE in order to co…
▽ More
Numerical methods for approximately solving partial differential equations (PDE) are at the core of scientific computing. Often, this requires high-resolution or adaptive discretization grids to capture relevant spatio-temporal features in the PDE solution, e.g., in applications like turbulence, combustion, and shock propagation. Numerical approximation also requires knowing the PDE in order to construct problem-specific discretizations. Systematically deriving such solution-adaptive discrete operators, however, is a current challenge. Here we present STENCIL-NET, an artificial neural network architecture for data-driven learning of problem- and resolution-specific local discretizations of nonlinear PDEs. STENCIL-NET achieves numerically stable discretization of the operators in an unknown nonlinear PDE by spatially and temporally adaptive parametric pooling on regular Cartesian grids, and by incorporating knowledge about discrete time integration. Knowing the actual PDE is not necessary, as solution data is sufficient to train the network to learn the discrete operators. A once-trained STENCIL-NET model can be used to predict solutions of the PDE on larger spatial domains and for longer times than it was trained for, hence addressing the problem of PDE-constrained extrapolation from data. To support this claim, we present numerical experiments on long-term forecasting of chaotic PDE solutions on coarse spatio-temporal grids. We also quantify the speed-up achieved by substituting base-line numerical methods with equation-free STENCIL-NET predictions on coarser grids with little compromise on accuracy.
△ Less
Submitted 18 January, 2021; v1 submitted 15 January, 2021;
originally announced January 2021.
-
Learning physically consistent mathematical models from data using group sparsity
Authors:
Suryanarayana Maddu,
Bevan L. Cheeseman,
Christian L. Müller,
Ivo F. Sbalzarini
Abstract:
We propose a statistical learning framework based on group-sparse regression that can be used to 1) enforce conservation laws, 2) ensure model equivalence, and 3) guarantee symmetries when learning or inferring differential-equation models from measurement data. Directly learning $\textit{interpretable}$ mathematical models from data has emerged as a valuable modeling approach. However, in areas l…
▽ More
We propose a statistical learning framework based on group-sparse regression that can be used to 1) enforce conservation laws, 2) ensure model equivalence, and 3) guarantee symmetries when learning or inferring differential-equation models from measurement data. Directly learning $\textit{interpretable}$ mathematical models from data has emerged as a valuable modeling approach. However, in areas like biology, high noise levels, sensor-induced correlations, and strong inter-system variability can render data-driven models nonsensical or physically inconsistent without additional constraints on the model structure. Hence, it is important to leverage $\textit{prior}$ knowledge from physical principles to learn "biologically plausible and physically consistent" models rather than models that simply fit the data best. We present a novel group Iterative Hard Thresholding (gIHT) algorithm and use stability selection to infer physically consistent models with minimal parameter tuning. We show several applications from systems biology that demonstrate the benefits of enforcing $\textit{priors}$ in data-driven modeling.
△ Less
Submitted 11 December, 2020;
originally announced December 2020.
-
c-lasso -- a Python package for constrained sparse and robust regression and classification
Authors:
Léo Simpson,
Patrick L. Combettes,
Christian L. Müller
Abstract:
We introduce c-lasso, a Python package that enables sparse and robust linear regression and classification with linear equality constraints. The underlying statistical forward model is assumed to be of the following form: \[ y = X β+ σε\qquad \textrm{subject to} \qquad Cβ=0 \] Here, $X \in \mathbb{R}^{n\times d}$is a given design matrix and the vector $y \in \mathbb{R}^{n}$ is a continuous or bina…
▽ More
We introduce c-lasso, a Python package that enables sparse and robust linear regression and classification with linear equality constraints. The underlying statistical forward model is assumed to be of the following form: \[ y = X β+ σε\qquad \textrm{subject to} \qquad Cβ=0 \] Here, $X \in \mathbb{R}^{n\times d}$is a given design matrix and the vector $y \in \mathbb{R}^{n}$ is a continuous or binary response vector. The matrix $C$ is a general constraint matrix. The vector $β\in \mathbb{R}^{d}$ contains the unknown coefficients and $σ$ an unknown scale. Prominent use cases are (sparse) log-contrast regression with compositional data $X$, requiring the constraint $1_d^T β= 0$ (Aitchion and Bacon-Shone 1984) and the Generalized Lasso which is a special case of the described problem (see, e.g, (James, Paulson, and Rusmevichientong 2020), Example 3). The c-lasso package provides estimators for inferring unknown coefficients and scale (i.e., perspective M-estimators (Combettes and Müller 2020a)) of the form \[ \min_{β\in \mathbb{R}^d, σ\in \mathbb{R}_{0}} f\left(Xβ- y,σ \right) + λ\left\lVert β\right\rVert_1 \qquad \textrm{subject to} \qquad Cβ= 0 \] for several convex loss functions $f(\cdot,\cdot)$. This includes the constrained Lasso, the constrained scaled Lasso, and sparse Huber M-estimators with linear equality constraints.
△ Less
Submitted 2 November, 2020;
originally announced November 2020.
-
Fast computation of latent correlations
Authors:
Grace Yoon,
Christian L. Müller,
Irina Gaynanova
Abstract:
Latent Gaussian copula models provide a powerful means to perform multi-view data integration since these models can seamlessly express dependencies between mixed variable types (binary, continuous, zero-inflated) via latent Gaussian correlations. The estimation of these latent correlations, however, comes at considerable computational cost, having prevented the routine use of these models on high…
▽ More
Latent Gaussian copula models provide a powerful means to perform multi-view data integration since these models can seamlessly express dependencies between mixed variable types (binary, continuous, zero-inflated) via latent Gaussian correlations. The estimation of these latent correlations, however, comes at considerable computational cost, having prevented the routine use of these models on high-dimensional data. Here, we propose a new computational approach for estimating latent correlations via a hybrid multi-linear interpolation and optimization scheme. Our approach speeds up the current state of the art computation by several orders of magnitude, thus allowing fast computation of latent Gaussian copula models even when the number of variables $p$ is large. We provide theoretical guarantees for the approximation error of our numerical scheme and support its excellent performance on simulated and real-world data. We illustrate the practical advantages of our method on high-dimensional sparse quantitative and relative abundance microbiome data as well as multi-view data from The Cancer Genome Atlas Project. Our method is implemented in the R package mixedCCA, available at https://github.com/irinagain/mixedCCA.
△ Less
Submitted 3 November, 2020; v1 submitted 24 June, 2020;
originally announced June 2020.
-
Stability selection enables robust learning of partial differential equations from limited noisy data
Authors:
Suryanarayana Maddu,
Bevan L. Cheeseman,
Ivo F. Sbalzarini,
Christian L. Müller
Abstract:
We present a statistical learning framework for robust identification of partial differential equations from noisy spatiotemporal data. Extending previous sparse regression approaches for inferring PDE models from simulated data, we address key issues that have thus far limited the application of these methods to noisy experimental data, namely their robustness against noise and the need for manua…
▽ More
We present a statistical learning framework for robust identification of partial differential equations from noisy spatiotemporal data. Extending previous sparse regression approaches for inferring PDE models from simulated data, we address key issues that have thus far limited the application of these methods to noisy experimental data, namely their robustness against noise and the need for manual parameter tuning. We address both points by proposing a stability-based model selection scheme to determine the level of regularization required for reproducible recovery of the underlying PDE. This avoids manual parameter tuning and provides a principled way to improve the method's robustness against noise in the data. Our stability selection approach, termed PDE-STRIDE, can be combined with any sparsity-promoting penalized regression model and provides an interpretable criterion for model component importance. We show that in particular the combination of stability selection with the iterative hard-thresholding algorithm from compressed sensing provides a fast, parameter-free, and robust computational framework for PDE inference that outperforms previous algorithmic approaches with respect to recovery accuracy, amount of data required, and robustness to noise. We illustrate the performance of our approach on a wide range of noise-corrupted simulated benchmark problems, including 1D Burgers, 2D vorticity-transport, and 3D reaction-diffusion problems. We demonstrate the practical applicability of our method on real-world data by considering a purely data-driven re-evaluation of the advective triggering hypothesis for an embryonic polarization system in C.~elegans. Using fluorescence microscopy images of C.~elegans zygotes as input data, our framework is able to recover the PDE model for the regulatory reaction-diffusion-flow network of the associated proteins.
△ Less
Submitted 17 July, 2019;
originally announced July 2019.
-
Regression models for compositional data: General log-contrast formulations, proximal optimization, and microbiome data applications
Authors:
Patrick L. Combettes,
Christian L. Müller
Abstract:
Compositional data sets are ubiquitous in science, including geology, ecology, and microbiology. In microbiome research, compositional data primarily arise from high-throughput sequence-based profiling experiments. These data comprise microbial compositions in their natural habitat and are often paired with covariate measurements that characterize physicochemical habitat properties or the physiolo…
▽ More
Compositional data sets are ubiquitous in science, including geology, ecology, and microbiology. In microbiome research, compositional data primarily arise from high-throughput sequence-based profiling experiments. These data comprise microbial compositions in their natural habitat and are often paired with covariate measurements that characterize physicochemical habitat properties or the physiology of the host. Inferring parsimonious statistical associations between microbial compositions and habitat- or host-specific covariate data is an important step in exploratory data analysis. A standard statistical model linking compositional covariates to continuous outcomes is the linear log-contrast model. This model describes the response as a linear combination of log-ratios of the original compositions and has been extended to the high-dimensional setting via regularization. In this contribution, we propose a general convex optimization model for linear log-contrast regression which includes many previous proposals as special cases. We introduce a proximal algorithm that solves the resulting constrained optimization problem exactly with rigorous convergence guarantees. We illustrate the versatility of our approach by investigating the performance of several model instances on soil and gut microbiome data analysis tasks.
△ Less
Submitted 3 March, 2019;
originally announced March 2019.
-
Perspective Maximum Likelihood-Type Estimation via Proximal Decomposition
Authors:
Patrick L. Combettes,
Christian L. Müller
Abstract:
We introduce an optimization model for maximum likelihood-type estimation (M-estimation) that generalizes a large class of existing statistical models, including Huber's concomitant M-estimator, Owen's Huber/Berhu concomitant estimator, the scaled lasso, support vector machine regression, and penalized estimation with structured sparsity. The model, termed perspective M-estimation, leverages the o…
▽ More
We introduce an optimization model for maximum likelihood-type estimation (M-estimation) that generalizes a large class of existing statistical models, including Huber's concomitant M-estimator, Owen's Huber/Berhu concomitant estimator, the scaled lasso, support vector machine regression, and penalized estimation with structured sparsity. The model, termed perspective M-estimation, leverages the observation that convex M-estimators with concomitant scale as well as various regularizers are instances of perspective functions. Such functions are amenable to proximal analysis, which leads to principled and provably convergent optimization algorithms via proximal splitting. Using a geometrical approach based on duality, we derive novel proximity operators for several perspective functions of interest. Numerical experiments on synthetic and real-world data illustrate the broad applicability of the proposed framework.
△ Less
Submitted 5 October, 2018; v1 submitted 15 May, 2018;
originally announced May 2018.
-
Construction of KAGRA: an Underground Gravitational Wave Observatory
Authors:
T. Akutsu,
M. Ando,
S. Araki,
A. Araya,
T. Arima,
N. Aritomi,
H. Asada,
Y. Aso,
S. Atsuta,
K. Awai,
L. Baiotti,
M. A. Barton,
D. Chen,
K. Cho,
K. Craig,
R. DeSalvo,
K. Doi,
K. Eda,
Y. Enomoto,
R. Flaminio,
S. Fujibayashi,
Y. Fujii,
M. -K. Fujimoto,
M. Fukushima,
T. Furuhata
, et al. (202 additional authors not shown)
Abstract:
Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferomet…
▽ More
Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.
△ Less
Submitted 11 December, 2017; v1 submitted 30 November, 2017;
originally announced December 2017.
-
Perspective Functions: Proximal Calculus and Applications in High-Dimensional Statistics
Authors:
Patrick L. Combettes,
Christian L. Müller
Abstract:
Perspective functions arise explicitly or implicitly in various forms in applied mathematics and in statistical data analysis. To date, no systematic strategy is available to solve the associated, typically nonsmooth, optimization problems. In this paper, we fill this gap by showing that proximal methods provide an efficient framework to model and solve problems involving perspective functions. We…
▽ More
Perspective functions arise explicitly or implicitly in various forms in applied mathematics and in statistical data analysis. To date, no systematic strategy is available to solve the associated, typically nonsmooth, optimization problems. In this paper, we fill this gap by showing that proximal methods provide an efficient framework to model and solve problems involving perspective functions. We study the construction of the proximity operator of a perspective function under general assumptions and present important instances in which the proximity operator can be computed explicitly or via straightforward numerical operations. These results constitute central building blocks in the design of proximal optimization algorithms. We showcase the versatility of the framework by designing novel proximal algorithms for state-of-the-art regression and variable selection schemes in high-dimensional statistics.
△ Less
Submitted 13 December, 2016; v1 submitted 5 October, 2016;
originally announced October 2016.
-
Small optic suspensions for Advanced LIGO input optics and other precision optical experiments
Authors:
G. Ciani,
M. A. Arain,
S. M. Aston,
D. Feldbaum,
P. Fulda,
J. Gleason,
M. Heintze,
R. M. Martin,
C. L. Mueller,
D. M. Nanda Kumar,
A. Pele,
D. H. Reitze,
P. Sainathan,
D. B. Tanner,
L. F. Williams,
G. Mueller
Abstract:
We report on the design and performance of small optic suspensions developed to suppress seismic motion of out-of-cavity optics in the Input Optics subsystem of the Advanced LIGO interferometric gravitational wave detector. These compact single stage suspensions provide isolation in all six degrees of freedom of the optic, local sensing and actuation in three of them, and passive damping for the o…
▽ More
We report on the design and performance of small optic suspensions developed to suppress seismic motion of out-of-cavity optics in the Input Optics subsystem of the Advanced LIGO interferometric gravitational wave detector. These compact single stage suspensions provide isolation in all six degrees of freedom of the optic, local sensing and actuation in three of them, and passive damping for the other three.
△ Less
Submitted 29 June, 2016;
originally announced June 2016.
-
Supplement: The Rate of Binary Black Hole Mergers Inferred from Advanced LIGO Observations Surrounding GW150914
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
P. A. Altin,
S. B. Anderson
, et al. (942 additional authors not shown)
Abstract:
Supplemental information for a Letter reporting the rate of binary black hole (BBH) coalescences inferred from 16 days of coincident Advanced LIGO observations surrounding the transient gravitational wave signal GW150914. In that work we reported various rate estimates whose 90\% credible intervals fell in the range $2$--$600 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$. Here we give details of our meth…
▽ More
Supplemental information for a Letter reporting the rate of binary black hole (BBH) coalescences inferred from 16 days of coincident Advanced LIGO observations surrounding the transient gravitational wave signal GW150914. In that work we reported various rate estimates whose 90\% credible intervals fell in the range $2$--$600 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$. Here we give details of our method and computations, including information about our search pipelines, a derivation of our likelihood function for the analysis, a description of the astrophysical search trigger distribution expected from merging BBHs, details on our computational methods, a description of the effects and our model for calibration uncertainty, and an analytic method of estimating our detector sensitivity that is calibrated to our measurements.
△ Less
Submitted 20 September, 2016; v1 submitted 13 June, 2016;
originally announced June 2016.
-
Generalized Stability Approach for Regularized Graphical Models
Authors:
Christian L. Müller,
Richard Bonneau,
Zachary Kurtz
Abstract:
Selecting regularization parameters in penalized high-dimensional graphical models in a principled, data-driven, and computationally efficient manner continues to be one of the key challenges in high-dimensional statistics. We present substantial computational gains and conceptual generalizations of the Stability Approach to Regularization Selection (StARS), a state-of-the-art graphical model sele…
▽ More
Selecting regularization parameters in penalized high-dimensional graphical models in a principled, data-driven, and computationally efficient manner continues to be one of the key challenges in high-dimensional statistics. We present substantial computational gains and conceptual generalizations of the Stability Approach to Regularization Selection (StARS), a state-of-the-art graphical model selection scheme. Using properties of the Poisson-Binomial distribution and convex non-asymptotic distributional modeling we propose lower and upper bounds on the StARS graph regularization path which results in greatly reduced computational cost without compromising regularization selection. We also generalize the StARS criterion from single edge to induced subgraph (graphlet) stability. We show that simultaneously requiring edge and graphlet stability leads to superior graph recovery performance independent of graph topology. These novel insights render Gaussian graphical model selection a routine task on standard multi-core computers.
△ Less
Submitted 23 May, 2016;
originally announced May 2016.
-
A First Targeted Search for Gravitational-Wave Bursts from Core-Collapse Supernovae in Data of First-Generation Laser Interferometer Detectors
Authors:
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
P. A. Altin,
S. B. Anderson,
W. G. Anderson,
K. Arai
, et al. (942 additional authors not shown)
Abstract:
We present results from a search for gravitational-wave bursts coincident with a set of two core-collapse supernovae observed between 2007 and 2011. We employ data from the Laser Interferometer Gravitational-wave Observatory (LIGO), the Virgo gravitational-wave observatory, and the GEO 600 gravitational-wave observatory. The targeted core-collapse supernovae were selected on the basis of (1) proxi…
▽ More
We present results from a search for gravitational-wave bursts coincident with a set of two core-collapse supernovae observed between 2007 and 2011. We employ data from the Laser Interferometer Gravitational-wave Observatory (LIGO), the Virgo gravitational-wave observatory, and the GEO 600 gravitational-wave observatory. The targeted core-collapse supernovae were selected on the basis of (1) proximity (within approximately 15 Mpc), (2) tightness of observational constraints on the time of core collapse that defines the gravitational-wave search window, and (3) coincident operation of at least two interferometers at the time of core collapse. We find no plausible gravitational-wave candidates. We present the probability of detecting signals from both astrophysically well-motivated and more speculative gravitational-wave emission mechanisms as a function of distance from Earth, and discuss the implications for the detection of gravitational waves from core-collapse supernovae by the upgraded Advanced LIGO and Virgo detectors.
△ Less
Submitted 19 May, 2016; v1 submitted 5 May, 2016;
originally announced May 2016.
-
Search for transient gravitational waves in coincidence with short duration radio transients during 2007-2013
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
others,
:,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca
, et al. (977 additional authors not shown)
Abstract:
We present an archival search for transient gravitational-wave bursts in coincidence with 27 single pulse triggers from Green Bank Telescope pulsar surveys, using the LIGO, Virgo and GEO interferometer network. We also discuss a check for gravitational-wave signals in coincidence with Parkes Fast Radio Bursts using similar methods. Data analyzed in these searches were collected between 2007 and 20…
▽ More
We present an archival search for transient gravitational-wave bursts in coincidence with 27 single pulse triggers from Green Bank Telescope pulsar surveys, using the LIGO, Virgo and GEO interferometer network. We also discuss a check for gravitational-wave signals in coincidence with Parkes Fast Radio Bursts using similar methods. Data analyzed in these searches were collected between 2007 and 2013. Possible sources of emission of both short-duration radio signals and transient gravitational-wave emission include starquakes on neutron stars, binary coalescence of neutron stars, and cosmic string cusps. While no evidence for gravitational-wave emission in coincidence with these radio transients was found, the current analysis serves as a prototype for similar future searches using more sensitive second-generation interferometers.
△ Less
Submitted 21 June, 2016; v1 submitted 5 May, 2016;
originally announced May 2016.
-
Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914
Authors:
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
P. A. Altin,
S. B. Anderson,
W. G. Anderson,
K. Arai
, et al. (1522 additional authors not shown)
Abstract:
This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the dif…
▽ More
This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.
△ Less
Submitted 21 July, 2016; v1 submitted 26 April, 2016;
originally announced April 2016.
-
The Sensitivity of the Advanced LIGO Detectors at the Beginning of Gravitational Wave Astronomy
Authors:
D. V. Martynov,
E. D. Hall,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
C. Adams,
R. X. Adhikari,
R. A. Anderson,
S. B. Anderson,
K. Arai,
M. A. Arain,
S. M. Aston,
L. Austin,
S. W. Ballmer,
M. Barbet,
D. Barker,
B. Barr,
L. Barsotti,
J. Bartlett,
M. A. Barton,
I. Bartos,
J. C. Batch,
A. S. Bell,
I. Belopolski,
J. Bergman
, et al. (239 additional authors not shown)
Abstract:
The Laser Interferometer Gravitational Wave Observatory (LIGO) consists of two widely separated 4 km laser interferometers designed to detect gravitational waves from distant astrophysical sources in the frequency range from 10 Hz to 10 kHz. The first observation run of the Advanced LIGO detectors started in September 2015 and ended in January 2016. A strain sensitivity of better than…
▽ More
The Laser Interferometer Gravitational Wave Observatory (LIGO) consists of two widely separated 4 km laser interferometers designed to detect gravitational waves from distant astrophysical sources in the frequency range from 10 Hz to 10 kHz. The first observation run of the Advanced LIGO detectors started in September 2015 and ended in January 2016. A strain sensitivity of better than $10^{-23}/\sqrt{\text{Hz}}$ was achieved around 100 Hz. Understanding both the fundamental and the technical noise sources was critical for increasing the observable volume in the universe. The average distance at which coalescing binary black hole systems with individual masses of 30 $M_\odot$ could be detected was 1.3 Gpc. Similarly, the range for binary neutron star inspirals was about 75 Mpc. With respect to the initial detectors, the observable volume of Universe increased respectively by a factor 69 and 43. These improvements allowed Advanced LIGO to detect the gravitational wave signal from the binary black hole coalescence, known as GW150914.
△ Less
Submitted 10 February, 2018; v1 submitted 1 April, 2016;
originally announced April 2016.
-
Localization and broadband follow-up of the gravitational-wave transient GW150914
Authors:
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
P. A. Altin,
S. B. Anderson,
W. G. Anderson,
K. Arai
, et al. (1522 additional authors not shown)
Abstract:
A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared wit…
▽ More
A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.
△ Less
Submitted 21 July, 2016; v1 submitted 26 February, 2016;
originally announced February 2016.
-
High-energy Neutrino follow-up search of Gravitational Wave Event GW150914 with ANTARES and IceCube
Authors:
S. Adrián-Martínez,
A. Albert,
M. André,
G. Anton,
M. Ardid,
J. -J. Aubert,
T. Avgitas,
B. Baret,
J. Barrios-Martí,
S. Basa,
V. Bertin,
S. Biagi,
R. Bormuth,
M. C. Bouwhuis,
R. Bruijn,
J. Brunner,
J. Busto,
A. Capone,
L. Caramete,
J. Carr,
S. Celli,
T. Chiarusi,
M. Circella,
A. Coleiro,
R. Coniglione
, et al. (1369 additional authors not shown)
Abstract:
We present the high-energy-neutrino follow-up observations of the first gravitational wave transient GW150914 observed by the Advanced LIGO detectors on Sept. 14th, 2015. We search for coincident neutrino candidates within the data recorded by the IceCube and ANTARES neutrino detectors. A possible joint detection could be used in targeted electromagnetic follow-up observations, given the significa…
▽ More
We present the high-energy-neutrino follow-up observations of the first gravitational wave transient GW150914 observed by the Advanced LIGO detectors on Sept. 14th, 2015. We search for coincident neutrino candidates within the data recorded by the IceCube and ANTARES neutrino detectors. A possible joint detection could be used in targeted electromagnetic follow-up observations, given the significantly better angular resolution of neutrino events compared to gravitational waves. We find no neutrino candidates in both temporal and spatial coincidence with the gravitational wave event. Within 500 s of the gravitational wave event, the number of neutrino candidates detected by IceCube and ANTARES were three and zero, respectively. This is consistent with the expected atmospheric background, and none of the neutrino candidates were directionally coincident with GW150914. We use this non-detection to constrain neutrino emission from the gravitational-wave event.
△ Less
Submitted 22 April, 2016; v1 submitted 17 February, 2016;
originally announced February 2016.
-
Calibration of the Advanced LIGO detectors for the discovery of the binary black-hole merger GW150914
Authors:
The LIGO Scientific Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
K. Ackley,
C. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
B. Allen,
P. A. Altin,
D. V. Amariutei,
S. B. Anderson,
W. G. Anderson,
K. Arai,
M. C. Araya,
C. C. Arceneaux,
J. S. Areeda,
K. G. Arun
, et al. (702 additional authors not shown)
Abstract:
In Advanced LIGO, detection and astrophysical source parameter estimation of the binary black hole merger GW150914 requires a calibrated estimate of the gravitational-wave strain sensed by the detectors. Producing an estimate from each detector's differential arm length control loop readout signals requires applying time domain filters, which are designed from a frequency domain model of the detec…
▽ More
In Advanced LIGO, detection and astrophysical source parameter estimation of the binary black hole merger GW150914 requires a calibrated estimate of the gravitational-wave strain sensed by the detectors. Producing an estimate from each detector's differential arm length control loop readout signals requires applying time domain filters, which are designed from a frequency domain model of the detector's gravitational-wave response. The gravitational-wave response model is determined by the detector's opto-mechanical response and the properties of its feedback control system. The measurements used to validate the model and characterize its uncertainty are derived primarily from a dedicated photon radiation pressure actuator, with cross-checks provided by optical and radio frequency references. We describe how the gravitational-wave readout signal is calibrated into equivalent gravitational-wave-induced strain and how the statistical uncertainties and systematic errors are assessed. Detector data collected over 38 calendar days, from September 12 to October 20, 2015, contain the event GW150914 and approximately 16 of coincident data used to estimate the event false alarm probability. The calibration uncertainty is less than 10% in magnitude and 10 degrees in phase across the relevant frequency band 20 Hz to 1 kHz.
△ Less
Submitted 28 February, 2017; v1 submitted 11 February, 2016;
originally announced February 2016.
-
The Rate of Binary Black Hole Mergers Inferred from Advanced LIGO Observations Surrounding GW150914
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
P. A. Altin,
S. B. Anderson
, et al. (942 additional authors not shown)
Abstract:
A transient gravitational-wave signal, GW150914, was identified in the twin Advanced LIGO detectors on September 14, 2015 at 09:50:45 UTC. To assess the implications of this discovery, the detectors remained in operation with unchanged configurations over a period of 39 d around the time of the signal. At the detection statistic threshold corresponding to that observed for GW150914, our search of…
▽ More
A transient gravitational-wave signal, GW150914, was identified in the twin Advanced LIGO detectors on September 14, 2015 at 09:50:45 UTC. To assess the implications of this discovery, the detectors remained in operation with unchanged configurations over a period of 39 d around the time of the signal. At the detection statistic threshold corresponding to that observed for GW150914, our search of the 16 days of simultaneous two-detector observational data is estimated to have a false alarm rate (FAR) of $< 4.9 \times 10^{-6} \, \mathrm{yr}^{-1}$, yielding a $p$-value for GW150914 of $< 2 \times 10^{-7}$. Parameter estimation followup on this trigger identifies its source as a binary black hole (BBH) merger with component masses $(m_1, m_2) = \left(36^{+5}_{-4},29^{+4}_{-4}\right) \, M_\odot$ at redshift $z = 0.09^{+0.03}_{-0.04}$ (median and 90\% credible range). Here we report on the constraints these observations place on the rate of BBH coalescences. Considering only GW150914, assuming that all BBHs in the Universe have the same masses and spins as this event, imposing a search FAR threshold of 1 per 100 years, and assuming that the BBH merger rate is constant in the comoving frame, we infer a 90% credible range of merger rates between $2$--$53 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ (comoving frame). Incorporating all search triggers that pass a much lower threshold while accounting for the uncertainty in the astrophysical origin of each trigger, we estimate a higher rate, ranging from $13$--$600 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ depending on assumptions about the BBH mass distribution. All together, our various rate estimates fall in the conservative range $2$--$600 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$.
△ Less
Submitted 20 September, 2016; v1 submitted 11 February, 2016;
originally announced February 2016.
-
Tests of general relativity with GW150914
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
P. A. Altin,
S. B. Anderson
, et al. (963 additional authors not shown)
Abstract:
The LIGO detection of GW150914 provides an unprecedented opportunity to study the two-body motion of a compact-object binary in the large velocity, highly nonlinear regime, and to witness the final merger of the binary and the excitation of uniquely relativistic modes of the gravitational field. We carry out several investigations to determine whether GW150914 is consistent with a binary black-hol…
▽ More
The LIGO detection of GW150914 provides an unprecedented opportunity to study the two-body motion of a compact-object binary in the large velocity, highly nonlinear regime, and to witness the final merger of the binary and the excitation of uniquely relativistic modes of the gravitational field. We carry out several investigations to determine whether GW150914 is consistent with a binary black-hole merger in general relativity. We find that the final remnant's mass and spin, as determined from the low-frequency (inspiral) and high-frequency (post-inspiral) phases of the signal, are mutually consistent with the binary black-hole solution in general relativity. Furthermore, the data following the peak of GW150914 are consistent with the least-damped quasi-normal mode inferred from the mass and spin of the remnant black hole. By using waveform models that allow for parameterized general-relativity violations during the inspiral and merger phases, we perform quantitative tests on the gravitational-wave phase in the dynamical regime and we determine the first empirical bounds on several high-order post-Newtonian coefficients. We constrain the graviton Compton wavelength, assuming that gravitons are dispersed in vacuum in the same way as particles with mass, obtaining a $90\%$-confidence lower bound of $10^{13}$ km. In conclusion, within our statistical uncertainties, we find no evidence for violations of general relativity in the genuinely strong-field regime of gravity.
△ Less
Submitted 23 October, 2018; v1 submitted 11 February, 2016;
originally announced February 2016.
-
Properties of the Binary Black Hole Merger GW150914
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca
, et al. (980 additional authors not shown)
Abstract:
On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in ge…
▽ More
On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in general relativity. GW150914 was produced by a nearly equal mass binary black hole of $36^{+5}_{-4} M_\odot$ and $29^{+4}_{-4} M_\odot$; for each parameter we report the median value and the range of the 90% credible interval. The dimensionless spin magnitude of the more massive black hole is bound to be $<0.7$ (at 90% probability). The luminosity distance to the source is $410^{+160}_{-180}$ Mpc, corresponding to a redshift $0.09^{+0.03}_{-0.04}$ assuming standard cosmology. The source location is constrained to an annulus section of $610$ deg$^2$, primarily in the southern hemisphere. The binary merges into a black hole of $62^{+4}_{-4} M_\odot$ and spin $0.67^{+0.05}_{-0.07}$. This black hole is significantly more massive than any other inferred from electromagnetic observations in the stellar-mass regime.
△ Less
Submitted 20 June, 2016; v1 submitted 11 February, 2016;
originally announced February 2016.
-
GW150914: First results from the search for binary black hole coalescence with Advanced LIGO
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
P. A. Altin,
S. B. Anderson
, et al. (957 additional authors not shown)
Abstract:
On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-wave Observatory (LIGO) simultaneously observed the binary black hole merger GW150914. We report the results of a matched-filter search using relativistic models of compact-object binaries that recovered GW150914 as the most significant event during the coincident observations between the two LIGO det…
▽ More
On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-wave Observatory (LIGO) simultaneously observed the binary black hole merger GW150914. We report the results of a matched-filter search using relativistic models of compact-object binaries that recovered GW150914 as the most significant event during the coincident observations between the two LIGO detectors from September 12 to October 20, 2015. GW150914 was observed with a matched filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203 000 years, equivalent to a significance greater than 5.1 σ.
△ Less
Submitted 27 April, 2016; v1 submitted 11 February, 2016;
originally announced February 2016.
-
All-sky search for long-duration gravitational wave transients with LIGO
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
D. V. Amariutei,
S. B. Anderson,
W. G. Anderson
, et al. (913 additional authors not shown)
Abstract:
We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10 - 500 s i…
▽ More
We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10 - 500 s in a frequency band of 40 - 1000 Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. All candidate triggers were consistent with the expected background; as a result we set 90% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. For signals from black hole accretion disk instabilities, we set upper limits on the source rate density between $3.4 \times 10^{-5}$ - $9.4 \times 10^{-4}$ Mpc$^{-3}$ yr$^{-1}$ at 90% confidence. These are the first results from an all-sky search for unmodeled long-duration transient gravitational waves.
△ Less
Submitted 18 February, 2016; v1 submitted 13 November, 2015;
originally announced November 2015.
-
First low frequency all-sky search for continuous gravitational wave signals
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
D. V. Amariutei,
M. Andersen
, et al. (899 additional authors not shown)
Abstract:
In this paper we present the results of the first low frequency all-sky search of continuous gravitational wave signals conducted on Virgo VSR2 and VSR4 data. The search covered the full sky, a frequency range between 20 Hz and 128 Hz with a range of spin-down between $-1.0 \times 10^{-10}$ Hz/s and $+1.5 \times 10^{-11}$ Hz/s, and was based on a hierarchical approach. The starting point was a set…
▽ More
In this paper we present the results of the first low frequency all-sky search of continuous gravitational wave signals conducted on Virgo VSR2 and VSR4 data. The search covered the full sky, a frequency range between 20 Hz and 128 Hz with a range of spin-down between $-1.0 \times 10^{-10}$ Hz/s and $+1.5 \times 10^{-11}$ Hz/s, and was based on a hierarchical approach. The starting point was a set of short Fast Fourier Transforms (FFT), of length 8192 seconds, built from the calibrated strain data. Aggressive data cleaning, both in the time and frequency domains, has been done in order to remove, as much as possible, the effect of disturbances of instrumental origin. On each dataset a number of candidates has been selected, using the FrequencyHough transform in an incoherent step. Only coincident candidates among VSR2 and VSR4 have been examined in order to strongly reduce the false alarm probability, and the most significant candidates have been selected. Selected candidates have been subject to a follow-up by constructing a new set of longer FFTs followed by a further incoherent analysis, still based on the FrequencyHough transform. No evidence for continuous gravitational wave signals was found, therefore we have set a population-based joint VSR2-VSR4 90$\%$ confidence level upper limit on the dimensionless gravitational wave strain in the frequency range between 20 Hz and 128 Hz. This is the first all-sky search for continuous gravitational waves conducted, on data of ground-based interferometric detectors, at frequencies below 50 Hz. We set upper limits in the range between about $10^{-24}$ and $2\times 10^{-23}$ at most frequencies. Our upper limits on signal strain show an improvement of up to a factor of $\sim$2 with respect to the results of previous all-sky searches at frequencies below $80~\mathrm{Hz}$.
△ Less
Submitted 17 October, 2019; v1 submitted 13 October, 2015;
originally announced October 2015.
-
A search of the Orion spur for continuous gravitational waves using a "loosely coherent" algorithm on data from LIGO interferometers
Authors:
J. Aasi,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
D. V. Amariutei,
M. Andersen,
S. B. Anderson,
W. G. Anderson
, et al. (898 additional authors not shown)
Abstract:
We report results of a wideband search for periodic gravitational waves from isolated neutron stars within the Orion spur towards both the inner and outer regions of our Galaxy. As gravitational waves interact very weakly with matter, the search is unimpeded by dust and concentrations of stars. One search disk (A) is $6.87^\circ$ in diameter and centered on…
▽ More
We report results of a wideband search for periodic gravitational waves from isolated neutron stars within the Orion spur towards both the inner and outer regions of our Galaxy. As gravitational waves interact very weakly with matter, the search is unimpeded by dust and concentrations of stars. One search disk (A) is $6.87^\circ$ in diameter and centered on $20^\textrm{h}10^\textrm{m}54.71^\textrm{s}+33^\circ33'25.29"$, and the other (B) is $7.45^\circ$ in diameter and centered on $8^\textrm{h}35^\textrm{m}20.61^\textrm{s}-46^\circ49'25.151"$. We explored the frequency range of 50-1500 Hz and frequency derivative from $0$ to $-5\times 10^{-9}$ Hz/s. A multi-stage, loosely coherent search program allowed probing more deeply than before in these two regions, while increasing coherence length with every stage.
Rigorous followup parameters have winnowed initial coincidence set to only 70 candidates, to be examined manually. None of those 70 candidates proved to be consistent with an isolated gravitational wave emitter, and 95% confidence level upper limits were placed on continuous-wave strain amplitudes. Near $169$ Hz we achieve our lowest 95% CL upper limit on worst-case linearly polarized strain amplitude $h_0$ of $6.3\times 10^{-25}$, while at the high end of our frequency range we achieve a worst-case upper limit of $3.4\times 10^{-24}$ for all polarizations and sky locations.
△ Less
Submitted 14 October, 2015; v1 submitted 12 October, 2015;
originally announced October 2015.
-
In-situ characterization of the thermal state of resonant optical interferometers via tracking of their higher-order mode resonances
Authors:
Chris L. Mueller,
Paul Fulda,
Rana X. Adhikari,
Koji Arai,
Aidan F. Brooks,
Rijuparna Chakraborty,
Valery V. Frolov,
Peter Fritschel,
Eleanor J. King,
David B. Tanner,
Hiroaki Yamamoto,
Guido Mueller
Abstract:
Thermal lensing in resonant optical interferometers such as those used for gravitational wave detection is a concern due to the negative impact on control signals and instrument sensitivity. In this paper we describe a method for monitoring the thermal state of such interferometers by probing the higher-order spatial mode resonances of the cavities within them. We demonstrate the use of this techn…
▽ More
Thermal lensing in resonant optical interferometers such as those used for gravitational wave detection is a concern due to the negative impact on control signals and instrument sensitivity. In this paper we describe a method for monitoring the thermal state of such interferometers by probing the higher-order spatial mode resonances of the cavities within them. We demonstrate the use of this technique to measure changes in the Advanced LIGO input mode cleaner cavity geometry as a function of input power, and subsequently infer the optical absorption at the mirror surfaces at the level of 1 ppm per mirror. We also demonstrate the generation of a useful error signal for thermal state of the Advanced LIGO power recycling cavity by continuously tracking the first order spatial mode resonance frequency. Such an error signal could be used as an input to thermal compensation systems to maintain the interferometer cavity geometries in the presence of transients in circulating light power levels, thereby maintaining optimal sensitivity and maximizing the duty-cycle of the detectors.
△ Less
Submitted 8 February, 2015;
originally announced February 2015.
-
Searches for continuous gravitational waves from nine young supernova remnants
Authors:
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca,
D. Amariutei,
S. B. Anderson
, et al. (871 additional authors not shown)
Abstract:
We describe directed searches for continuous gravitational waves in data from the sixth LIGO science data run. The targets were nine young supernova remnants not associated with pulsars; eight of the remnants are associated with non-pulsing suspected neutron stars. One target's parameters are uncertain enough to warrant two searches, for a total of ten. Each search covered a broad band of frequenc…
▽ More
We describe directed searches for continuous gravitational waves in data from the sixth LIGO science data run. The targets were nine young supernova remnants not associated with pulsars; eight of the remnants are associated with non-pulsing suspected neutron stars. One target's parameters are uncertain enough to warrant two searches, for a total of ten. Each search covered a broad band of frequencies and first and second frequency derivatives for a fixed sky direction. The searches coherently integrated data from the two LIGO interferometers over time spans from 5.3-25.3 days using the matched-filtering F-statistic. We found no credible gravitational-wave signals. We set 95% confidence upper limits as strong (low) as $4\times10^{-25}$ on intrinsic strain, $2\times10^{-7}$ on fiducial ellipticity, and $3\times10^{-6}$ on r-mode amplitude. These beat the indirect limits from energy conservation and are within the range of theoretical predictions for neutron-star ellipticities and r-mode amplitudes.
△ Less
Submitted 12 August, 2021; v1 submitted 18 December, 2014;
originally announced December 2014.
-
A directed search for gravitational waves from Scorpius X-1 with initial LIGO
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca
, et al. (875 additional authors not shown)
Abstract:
We present results of a search for continuously-emitted gravitational radiation, directed at the brightest low-mass X-ray binary, Scorpius X-1. Our semi-coherent analysis covers 10 days of LIGO S5 data ranging from 50-550 Hz, and performs an incoherent sum of coherent $\mathcal{F}$-statistic power distributed amongst frequency-modulated orbital sidebands. All candidates not removed at the veto sta…
▽ More
We present results of a search for continuously-emitted gravitational radiation, directed at the brightest low-mass X-ray binary, Scorpius X-1. Our semi-coherent analysis covers 10 days of LIGO S5 data ranging from 50-550 Hz, and performs an incoherent sum of coherent $\mathcal{F}$-statistic power distributed amongst frequency-modulated orbital sidebands. All candidates not removed at the veto stage were found to be consistent with noise at a 1% false alarm rate. We present Bayesian 95% confidence upper limits on gravitational-wave strain amplitude using two different prior distributions: a standard one, with no a priori assumptions about the orientation of Scorpius X-1; and an angle-restricted one, using a prior derived from electromagnetic observations. Median strain upper limits of 1.3e-24 and 8e-25 are reported at 150 Hz for the standard and angle-restricted searches respectively. This proof of principle analysis was limited to a short observation time by unknown effects of accretion on the intrinsic spin frequency of the neutron star, but improves upon previous upper limits by factors of ~1.4 for the standard, and 2.3 for the angle-restricted search at the sensitive region of the detector.
△ Less
Submitted 1 December, 2014;
originally announced December 2014.
-
Narrow-band search of continuous gravitational-wave signals from Crab and Vela pulsars in Virgo VSR4 data
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca
, et al. (878 additional authors not shown)
Abstract:
In this paper we present the results of a coherent narrow-band search for continuous gravitational-wave signals from the Crab and Vela pulsars conducted on Virgo VSR4 data. In order to take into account a possible small mismatch between the gravitational wave frequency and two times the star rotation frequency, inferred from measurement of the electromagnetic pulse rate, a range of 0.02 Hz around…
▽ More
In this paper we present the results of a coherent narrow-band search for continuous gravitational-wave signals from the Crab and Vela pulsars conducted on Virgo VSR4 data. In order to take into account a possible small mismatch between the gravitational wave frequency and two times the star rotation frequency, inferred from measurement of the electromagnetic pulse rate, a range of 0.02 Hz around two times the star rotational frequency has been searched for both the pulsars. No evidence for a signal has been found and 95$\%$ confidence level upper limits have been computed both assuming polarization parameters are completely unknown and that they are known with some uncertainty, as derived from X-ray observations of the pulsar wind torii. For Vela the upper limits are comparable to the spin-down limit, computed assuming that all the observed spin-down is due to the emission of gravitational waves. For Crab the upper limits are about a factor of two below the spin-down limit, and represent a significant improvement with respect to past analysis. This is the first time the spin-down limit is significantly overcome in a narrow-band search.
△ Less
Submitted 30 October, 2014;
originally announced October 2014.
-
Characterization of the LIGO detectors during their sixth science run
Authors:
The LIGO Scientific Collaboration,
The Virgo Collaboration,
J. Aasi,
J. Abadie,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
C. Adams,
T. Adams,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
B. Allen,
A. Allocca,
E. Amador. Ceron,
D. Amariutei,
R. A. Anderson,
S. B. Anderson,
W. G. Anderson
, et al. (846 additional authors not shown)
Abstract:
In 2009-2010, the Laser Interferometer Gravitational-wave Observa- tory (LIGO) operated together with international partners Virgo and GEO600 as a network to search for gravitational waves of astrophysical origin. The sensitiv- ity of these detectors was limited by a combination of noise sources inherent to the instrumental design and its environment, often localized in time or frequency, that cou…
▽ More
In 2009-2010, the Laser Interferometer Gravitational-wave Observa- tory (LIGO) operated together with international partners Virgo and GEO600 as a network to search for gravitational waves of astrophysical origin. The sensitiv- ity of these detectors was limited by a combination of noise sources inherent to the instrumental design and its environment, often localized in time or frequency, that couple into the gravitational-wave readout. Here we review the performance of the LIGO instruments during this epoch, the work done to characterize the de- tectors and their data, and the effect that transient and continuous noise artefacts have on the sensitivity of LIGO to a variety of astrophysical sources.
△ Less
Submitted 18 November, 2014; v1 submitted 28 October, 2014;
originally announced October 2014.
-
Searching for stochastic gravitational waves using data from the two co-located LIGO Hanford detectors
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
J. Abadie,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
B. Allen,
A. Allocca,
E. Amado. Ceron,
D. Amariutei,
R. A. Anderson,
S. B. Anderson
, et al. (852 additional authors not shown)
Abstract:
Searches for a stochastic gravitational-wave background (SGWB) using terrestrial detectors typically involve cross-correlating data from pairs of detectors. The sensitivity of such cross-correlation analyses depends, among other things, on the separation between the two detectors: the smaller the separation, the better the sensitivity. Hence, a co-located detector pair is more sensitive to a gravi…
▽ More
Searches for a stochastic gravitational-wave background (SGWB) using terrestrial detectors typically involve cross-correlating data from pairs of detectors. The sensitivity of such cross-correlation analyses depends, among other things, on the separation between the two detectors: the smaller the separation, the better the sensitivity. Hence, a co-located detector pair is more sensitive to a gravitational-wave background than a non-co-located detector pair. However, co-located detectors are also expected to suffer from correlated noise from instrumental and environmental effects that could contaminate the measurement of the background. Hence, methods to identify and mitigate the effects of correlated noise are necessary to achieve the potential increase in sensitivity of co-located detectors. Here we report on the first SGWB analysis using the two LIGO Hanford detectors and address the complications arising from correlated environmental noise. We apply correlated noise identification and mitigation techniques to data taken by the two LIGO Hanford detectors, H1 and H2, during LIGO's fifth science run. At low frequencies, 40 - 460 Hz, we are unable to sufficiently mitigate the correlated noise to a level where we may confidently measure or bound the stochastic gravitational-wave signal. However, at high frequencies, 460-1000 Hz, these techniques are sufficient to set a $95%$ confidence level (C.L.) upper limit on the gravitational-wave energy density of Ω(f)<7.7 x 10^{-4} (f/ 900 Hz)^3, which improves on the previous upper limit by a factor of $\sim 180$. In doing so, we demonstrate techniques that will be useful for future searches using advanced detectors, where correlated noise (e.g., from global magnetic fields) may affect even widely separated detectors.
△ Less
Submitted 2 December, 2014; v1 submitted 22 October, 2014;
originally announced October 2014.
-
Sparse and compositionally robust inference of microbial ecological networks
Authors:
Zachary D. Kurtz,
Christian L. Mueller,
Emily R. Miraldi,
Dan R. Littman,
Martin J. Blaser,
Richard A. Bonneau
Abstract:
16S-ribosomal sequencing and other metagonomic techniques provide snapshots of microbial communities, revealing phylogeny and the abundances of microbial populations across diverse ecosystems. While changes in microbial community structure are demonstrably associated with certain environmental conditions, identification of underlying mechanisms requires new statistical tools, as these datasets pre…
▽ More
16S-ribosomal sequencing and other metagonomic techniques provide snapshots of microbial communities, revealing phylogeny and the abundances of microbial populations across diverse ecosystems. While changes in microbial community structure are demonstrably associated with certain environmental conditions, identification of underlying mechanisms requires new statistical tools, as these datasets present several technical challenges. First, the abundances of microbial operational taxonomic units (OTUs) from 16S datasets are compositional, and thus, microbial abundances are not independent. Secondly, microbial sequencing-based studies typically measure hundreds of OTUs on only tens to hundreds of samples; thus, inference of OTU-OTU interaction networks is severely under-powered, and additional assumptions are required for accurate inference. Here, we present SPIEC-EASI (SParse InversE Covariance Estimation for Ecological Association Inference), a statistical method for the inference of microbial ecological interactions from metagenomic datasets that addresses both of these issues. SPIEC-EASI combines data transformations developed for compositional data analysis with a graphical model inference framework that assumes the underlying ecological interaction network is sparse. To reconstruct the interaction network, SPIEC-EASI relies on algorithms for sparse neighborhood and inverse covariance selection. Because no large-scale microbial ecological networks have been experimentally validated, SPIEC-EASI comprises computational tools to generate realistic OTU count data from a set of diverse underlying network topologies. SPIEC-EASI outperforms state-of-the-art methods in terms of edge recovery and network properties on realistic synthetic data under a variety of scenarios. SPIEC-EASI also reproducibly predicts previously unknown microbial interactions using data from the American Gut project.
△ Less
Submitted 13 February, 2015; v1 submitted 18 August, 2014;
originally announced August 2014.
-
Multimessenger Search for Sources of Gravitational Waves and High-Energy Neutrinos: Results for Initial LIGO-Virgo and IceCube
Authors:
The IceCube Collaboration,
The LIGO Scientific Collaboration,
The Virgo Collaboration,
M. G. Aartsen,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
D. Altmann,
T. Anderson,
C. Arguelles,
T. C. Arlen,
J. Auffenberg,
X. Bai,
S. W. Barwick,
V. Baum,
J. J. Beatty,
J. Becker Tjus,
K. -H. Becker,
S. BenZvi,
P. Berghaus,
D. Berley,
E. Bernardini,
A. Bernhard
, et al. (1166 additional authors not shown)
Abstract:
We report the results of a multimessenger search for coincident signals from the LIGO and Virgo gravitational-wave observatories and the partially completed IceCube high-energy neutrino detector, including periods of joint operation between 2007-2010. These include parts of the 2005-2007 run and the 2009-2010 run for LIGO-Virgo, and IceCube's observation periods with 22, 59 and 79 strings. We find…
▽ More
We report the results of a multimessenger search for coincident signals from the LIGO and Virgo gravitational-wave observatories and the partially completed IceCube high-energy neutrino detector, including periods of joint operation between 2007-2010. These include parts of the 2005-2007 run and the 2009-2010 run for LIGO-Virgo, and IceCube's observation periods with 22, 59 and 79 strings. We find no significant coincident events, and use the search results to derive upper limits on the rate of joint sources for a range of source emission parameters. For the optimistic assumption of gravitational-wave emission energy of $10^{-2}$\,M$_\odot$c$^2$ at $\sim 150$\,Hz with $\sim 60$\,ms duration, and high-energy neutrino emission of $10^{51}$\,erg comparable to the isotropic gamma-ray energy of gamma-ray bursts, we limit the source rate below $1.6 \times 10^{-2}$\,Mpc$^{-3}$yr$^{-1}$. We also examine how combining information from gravitational waves and neutrinos will aid discovery in the advanced gravitational-wave detector era.
△ Less
Submitted 9 October, 2014; v1 submitted 3 July, 2014;
originally announced July 2014.
-
Improved Upper Limits on the Stochastic Gravitational-Wave Background from 2009-2010 LIGO and Virgo Data
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca,
D. Amariutei,
M. Andersen
, et al. (824 additional authors not shown)
Abstract:
Gravitational waves from a variety of sources are predicted to superpose to create a stochastic background. This background is expected to contain unique information from throughout the history of the universe that is unavailable through standard electromagnetic observations, making its study of fundamental importance to understanding the evolution of the universe. We carry out a search for the st…
▽ More
Gravitational waves from a variety of sources are predicted to superpose to create a stochastic background. This background is expected to contain unique information from throughout the history of the universe that is unavailable through standard electromagnetic observations, making its study of fundamental importance to understanding the evolution of the universe. We carry out a search for the stochastic background with the latest data from LIGO and Virgo. Consistent with predictions from most stochastic gravitational-wave background models, the data display no evidence of a stochastic gravitational-wave signal. Assuming a gravitational-wave spectrum of Omega_GW(f)=Omega_alpha*(f/f_ref)^alpha, we place 95% confidence level upper limits on the energy density of the background in each of four frequency bands spanning 41.5-1726 Hz. In the frequency band of 41.5-169.25 Hz for a spectral index of alpha=0, we constrain the energy density of the stochastic background to be Omega_GW(f)<5.6x10^-6. For the 600-1000 Hz band, Omega_GW(f)<0.14*(f/900 Hz)^3, a factor of 2.5 lower than the best previously reported upper limits. We find Omega_GW(f)<1.8x10^-4 using a spectral index of zero for 170-600 Hz and Omega_GW(f)<1.0*(f/1300 Hz)^3 for 1000-1726 Hz, bands in which no previous direct limits have been placed. The limits in these four bands are the lowest direct measurements to date on the stochastic background. We discuss the implications of these results in light of the recent claim by the BICEP2 experiment of the possible evidence for inflationary gravitational waves.
△ Less
Submitted 12 August, 2014; v1 submitted 17 June, 2014;
originally announced June 2014.
-
First all-sky search for continuous gravitational waves from unknown sources in binary systems
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca,
D. Amariutei,
M. Andersen
, et al. (827 additional authors not shown)
Abstract:
We present the first results of an all-sky search for continuous gravitational waves from unknown spinning neutron stars in binary systems using LIGO and Virgo data. Using a specially developed analysis program, the TwoSpect algorithm, the search was carried out on data from the sixth LIGO Science Run and the second and third Virgo Science Runs. The search covers a range of frequencies from 20 Hz…
▽ More
We present the first results of an all-sky search for continuous gravitational waves from unknown spinning neutron stars in binary systems using LIGO and Virgo data. Using a specially developed analysis program, the TwoSpect algorithm, the search was carried out on data from the sixth LIGO Science Run and the second and third Virgo Science Runs. The search covers a range of frequencies from 20 Hz to 520 Hz, a range of orbital periods from 2 to ~2,254 h and a frequency- and period-dependent range of frequency modulation depths from 0.277 to 100 mHz. This corresponds to a range of projected semi-major axes of the orbit from ~0.6e-3 ls to ~6,500 ls assuming the orbit of the binary is circular. While no plausible candidate gravitational wave events survive the pipeline, upper limits are set on the analyzed data. The most sensitive 95% confidence upper limit obtained on gravitational wave strain is 2.3e-24 at 217 Hz, assuming the source waves are circularly polarized. Although this search has been optimized for circular binary orbits, the upper limits obtained remain valid for orbital eccentricities as large as 0.9. In addition, upper limits are placed on continuous gravitational wave emission from the low-mass x-ray binary Scorpius X-1 between 20 Hz and 57.25 Hz.
△ Less
Submitted 17 September, 2014; v1 submitted 30 May, 2014;
originally announced May 2014.
-
Methods and results of a search for gravitational waves associated with gamma-ray bursts using the GEO600, LIGO, and Virgo detectors
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca,
D. Amariutei,
M. Andersen,
R. A. Anderson,
S. B. Anderson
, et al. (868 additional authors not shown)
Abstract:
In this paper we report on a search for short-duration gravitational wave bursts in the frequency range 64 Hz-1792 Hz associated with gamma-ray bursts (GRBs), using data from GEO600 and one of the LIGO or Virgo detectors. We introduce the method of a linear search grid to analyse GRB events with large sky localisation uncertainties such as the localisations provided by the Fermi Gamma-ray Burst Mo…
▽ More
In this paper we report on a search for short-duration gravitational wave bursts in the frequency range 64 Hz-1792 Hz associated with gamma-ray bursts (GRBs), using data from GEO600 and one of the LIGO or Virgo detectors. We introduce the method of a linear search grid to analyse GRB events with large sky localisation uncertainties such as the localisations provided by the Fermi Gamma-ray Burst Monitor (GBM). Coherent searches for gravitational waves (GWs) can be computationally intensive when the GRB sky position is not well-localised, due to the corrections required for the difference in arrival time between detectors. Using a linear search grid we are able to reduce the computational cost of the analysis by a factor of O(10) for GBM events. Furthermore, we demonstrate that our analysis pipeline can improve upon the sky localisation of GRBs detected by the GBM, if a high-frequency GW signal is observed in coincidence. We use the linear search grid method in a search for GWs associated with 129 GRBs observed satellite-based gamma-ray experiments between 2006 and 2011. The GRBs in our sample had not been previously analysed for GW counterparts. A fraction of our GRB events are analysed using data from GEO600 while the detector was using squeezed-light states to improve its sensitivity; this is the first search for GWs using data from a squeezed-light interferometric observatory. We find no evidence for GW signals, either with any individual GRB in this sample or with the population as a whole. For each GRB we place lower bounds on the distance to the progenitor, assuming a fixed GW emission energy of $10^{-2} M_{\odot}c^{2}$, with a median exclusion distance of 0.8 Mpc for emission at 500 Hz and 0.3 Mpc at 1 kHz. The reduced computational cost associated with a linear search grid will enable rapid searches for GWs associated with Fermi GBM events in the Advanced detector era.
△ Less
Submitted 1 July, 2014; v1 submitted 5 May, 2014;
originally announced May 2014.
-
Search for gravitational radiation from intermediate mass black hole binaries in data from the second LIGO-Virgo joint science run
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca,
D. Amariutei,
M. Andersen
, et al. (825 additional authors not shown)
Abstract:
This paper reports on an unmodeled, all-sky search for gravitational waves from merging intermediate mass black hole binaries (IMBHB). The search was performed on data from the second joint science run of the LIGO and Virgo detectors (July 2009 - October 2010) and was sensitive to IMBHBs with a range up to $\sim 200$ Mpc, averaged over the possible sky positions and inclinations of the binaries wi…
▽ More
This paper reports on an unmodeled, all-sky search for gravitational waves from merging intermediate mass black hole binaries (IMBHB). The search was performed on data from the second joint science run of the LIGO and Virgo detectors (July 2009 - October 2010) and was sensitive to IMBHBs with a range up to $\sim 200$ Mpc, averaged over the possible sky positions and inclinations of the binaries with respect to the line of sight. No significant candidate was found. Upper limits on the coalescence-rate density of nonspinning IMBHBs with total masses between 100 and $450 \ \mbox{M}_{\odot}$ and mass ratios between $0.25$ and $1\,$ were placed by combining this analysis with an analogous search performed on data from the first LIGO-Virgo joint science run (November 2005 - October 2007). The most stringent limit was set for systems consisting of two $88 \ \mbox{M}_{\odot}$ black holes and is equal to $0.12 \ \mbox{Mpc}^{-3} \ \mbox{Myr}^{-1}$ at the $90\%$ confidence level. This paper also presents the first estimate, for the case of an unmodeled analysis, of the impact on the search range of IMBHB spin configurations: the visible volume for IMBHBs with nonspinning components is roughly doubled for a population of IMBHBs with spins aligned with the binary's orbital angular momentum and uniformly distributed in the dimensionless spin parameter up to 0.8, whereas an analogous population with antialigned spins decreases the visible volume by $\sim 20\%\,$.
△ Less
Submitted 19 June, 2014; v1 submitted 8 April, 2014;
originally announced April 2014.
-
Search for gravitational waves associated with gamma-ray bursts detected by the InterPlanetary Network
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca,
D. Amariutei,
M. Andersen,
R. A. Anderson,
S. B. Anderson
, et al. (879 additional authors not shown)
Abstract:
We present the results of a search for gravitational waves associated with 223 gamma-ray bursts (GRBs) detected by the InterPlanetary Network (IPN) in 2005-2010 during LIGO's fifth and sixth science runs and Virgo's first, second and third science runs. The IPN satellites provide accurate times of the bursts and sky localizations that vary significantly from degree scale to hundreds of square degr…
▽ More
We present the results of a search for gravitational waves associated with 223 gamma-ray bursts (GRBs) detected by the InterPlanetary Network (IPN) in 2005-2010 during LIGO's fifth and sixth science runs and Virgo's first, second and third science runs. The IPN satellites provide accurate times of the bursts and sky localizations that vary significantly from degree scale to hundreds of square degrees. We search for both a well-modeled binary coalescence signal, the favored progenitor model for short GRBs, and for generic, unmodeled gravitational wave bursts. Both searches use the event time and sky localization to improve the gravitational-wave search sensitivity as compared to corresponding all-time, all-sky searches. We find no evidence of a gravitational-wave signal associated with any of the IPN GRBs in the sample, nor do we find evidence for a population of weak gravitational-wave signals associated with the GRBs. For all IPN-detected GRBs, for which a sufficient duration of quality gravitational-wave data is available, we place lower bounds on the distance to the source in accordance with an optimistic assumption of gravitational-wave emission energy of $10^{-2}M_{\odot}c^2$ at 150 Hz, and find a median of 13 Mpc. For the 27 short-hard GRBs we place 90% confidence exclusion distances to two source models: a binary neutron star coalescence, with a median distance of 12Mpc, or the coalescence of a neutron star and black hole, with a median distance of 22 Mpc. Finally, we combine this search with previously published results to provide a population statement for GRB searches in first-generation LIGO and Virgo gravitational-wave detectors, and a resulting examination of prospects for the advanced gravitational-wave detectors.
△ Less
Submitted 17 April, 2014; v1 submitted 26 March, 2014;
originally announced March 2014.
-
Search for gravitational wave ringdowns from perturbed intermediate mass black holes in LIGO-Virgo data from 2005-2010
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca,
D. Amariutei,
M. Andersen,
R. Anderson
, et al. (827 additional authors not shown)
Abstract:
We report results from a search for gravitational waves produced by perturbed intermediate mass black holes (IMBH) in data collected by LIGO and Virgo between 2005 and 2010. The search was sensitive to astrophysical sources that produced damped sinusoid gravitational wave signals, also known as ringdowns, with frequency $50\le f_{0}/\mathrm{Hz} \le 2000$ and decay timescale…
▽ More
We report results from a search for gravitational waves produced by perturbed intermediate mass black holes (IMBH) in data collected by LIGO and Virgo between 2005 and 2010. The search was sensitive to astrophysical sources that produced damped sinusoid gravitational wave signals, also known as ringdowns, with frequency $50\le f_{0}/\mathrm{Hz} \le 2000$ and decay timescale $0.0001\lesssim τ/\mathrm{s} \lesssim 0.1$ characteristic of those produced in mergers of IMBH pairs. No significant gravitational wave candidate was detected. We report upper limits on the astrophysical coalescence rates of IMBHs with total binary mass $50 \le M/\mathrm{M}_\odot \le 450$ and component mass ratios of either 1:1 or 4:1. For systems with total mass $100 \le M/\mathrm{M}_\odot \le 150$, we report a 90%-confidence upper limit on the rate of binary IMBH mergers with non-spinning and equal mass components of $6.9\times10^{-8}\,$Mpc$^{-3}$yr$^{-1}$. We also report a rate upper limit for ringdown waveforms from perturbed IMBHs, radiating 1% of their mass as gravitational waves in the fundamental, $\ell=m=2$, oscillation mode, that is nearly three orders of magnitude more stringent than previous results.
△ Less
Submitted 21 May, 2014; v1 submitted 20 March, 2014;
originally announced March 2014.
-
Implementation of an F-statistic all-sky search for continuous gravitational waves in Virgo VSR1 data
Authors:
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca,
D. Amariutei,
M. Andersen,
R. Anderson,
S. B. Anderson
, et al. (826 additional authors not shown)
Abstract:
We present an implementation of the $\mathcal{F}$-statistic to carry out the first search in data from the Virgo laser interferometric gravitational wave detector for periodic gravitational waves from a priori unknown, isolated rotating neutron stars. We searched a frequency $f_0$ range from 100 Hz to 1 kHz and the frequency dependent spindown $f_1$ range from…
▽ More
We present an implementation of the $\mathcal{F}$-statistic to carry out the first search in data from the Virgo laser interferometric gravitational wave detector for periodic gravitational waves from a priori unknown, isolated rotating neutron stars. We searched a frequency $f_0$ range from 100 Hz to 1 kHz and the frequency dependent spindown $f_1$ range from $-1.6\,(f_0/100\,{\rm Hz}) \times 10^{-9}\,$ Hz/s to zero. A large part of this frequency - spindown space was unexplored by any of the all-sky searches published so far. Our method consisted of a coherent search over two-day periods using the $\mathcal{F}$-statistic, followed by a search for coincidences among the candidates from the two-day segments. We have introduced a number of novel techniques and algorithms that allow the use of the Fast Fourier Transform (FFT) algorithm in the coherent part of the search resulting in a fifty-fold speed-up in computation of the $\mathcal{F}$-statistic with respect to the algorithm used in the other pipelines. No significant gravitational wave signal was found. The sensitivity of the search was estimated by injecting signals into the data. In the most sensitive parts of the detector band more than 90% of signals would have been detected with dimensionless gravitational-wave amplitude greater than $5 \times 10^{-24}$.
△ Less
Submitted 10 April, 2014; v1 submitted 20 February, 2014;
originally announced February 2014.
-
The NINJA-2 project: Detecting and characterizing gravitational waveforms modelled using numerical binary black hole simulations
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the NINJA-2 Collaboration,
:,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca
, et al. (867 additional authors not shown)
Abstract:
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of…
▽ More
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF]
△ Less
Submitted 5 January, 2014;
originally announced January 2014.