ANNA: a toolbox for Newtonian Noise Analysis

This paper introduces ANNA, a MATLAB and Python toolbox that efficiently computes Newtonian noise for the Einstein Telescope by integrating seismic wave fields defined on 3D finite element meshes, with its accuracy verified against analytical solutions for both full-space and half-space scenarios.

The Gist

Imagine you are trying to listen to a whisper in a library, but the floor beneath you is shaking because someone is walking outside. Even though you can't see the person, their footsteps create tiny vibrations in the ground. These vibrations change the density of the soil just a little bit, and because mass attracts mass, that shifting soil pulls on the sensitive equipment in the library, making it wobble.

This is the core problem facing the Einstein Telescope (ET), a futuristic machine designed to listen to the "whispers" of the universe—gravitational waves from black holes colliding billions of light-years away. To hear these faint signals, the telescope needs to be incredibly quiet, sensitive enough to detect movements smaller than a proton. But the Earth itself is noisy.

The Invisible "Gravity Wind"

The paper introduces a tool called ANNA (a toolbox for Newtonian Noise Analysis) to solve a specific type of noise called Newtonian Noise.

Think of seismic waves (vibrations in the ground) like ripples in a pond. When these ripples pass through the soil, they squeeze and stretch the dirt. Since gravity depends on how much mass is in a specific spot, these squeezes create tiny, invisible "gravity winds" that blow on the telescope's mirrors. The mirrors move not because they are being hit, but because the gravity around them is shifting. This is the "noise" that drowns out the cosmic whispers.

What ANNA Does: The "Gravity Calculator"

ANNA is essentially a sophisticated calculator and simulator for these gravity winds.

1. The Map (The Mesh): Imagine the ground around the telescope is a giant 3D jigsaw puzzle made of tiny blocks (called "finite elements"). Some blocks are simple pyramids, others are complex bricks. ANNA can handle all these shapes.
2. The Input: Scientists feed ANNA a map of how the ground is shaking (the seismic wave field).
3. The Magic: ANNA uses a mathematical technique called "Gaussian quadrature" (think of it as a super-precise way of adding up millions of tiny numbers) to calculate exactly how much the shifting dirt will pull on the telescope's mirrors.
4. The Output: It tells the scientists two things:
   • The Total Pull: How much the mirrors will move.
   • The Breakdown: How much of that pull comes from the ground right under the mirror (bulk) versus the ground further away or at the surface.

Why It's a Big Deal

Before ANNA, figuring out this noise in complex, rocky, or uneven ground was like trying to predict the weather by looking at a single cloud. It was hard and often inaccurate.

ANNA is like a weather forecast model for gravity. It allows scientists to simulate how the telescope will behave in different types of soil, with different vibrations, before they even build it.

• The Test Drive: The authors tested ANNA by simulating simple scenarios (like a flat, uniform ground with standard waves) and compared the results to known mathematical formulas. ANNA matched the formulas perfectly.
• The Real World: They also tested it on more complex scenarios, like waves hitting the surface of the ground (Rayleigh waves), and it still worked beautifully.

The Bottom Line

ANNA is a software toolkit (available for MATLAB, Octave, and Python) that helps scientists design the Einstein Telescope. It ensures that when the telescope is finally built deep underground, the engineers know exactly how to shield it from the "gravity wind" caused by the Earth's natural vibrations. Without tools like ANNA, the telescope might be too noisy to hear the universe's most distant secrets.

In short: ANNA is the tool that helps us build a better ear for the universe by teaching us how to ignore the noise of our own planet.

Technical Summary

1. The Problem: Newtonian Noise in Gravitational Wave Detection

The Einstein Telescope (ET), a proposed third-generation underground gravitational wave observatory, aims for unprecedented sensitivity down to 3 Hz. At these low frequencies, the observatory faces a fundamental limitation known as Newtonian Noise (NN).

• Mechanism: Seismic waves (caused by anthropogenic or natural vibrations) propagate through the soil, creating density fluctuations.
• Effect: These fluctuations generate varying gravitational attractions that directly perturb the position of the laser interferometer's mirrors (test masses).
• Challenge: Unlike seismic noise that can be mechanically isolated, Newtonian noise acts directly via gravity and cannot be shielded. Accurately modeling this coupling is critical for designing mitigation strategies (e.g., sensor arrays and feed-forward subtraction) for the ET.

2. Methodology: The ANNA Toolbox

The paper introduces ANNA (Newtonian Noise Analysis), a computational toolbox designed to calculate Newtonian noise from complex seismic wave fields.

• Core Approach: The tool computes NN by integrating density fluctuations derived from seismic wave fields over the soil domain surrounding a test mass.
• Numerical Technique:
  • It utilizes Gaussian quadrature for integration.
  • It operates on 3D finite element meshes (FEM).
• Mesh Compatibility: The toolbox supports a wide variety of element types to accommodate complex geometries:
  • Tetrahedral elements: Linear (4-node) and Quadratic (10-node).
  • Brick (Hexahedral) elements: Linear (8-node) and Quadratic (20-node).
• Workflow:
  1. The user defines or interpolates a seismic wave field onto a finite element mesh.
  2. ANNA computes the total Newtonian noise.
  3. It decomposes the results into bulk (volume) and surface contributions.
• Software Environment:
  • Primary implementation: MATLAB.
  • Compatibility: GNU Octave (open-source alternative).
  • Availability: A Python version is also provided.

3. Key Contributions

• Versatile FEM Framework: Unlike previous methods that may rely on simplified analytical models or specific grid types, ANNA supports heterogeneous media through flexible 3D meshing (linear and quadratic elements).
• Decomposition Capability: The ability to separate bulk and surface contributions allows researchers to pinpoint the specific sources of noise within the soil volume versus the surface, aiding in targeted mitigation.
• Open Accessibility: By providing versions for MATLAB, Octave, and Python, the tool lowers the barrier to entry for the gravitational wave community.
• Physical Consistency: The framework is designed to handle complex, real-world geological heterogeneities that analytical solutions cannot easily address.

4. Results and Verification

The authors verified the toolbox against known analytical solutions in two distinct scenarios, demonstrating "excellent agreement":

1. Homogeneous Full Space:

   • Scenario: Plane P- and S-waves propagating in an elastic, homogeneous full space.
   • Setup: A mirror suspended in a spherical cavity.
   • Assumption: The wavelength is much larger than the cavity radius, allowing wave scattering to be ignored.
   • Outcome: ANNA's numerical results matched the analytical solutions perfectly.

2. Homogeneous Half-Space:

   • Scenario: A Rayleigh wave propagating on the free surface of a homogeneous elastic half-space.
   • Setup: A test mass suspended at a finite distance above the surface.
   • Outcome: The toolbox reported similar high-level agreement with analytical solutions.

5. Significance

The ANNA toolbox represents a significant step forward in the preparation for the Einstein Telescope and future gravitational wave observatories.

• Enabling Complex Modeling: It provides a computationally efficient and physically consistent method to model gravitational-seismic coupling in heterogeneous media, which is essential for realistic underground site characterization.
• Design Optimization: By accurately predicting Newtonian noise in complex geological settings, engineers can better design sensor arrays and subtraction algorithms to mitigate this noise floor, ensuring the ET can achieve its 3 Hz sensitivity goal.
• Community Resource: As an open-source, multi-platform tool, it standardizes the approach to NN analysis, facilitating collaboration and reproducibility within the scientific community.