{"id":"4300e8fa-eaee-4953-96f2-8a94f91aef84","arxiv_id":"2502.08890","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"νSpaceSim is an end-to-end Monte Carlo package that simulates tau-neutrino earth-skimming events, from underground interaction through lepton propagation and decay to optical and radio signals at a user-defined detector.","lead":"The paper describes νSpaceSim, an open-source simulation package that models how cosmic neutrinos passing through the Earth can create detectable upward air showers. It is a design and sensitivity tool for next-generation balloon and space neutrino observatories.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The released version's EAS module uses the Greisen parametrization rather than a full composite-shower model, so the quantitative sensitivity curves in Figures 7 and 8 may carry an unquantified systematic bias that the paper does not bound.","rationale":"The reader's weakest-assumption statement already identifies the chained approximations, including the Greisen parametrization, as not validated end-to-end. My pass converges on the same point but isolates the Greisen-based EAS module as the single most load-bearing element. This is because the other major links in the chain are externally benchmarked codes: NuPyProp/NuTauSim are published and validated packages, ZHAireS is a well-established radio Monte Carlo, and MERRA-2 is a standard atmospheric reanalysis. The Greisen parametrization, by contrast, is a simple analytic fit that the paper itself is in the process of replacing with CONEX-based composite showers, and the paper's own Figure 6 shows that the replacement accounts for mean and variance in a way that the current default does not. Since the headline quantitative outputs (Figures 7 and 8) are sensitivity curves, and since trigger rates depend on the detailed shower profile, the current public release's sensitivity numbers are not yet demonstrated to be reliable. The proposed check is feasible because the code is open source and the CONEX-based module already exists in beta, so the two configurations can be compared directly. This does not change the reader's conditional verdict: the paper is honest about the in-progress status, the package is real and publicly available, and the concern is addressable by validation rather than fatal. I therefore recommend keeping the verdict unchanged rather than moving it to accept or reject.","tokens_in":7383,"tokens_out":4231,"duration_ms":48418,"concrete_test":"Run the identical detector configuration used for Figure 7 twice with νSpaceSim: once with the current default Greisen EAS module and once with the CONEX-parametrized composite-shower module shown in beta in Figure 6, keeping all other inputs (Earth-propagation lookup tables, atmosphere, clouds, detector geometry) fixed. Compare the resulting 90% CL sensitivity curves and the differential event-rate distributions. If the difference in the sensitivity limit at the peak energy exceeds roughly 20%, the current public release does not yet support the quantitative sensitivity claims, and the paper should describe Figures 7 and 8 as preliminary pending validation of the EAS module. If the difference is below roughly 20%, the Greisen approximation is adequate for the claimed sensitivity-estimation purpose and the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing approximation in the released version is the Greisen parametrization for EAS longitudinal development, flagged in Section 3 and the caption to Figure 6: the paper states that this parametrization is 'currently employed in νSpaceSim', while the CONEX-parametrized composite-shower implementation is 'currently in beta' and 'future additions will include ... signals from the muon component of regular EAS'. The Greisen function provides only a mean electromagnetic cascade profile; it does not capture the composite electromagnetic-plus-hadronic character of tau-decay-initiated showers, the event-by-event fluctuations in Xmax and Nmax, or the muon content. The optical Cherenkov and radio signals that drive triggering depend directly on Xmax, Nmax, and the full longitudinal profile, so any mis-modeling of the shower development propagates directly into the predicted event rates and the 90% CL sensitivity curves in Figures 7 and 8. The paper does not validate the end-to-end chain against a real detected upward shower or against an independent full Monte Carlo, so there is currently no quantitative bound on this systematic error. This is not an allegation of fraud; it is a request to quantify the effect of the default EAS approximation before the quoted sensitivities are used to guide experiment design.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents νSpaceSim, an open-source Python package that simulates the full chain from cosmic neutrinos interacting in the Earth, through tau-lepton production and propagation, tau decay and extensive air shower development, to optical Cherenkov and radio emission, atmospheric propagation, and detector response. The authors describe the modular architecture, the use of precomputed lookup tables from external packages (NuPyProp, NuTauSim/NuLeptonSim, Pythia8, ZHAireS, MERRA-2, CONEX), and the package's speed and public availability. The manuscript reports example outputs, including a simulated sensitivity curve for an SPB-2-like balloon instrument and a comparison with the GW170817 sensitivities of ANTARES, IceCube, and Auger. The stated goal is to provide a user-friendly, end-to-end tool for designing and estimating the sensitivity of balloon-borne, space-based, and ground-based detectors to Earth-emergent neutrino-induced showers.","tokens_in":7771,"tokens_out":4649,"duration_ms":44677,"significance":"If the package performs as claimed, it would be a valuable community resource for the UHECR and high-energy neutrino community. The authors explicitly credit the package's public availability through pip, GitHub, and HEASARC, and the modular design with independently published physics packages is a genuine strength. The paper also demonstrates a fast workflow and a direct interface with experimental frameworks such as EUSO-Offline. However, the central claim of being a 'comprehensive' end-to-end simulation that can be used to calculate reliable sensitivity estimates is currently not supported by quantitative validation. The manuscript's own text acknowledges that the default shower model is the Greisen parametrization, that the improved composite-shower implementation is in beta, and that several atmospheric and detector effects are treated with simplifying assumptions. These limitations directly affect the quoted sensitivity projections, so the significance of the numerical results in Figures 7 and 8 is contingent on added validation.","major_comments":[{"comment":"The default EAS module in the released version uses the Greisen parametrization, which the paper itself describes as approximate in comparison with the CONEX-parametrized composite-shower implementation currently in beta. The Greisen profile is a mean electromagnetic cascade function; it does not describe event-by-event fluctuations in Xmax and Nmax, the muon content, or the composite electromagnetic-plus-hadronic structure of tau-decay-initiated showers. Since the optical Cherenkov and radio signals that drive triggering depend directly on the longitudinal profile, this approximation propagates into the sensitivity curves in Figures 7 and 8. The manuscript does not provide a quantitative bound on this systematic error, and the end-to-end chain is not benchmarked against an independent full Monte Carlo or a real detected event. Please add a validation study comparing the Greisen-based default with a CONEX/ZHAireS-based chain for representative energies and emergence angles, or explicitly re-label the sensitivity curves as preliminary projections that do not yet incorporate the improved shower model.","section":"Section 3 and Figure 6"},{"comment":"The comparison of the simulated SPB-2 sensitivity with the measured 14-day 90% confidence sensitivities of ANTARES, IceCube, and Auger to GW170817 lacks the assumptions needed for reproducibility. The manuscript does not specify the input neutrino flux model and spectrum, the source time profile used for GW170817, the detector trigger and background models, or the exact definition of the 14-day window. Without these details, the SPB-2 curve in Figure 8 cannot be reproduced or interpreted quantitatively. Please provide the full configuration, ideally as a TOML file in the public repository, and state the assumed transient source model.","section":"Section 4 and Figure 8"},{"comment":"The treatment of clouds as opaque to Cherenkov light is a strong assumption that can bias the optical sensitivity in either direction, and the manuscript gives no estimate of the size of this effect. Since the MERRA-2 database option is available, the authors should quantify the impact of the opaque-cloud assumption on the projected sensitivity in Figures 7 and 8, or at minimum state clearly whether the curves assume no clouds, a constant cloud layer, or MERRA-2-derived clouds. A sensitivity test varying cloud altitude and opacity would make the limitation concrete.","section":"Section 3, cloud and atmospheric modeling"}],"minor_comments":[{"comment":"The affiliation line contains a typo: 'Department of Phisics' should be 'Department of Physics'.","section":"Affiliation"},{"comment":"Reference [36] is listed with a garbled author string ('Antares, I., Auger, P., Scientific, L.I.G.O., ...') and should be corrected to the published collaboration author list; the reference style is also inconsistent with the rest of the bibliography.","section":"References"},{"comment":"The performance claim '10^6 events in ~5 minutes using a MacBook Pro with an M2Max processor' should specify the software version, the number of threads used, and whether this runtime includes the full end-to-end pipeline or only a subset of modules, so that the claim can be reproduced.","section":"Section 2"},{"comment":"The capitalization of package and collaboration names is inconsistent: 'nuPyProp' and 'NuPyProp' are both used, as are 'nuspacesim' and 'NuSpaceSim', and 'Euso Offline' appears in place of 'EUSO-Offline'.","section":"Throughout"},{"comment":"The phrase 'counts with inherent multi-core processing via Dask' is awkward; 'includes inherent multi-core processing via Dask' or 'uses Dask for multi-core processing' would be clearer.","section":"Section 2"},{"comment":"The title 'Nuspacesim Results Dashboard' in Figure 4 should be 'νSpaceSim Results Dashboard' for consistency with the paper's notation.","section":"Figure 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a proceedings contribution, and the package itself appears real and useful. My main concern is that the title and abstract claim a 'comprehensive' end-to-end simulation, while the manuscript's own text acknowledges that the default shower model is approximate and that key improvements are in beta. The burden is on the authors to either validate the current default against a full Monte Carlo or to present the sensitivity curves as clearly labeled preliminary results. This is fixable within the scope of the manuscript, so I recommend major revision rather than rejection. The reader's stress-test concern about the Greisen parametrization is substantiated by Section 3 and Figure 6, and it is the load-bearing issue for the numerical results."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know that νSpaceSim is a real, public, pip-installable code that ties together NuPyProp, Pythia8, ZHAireS, CONEX, and MERRA-2 into a fast end-to-end sensitivity estimator for upward tau air showers. The paper's best feature is its transparency: it clearly labels the composite EAS with CONEX as beta, says the Greisen parametrization is currently used, and lists future upgrades. That honesty is worth something.\n\nWhat's new is the integration itself: a modular Python wrapper, cloud opacity treatment, transient source pointing with a GW170817 comparison, and user-customizable detector geometry. The physics components come from prior published codes, but the public release and the ability to swap libraries (NuLeptonSim, TauRunner) are a real contribution to the field.\n\nThe main soft spot is the lack of end-to-end validation. There is no benchmark against a real detected event or an independent full Monte Carlo. The Greisen parametrization gives only a mean electromagnetic profile; it misses Xmax and Nmax fluctuations and the composite hadronic-muonic character that directly affect Cherenkov and radio triggers. So the 90% CL sensitivity curves in Figures 7 and 8 carry an unquantified systematic. The paper acknowledges the approximation but does not bound its effect. That is the load-bearing issue, and it is fixable: a short validation section or a clear caveat would change the verdict quickly.\n\nFigure 8 is also a bit overreaching. It compares a simulated SPB-2 sensitivity to measured limits from Auger, ANTARES, and IceCube for GW170817. Those are different objects: a simulation of a detector that was not there versus actual limits. The comparison is suggestive, not a direct apples-to-apples measurement, and the paper should say so more carefully.\n\nMinor points: the self-cited exit probability figures come from externally benchmarked codes, so that is fine. Some features are labeled beta, which is honest but means the paper is as much a status report as a physics result.\n\nWho is this for? Anyone working on EUSO-SPB2, POEMMA, PUEO, or ground-based detectors who needs a fast first-pass sensitivity tool. It is worth serious refereeing because the code is real, public, and already in use. The referee should ask for a sensitivity study of the default shower model, or at least a prominent statement that the quoted numbers are provisional until the CONEX upgrade is default.\n\nMy recommendation: send it to peer review, with the expectation of a revision that either validates the end-to-end chain or clearly quantifies the shower-model systematic.","headline":"A genuinely useful, openly released end-to-end simulation package for earth-skimming tau neutrinos; the paper is honest about its approximations, but the unvalidated Greisen shower model means the quoted sensitivity curves are indicative, not definitive.","tokens_in":8206,"tokens_out":1911,"would_cite":true,"duration_ms":20496,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"νSpaceSim claims one modular pipeline can model the whole neutrino-to-detector chain and produce sensitivity estimates for any optical or radio instrument.","keywords":["Earth-skimming neutrinos","tau neutrinos","upward extensive air showers","Cherenkov light","radio emission","detector sensitivity","balloon experiments","space-based detectors"],"falsifier":"Take a fixed tau-emergence geometry and energy, run νSpaceSim and an independent full air-shower Monte Carlo with the same atmosphere, and compare the predicted Cherenkov photon counts and radio pulse amplitudes at the detector; a disagreement beyond the quoted statistical uncertainties would show that the current sensitivity curves are not reliable.","tokens_in":7165,"feed_emoji":"🌍","tokens_out":10110,"duration_ms":98609,"temperature":0.7,"pith_summary":"νSpaceSim is an open-source simulation package whose goal is to model the complete chain by which an ultra-high-energy tau neutrino crosses the Earth, emerges as a tau lepton, decays in the atmosphere, and produces an upward-moving air shower whose optical Cherenkov and radio signals reach a user-defined detector. The paper's central claim is that this one package can calculate the sensitivity of any optical and/or radio instrument in a straightforward and user-friendly way, whether the detector sits on the ground, on a balloon, or in space. This matters because PeV-to-EeV neutrinos are so rare and weakly interacting that direct detection is impractical; treating the entire Earth as the neutrino target and the atmosphere as the signal generator gives an effective detector mass far beyond a gigaton, so accurate and fast sensitivity modeling is essential for designing the next experiments and interpreting data from recent ones.","feed_headline":"One simulation traces Earth-skimming neutrinos to detector signals","feed_subtitle":"νSpaceSim models tau-neutrino Earth propagation, air-shower light, and instrument response in a single user-configurable pipeline.","key_machinery":"The carrying mechanism is a sampled-library, modular pipeline: a vectorized Python wrapper that chains together pre-built lookup tables and external modules for the computationally expensive physical stages, then folds in a user-supplied detector configuration (altitude, effective area, frequency band, pointing) at the final stage. Lookup tables provide the probability that a tau exits the Earth and its emerging energy; a Monte-Carlo event generator produces the tau decay products; a radio-emission code computes the shower's radio signal; atmospheric reanalysis data supply cloud and weather conditions; and parameterized shower profiles are being integrated to replace the current Greisen parametrization. Because the wrapper is modular, a user can swap in alternative propagation tables or detector geometries through a TOML parameter file, and the same pipeline yields sensitivity curves for ground-based, balloon-borne, and space-based instruments.","core_discovery":"The paper claims that the relevant physics—neutrino charged- and neutral-current interactions inside the Earth, tau regeneration, tau decay, air-shower development, optical and radio emission, atmospheric attenuation, and instrument response—can be assembled into one modular, vectorized simulation, and that the result is fast enough to generate a million events in about five minutes on a laptop. The package's workflow is demonstrated with a simulated balloon mission at 33 km altitude pointing between the Earth's limb and 6.4 degrees below it, and with a target-of-opportunity comparison showing that the simulated 14-day sensitivity of a balloon detector to a binary-neutron-star merger would have been competitive with the sensitivities of the large observatories that actually observed it.","pith_inferences":["A natural next validation, not performed in the paper, would be to compare νSpaceSim's predicted Cherenkov and radio event rates with an independent full air-shower Monte Carlo for identical tau-emergence conditions; the paper's own note that the shower-shape parametrization is still being upgraded suggests this stage is the least settled.","The treatment of clouds as fully opaque to Cherenkov light means the optical sensitivity estimates are conservative in cloudy conditions; modeling partial transmission and scattering would probably increase predicted photon counts on partially cloudy lines of sight.","If the end-to-end chain is accurate, the same modular wrapper could be extended to the muon channel, which the paper mentions as possible but does not model in detail.","The modular design invites a systematic-error study that the paper does not report: swapping the built-in Earth-propagation lookup tables for an alternative public code would show how much the final sensitivity curves move, and that spread could be quoted as a model uncertainty."],"forward_implications":["A detector team can obtain a first sensitivity estimate by editing a parameter file, without writing custom simulation code for each physical stage.","Balloon- and space-mission designs can be optimized for neutrino-triggered upward showers while the package is still under active development.","A large ground-based fluorescence observatory can use νSpaceSim output to estimate its exposure to neutrino-induced upward extensive air showers.","Transient-source alerts can be paired with simulated 14-day sensitivities, giving a direct comparison of balloon and ground experiments for target-of-opportunity follow-up.","Users can exchange the Earth-propagation lookup tables for alternative public codes to quantify how much the resulting sensitivity depends on that modeling choice."],"supporting_citations":[{"why":"Supplies the lookup-table model of tau propagation and Earth-exit probability that feeds the simulated lepton chain.","marker":"[21]"},{"why":"Provides an alternative tau-propagation code used to cross-check the exit-probability stage.","marker":"[22]"},{"why":"Generates the decay products of the emerging tau that initiate the air shower.","marker":"[23]"},{"why":"Supplies atmospheric data used to model clouds, aerosols, ozone, and scattering along the signal path.","marker":"[24]"},{"why":"Computes the radio emission from the simulated extensive air showers.","marker":"[25]"},{"why":"Defines the CONEX-compatible output format that lets νSpaceSim showers be fed into external experiment frameworks.","marker":"[30]"},{"why":"Provides the planned universal Cherenkov-yield calculation that will replace the current Cherenkov modeling.","marker":"[35]"},{"why":"Gives the expected energy and angular distributions of Earth-emerging taus used to characterize the signal.","marker":"[20]"}],"fun_headline_variants":["νSpaceSim: one code from neutrino to air-shower signal","Earth-skimming neutrinos simulated end-to-end in minutes","Fast full-chain neutrino shower simulation for UHECR hunters","From Earth-bound tau to detector: νSpaceSim does it all","Million cosmic neutrino events on a laptop in five minutes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the chained approximations (lookup-table Earth-exit probabilities, the Greisen shower parametrization, the radio-emission model, and the treatment of clouds as opaque to Cherenkov light) each stay accurate enough that their errors do not accumulate; the paper does not validate the full chain against a real detected event.","fun_headline_variants_meta":{"raw":{"variants":["νSpaceSim: one code from neutrino to air-shower signal","Earth-skimming neutrinos simulated end-to-end in minutes","Fast full-chain neutrino shower simulation for UHECR hunters","From Earth-bound tau to detector: νSpaceSim does it all","Million cosmic neutrino events on a laptop in five minutes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000773,"raw_usage":{"total_tokens":3444,"prompt_tokens":989,"completion_tokens":2455,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":605,"completion_tokens_details":{"reasoning_tokens":2368}},"tokens_in":605,"tokens_out":2455,"duration_ms":17329,"temperature":1.0,"reasoning_tokens":2368,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T23:17:42.372686+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a fixed tau-emergence geometry and energy, run νSpaceSim and an independent full air-shower Monte Carlo with the same atmosphere, and compare the predicted Cherenkov photon counts and radio pulse amplitudes at the detector; a disagreement beyond the quoted statistical uncertainties would show that the current sensitivity curves are not reliable.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides an alternative tau-propagation code used to cross-check the exit-probability stage."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Generates the decay products of the emerging tau that initiate the air shower."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies atmospheric data used to model clouds, aerosols, ozone, and scattering along the signal path."},{"cited_title":"R., & Zas, E","cited_arxiv_id":null,"evidence_quote":"Computes the radio emission from the simulated extensive air showers."},{"cited_title":"J., & Bergman, D","cited_arxiv_id":null,"evidence_quote":"Provides the planned universal Cherenkov-yield calculation that will replace the current Cherenkov modeling."},{"cited_title":"H., Krizmanic, J","cited_arxiv_id":null,"evidence_quote":"Gives the expected energy and angular distributions of Earth-emerging taus used to characterize the signal."}],"review_version":1}