warp-sensitivity-analysis

Multi-scale theoretical physics testing framework for warp drive signatures.

Overview

This package provides a comprehensive analysis pipeline spanning three regimes:

1. Signal vs. Noise Analysis: Compare synthetic warp signals against detector sensitivity curves
2. Semi-Classical Testing: Generate Post-Newtonian corrections and compare with precision experiments
3. Strong-Curvature Models: Explore Planck-scale physics with toy black hole and cosmological models

Repository Structure

. ├── README.md ├── EXTENDED_FRAMEWORK.md # Detailed framework documentation
├── mock_data.ndjson # Synthetic signal time-series/spectra JSON-lines ├── mock_data.am # AsciiMath metadata for mock_data ├── sensitivity_curve.csv # Generated noise curve (frequency, noise PSD) ├── sensitivity_curve.am # AsciiMath metadata for sensitivity curve ├── generate_sensitivity_curve.py # Noise curve generator
├── analyze_sensitivity.py # Original signal vs noise analysis
├── ├── semi_classical/ # Post-Newtonian analysis
│ ├── compute_pn_corrections.py # PN expansion generator
│ ├── analyze_pn_tests.py # Experimental comparison
│ ├── theory_params.am # Theory configuration
│ ├── pn_config.am # PN expansion settings
│ └── pn_data/ # Experimental datasets
│ ├── ligo_data.csv # LIGO sensitivity curves
│ ├── ligo_data.am # LIGO metadata
│ └── atomic_interf.am # Atomic interferometry data
│
└── strong_curvature/ # Planck-scale models
├── generate_2d_blackhole.py # 2D black hole toy models
├── minisuperspace_cosmo.py # FRW minisuperspace cosmology
├── compare_strong_models.py # Model unification
├── blackhole_config.am # Black hole parameters
└── cosmo_config.am # Cosmology parameters
```markdown
## Prerequisites
- Python 3.7+ - NumPy - SciPy
- SymPy (for symbolic PN calculations)
- ndjson - tqdm (optional, for progress bars)
Install via:
```bash
pip install numpy scipy sympy python-ndjson tqdm

Usage

1. Generate Noise Curve

python generate_sensitivity_curve.py

This produces:

• sensitivity_curve.csv

  • Columns: frequency_Hz, noise_strain_per_sqrtHz

• sensitivity_curve.am

  • One-line AsciiMath metadata: model name, parameters, point count

2. Analyze Signal Detectability

python analyze_sensitivity.py ` --mock mock_data.ndjson `
--meta mock_data.am ` --noise sensitivity_curve.csv `
--nmeta sensitivity_curve.am ` --out sensitivity_comparison.ndjson `
--oam sensitivity_comparison.am

Outputs:

• sensitivity_comparison.ndjson
  JSON-line records:

{
"label": "signal1",
"detectable": true,
"snr": 12.5
}

• sensitivity_comparison.am
  AsciiMath summary: detector name, injection count, detection thresholds

3. Post-Newtonian Analysis

Generate PN corrections for warp drive theories:

cd semi_classical/
python compute_pn_corrections.py \
--theory theory_params.am \
--pn-config pn_config.am \
--out pn_waveforms.ndjson \
--oam pn_summary.am

Configuration files:

• theory_params.am: Contains WarpVelocity, MassParameter, TheoryType
• pn_config.am: Contains MaxPNOrder, FrequencyRange, ExpansionParameter

Outputs:

• pn_waveforms.ndjson: PN corrections with observational signatures
• pn_summary.am: Analysis metadata

Analyze against experimental data:

python analyze_pn_tests.py \
--pn-data pn_waveforms.ndjson \
--pn-meta pn_summary.am \
--exp-data pn_data/ligo_data.csv \
--exp-meta pn_data/ligo_data.am \
--out pn_analysis.ndjson \
--oam pn_analysis.am \
--snr-threshold 5.0

Flags and Input Specifications:

• --pn-data: PN corrections from compute_pn_corrections.py (.ndjson format)
  • Each entry contains: pn_order, correction (symbolic expressions), signature (frequency-domain observational signatures), v_over_c, mass_parameter
