{"id":"be698ec0-0e64-470d-a3c5-a951d8d6587f","arxiv_id":"2509.07822","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Ionization cross sections of biomolecules hit by protons collapse onto a near-universal curve when divided by a weighted atom count, and a simple per-family fitted formula reproduces it.","lead":"This paper reviews and extends a computational shortcut, the IAM-PCM, for predicting how protons knock electrons off molecules, including large biomolecules. It reports that ionization cross sections for many biomolecules collapse onto one curve when scaled by a molecule-specific number, and that this curve can be captured by a simple fitted formula with two parameters per molecule family.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Scaling claim rests on single-scattering assumption that authors themselves flag as possibly overcorrecting; no test yet distinguishes model artifact from physical prediction.","rationale":"The paper is internally consistent as a model review: the weighted-sum definitions, IEM multinomial probabilities, and orientation averaging in Sec. II check out, and the scaling collapse is convincingly demonstrated on the model's own output, with error bands under 10% as quoted. The parametrization Eq. (25) is an economical representation of those output curves. However, the model's foundational premise—single scattering—is exactly the load-bearing assumption for every cross section and scaling curve presented. The authors themselves flag in Sec. IV that IAM-PCM 'appears to overcorrect additivity rule predictions' and that the single-scattering condition may need relaxation. If that condition fails, the central claim that net ionization cross sections follow from atomic composition alone would be an artifact of the model rather than a physical prediction. The proposed computational test—allowing a second scattering event and checking whether the reduced curves shift beyond the ~10% band—would directly settle this. The reader's CONDITIONAL verdict is appropriate; this concern does not change it, but it sharpens the condition: the scaling's physical validity must be tested against the single-scattering assumption, not only against the model's own output.","tokens_in":14345,"tokens_out":10654,"duration_ms":119723,"concrete_test":"Modify the IAM-PCM to allow a second scattering event along the projectile path (e.g., via a mean-free-path criterion or a two-event weight model) and recompute the orientation-averaged net ionization cross sections for a representative pyrimidine, purine, amino acid, and nucleotide over 10 keV–10 MeV. If the reduced cross sections shift by more than the quoted ~10% band or no longer collapse onto a common curve, the scaling/parametrization in Sec. IIIC is contingent on the single-scattering condition, not a robust molecular property. Alternatively, compare the single-scattering IAM-PCM net ionization for uracil with a full TD-DFT/molecular calculation at 100 keV; a >20% deviation would indicate the assumption matters.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central scaling claim (Sec. IIIC, Eqs. (20)-(25)) is built entirely on the IAM-PCM single-scattering condition: net cross sections are unions of atomic disks, formalized by weight factors s_x|j in Eqs. (2)-(6) that average over m pixels along a projectile path with weight 1/m. This is not a harmless bookkeeping detail; it is what turns the additivity-rule sum into the IAM-PCM union area. Every reduced ionization curve that collapses in Figs. 6-7 and every parameter in Table I inherits this assumption. The paper itself, in Sec. IV, states that IAM-PCM 'appears to overcorrect additivity rule predictions of net cross sections' and explicitly questions 'whether the single-scattering condition that is used as a model assumption needs to be relaxed.' Thus the authors acknowledge that the central mechanism may be quantitatively wrong for the very class of targets (complex biomolecules) to which the parametrization is applied. If the single-scattering condition is violated, the 'remarkable accuracy' of Eq. (25) is an accuracy for an incorrect model, and the claim that cross sections follow from atomic composition alone is not physically justified. The paper provides no independent test of this assumption; the only validation points are a private communication (Ref. [15]) and experiments with factor-of-two discrepancies (Fig. 4 right panel).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reviews the independent atom model with pixel counting (IAM-PCM) and presents a new conceptualization in which pixels in an impact-parameter area decomposition represent scattering events. The method is applied to proton collisions with ten-electron systems and a range of biomolecules (pyrimidines, purines, amino acids, nucleotides). The central new claim is that net ionization cross sections, divided by a composition-dependent parameter N_IAM (Eq. (20)), collapse onto a common reduced curve, and that this curve is accurately reproduced by a simple parametrized Bethe-type formula, Eq. (25), with per-family fitted parameters alpha and beta (Table I) and Bethe constants A=100, B=120 read from the authors' TC-BGM atomic cross sections. The paper also discusses charge-state correlated capture and ionization cross sections, comparing with experiments for H2O, CH4, NH3, and HF.","tokens_in":14676,"tokens_out":5892,"duration_ms":69148,"significance":"If fully established, the paper would provide an extremely economical, composition-only description of net ionization cross sections for biologically relevant molecules under proton impact over 10 keV to 10 MeV, which is directly relevant to hadron therapy and radiation dosimetry. The formal derivation of the IEM multiple-capture/ionization probabilities (Eqs. (7)-(11)) is standard and internally consistent, and the scaling collapse is demonstrated with quoted error bands in Figs. 6-7. The authors are also candid about model assumptions, explicitly flagging the single-scattering condition and the ad hoc renormalization of multiple capture in Sec. IV. However, the paper's headline claim of 'remarkable accuracy' for Eq. (25) is currently a statement about the quality of a fit to the authors' own calculations, not an independently validated prediction. The physical status of the scaling law rests on the contested single-scattering assumption, and the available experimental validation is mixed, with factor-of-two discrepancies in some ratios. The manuscript would be strengthened by a predictive test that separates the model's internal consistency from a robust physics claim.","major_comments":[{"comment":"The central claim that Eq. (25) captures the scaling 'with remarkable accuracy' is a statement about a fit to the authors' own IAM-PCM calculations. The parameters alpha and beta are fitted per molecular family (Table I), and A=100, B=120 are read off the TC-BGM atomic curves. The reduced cross sections in Figs. 6-7 are the very data used to determine the family parameters and N_IAM. No independent test is presented (e.g., a molecule held out from the fitting, or comparison with experimental net-ionization data for a biomolecule not already used to set parameters). Please either provide such a test or explicitly characterize Eq. (25) as an interpolation formula for the IAM-PCM results rather than a validated prediction.","section":"Sec. IIIC, Eqs. (20)-(25), Table I"},{"comment":"The scaling collapse and all molecular cross sections inherit the single-scattering condition formalized in Eqs. (2)-(6), where each of m pixels along a projectile path contributes with weight 1/m and multiple scattering is excluded. In Sec. IV the authors state that IAM-PCM 'appears to overcorrect additivity rule predictions of net cross sections' and question whether the single-scattering condition needs to be relaxed. Since the paper provides no diagnostic that distinguishes a physical scaling law from a property of this assumption, the central claim remains conditional. The only external checks cited are a private communication (Ref. [15]) and data with factor-of-two discrepancies (Fig. 4, right panel). A quantitative sensitivity study (e.g., comparing the 1/m weighting with a first-hit alternative, or estimating the magnitude of multiple scattering) would be needed to support the ph","section":"Sec. IIA and Sec. IV"},{"comment":"The coefficients in N_IAM (0.25 for H, 0.9 for F, 1.5 for P, 1 for C/N/O) are said to 'follow from the atomic cross sections,' but no derivation or fitting procedure is supplied. Together with the per-family alpha and beta, these coefficients form a five-parameter atomic weighting plus two-parameter family fit. The paper should state how the weights were obtained (e.g., from high-energy ratios of the atomic cross sections in Fig. 1) and ideally provide a table of those ratios. Without this, the 'single parameter' characterization is misleading and the parametrization is not independently reproducible.","section":"Sec. IIIB, Eq. (20)"}],"minor_comments":[{"comment":"The expression (100 ln E + 120)/E is dimensionally awkward if E is in keV. Please specify that E is to be taken in keV inside the logarithm, or introduce a dimensionless ratio E/E0 with E0=1 keV.","section":"Eq. (25)"},{"comment":"The reported 'average err' values (10%, 6%, 2%) are not defined. Please state how the average is calculated (e.g., over energy points, over molecules) and whether it is a maximum deviation or a root-mean-square band width.","section":"Figs. 6-7"},{"comment":"Several experimental data points appear without error bars. Adding error bars would clarify the claimed factor-of-two discrepancies and the 'complete coincidence count' argument for Werner et al.","section":"Fig. 4, right panel"},{"comment":"The ad hoc renormalization of multiple capture into single capture is clearly labeled, but it may help to note explicitly that this procedure preserves net capture and does not alter net ionization cross sections, so that readers understand which results are affected.","section":"Sec. IIB, Eqs. (17)-(18)"},{"comment":"There are a few typographical errors: 'endavours' (p. 2), 'reasonabe' (p. 7), 'purins' (Sec. IIIC) and 'multply-charged' (Sec. IV). Also, the phrase 'This book' in the Introduction indicates a book-chapter origin; the journal submission may need adjustments in style.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reads more like a book chapter than a standalone research article, with extensive review material and long figure descriptions. The central new result—the scaling parametrization—would be acceptable as a practical engineering model, but the abstract's 'remarkable accuracy' overstates its status. The major issue is that Eq. (25) is a fit to the same IAM-PCM data it claims to capture, and the key underlying assumption (single-scattering) is openly questioned by the authors. A major revision that includes a genuine predictive test (e.g., leave-one-molecule-out validation or a comparison with a not-yet-used experimental dataset) and a quantitative sensitivity analysis of the single-scattering condition would substantially increase confidence. I do not see this as a reject, because the model is coherent and the limitations are acknowledged; it is a matter of calibrating the claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a book-chapter-style review of the authors' own IAM-PCM program. Most of it summarizes their previous papers; the genuinely new items are the HF target results, the per-family alpha/beta fit in Table I, and the four-family reduced-curve compilation. If you need a fast composition-only estimate of ionization yields for biomolecules in the hadron-therapy range, this is a handy source.\n\nThe paper does two things well. It lays out the formal machinery (weighted sum, multinomial probabilities, orientation averaging) carefully, and it convincingly demonstrates that the reduced cross-section collapse is a robust property of the model's own output—error bands under 10% across pyrimidines, purines, amino acids, and nucleotides. The pixel-as-scattering-event reinterpretation is a clearer way to present the method than the earlier geometric language. The authors also deserve credit for flagging their own uncertainties: they openly call the capture renormalization 'ad hoc' and concede in Sec. IV that IAM-PCM may overcorrect the additivity rule and that the single-scattering condition might need to be relaxed.\n\nThe main soft spot is the predictive framing. Equation (25) is fit to the same data it supposedly reproduces: alpha and beta come from Table I, and A and B are read off the atomic curves. So the 'remarkable accuracy' in the abstract is really just interpolation error. That is fine as a parametrization, but it should not be sold as a prediction. Second, the scaling itself is built entirely on the single-scattering assumption. The paper's own admission that this assumption may be wrong means the universal curve could be an artifact of the model, not a physical law. There is no out-of-sample test to settle that. Third, the sigma_q/sigma_1 universality rests on experiments that disagree by factors of two and on the ad hoc capture renormalization. And one key validation point is a private communication.\n\nNone of this is fatal for a methods summary, but it means the strong claims should be treated as hypotheses, not established results. The paper deserves a serious referee, but the referee should require a reframing of Eq. (25) as interpolation, an out-of-sample test, and a citable version of the corrected OCTOPUS comparison.","headline":"A solid methods summary and compilation of the authors' own IAM-PCM results, with one genuinely new parametrization; but the 'prediction' in Eq. (25) is really an interpolation, and the scaling rests on an assumption the authors themselves question.","tokens_in":15189,"tokens_out":2310,"would_cite":false,"duration_ms":26150,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Proton-impact ionization of biomolecules collapses onto one atomic-composition curve.","keywords":["independent atom model","pixel counting method","ion-molecule collisions","proton impact","net ionization cross sections","biomolecules","additivity rule","charge-state correlated cross sections"],"falsifier":"Measure the net ionization cross section of a nucleotide not in the fit set, for example deoxyguanosine triphosphate (dGTP, C10H16N5O13P3, N_IAM = 36.5 by Eq. (20)), for 10 keV to 10 MeV protons and compare with Eq. (25) using the nucleotide alpha and beta from Table I. The paper's claim is that the reduced cross section sigma/N_IAM stays on the same narrow band as the six fitted nucleotides; a deviation outside the claimed ~2% would falsify the parametrization. A second decisive test is complete-coincidence proton-on-methane and proton-on-ammonia electron-production measurements at about 1 Me","tokens_in":2007,"feed_emoji":"🧬","tokens_out":2595,"duration_ms":102012,"temperature":0.7,"pith_summary":"This paper argues that the independent atom model with pixel counting (IAM-PCM), which builds molecular collision cross sections from atomic ones by treating each atom's cross section as an overlapping pixelated area, can be re-read as a scattering-event picture: every pixel is a potential event, and along a projectile path crossing m pixels only one event counts, with weight 1/m. On this basis the authors compute net and charge-state-resolved capture and ionization cross sections for 10 keV to 10 MeV protons on systems ranging from neon to nucleotides. Their central quantitative claim is that net ionization cross sections for pyrimidines, purines, amino acids, and nucleotides, divided by a single composition number N_IAM, collapse onto a common reduced curve, and that a two-term effective-charge formula reproduces each family's curve. If this holds, ionization cross sections for any molecule in these families would follow from its atomic composition alone, without solving molecular dynamics. That matters for radiation-dose modelling and for estimating cross sections of unmeasured biomolecules.","feed_headline":"Atomic count fixes proton-ionization cross sections","feed_subtitle":"Net electron emission of four biomolecule families collapses onto one curve set by composition alone.","key_machinery":"The central object is the IAM-PCM weight factor s_{x|j}, the fractional contribution of atom j to the union of overlapping atomic cross-section disks, computed by pixel counting; it converts the simple additivity-rule sum into a molecular cross section that depends on molecular orientation and projectile energy. The scaling claim rides on N_IAM, a linear composition number derived from the relative sizes of atomic net ionization cross sections, and on the effective projectile charge Q_eff(E) = 1 - alpha exp(-beta E^0.4), which bends the high-energy Bethe form down to capture-dominated energies. Together these define Eq. (25), the parametrization that is the paper's headline result.","core_discovery":"The paper's central claim is that net electron-production (ionization) cross sections for proton collisions with biomolecules obey a universal reduced scaling: when the molecular cross section is divided by N_IAM = 0.25 n_H + n_C + n_N + n_O + 0.9 n_F + 1.5 n_P, all molecules within a family fall on one curve at high energies, and the modified Bethe form sigma_ion^mod(E) = N_IAM [1 - alpha exp(-beta E^0.4)]^2 (100 ln E + 120)/E reproduces the per-family curves from 10 keV to 10 MeV. The authors also argue that the IAM-PCM can be conceptualized equivalently by associating each pixel in the pixelized atomic cross-section disks with a scattering event, taking only one event per projectile path","pith_inferences":["If the N_IAM scaling is real, it likely extends to other planar bio-relevant heterocycles such as imidazoles or thiophenes because their valence structure is similar; computing the reduced curve for one such molecule would test this extension cheaply.","The formula can be inverted: measured cross sections for a few molecules in a family could fix alpha and beta independently of any particular atomic-collision code, making Eq. (25) an empirical predictive tool.","Methane's failure to join the reduced biomolecular curve suggests the scaling applies to molecules with a conjugated or at least multi-heavy-atom core; finding the smallest molecule that still collapses would map the boundary of the effect.","The pixel-scattering reinterpretation ties the weight factor to a mean number of atoms encountered per projectile path; that number could be compared directly with molecular geometry to check the single-scattering condition without new scattering data."],"forward_implications":["Net ionization cross sections for any molecule in the four families—pyrimidines, purines, amino acids, and nucleotides—can be estimated directly from its chemical formula and the family's two fitted parameters, with no molecular scattering calculation.","Composition rather than electronic structure orders electron production: a molecule with larger N_IAM has proportionally larger ionization cross section, with the reported spread kept to under 10% (pyrimidines), 6% (purines, amino acids), and 2% (nucleotides).","The ratios of q-fold to single electron production, sigma_q/sigma_1, are predicted to be nearly universal for the ten-electron molecules HF, H2O, NH3, CH4 above about 100 keV, independent of N_IAM; this is directly testable with complete coincidence experiments.","Relaxing the single-scattering assumption—the paper's stated next step—would change the weight factors and could remove the apparent overcorrection seen for multiply-charged projectiles.","For radiation-damage and hadron-therapy modelling, the parametrization turns a library of molecular targets into a closed-form input from atomic composition alone."],"supporting_citations":[{"why":"Introduces the IAM-PCM geometric picture of molecular cross sections as overlapping atomic areas, the foundation of the weight-factor sum.","marker":"[6]"},{"why":"Lays out the pixelization technique and the full IAM-PCM theory, establishing the operational machinery used throughout.","marker":"[7]"},{"why":"Applies IAM-PCM to biomolecules and derives the scaling of net cross sections with molecule size; the direct predecessor of the present N_IAM collapse.","marker":"[8]"},{"why":"Develops the charge-state-correlated independent-electron-model analysis and the renormalization of multiple capture into single capture.","marker":"[13]"},{"why":"Provides converged atomic cross sections for hydrogen used to benchmark the underlying atomic collision data.","marker":"[22]"},{"why":"Supplies the Bethe-Born form and generalized oscillator-strength treatment that Eq. (25) modifies with the alpha-beta suppression factor.","marker":"[25]"},{"why":"Provides complete-coincidence experimental q-fold electron-removal data for proton-water collisions used to test the predicted universal ratios.","marker":"[26]"},{"why":"Provides proton-water ratio data that differ from [26]; the discrepancy motivates the call for new complete-coincidence measurements.","marker":"[27]"},{"why":"An independent implementation of IAM-PCM with time-dependent density functional theory atomic cross sections; its corrected results match the present net capture curves.","marker":"[14]"}],"fun_headline_variants":["Biomolecule ionization scales with atomic makeup","One formula predicts proton ionization cross sections","Weighted atom count collapses ionization curves","Pixel counting reveals universal ionization scaling","Atomic weight factors unify proton impact ionization"],"cache_read_input_tokens":16768,"weakest_assumption_plain":"The load-bearing premise is the single-scattering condition—that a proton whose path crosses m atomic pixels counts exactly one scattering event with weight 1/m—and the paper itself flags this condition, together with the ad hoc renormalization of multiple capture into single capture, as the spot to relax.","fun_headline_variants_meta":{"raw":{"variants":["Biomolecule ionization scales with atomic makeup","One formula predicts proton ionization cross sections","Weighted atom count collapses ionization curves","Pixel counting reveals universal ionization scaling","Atomic weight factors unify proton impact ionization"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000533,"raw_usage":{"total_tokens":2403,"prompt_tokens":750,"completion_tokens":1653,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":494,"completion_tokens_details":{"reasoning_tokens":1591}},"tokens_in":494,"tokens_out":1653,"duration_ms":15295,"temperature":1.0,"reasoning_tokens":1591,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T21:40:50.187074+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the net ionization cross section of a nucleotide not in the fit set, for example deoxyguanosine triphosphate (dGTP, C10H16N5O13P3, N_IAM = 36.5 by Eq. (20)), for 10 keV to 10 MeV protons and compare with Eq. (25) using the nucleotide alpha and beta from Table I. The paper's claim is that the reduced cross section sigma/N_IAM stays on the same narrow band as the six fitted nucleotides; a deviation outside the claimed ~2% would falsify the parametrization. A second decisive test is complete-coincidence proton-on-methane and proton-on-ammonia electron-production measurements at about 1 Me","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the IAM-PCM geometric picture of molecular cross sections as overlapping atomic areas, the foundation of the weight-factor sum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Lays out the pixelization technique and the full IAM-PCM theory, establishing the operational machinery used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Applies IAM-PCM to biomolecules and derives the scaling of net cross sections with molecule size; the direct predecessor of the present N_IAM collapse."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Develops the charge-state-correlated independent-electron-model analysis and the renormalization of multiple capture into single capture."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides converged atomic cross sections for hydrogen used to benchmark the underlying atomic collision data."},{"cited_title":"Fedus and G","cited_arxiv_id":null,"evidence_quote":"Supplies the Bethe-Born form and generalized oscillator-strength treatment that Eq. (25) modifies with the alpha-beta suppression factor."},{"cited_title":"Werner, K","cited_arxiv_id":null,"evidence_quote":"Provides complete-coincidence experimental q-fold electron-removal data for proton-water collisions used to test the predicted universal ratios."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides proton-water ratio data that differ from [26]; the discrepancy motivates the call for new complete-coincidence measurements."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"An independent implementation of IAM-PCM with time-dependent density functional theory atomic cross sections; its corrected results match the present net capture curves."}],"review_version":1}