{"id":"b52e701a-4081-45ef-ac2d-ddcfb58b30eb","arxiv_id":"1908.08263","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"nOOn is a ROOT-based Monte Carlo generator that produces event-by-event forward neutrons from electromagnetic dissociation in ultra-peripheral heavy-ion collisions, using measured photoneutron cross sections and a logarithmic extrapolation to LHC energies.","lead":"This paper introduces nOOn, a Monte Carlo program that generates forward neutrons from electromagnetic dissociation in ultra-peripheral heavy-ion collisions. It lets experimenters add realistic neutron multiplicities and energies to existing vector meson generators such as STARlight, which currently do not produce them.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Validation gap: the neutron-multiplicity model is extrapolated from 140 MeV to LHC energies, and no end-to-end comparison of nOOn output to ALICE EMD data is shown, so the central physics claim is conditional.","rationale":"The reader's weakest assumption identifies the same load-bearing concern I would stress: the multiplicity model is extrapolated from photon energies below 140 MeV to the full LHC range, with no end-to-end comparison of nOOn output to measured EMD data. The paper's independent supports, namely the available code, the ROOT-based interface, and the worked examples with STARlight and the hot-spot model, are real but do not test whether the generated neutron multiplicities are physically accurate. I see no internal inconsistency that would overturn the conditional verdict, though Eq. (19) appears to have a combinatorial factor of 2 in the P_{1n}P_{2n} term that the compound-Poisson expansion would not produce; this should be checked in the code, but it is secondary to the missing validation of the high-energy extrapolation. If the proposed ALICE comparison reproduces the measured EMD rates, the concern would be resolved and the verdict could move to ACCEPT; if it fails, the generator's central physics claim would not be supported.","tokens_in":13441,"tokens_out":14541,"duration_ms":159322,"concrete_test":"Run nOOn in standalone mode for Pb-Pb at sqrt(s_NN) = 2.76 TeV and compare the predicted inclusive EMD cross section and the 1n, 2n, and Xn multiplicity fractions directly with the ALICE measurement [23], and also with RELDIS [24,25]. If the ALICE Xn cross section and the neutron-multiplicity fractions are reproduced within the quoted uncertainties, the Sec. 3.1 logarithmic extrapolation is validated for the LHC regime; if they are not, the central claim of reliable event-by-event neutron generation is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that nOOn reliably generates event-by-event forward neutrons for RHIC/LHC depends on the multiplicity model of Sec. 3.1. The average and dispersion of the neutron multiplicity are fitted to photoneutron data only up to 140 MeV (Fig. 7) and then extrapolated logarithmically to ~10^9 MeV; the branching-ratio map of Fig. 9 is built from that fit with a Gaussian approximation. At LHC energies the equivalent-photon spectrum extends to tens of GeV, so a non-negligible part of the EMD cross section, especially for multi-neutron final states and for Xn tagging, originates from photons above the fitted region. The paper does not validate the final generated events: it compares only the average/dispersion fit to RELDIS in Fig. 7, not the nOOn output to measured EMD cross sections or multiplicity distributions such as the ALICE measurement [23]. The generator could be implemented exactly as described and still emit neutrons with incorrect multiplicities if this extrapolation is biased, so the central claim is conditional on an unvalidated physics input.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents nOOn, a ROOT-based Monte Carlo generator that produces forward neutrons from electromagnetic dissociation (EMD) in ultra-peripheral collisions (UPC) of heavy ions. The generator computes nuclear break-up probabilities using measured photoneutron cross sections, a logarithmic extrapolation of the mean and dispersion of neutron multiplicity from 140 MeV to 1e9 MeV, a Gaussian approximation for the multiplicity shape, and a Poisson-based convolution of multiple excitations with truncation and renormalisation. It can be used either as an afterburner for the STARlight generator or with theoretical photonuclear cross sections as input, and it outputs final-state neutrons in a TTree format. The paper includes examples for coherent rho0 and J/psi production and discusses the implementation details and program flow.","tokens_in":13715,"tokens_out":3635,"duration_ms":39238,"significance":"If the generator is validated, it fills a practical gap: STARlight, the most commonly used UPC Monte Carlo, does not produce final-state neutrons, while nOOn does, and it is designed for straightforward integration with the simulation frameworks of RHIC and LHC experiments. The paper's strengths are that it builds on established measured data for the GDR region, uses evaluated nuclear data from ENDF, and provides an open-source, GPL-licensed program with reproducible event-by-event output. The central physics risk is the large extrapolation of neutron-multiplicity parameters from 140 MeV to LHC energies and the absence of a direct end-to-end comparison of the generator output to measured Pb-Pb EMD data, which leaves the quantitative reliability of the generated neutron multiplicities at LHC conditions conditional.","major_comments":[{"comment":"The mean and dispersion of the neutron multiplicity are fitted to data only up to 140 MeV and then extrapolated logarithmically to 1e9 MeV, while at LHC energies a non-negligible fraction of the EMD cross section comes from photons above the fitted region. The only validation shown is the comparison to the RELDIS model in Fig. 7, not a comparison of nOOn output to a measured Pb-Pb observable. The paper should add a quantitative validation against the ALICE EMD cross-section measurement (Ref. [23]) or, at minimum, provide an estimate of the systematic uncertainty propagated from the extrapolation to the generated multiplicities.","section":"Section 3.1, Fig. 7"},{"comment":"The Poisson assumption for the number of Coulomb excitations (Eq. (12)) is introduced without justification, and the renormalisation procedure shown in Fig. 5 indicates that about 20% of the break-up probability is lost at small impact parameters before renormalisation. The claim that increasing the number of excitations to six would not significantly change predictions is not demonstrated quantitatively. Since the renormalisation reshapes the multiplicity distribution and could affect the impact-parameter dependence used in Eq. (10), the paper should justify the Poisson ansatz and quantify the effect of the truncation and renormalisation on, for example, the 0n0n and XnXn fractions.","section":"Section 2.3, Eqs. (12) and (20), Fig. 5"},{"comment":"The Gaussian shape of the neutron multiplicity distribution is inferred from deconvoluted data at only two photon energies, 199 and 390 MeV, as shown in Fig. 8. This Gaussian shape is then applied over the entire extrapolated range up to 1e9 MeV in constructing the branching-ratio map of Fig. 9. The paper should discuss whether the Gaussian approximation remains reasonable at higher photon energies where additional reaction channels open, and ideally compare the branching-ratio map with the predictions of a model such as RELDIS over a wider energy range.","section":"Section 3.1, Fig. 8"}],"minor_comments":[{"comment":"The fit function for the logarithmic extrapolation of the mean and dispersion is described only verbally; the explicit functional form and the fit parameters should be stated in the text or in an appendix for reproducibility.","section":"Section 3.1"},{"comment":"The program is currently restricted to Pb only, yet the abstract and introduction refer to RHIC as a target application. The restriction should be stated more prominently in the abstract, or an indication of planned support for other nuclei (e.g., Au) should be given.","section":"Program Summary and Abstract"},{"comment":"For photon energies above 140 MeV, the last bin of the ENDF emission spectrum is reused; the potential effect of this approximation on the generated neutron energies and on the response of zero-degree calorimeters should be commented on.","section":"Section 4, 'Particle generation'"},{"comment":"The numerical labels in the heatmap of Fig. 13 are difficult to read at the published size; increasing the font size or using a different colour scale would improve readability.","section":"Fig. 13"},{"comment":"The photon flux expression contains a term proportional to 1/gamma^2 K0^2; the text says K1 is a Bessel function but does not state the order of K0, which is clear from context but should be specified for completeness.","section":"Section 2, Eq. (6)"}],"recommendation":"major_revision","confidential_remarks":"This is a useful code paper for CPC, and the authors have made a genuine effort to base the generator on measured and evaluated data. The main concern is that the central physics input—the high-energy extrapolation of the neutron multiplicity—is not validated against any measured Pb-Pb observable, while such data exist (e.g., ALICE Ref. [23]). The authors should be asked to provide a direct comparison of nOOn predictions to measured EMD cross sections or neutron multiplicity distributions, or to explicitly bound the resulting systematic uncertainty. This is a load-bearing point for the claims made in the conclusions, so it should not be left as a purely cosmetic improvement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper to know about: nOOn is a Monte Carlo generator that finally puts final-state forward neutrons into UPC events. STARlight gives you neutron-tagged cross sections but no neutron particles; this code takes either STARlight events or theory rapidity distributions, samples neutron multiplicities from measured photoneutron cross sections, and writes boosted neutrons into a TTree. That is the genuinely new part, and it is useful.\n\nWhat it does well: the input data are mostly measured or evaluated nuclear data (photoneutron cross sections, ENDF spectra), the Poisson combinatorics for multiple excitations are standard and clearly presented, and the renormalisation of the truncated sum is honest. The code is public and GPL, and the interface to ROOT/STARlight is sensible. The example applications with rho0 and J/psi show how the tool can separate 0n0n, 0nXn, XnXn classes.\n\nThe soft spot is exactly where the reader puts it. The mean and dispersion of neutron multiplicity are fitted to data below 140 MeV and extrapolated logarithmically to 1e9 MeV, with a Gaussian shape for the branching ratios. Since the photon energy of the accompanying vector meson can be several GeV at the LHC, a non-negligible part of the tagged cross section comes from this extrapolation. The authors compare the fit to RELDIS in Fig. 7, but they do not compare the final generated multiplicities or EMD cross sections to a measured Pb-Pb data set such as ALICE's 2.76 TeV measurement. That is a real validation gap, and it makes the absolute normalisation of high-multiplicity classes less certain than the paper implies. It is not a fatal flaw: the authors state the approximation openly, and for GDR-dominated or low-energy applications the tool is on solid ground. But a benchmark against measured forward-neutron data would materially increase confidence.\n\nMinor: the GitHub repository is not pinned to a commit in the paper, which makes reproducibility harder for a code paper.\n\nWho is this for: experimental UPC people at RHIC and the LHC, especially for trigger studies and for neutron-tagging based photon-energy separation. It is not a new physics result.\n\nI would send it to peer review. A serious referee can check the code and ask for the missing validation; the tool fills a real gap and deserves to be in the literature. My verdict would be conditional accept after a comparison to existing EMD measurements and a pinned repository version.","headline":"A useful UPC neutron generator that fills a real gap, but with an unvalidated high-energy extrapolation and no end-to-end benchmark against measured EMD data; worth peer review after that gap is addressed.","tokens_in":14179,"tokens_out":2881,"would_cite":true,"duration_ms":31475,"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":"The nOOn Monte Carlo program generates, event by event, the forward neutrons produced by electromagnetic dissociation in ultra-peripheral heavy-ion collisions.","keywords":["ultra-peripheral collisions","forward neutrons","electromagnetic dissociation","Monte Carlo generator","photonuclear cross sections","vector meson photoproduction","nuclear break-up","heavy-ion collisions"],"falsifier":"Compare the generator's predicted fractions of events with 0, 1, or many forward neutrons against measured forward-neutron multiplicities in Pb-Pb ultra-peripheral collisions at LHC energies; a disagreement that grows with photon energy would falsify the logarithmic extrapolation from 140 MeV to $10^9$ MeV.","tokens_in":13238,"feed_emoji":"⚛️","tokens_out":9156,"duration_ms":85473,"temperature":0.7,"pith_summary":"The paper introduces a Monte Carlo program called nOOn (noon) for ultra-peripheral collisions, in which two heavy nuclei pass without touching and interact through photon exchange. It generates, event by event, the forward neutrons that are produced when such an interaction excites one or both nuclei and they break up. The motivation is that experiments use these neutrons as triggers and as tags to separate the two photon-energy solutions in vector-meson photoproduction, so simulations need them included. nOOn builds neutron emission from measured photonuclear cross sections and a few modelling assumptions, and can be attached to existing vector-meson event generators or to theoretical photonuclear cross-section predictions.","feed_headline":"nOOn generates the forward neutrons that UPC experiments trigger on","feed_subtitle":"Neutron-tagged vector-meson events at RHIC and the LHC become simulatable for the first time.","key_machinery":"The central object is the nOOn program itself, a C++ Monte Carlo generator for events with forward neutrons. Its load-bearing mechanism is the combination of three impact-parameter-dependent probabilities: the hard photoproduction probability $P_P(b)$, the nuclear break-up probability $P_{ij}(b)$, factorised as $P_i(b) P_j(b)$ and constructed from Poisson-distributed, independent Coulomb excitations, and the no-hadronic-interaction factor $\\exp(-P_H(b))$; these are weighted by the semi-classical photon flux to assign neutron multiplicities for a given photon energy. The energy behaviour is carried by a branching-ratio map obtained by extrapolating the measured mean and dispersion of neutron multiplicity logarithmically up to photon energies of $10^9$ MeV, supplemented by a Gaussian approximation for the multiplicity shape and evaluated nuclear data tables for neutron energies.","core_discovery":"nOOn claims to supply complete per-event simulation of the forward neutrons from Coulomb break-up in ultra-peripheral heavy-ion collisions. It starts from measured total and partial photonuclear cross sections for 208Pb, builds impact-parameter-dependent probabilities for each nucleus to emit any number of neutrons, treats multiple photon exchanges as independent Poisson excitations, and factorises the two sides. The photon-energy dependence is handled by a branching-ratio map built from a logarithmic extrapolation of the measured mean and dispersion of neutron multiplicity, with a Gaussian shape for the multiplicity distribution; neutron energies are sampled from evaluated nuclear data and then boosted to the laboratory frame. The output is a standard list of particles per event, ready for detector simulation.","pith_inferences":["Because the neutron tag effectively selects the photon energy, this generator could also be used to correct the acceptance of forward J/psi measurements, not just to build triggers; the paper does not develop that use.","The Poisson and independence assumptions could be tested directly against measured forward-neutron multiplicity distributions from LHC Pb-Pb runs, a comparison the paper does not include.","The Gaussian shape used above 140 MeV is an interpolation; data at intermediate photon energies would either validate it or force a more detailed shape model."],"forward_implications":["Simulated ultra-peripheral collision events can now include realistic forward-neutron signals, allowing trigger and acceptance studies for neutron-based selectors.","Vector-meson measurements can be separated into 0n0n, 0nXn, and XnXn neutron-tag classes, which helps pin down the photon energy when rapidity alone is ambiguous.","Theoretical photonuclear cross sections can be expanded into rapidity-dependent, neutron-tagged cross-section predictions, as the paper demonstrates for rho0 and J/psi.","The same program structure applies to other photon-induced processes in ultra-peripheral collisions, such as jet production and light-by-light scattering."],"supporting_citations":[{"why":"Supplies the photoneutron cross-section data for 208Pb that anchor the total break-up probability at GDR energies.","marker":"[30]"},{"why":"Provides the partial cross sections for neutron multiplicities up to 10 and photon energies up to 140 MeV used for the low-energy region.","marker":"[31]"},{"why":"Gives the total photonuclear absorption cross section in the delta-resonance region and the multiplicity-shape data at 199 and 390 MeV used to justify the Gaussian approximation.","marker":"[32]"},{"why":"Provides the Regge parametrisation used to extend the photonuclear cross section to the highest LHC photon energies.","marker":"[39, 40]"},{"why":"Reports the measured mean and dispersion of neutron multiplicity that are logarithmically extrapolated to build the branching-ratio map.","marker":"[41]"},{"why":"Serves as the dissociation model against which the extrapolated neutron multiplicity is compared for consistency.","marker":"[25]"},{"why":"Supplies the factorised cross-section formula with nuclear break-up that the generator's probability construction implements.","marker":"[26]"},{"why":"Provides the external vector-meson event sample that nOOn reads as input in the example application.","marker":"[29]"},{"why":"Supplies the evaluated photo-neutron emission spectra used to sample neutron energies in the nucleus rest frame.","marker":"[43]"}],"fun_headline_variants":["nOOn: first Monte Carlo for UPC forward neutrons","Simulate neutron-tagged UPC events with nOOn","UPC forward neutrons enter simulation with nOOn","nOOn interfaces with vector-meson generators for UPC","Generate forward neutrons in UPC with nOOn"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The generator assumes that the average and spread of neutron multiplicity, measured only up to photon energies of 140 MeV, continue to grow logarithmically all the way to about $10^9$ MeV; if that extrapolation is wrong, the produced neutron multiplicities at LHC energies will be biased.","fun_headline_variants_meta":{"raw":{"variants":["nOOn: first Monte Carlo for UPC forward neutrons","Simulate neutron-tagged UPC events with nOOn","UPC forward neutrons enter simulation with nOOn","nOOn interfaces with vector-meson generators for UPC","Generate forward neutrons in UPC with nOOn"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000263,"raw_usage":{"total_tokens":1558,"prompt_tokens":861,"completion_tokens":697,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":477,"completion_tokens_details":{"reasoning_tokens":619}},"tokens_in":477,"tokens_out":697,"duration_ms":6626,"temperature":1.0,"reasoning_tokens":619,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:44:18.256344+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the generator's predicted fractions of events with 0, 1, or many forward neutrons against measured forward-neutron multiplicities in Pb-Pb ultra-peripheral collisions at LHC energies; a disagreement that grows with photon energy would falsify the logarithmic extrapolation from 140 MeV to $10^9$ MeV.","supporting_citations":[{"cited_title":"Veyssiere, H","cited_arxiv_id":null,"evidence_quote":"Supplies the photoneutron cross-section data for 208Pb that anchor the total break-up probability at GDR energies."},{"cited_title":"Lepretre, H","cited_arxiv_id":null,"evidence_quote":"Provides the partial cross sections for neutron multiplicities up to 10 and photon energies up to 140 MeV used for the low-energy region."},{"cited_title":"Carlos, H","cited_arxiv_id":null,"evidence_quote":"Gives the total photonuclear absorption cross section in the delta-resonance region and the multiplicity-shape data at 199 and 390 MeV used to justify the Gaussian approximation."},{"cited_title":"Lepretre, H","cited_arxiv_id":null,"evidence_quote":"Reports the measured mean and dispersion of neutron multiplicity that are logarithmically extrapolated to build the branching-ratio map."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the evaluated photo-neutron emission spectra used to sample neutron energies in the nucleus rest frame."}],"review_version":1}