REVIEW 4 major objections 6 minor 13 references
Multi-IMPT: a biologically equivalent approach to proton ARC therapy
T0 review · 4 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read A scheme that rotates through six four-field IMPT plans can match proton arc therapy's plan quality using biologically effective dose as the comparison metric.
desk verdict Useful planning idea, but the reported plans violate the paper's own BED constraints, so the ARC-equivalence claim is not yet demonstrated. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing identity is the linear-quadratic biologically effective dose to an OAR voxel, $BED_j^m = \sum_{t=1}^{T} d_{jt}^m + \rho_m \sum_{t=1}^{T} (d_{jt}^m)^2$ with $\rho_m = 1/(\alpha_m/\beta_m)$, evaluated with $\alpha/\beta = 2$ Gy for every OAR. Because BED is nonlinear in per-fraction dose, changing which four beams are used in a fraction changes the OAR cost, letting the optimizer spread dose unevenly across fractions to reproduce an arc plan's cumulative BED. The optimization minimizes physical target dose deviation subject to BED-max, BED-mean, and BED-DVH constraints for OAR, DVH-min and max-dose constraints for the target, and a minimum-monitor-unit constraint, solved by iterative convex relaxation and the alternating direction method of multipliers.
What would settle it
Recompute the four planning comparisons with OAR $\alpha/\beta$ set to, say, 3 Gy or 10 Gy, or measure normal-tissue damage in an animal model treated with both schedules, and check whether multi-IMPT's OAR BED stays within the paper's stated clinically insignificant margin of ARC; a reversal of the prostate bladder or rectum BED50 advantage or a widening lung difference would falsify the equivalence claim.
Extended reading notes
Core claim
The central claim is that multi-IMPT—delivering a different subset of four beam angles in each fraction and optimizing BED for OAR while targeting physical dose to the target—produces plan quality equivalent to spot-scanning proton arc therapy. In the prostate case, multi-IMPT lowered bladder BED50 from 44.60 Gy to 24.49 Gy and rectum BED50 from 45.02 Gy to 24.69 Gy compared with ARC, while target conformity was similar. Lung and head-and-neck cases were dosimetrically similar; the brain case was slightly worse with conformity index 0.863 versus 0.904 for ARC. The paper concludes that the dosimetric differences are not clinically significant and that multi-IMPT is an ARC-equivalent delivery scheme.
Load-bearing premise
The equivalence holds only if the linear-quadratic BED model with $\alpha/\beta = 2$ Gy for all organs at risk describes real normal-tissue response, and only if modeling proton ARC as the same plan repeated identically in every fraction is faithful; if either fails, the biological equivalence result may not carry to patients.
Editorial extensions
If this is right
- Clinics with fixed-gantry or limited-angle proton systems could deliver plans with ARC-like OAR sparing by rotating through a handful of fixed-beam IMPT plans across fractions.
- Standard IMPT delivery, quality assurance, and delivery workflows could be used, avoiding the energy-switching overhead and continuous gantry rotation of ARC delivery.
- Because the optimization separates by fraction, the multi-IMPT inverse problem is computationally cheaper than a single full-arc optimization.
- BED-based objectives give a principled way to evaluate nonuniform fractionation in which each fraction's dose distribution differs.
- For prostate, multi-IMPT may improve on ARC for OAR BED while keeping target coverage; for brain, ARC retains a modest advantage.
Reading between the lines
- The paper leaves implicit that its six fixed angle combinations are only one possible schedule; an optimizer that selects per-fraction angle subsets or plan order could push OAR BED even lower, and may be needed in anatomies not tested here.
- Because BED is nonlinear in per-fraction dose, the method suggests a general principle that spatiotemporal fractionation can buy normal-tissue sparing without arc hardware, a principle that could extend to photon therapy or to combining IMPT with ultra-high-dose-rate delivery.
- A testable extension is a robustness study: if OAR $\alpha/\beta$ varies by patient or endpoint, one could optimize multi-IMPT across a range of $\alpha/\beta$ values and compare with equivalent ARC plans rather than the fixed 2 Gy used here.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes multi-IMPT, a delivery scheme in which a different subset of four beam angles is used in each treatment fraction, and optimizes the biologically effective dose (BED) to organs at risk while maintaining physical dose to the target, with the goal of approximating proton arc therapy (ARC) plan quality. The authors formulate a non-convex BED-constrained optimization problem, solve it with iterative convex relaxation and ADMM, and present planning comparisons against an ARC model (equal dose per fraction, beams spaced at 15° over 360°) for prostate, lung, brain, and head-and-neck cases. They report similar or slightly better plan quality for multi-IMPT in most cases and conclude that multi-IMPT is biologically equivalent to ARC.
Significance. If substantiated, the approach would be clinically valuable: proton centers without arc delivery capabilities could approximate ARC dosimetric benefits using standard IMPT delivery. The paper has several strengths: it formulates a concrete, reproducible optimization model; the method is computationally tractable; and the planning comparisons cover four distinct anatomical sites. The central claim, however, rests on the reported plans satisfying the stated BED constraints and on the ARC comparator faithfully representing clinical proton arc therapy. The manuscript also makes a falsifiable prediction that multi-IMPT can match ARC's OAR BED, which is a useful contribution. The main weaknesses are that several reported plans violate the paper's own constraints, the ARC model is simplified, and the results are descriptive rather than statistically quantified.
major comments (4)
- [Tables 1 and 3] The reported plans violate the BED constraints that the paper defines as hard constraints in Eq. (1). In Table 1, ARC bladder BED50 is 44.60 Gy against an upper bound of 40 Gy, bladder BED20 is 67.75 Gy against 63 Gy, rectum BED50 is 45.02 Gy against 40 Gy, and rectum BED20 is 65.70 Gy against 63 Gy. In Table 3, brainstem BEDmax is 85.25 Gy (ARC) and 87.91 Gy (multi-IMPT) against 83.7 Gy, and brain BEDmax is 104.31 Gy (ARC) and 109.37 Gy (multi-IMPT) against 96 Gy. Since the paper states that BEDp means at most p% of OAR voxels should receive BED greater than the upper bound, these values are direct violations, not reporting artifacts. The equivalence comparison is therefore between infeasible, non-clinically-deliverable plans, and the conclusion that multi-IMPT can replace ARC is not supported by the presented data. The authors should re-optimize with constraints enforced or provide a clear explanation of why these violations are acceptable.
- [Section 2.1 and Eq. (1)] There is an internal inconsistency between the stated target DVH-min constraint and the implemented constraint. Section 2.1, constraint 4, defines a physical-dose constraint: at least p fraction of target voxels receive physical dose d_jt^T >= d_min. However, Eq. (1) and Eq. (2) replace this with a BED constraint of the form sum_t z_jt^T + rho_T sum_t (z_jt^T)^2 >= BED_DVH^T, which is a different object and is not defined in the text. This means the optimization that produced the reported plans may not match the problem formulation in the paper. Please clarify which constraint was actually used, define BED_DVH^T and rho_T, and correct the equations accordingly.
- [Section 2.2] The ARC comparator is modeled as u_t = u for all fractions, i.e., the same plan is delivered in every fraction with all fields active. This equal-dose-per-fraction assumption is a simplification of clinical proton arc therapy, where the gantry rotates continuously and the delivered dose per control point may vary. The equivalence claim is therefore only with respect to this simplified ARC model. The authors should justify that this is a representative ARC baseline or investigate how the comparison changes under a more general ARC delivery model.
- [Section 3] The results are based on four single-patient cases (one per site) and are compared descriptively, with no error bars, uncertainty quantification, or statistical tests. The statement in the conclusion that the dosimetric differences are 'not clinically significant' is consequently not substantiated. The authors should either add a quantitative analysis (e.g., multiple patients per site or at least explicit thresholds of clinical significance) or temper the claims to reflect the limited evidence.
minor comments (6)
- [Eq. (1)] The symbol 'px' is used for the prescription dose but is never defined; please define it explicitly.
- [Section 2.1] The active index set uses the condition j >= p*n, which is ambiguous when p*n is not an integer; please clarify the rounding convention and whether the constraint should apply to the floor or ceiling of p*n.
- [Appendix A] The projection formulas for z_jt contain expressions for q that are not fully legible due to typesetting; please rewrite them explicitly so the algorithm is reproducible.
- [Tables 2 and 4] Table 2 lists Heart BEDmean and Table 4 lists R Parotid BEDmax without upper bounds; please state explicitly whether these structures were unconstrained or whether the bounds were omitted from the table.
- [Section 2.1, BED-DVH constraint] The definition of the active index set uses j >= p*n, but for a constraint that at most p% of voxels exceed the bound, the relevant index should be the first voxel above the allowed fraction; please reconcile the indexing with the stated clinical meaning.
- [Algorithm 1] In the algorithm outline, Step 4b says 'update primal variables' and Step 4c says 'update dual variables', but in the Appendix the order of presentation is reversed; please make the algorithm listing and the detailed description consistent.
Circularity Check
No significant circularity: the ARC-equivalence claim is an empirical planning comparison with independent beam geometries; self-citations are algorithmic and non-load-bearing.
full rationale
The paper's central claim is a treatment-planning comparison, not a parameter-fit or a derived theorem. The multi-IMPT plan and the ARC plan are generated from different beam-angle subsets and fractionation patterns (Eq. (1) with fraction-dependent u_t for multi-IMPT versus u_t = u for ARC), and their BED-to-OAR values are then compared as an endpoint. The BED formula with alpha/beta = 2 Gy is a stated biological assumption, and the fact that the same BED model appears both in the constraints and in the comparison metric makes the equivalence claim model-internal, but it does not make the comparison equal by construction: no constant is fitted to the reported BED differences, and no equation forces the two plans to coincide. The self-citations (ICR/ADMM/MMU, refs. [29,31,32,33,34,36]) are algorithmic provenance; Algorithm 1 and the appendix specify the variable updates, so the equivalence conclusion does not rest on an unverified self-cited theorem. The reported BED upper-bound violations (e.g., Table 1 prostate ARC bladder BED50 44.60 Gy versus the 40 Gy upper bound; Table 3 brain multi-IMPT BEDmax 109.37 Gy versus the 96 Gy upper bound) are a correctness and deliverability risk, not a circularity, because they do not make the comparison an identity. Overall circularity score is 2 due solely to minor non-load-bearing self-citations; the core comparison is self-contained against the stated ARC model.
Assumptions & free parameters
free parameters (5)
- alpha_over_beta_OAR =
2 Gy (all OAR)
- MMU_threshold_g =
not reported
- ARC_angle_spacing =
15 degrees over 360
- multi_IMPT_beam_subsets =
six rotated 4-beam sets at 0,15,30,45,60,75 deg offsets
- optimization_weights_mu_w =
not reported
assumptions (5)
- domain assumption Linear-quadratic BED model correctly predicts normal tissue complication for fractionated proton therapy.
- domain assumption alpha/beta = 2 Gy is representative for all organs at risk.
- domain assumption The ARC model with equal dose per fraction (u_t = u) represents clinical spot-scanning proton arc therapy.
- domain assumption MatRad dose influence matrices are accurate enough for both plans.
- ad hoc to paper ADMM applied to the non-convex problem (2) with random initialization converges to a good local solution.
Cite this review
Pith. "Pith review of Multi-IMPT: a biologically equivalent approach to proton ARC therapy." pith.science (2026). https://pith.science/paper/BUCQUDUF
@misc{pith2026241117578,
author = {Pith},
title = {Pith review of: Multi-IMPT: a biologically equivalent approach to proton ARC therapy},
year = {2026},
howpublished = {\url{https://pith.science/paper/BUCQUDUF}},
note = {Machine review of arXiv:2411.17578}
}
read the original abstract
Objective: Proton spot-scanning arc therapy (ARC) is an emerging modality that can improve the high-dose conformity to targets compared with standard intensity-modulated proton therapy (IMPT). However, the efficient treatment delivery of ARC is challenging due to the required frequent energy changes during the continuous gantry rotation. This work proposes a novel method that delivers a multiple IMPT (multi-IMPT) plan that is equivalent to ARC in terms of biologically effective dose (BED). Approach: The proposed multi-IMPT method utilizes a different subset of limited number of beam angles in each fraction for dose delivery. Due to the different dose delivered to organs at risk (OAR) in each fraction, we optimize biologically effective dose (BED) for OAR and the physical dose delivered for target in each fraction. The BED-based multi-IMPT inverse optimization problem is solved via the iterative convex relaxation method and the alternating direction method of multipliers. The effectiveness of the proposed multi-IMPT method is evaluated in terms of dose objectives in comparison with ARC. Main results: Multi-IMPT provided similar plan quality with ARC. For example, multi-IMPT provided better OAR sparing and slightly better target dose coverage for the prostate case; similar dose distribution for the lung case; slightly worse dose coverage for the brain case; better dose coverage but slightly higher BED in OAR for the head-and-neck case. Significance: We have proposed a multi-IMPT approach to deliver ARC-equivalent plan quality. Keywords: biologically effective dose (BED), proton arc therapy
Reference graph
Works this paper leans on
-
[1]
Introduction Sandison et al. [1] proposed passive-scattering-based proton arc radiotherapy (PS-ARC) and demonstrated an approach could improve the target dose conformity. However, the implementation of PS-ARC [2, 3, 4] faced several limitations including the need to change beam compensator and range modulation wheel during gantry rotation. The technologic...
-
[2]
!∈𝑅#!×%": dose influence matrix of 𝑚-th OAR during fraction 𝑡; 𝑘
Problem Formulation In this section, we start by defining the parameters for the optimization problem (including the definition of the BED), followed by defining the constraints to minimize the BED in OAR, and finally introducing the complete optimization model for multi-IMPT. 2.1. Defining parameters, decision variables and constraints in the multi-IMPT ...
-
[3]
BED-mean constraint for OAR [21, 22, 23, 24]: Let 𝑀+ be the set of OAR whose small portion can be damaged without affecting their function. For such OAR, the BED-mean constraint bounds the mean BED (𝐵𝐸𝐷!.,#!) delivered to all voxels in OAR 𝑚. Thus, we define CC𝐴&"!𝑢"( ")*+ CC𝜌!(𝐴&"!𝑢")+( ")* ≤#! &)* 𝑛!×𝐵𝐸𝐷!.,#! ∀ 𝑚∈𝑀+.#! &)*
-
[4]
BED-DVH max constraint for OAR [21, 22, 23, 24]: Consider the set of OAR 𝑀/. The BED-DVH constraints states that for any OAR 𝑚∈𝑀/, at most 𝑝 fraction of voxels should receive BED larger than 𝐵𝐸𝐷01!, i.e., 𝐵𝐸𝐷&!=∑𝐴&"!𝑢"(")*+ ∑𝜌!(𝐴&"!𝑢")+≥(")* 𝐵𝐸𝐷01! for at most 𝑝×𝑛! voxels. One of the commonly used techniques to define the DVH max constraint is to first de...
-
[5]
DVH min constraint for target [25, 26]: DVH min constraint ensures that at least 𝑝 fraction of the target voxels receive physical dose larger than 𝑑01' in each fraction 𝑡, i.e., 𝑑&"'≥𝑑01', for at least 𝑝×𝑛' voxels. To define the DVH min constraint, we first define the active index set for the target as Ω'={𝑗∈[𝑛'2] | 𝑗≤𝑝×𝑛!,𝑑&"'≤𝑑01'}, 6 where [𝑛'2] is the...
-
[6]
Conclusion and Discussion In this work, we propose an ARC-equivalent IMPT method, termed multi-IMPT, which consists of multiple IMPT plans. The multi-IMPT method utilizes different combinations of a small subset of beams in each fraction, which sums up to a large set of beams as in ARC. It was shown that the multi-IMPT can deliver the dose coverage equiva...
-
[7]
! for all 𝑚∈𝑀+: For each 𝑚∈𝑀+, we fix all the variables except 𝑧&
Updating 𝑧"! for all 𝑚∈𝑀+: For each 𝑚∈𝑀+, we fix all the variables except 𝑧&"! in Eq. (3). The resulting minimization problem over 𝑧&"! is 𝑚𝑖𝑛CC;𝐴&"!𝑢"+𝜆&"!−𝑧&"!<+#! &)* ( ")* s.t.CCk𝑧&"!+12𝜌!l+#! &)* ≤( ")* 𝑛!𝐵𝐸𝐷!.,#!𝜌!+𝑇𝑛!4𝜌!+. We observe that the optimal solution to this problem is the projection of 𝐴&"!𝑢"+𝜆&"! onto the inequality constraint. Thus, the...
-
[8]
! for all 𝑚∈𝑀/: For each 𝑚∈𝑀*, 𝑗∈ΩE, we use the same procedure as outlined for updating 𝑧
Updating 𝑧"! for all 𝑚∈𝑀/: For each 𝑚∈𝑀*, 𝑗∈ΩE, we use the same procedure as outlined for updating 𝑧"! for 𝑚∈𝑀*. 6. Updating 𝑧"': For each 𝑗∈[𝑛'], we fix all variables except 𝑧&"' in Eq. (3). The resulting minimization problem is 𝑚𝑖𝑛C;𝐴&"'𝑢"+𝜆&"'−𝑧&"'<+( ")*s.t.Ck𝑧&"'+12𝜌'l+≥( ")* 𝐵𝐸𝐷01'𝜌'+𝑇4𝜌'+. The optimal solution to this problem is the projection of 𝐴...
work page 1997
Show all 13 references
-
[9]
Treatment planning and verification of proton therapy using spot scanning: initial experiences,
A. J. Lomax, T. Böhringer, A. Bolsi, D. Coray, F. Emert, G. Goitein, M. Jermann, S. Lin, E. Pedroni, H. Rutz and others, "Treatment planning and verification of proton therapy using spot scanning: initial experiences," Medical physics, vol. 31, p. 3150–3157, 2004. [9] E. Pedro...
2004
-
[16]
Fast robust optimization of proton PBS arc therapy plans using early energy layer selection and spot assignment,
E. Engwall, C. Battinelli, V. Wase, O. Marthin, L. Glimelius, R. Bokrantz, B. Andersson and A. Fredriksson, "Fast robust optimization of proton PBS arc therapy plans using early energy layer selection and spot assignment," Physics in Medicine & Biology, vol. 67, p. 065010, 202...
2022
-
[25]
Clinically relevant intensity modulation opitimization using physical criteria,
T. Bortfeld, "Clinically relevant intensity modulation opitimization using physical criteria," in XII International Conference on the Use of Computers in Radiation Therapy, 1997, 1997. [26] Q. Wu and R. Mohan, "Algorithms and functionality of an intensity modulated radiotherap...
1997
-
[35]
Distributed optimization and statistical learning via the alternating direction method of multipliers,
S. Boyd, N. Parikh, E. Chu, B. Peleato, J. Eckstein and others, "Distributed optimization and statistical learning via the alternating direction method of multipliers," Foundations and Trends® in Machine learning, vol. 3, p. 1–122, 2011. [36] H. Gao, "Robust fluence map optimi...
2011
-
[44]
Improving dosimetric outcome for hippocampus and cochlea sparing whole brain radiotherapy using spot-scanning proton arc therapy,
X. Ding, J. Zhou, X. Li, K. Blas, G. Liu, Y. Wang, A. Qin, P. Chinnaiyan, D. Yan, C. Stevens and others, "Improving dosimetric outcome for hippocampus and cochlea sparing whole brain radiotherapy using spot-scanning proton arc therapy," Acta Oncologica, vol. 58, p. 483–490, 20...
2019
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.