• --pn-meta: PN metadata file (.am format) containing analysis parameters
  • Required keys: max_pn_order, warp_velocity, frequency_range, n_corrections
• --exp-data: Experimental sensitivity data (.csv with columns: frequency_hz, strain_sensitivity, experiment_type)
  • Compatible with LIGO, Virgo, atomic interferometry, and pulsar timing datasets
• --exp-meta: Experimental metadata (.am with keys: ExperimentType, InstrumentName, FrequencyRange, StrainSensitivity, DataSource)
• --snr-threshold: SNR threshold for detectability (default: 5.0, recommended range: 3.0-10.0)

Outputs:

• pn_analysis.ndjson: Detection analysis for each PN order with fields:
  • pn_order: PN correction order ("1PN", "2PN", "3PN", etc.)
  • theory_params: Contains v_over_c and mass_parameter values from theory
  • experimental_analysis: Per-experiment detectability results including SNR, parameter bounds, frequency ranges
  • overall_detectability: Boolean indicating if any experiment can detect this PN order
  • best_snr: Highest SNR achieved across all experimental datasets
• pn_analysis.am: Summary metadata including detection statistics, parameter constraints, and experiment coverage

4. Strong-Curvature Models

Generate 2D black hole curvature data:

cd strong_curvature/
python generate_2d_blackhole.py \
--model-config blackhole_config.am \
--out blackhole_data.ndjson \
--oam blackhole_summary.am

Configuration (blackhole_config.am):

• ModelType: "2d_schwarzschild" or "warp_bubble"
• MassParameter: Characteristic mass scale
• PlanckLength: Planck length (default: 1e-35)
• BubbleThickness: For warp bubble models

Outputs:

• blackhole_data.ndjson: Curvature invariants and quantum parameters
• blackhole_summary.am: Model summary and regime classification

Generate minisuperspace cosmology:

python minisuperspace_cosmo.py \
--cosmo-config cosmo_config.am \
--out cosmo_data.ndjson \
--oam cosmo_summary.am

Compare and unify strong-curvature models:

python compare_strong_models.py \
--models blackhole_data.ndjson cosmo_data.ndjson \
--meta blackhole_summary.am cosmo_summary.am \
--out unified_strong_models.ndjson \
--oam unified_summary.am

Flags and Input Specifications:

• --models: List of model data files (.ndjson) from toy model generators
  • Multiple files can be specified: --models file1.ndjson file2.ndjson file3.ndjson
  • Each .ndjson contains model results with curvature invariants, geodesic analysis, and Planck-scale physics indicators
  • Expected structure: model_type, curvature_analysis, geodesic_analysis, planck_scale_analysis
• --meta: Corresponding metadata files (.am) with model parameters and analysis settings
  • Must match order of --models files: --meta meta1.am meta2.am meta3.am
  • Required keys: model_type, analysis parameters specific to each model class

Outputs:

• unified_strong_models.ndjson: Unified comparison with regime classification, containing:
  • model_id: Unique identifier for each model in the comparison
  • model_type: Type of model ("2d_schwarzschild", "warp_bubble", "frw_minisuperspace", etc.)
  • regime_classification: Overall physical regime ("classical", "semi_classical", "quantum_gravity")
  • curvature_analysis: Extracted curvature scales and quantum gravity parameters
  • classical_gr_valid: Boolean indicating if classical General Relativity is adequate
  • quantum_correction_strength: Estimated strength of quantum gravity corrections (0.0-1.0)
  • requires_quantum_gravity: Boolean indicating if quantum effects dominate
  • parameter_ranges: Valid parameter ranges for classical and quantum regimes
• unified_summary.am: Summary statistics including:
  • regime_distribution: Count of models in each regime (classical/semi_classical/quantum_gravity)
  • n_quantum_gravity: Number of models requiring full quantum gravity treatment
  • n_classical_valid: Number of models where classical GR remains valid
  • average_quantum_correction: Average quantum correction strength across all models
  • model_types_analyzed: List of model types included in comparison
  • planck_scale_physics_important: Boolean indicating if any models reach Planck-scale physics

Input/Output Formats

AsciiMath Metadata Files (.am)

All configuration and metadata files use the AsciiMath format with consistent key-value pairs:

