skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, April 16 until 2:00 AM ET on Friday, April 17 due to maintenance. We apologize for the inconvenience.


Title: Detection of anomalies amongst LIGO’s glitch populations with autoencoders
Abstract Gravitational wave (GW) interferometers are able to detect a change in distance of ~1/10 000th the size of a proton. Such sensitivity leads to large rates of non-gaussian, transient bursts of noise, also known as glitches, which hinder the detection and parameter estimation of short- and long-lived GW signals in the main detector strain. Glitches, come in a wide range of frequency-amplitude-time morphologies and may be caused by environmental or instrumental processes, so a key step towards their mitigation is to understand their population. Current approaches for their identification use supervised models to learn their morphology in the main strain with a fixed set of classes, but do not consider relevant information provided by auxiliary channels that monitor the state of the interferometers. In this work, we present an unsupervised algorithm to find anomalous glitches. Firstly, we encode a subset of auxiliary channels from Laser Interferometer Gravitational-Wave Observatory Livingston in the fractal dimension (FD), which measures the complexity of the signal. For this aim, we speed up the fractal dimension calculation to encode 1 h of data in 11 s. Secondly, we learn the underlying distribution of the data using an autoencoder with cyclic periodic convolutions. In this way, we learn the underlying distribution of glitches and we uncover unknown glitch morphologies, and overlaps in time between different glitches and misclassifications. This led to the discovery of 6.6 % anomalies in the input data. The results of this investigation stress the learnable structure of auxiliary channels encoded in FD and provide a flexible framework for glitch discovery.  more » « less
Award ID(s):
2308693
PAR ID:
10510717
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Institute of Physics
Date Published:
Journal Name:
Classical and Quantum Gravity
Volume:
41
Issue:
5
ISSN:
0264-9381
Page Range / eLocation ID:
055004
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Objective.In vivoimaging assessments of skeletal muscle structure and function allow for longitudinal quantification of tissue health. Magnetic resonance elastography (MRE) non-invasively quantifies tissue mechanical properties, allowing for evaluation of skeletal muscle biomechanics in response to loading, creating a better understanding of muscle functional health.Approach. In this study, we analyze the anisotropic mechanical response of calf muscles using MRE with a transversely isotropic, nonlinear inversion algorithm (TI-NLI) to investigate the role of muscle fiber stiffening under load. We estimate anisotropic material parameters including fiber shear stiffness ( μ 1 ), substrate shear stiffness ( μ 2 ), shear anisotropy ( ϕ ), and tensile anisotropy ( ζ ) of the gastrocnemius muscle in response to both passive and active tension.Main results. In passive tension, we found a significant increase in μ 1 , ϕ , and ζ with increasing muscle length. While in active tension, we observed increasing μ 2 and decreasing ϕ and ζ during active dorsiflexion and plantarflexion—indicating less anisotropy—with greater effects when the muscles act as agonist.Significance. The study demonstrates the ability of this anisotropic MRE method to capture the multifaceted mechanical response of skeletal muscle to tissue loading from muscle lengthening and contraction. 
    more » « less
  2. Abstract Surface properties are critical to the capabilities of superconducting microwave devices. The native oxide of niobium-based devices is thought to consist of a thin normal conducting layer. To improve understanding on the importance of this layer, an attempt was made to replace it with a more easily controlled gold film. A niobium sample host microwave cavity was used to measure the surface resistance in continuous wave operation at 4.0 GHz and 5.2 GHz . Sample conditions studied include temperatures ranging from 1.6 K to 4.2 K with RF magnetic fields on the sample surface ranging from 1 mT to the maximum field before the superconducting properties were lost (quench field). The nominal film thickness of the gold layer was increased from 0.1 nm to 2.0 nm in five steps to study the impact of the normal layer thickness on surface resistance on a single niobium substrate. The 0.1 nm film was found to reduce the surface resistance of the sample and to enhance the quench field. With the exception of the final step from a 1.5 nm gold film to 2.0 nm , the magnitude of the surface resistance increased substantially with gold film thickness. The nature of the surface resistance field-dependence appeared to be roughly independent from the gold layer thickness. This initial study provides new perspectives and suggests avenues for optimizing and designing surfaces for resonant cavities in particle accelerators and quantum information applications. 
    more » « less
  3. Abstract This paper considers the Westervelt equation, one of the most widely used models in nonlinear acoustics, and seeks to recover two spatially-dependent parameters of physical importance from time-trace boundary measurements. Specifically, these are the nonlinearity parameter κ ( x ) often referred to as B / A in the acoustics literature and the wave speed c 0 ( x ) . The determination of the spatial change in these quantities can be used as a means of imaging. We consider identifiability from one or two boundary measurements as relevant in these applications. For a reformulation of the problem in terms of the squared slowness s = 1 / c 0 2 and the combined coefficient η = κ c 0 2 we devise a frozen Newton method and prove its convergence. The effectiveness (and limitations) of this iterative scheme are demonstrated by numerical examples. 
    more » « less
  4. Abstract Polyatomic molecules have been identified as sensitive probes of charge-parity violating and parity violating physics beyond the Standard Model (BSM). For example, many linear triatomic molecules are both laser-coolable and have parity doublets in the ground electronic X ˜ 2 Σ + ( 010 ) state arising from the bending vibration, both features that can greatly aid BSM searches. Understanding the X ˜ 2 Σ + ( 010 ) state is a crucial prerequisite to precision measurements with linear polyatomic molecules. Here, we characterize the fundamental bending vibration of 174 YbOH using high-resolution optical spectroscopy on the nominally forbidden X ˜ 2 Σ + ( 010 ) A ˜ 2 Π 1 / 2 ( 000 ) transition at 588 nm. We assign 39 transitions originating from the lowest rotational levels of the X ˜ 2 Σ + ( 010 ) state, and accurately model the state’s structure with an effective Hamiltonian using best-fit parameters. Additionally, we perform Stark and Zeeman spectroscopy on the X ˜ 2 Σ + ( 010 ) state and fit the molecule-frame dipole moment to D m o l = 2.16 ( 1 ) Dand the effective electrong-factor to g S = 2.07 ( 2 ) . Further, we use an empirical model to explain observed anomalous line intensities in terms of interference from spin–orbit and vibronic perturbations in the excited A ˜ 2 Π 1 / 2 ( 000 ) state. Our work is an essential step toward searches for BSM physics in YbOH and other linear polyatomic molecules. 
    more » « less
  5. Abstract We present a theory based on the conventional two-term (i.e. Lorentzian) approximation to the exact solution of the Boltzmann equation in non-magnetized weakly ionized plasma to efficiently obtain the electron rate and transport coefficients in a magnetized plasma for an arbitrary magnitude and direction of applied electric field E and magnetic field B . The proposed transcendental method does not require the two-term solution of the Boltzmann equation in magnetized plasma, based on which the transport parameters vary as a function of the reduced electric field E / N , reduced electron cyclotron frequency ω c e / N , and angle E , B between E and B vectors, whereNis the density of neutrals. Comparisons between the coefficients derived from BOLSIG+’s solution (obtained via the two-term expansion when B 0 ) and coefficients of the presented method are illustrated for air, a mixture of molecular hydrogen (H2) and helium (He) representing the giant gas planets of the Solar System, and pure carbon dioxide (CO2). The new approach may be used in the modeling of magnetized plasma encountered in the context of transient luminous events, e.g. sprite streamers in the atmosphere of Earth and Jupiter, in modeling the propagation of lightning’s electromagnetic pulses in Earth’s ionosphere, and in various laboratory and industrial applications of nonthermal plasmas. 
    more » « less