{
  "classification": "Original finite-source numerical benchmark using published optical equations; not a mission forecast or reproduction of published planet images",
  "sources": [
    {
      "doi": "10.1103/PhysRevD.100.084018",
      "url": "https://arxiv.org/pdf/1908.01948",
      "locators": "Eq9 PDF p5; Eqs81-82 PDF p18; Eq83 PDF p19"
    },
    {
      "doi": "10.1103/PhysRevD.103.124038",
      "url": "https://arxiv.org/abs/2012.05477",
      "locators": "Eqs1,4-6; forward and inverse image models"
    }
  ],
  "parameters": {
    "wavelength_m": 1e-06,
    "z_AU": 650,
    "rg_m": 2953.25,
    "AU_m": 149597870700.0,
    "aperture_diameter_m": 1.0,
    "alpha_per_m": 48.96949080575494,
    "peak_point_gain": 116589636790.06859,
    "first_zero_radius_m": 0.04910864944940688,
    "on_axis_aperture_gain": 3015856273.7201314,
    "grid_side": 24,
    "field_width_m": 6.0,
    "cell_width_m": 0.25,
    "radial_kernel_step_m": 0.001,
    "radial_kernel_max_m": 9,
    "cell_quadrature_order": 32,
    "pupil_radial_order": 96,
    "pupil_angular_samples": 192,
    "noise_seed": 271828
  },
  "operator": "A_ij = integral over square source cell j of aperture-area-average J0(alpha|x_i-u+r|)^2 d2u. Fixed observation centers; zero source outside 6m square; no circular convolution or row normalization. Source cell brightness is dimensionless; clean signal has relative brightness times m2 units, with common gain omitted.",
  "inverse": "Unconstrained Tikhonov: V diag(s/(s2+lambda2)) U^T y; lambda = control times largest singular value. No positivity clipping in computations or downloadable data. Exact same operator generates and reconstructs measurements; zero-noise recovery is a code consistency check, not physical validation.",
  "noise": "Fixed Gaussian sample shared across cases, sigma = selected fraction times maximum clean signal. It is analyst-selected noise, not a photon count, corona calculation or the published output SNR7.",
  "mismatch_check": "Additional no-noise signal generated with a 1.02m aperture and reconstructed using the nominal 1m aperture, regularization0.01. This tests one deliberately chosen PSF error, not a measured uncertainty or comprehensive robustness analysis.",
  "display": "True and recovered brightness use fixed [0,1] scale, error [-1,1]. Out-of-range values saturate only for display; extrema/counts and unmodified data are supplied. Measurement color scale fixed [0,max clean].",
  "limitations": [
    "Artificial small patch, zero exterior: omits contamination from the rest of a planet.",
    "Spherical static Sun and monochromatic source; no solar multipoles, pointing error or PSF mismatch.",
    "No telescope focal-plane coronagraph/detector model, spectral response, photon budget, real corona noise or exposure prediction.",
    "No rotating planet, cloud evolution, scan timing or navigation dynamics.",
    "No source-matched NIAC image reproduction or achieved planetary resolution."
  ],
  "checks": {
    "analytic_on_axis_relative_average": 0.025867275658045707,
    "quadrature_on_axis_absolute_error": 6.938893903907228e-17,
    "pupil_refinement_max_absolute_error_15_radii": 2.2898349882893854e-16,
    "cell_8_to_16_relative_frobenius": 3.315981836508818e-05,
    "cell_16_to_32_relative_frobenius": 1.0760454001415865e-07,
    "radial_interpolation_2mm_to_1mm_relative_frobenius": 1.3141569215677153e-06,
    "operator_symmetry_max_absolute": 0.0,
    "operator_min_entry": 4.998678378044934e-05,
    "same_model_no_noise_unregularized_rmse": 2.298373825090486e-12,
    "condition_number": 292846.89673757314,
    "different_quadrature_no_noise_unregularized_rmse": 1.8632097526519333e-06,
    "aperture_mismatch_1p02m_generated_1m_inverse_no_noise_r0p01_rmse": 0.06570859166937547,
    "cell_refinement_forward_pattern_relative_l2": 1.4789011047848709e-08
  },
  "software": {
    "python": "3.12.10",
    "numpy": "2.5.3",
    "scipy": "1.18.1",
    "matplotlib": "3.11.2"
  }
}