[ key1 = value1, key2 = "string_value", key3 = 1.23e-4, ... ]

Configuration Keys by Script:

Semi-Classical Analysis:

• theory_params.am:

  • WarpVelocity: Dimensionless warp velocity (v/c ratio)
  • MassParameter: Characteristic mass scale parameter
  • TheoryType: Theory model type ("alcubierre_warp", "van_den_broeck", etc.)
  • MetricSignature: Metric signature convention ("mostly_plus", "mostly_minus")

• pn_config.am:

  • MaxPNOrder: Maximum Post-Newtonian order to compute (integer: 1, 2, 3, ...)
  • FrequencyRange: Comma-separated frequency range in Hz ("10,1000")
  • ExpansionParameter: PN expansion parameter ("v_over_c", "frequency")
  • SymbolicComputation: Boolean for symbolic vs numerical computation

Strong-Curvature Models:

• blackhole_config.am:

  • ModelType: Model type ("2d_schwarzschild", "warp_bubble")
  • MassParameter: Characteristic mass scale in natural units
  • PlanckLength: Planck length in meters (default: 1e-35)
  • BubbleThickness: Warp bubble thickness parameter (for warp_bubble models)

• cosmo_config.am:

  • ModelType: Cosmology model ("frw_minisuperspace", "kasner")
  • HubbleParameter: Hubble parameter H₀ in units of 100 km/s/Mpc
  • OmegaMatter: Matter density parameter Ωₘ
  • PlanckLength: Planck length in meters
  • ScaleFactorRange: Comma-separated scale factor range ("1e-10,1e10")

Experimental Data:

• ligo_data.am, atomic_interf.am:
  • ExperimentType: Type of experiment ("gravitational_wave", "atomic_interferometry")
  • InstrumentName: Specific instrument name ("LIGO_H1", "LIGO_L1", "atom_interferometer")
  • FrequencyRange: Sensitive frequency range in Hz
  • StrainSensitivity: Characteristic strain sensitivity
  • DataSource: Source of experimental data

JSON-Lines Data Files (.ndjson)

• One JSON object per line for efficient streaming and processing

• For mock_data.ndjson:

  { "label": "signal1", "time_series": [0.1, 0.2, ...], "sampling_rate": 16384 }

• For pn_waveforms.ndjson:

  { "pn_order": "1PN", "correction": {...}, "signature": {...}, "v_over_c": 0.1, "mass_parameter": 1e-6 }

• For blackhole_data.ndjson:

  { "model_type": "2d_schwarzschild", "curvature_analysis": {...}, "geodesic_analysis": {...} }

CSV Data Files (.csv)

• Standard comma-separated with header row
• For sensitivity curves: frequency_Hz, noise_strain_per_sqrtHz
• For experimental data: frequency_hz, strain_sensitivity, experiment_type

Examples

# 1. Build the noise curve
python generate_sensitivity_curve.py
# 2. Run detection analysis
python analyze_sensitivity.py \
--mock mock_data.ndjson \
--meta mock_data.am \
--noise sensitivity_curve.csv \
--nmeta sensitivity_curve.am \
--out sensitivity_comparison.ndjson \
--oam sensitivity_comparison.am
# 3. Generate and analyze PN correctionscd semi_classical/
python compute_pn_corrections.py \
--theory theory_params.am \
--pn-config pn_config.am \
--out pn_waveforms.ndjson \
--oam pn_summary.am
python analyze_pn_tests.py \
--pn-data pn_waveforms.ndjson \
--pn-meta pn_summary.am \
--exp-data pn_data/ligo_data.csv \
--exp-meta pn_data/ligo_data.am \
--out pn_analysis.ndjson \
--oam pn_analysis.am
# 4. Generate strong-curvature modelscd ../strong_curvature/
python generate_2d_blackhole.py \
--model-config blackhole_config.am \
--out blackhole_data.ndjson \
--oam blackhole_summary.am
python compare_strong_models.py \
--models blackhole_data.ndjson \
--meta blackhole_summary.am \
--out unified_models.ndjson \
--oam unified_summary.am
# 5. Inspect results
head ../sensitivity_comparison.ndjson
cat ../sensitivity_comparison.am
head pn_analysis.ndjson
head unified_models.ndjson