{"id":"58072975-5665-49e0-80fd-165ed65a115f","arxiv_id":"2607.11006","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"AMIGO reconstructs all-sky radio maps from lunar-orbit array visibilities by combining mini-batch gradient descent with ALM priors, cutting memory and time costs enough for NSIDE=256 maps on mock DSL data.","lead":"AMIGO is a new radio-imaging algorithm that reconstructs all-sky maps from lunar-orbit interferometers by pairing mini-batch gradient descent with augmented Lagrangian multipliers. It makes high-resolution imaging computationally practical for missions like DSL that face huge data volumes and incomplete 3D baselines.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified beyond the reader's already-flagged mock-only limitation.","rationale":"The paper's strongest claim is algorithmic feasibility (reduced scaling + controllable priors) demonstrated on realistic-orbit mocks. The reader already isolates the correct weakest assumption—idealized mocks omit calibration, baseline and beam systematics that the authors themselves flag as future work (§6). No deeper internal flaw (e.g., incorrect gradient derivation in Appendix A, broken ALM residual control, or mis-reported scaling) is evident from the equations, tables or figures. Hierarchical initialization and batch-order robustness (§5.2–5.3) further support reliability under the tested conditions. Consequently the CONDITIONAL verdict and HIGH confidence remain appropriate; an ACCEPT would still require end-to-end tests with realistic systematics and public code. The concrete test above simply operationalizes the already-acknowledged gap without inventing a new objection.","tokens_in":20159,"tokens_out":488,"duration_ms":5527,"concrete_test":"Re-run the NSIDE=64, 3 MHz case of Table 3 after injecting the phase errors from baseline determination and time synchronization quantified in Zhou et al. (2026) at the levels expected for DSL; if wall-clock time or iteration count rises by more than a factor of ~3 or MSE exceeds ~0.15, the pure-mock feasibility claim would need quantitative qualification.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader's weakest_assumption correctly identifies the softest point: all validation uses idealized mocks (opaque Moon, perfect beams, white Gaussian noise, no calibration/baseline errors; §3, §6). That limitation is already stated by the authors and underpins the CONDITIONAL verdict. No additional load-bearing internal inconsistency appears. The claimed scalings (memory O(N) vs O(N^{2}), time <O(MN^{2}) vs O(N^{3}+MN^{2})) are supported by the measured wall-clock/memory numbers and fitted templates in §5.4/Table 3/Figure 7 for NSIDE 16–64; hierarchical initialization and ALM residual control (§2.2–2.3, §5.1–5.3) behave as described on the mocks. The central feasibility claim therefore holds under the paper's stated premises.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper introduces AMIGO, an imaging algorithm for lunar-orbiting interferometers such as DSL that combines Mini-Batch Gradient Descent with the Augmented Lagrangian Multiplier method. It formulates all-sky reconstruction as a constrained optimization problem (data fidelity plus non-negative temperature and angular-power-spectrum equality constraints), converts it via ALM into an unconstrained problem whose first sub-step is solved by MBGD, and controls prior strength through residual thresholds ΔH and ΔG. Mock visibilities generated under a realistic DSL orbit (breathing baselines, lunar shading, short-dipole beams) are used to reconstruct maps at 3, 10 and 30 MHz up to NSIDE = 256. The authors report MSE/SSIM/ρℓ metrics, demonstrate hierarchical initialization, and show that memory and wall-clock scaling improve relative to the Tikhonov-regularization matrix-inversion approach of Deng et al. (2026).","tokens_in":20399,"tokens_out":1045,"duration_ms":8203,"significance":"If the claimed scalings hold under more realistic systematics, AMIGO would remove a genuine computational bottleneck for all-sky imaging with dynamic 3-D baselines, making NSIDE ≥ 256 reconstructions feasible on modest hardware. The ALM residual-threshold mechanism supplies a transparent, adjustable way to inject physically motivated priors (non-negativity and Cℓ) without hard-wiring a full prior map. The hierarchical-initialization strategy and the explicit GPU/CPU benchmarks against D26 are concrete engineering contributions that other lunar-array concepts can reuse. The work is therefore a useful algorithmic step for the DSL and similar missions, even though it remains a proof-of-concept on idealized mocks.","major_comments":[{"comment":"§6 and §3 explicitly limit validation to mocks that omit calibration errors, baseline uncertainties and beam-model errors (opaque Moon, perfect antennas pointed at lunar center, white Gaussian noise only). The central feasibility claim for real DSL data therefore rests on an untested extrapolation. A quantitative stress test—e.g., injecting the phase errors already characterized in Zhou et al. (2026) and re-measuring MSE/SSIM/ρℓ—is needed before the claim can be regarded as demonstrated rather than projected.","section":null},{"comment":"§5.1 and Table 2: when the prior Cℓ is scaled by 1.5 the reconstruction develops artificial stripes and MSE rises from 0.053 to 0.303. The only diagnostic offered is “unacceptable iteration count.” Because a true Cℓ will never be known a priori, the paper should supply a practical, data-driven procedure (or at least a quantitative criterion) for detecting and correcting such mismatches without access to the input map.","section":null}],"minor_comments":[{"comment":"Table 3 / Figure 7: wall-clock times are hardware-specific (RTX 3090 + EPYC 9654). State the software stack (CUDA version, BLAS library) and report FLOPs or iteration-normalized timings so that the O(Nvis Npix) scaling can be verified independently.","section":null},{"comment":"Equation (9) and surrounding text: the learning-rate schedule (initial α = 1.0, decay 0.1, increase 1.2, floor 10−7) is described but never ablated. A short sensitivity check would strengthen the claim that imaging quality is insensitive to batch configuration.","section":null},{"comment":"Figures 1–3: color-bar ranges differ between frequencies; a common fractional-residual scale or an additional panel of absolute residuals would aid visual comparison.","section":null},{"comment":"Appendix A: the factor of 2 that appears after the spherical-harmonic chain rule is correct for real maps, but a one-sentence reminder that aℓ,−m = (−1)m a*ℓm would help readers less familiar with HEALPix conventions.","section":null},{"comment":"Typographical: “Augmented lagrangian” should be consistently capitalized; “breathing” strategy is introduced without a citation to the earlier DSL papers that define it.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a solid incremental contribution within the DSL series. The mock-only limitation is already flagged by the authors and does not constitute an internal inconsistency; a minor-revision request for a modest systematics injection and a clearer prior-mismatch protocol is proportionate. Fit for the journal is good."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a practical methods paper that does what it claims. The new piece is not inventing MBGD or ALM, but wiring them together for the DSL geometry—all-sky FOV, dynamic 3-D baselines, Moon blockage—and giving quantitative control of prior strength through residual thresholds on C_ℓ and non-negativity. That combination, plus hierarchical initialization, is what lets them reach NSIDE=256 where their earlier Tikhonov matrix approach (D26) stalls.\n\nWhat they do well is concrete. The optimization setup in §2 is clean, the gradient derivation in Appendix A is correct, and the scaling benchmarks in §5.4 (Table 3, Figure 7) show memory dropping from O(N^{2}) to O(N) and wall-clock scaling better than the matrix route. Imaging metrics (MSE, SSIM, multipole correlation) are reported across frequencies and resolutions, and they beat no-prior TR when the comparison is fair. Prior-mismatch tests and batch-order checks are sensible; the authors flag that large prior errors show up as excessive iteration counts. They also state the obvious limit: everything is on idealized mocks (opaque Moon, perfect beams, white noise, no calibration or baseline errors).\n\nSoft spots are real but proportionate. The free parameters (Δ_thres, batch size, learning-rate schedule, ρ) are tuned by hand; the paper shows the results are not wildly sensitive inside the tested range, but a real pipeline will still need a recipe. Code is “available on request,” which is weaker than a public release for a methods claim. None of that breaks the central feasibility argument under the premises they actually test.\n\nThis is for people building low-frequency space interferometers or writing the next map-maker for lunar-orbit arrays. A serious referee should see it; the math and the mock evidence are solid enough to deserve that time. I would cite the scaling and prior-control results if I were working on similar imaging, and I would accept it for peer review.","headline":"Solid engineering paper: AMIGO makes all-sky lunar-orbit imaging computationally tractable on mocks via MBGD+ALM, with clean scaling and honest limits.","tokens_in":21024,"tokens_out":497,"would_cite":true,"duration_ms":5668,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"AMIGO makes all-sky imaging from a lunar-orbit radio array computationally feasible by cutting memory and time costs while folding in controllable physical priors.","keywords":["radio interferometry","aperture synthesis","lunar orbit array","all-sky imaging","mini-batch gradient descent","augmented Lagrangian","angular power spectrum prior"],"falsifier":"Run AMIGO on the same mock data after injecting realistic baseline-determination phase errors or primary-beam model mismatches; if the residual maps and MSE degrade beyond the levels reported for pure thermal noise, the claim of practical feasibility fails.","tokens_in":21024,"feed_emoji":"🌕","tokens_out":660,"duration_ms":5070,"temperature":0.7,"pith_summary":"Low-frequency radio astronomy below 30 MHz is blocked from the ground by the ionosphere and Earth-based interference, so missions like the proposed DSL lunar-orbit array aim to open that window. The catch is that an orbiting array sees nearly the whole sky with constantly changing three-dimensional baselines, and classical image reconstruction becomes prohibitively expensive in memory and time. This paper introduces AMIGO, which replaces full matrix inversion with mini-batch gradient descent and wraps physical constraints (non-negative sky temperature and a prior angular power spectrum) inside an augmented Lagrangian multiplier. The result is memory that scales only with the number of sky pixels and wall-clock time that stays well below the cubic cost of traditional methods. Tests on mock DSL visibilities recover diffuse structure and bright sources from 3 to 30 MHz up to NSIDE=256, showing that high-resolution all-sky maps from lunar orbit are now within practical reach.","feed_headline":"Lunar-orbit radio maps become computable at high resolution","feed_subtitle":"AMIGO cuts memory and time so DSL-style arrays can image the whole sky below 30 MHz","key_machinery":"The AMIGO iteration: an outer ALM loop that updates Lagrange multipliers and penalty parameters for the equality constraint H_ℓ=0 (power-spectrum match) and inequality G_n≥0 (positive temperature), and an inner mini-batch gradient-descent step that updates sky temperatures from subsets of the visibility data without ever forming the full response matrix.","core_discovery":"AMIGO reconstructs all-sky maps from the dynamic, non-coplanar baselines of a lunar-orbit array by solving a constrained least-squares problem with mini-batch gradient descent and an augmented Lagrangian multiplier; memory falls from O(N^{2}) to O(N) and time falls from O(N^{3}+MN^{2}) to O(<MN^{2}), while adjustable residual thresholds let the user dial the strength of non-negativity and angular-power-spectrum priors.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["AMIGO reconstructs all-sky lunar maps with mini-batch gradient descent","Lunar-orbit arrays image whole sky as AMIGO drops memory to linear cost","AMIGO solves dynamic baselines for DSL-style all-sky radio maps","Constrained optimization via ALM lets lunar arrays map below 30 MHz","AMIGO turns non-coplanar lunar visibilities into full-sky images"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"That idealized mock visibilities—perfect orbits, opaque Moon, ideal beams, and pure thermal noise with no calibration or baseline errors—are realistic enough to prove the method will work on actual flight data.","fun_headline_variants_meta":{"raw":{"variants":["AMIGO reconstructs all-sky lunar maps with mini-batch gradient descent","Lunar-orbit arrays image whole sky as AMIGO drops memory to linear cost","AMIGO solves dynamic baselines for DSL-style all-sky radio maps","Constrained optimization via ALM lets lunar arrays map below 30 MHz","AMIGO turns non-coplanar lunar visibilities into full-sky images"]},"model":"grok-4.5","effort":"low","cost_usd":0.003836,"raw_usage":{"total_tokens":1256,"prompt_tokens":829,"num_sources_used":0,"completion_tokens":102,"cost_in_usd_ticks":38360000,"prompt_tokens_details":{"text_tokens":829,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":325,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":829,"tokens_out":102,"duration_ms":3355,"temperature":1.0,"reasoning_tokens":325,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T07:41:41.913904+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Run AMIGO on the same mock data after injecting realistic baseline-determination phase errors or primary-beam model mismatches; if the residual maps and MSE degrade beyond the levels reported for pure thermal noise, the claim of practical feasibility fails.","supporting_citations":[],"review_version":1}