<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Oncol.</journal-id>
<journal-title>Frontiers in Oncology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Oncol.</abbrev-journal-title>
<issn pub-type="epub">2234-943X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fonc.2024.1369065</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Oncology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Feasibility and constraints of Bragg peak FLASH proton therapy treatment planning</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>L&#xf6;vgren</surname>
<given-names>Nathalie</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2627745"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fagerstr&#xf6;m Kristensen</surname>
<given-names>Ingrid</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Petersson</surname>
<given-names>Kristoffer</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1445606"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Oncology, Oxford Institute for Radiation Oncology, University of Oxford</institution>, <addr-line>Oxford</addr-line>, <country>United Kingdom</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Clinical Oncology, Department of Clinical Sciences Lund, Lund University</institution>, <addr-line>Lund</addr-line>, <country>Sweden</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Department of Haematology, Oncology and Radiation Physics, Sk&#xe5;ne University Hospital</institution>, <addr-line>Lund</addr-line>, <country>Sweden</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Ranald Mackay, The University of Manchester, United Kingdom</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Yoichi Watanabe, University of Minnesota Twin Cities, United States</p>
<p>Yvonne Dzierma, Saarland University Hospital, Germany</p>
<p>Sunyoung Jang, The Pennsylvania State University, United States</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Nathalie L&#xf6;vgren, <email xlink:href="mailto:Nathalie.Lovgren@oncology.ox.ac.uk">Nathalie.Lovgren@oncology.ox.ac.uk</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>26</day>
<month>04</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>14</volume>
<elocation-id>1369065</elocation-id>
<history>
<date date-type="received">
<day>11</day>
<month>01</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>02</day>
<month>04</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 L&#xf6;vgren, Fagerstr&#xf6;m Kristensen and Petersson</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>L&#xf6;vgren, Fagerstr&#xf6;m Kristensen and Petersson</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>FLASH proton therapy (FLASH-PT) requires ultra-high dose rate (&#x2265;&#xa0;40 Gy/s) protons to be delivered in a short timescale whilst conforming to a patient-specific target. This study investigates the feasibility and constraints of Bragg peak FLASH-PT treatment planning, and compares the <italic>in silico</italic> results produced to plans for intensity modulated proton therapy (IMPT).</p>
</sec>
<sec>
<title>Materials and method</title>
<p>Bragg peak FLASH-PT and IMPT treatment plans were generated for bone (n=3), brain (n=3), and lung (n=4) targets using the MIROpt research treatment planning system and the Conformal FLASH library developed by Applications SA from the open-source version of UCLouvain. FLASH-PT beams were simulated using monoenergetic spot-scanned protons traversing through a conformal energy modulator, a range shifter, and an aperture. A dose rate constraint of &#x2265; 40 Gy/s was included in each FLASH-PT plan optimisation.</p>
</sec>
<sec>
<title>Results</title>
<p>Space limitations in the FLASH-PT adapted beam nozzle imposed a maximum target width constraint, excluding 4 cases from the study. FLASH-PT plans did not satisfy the imposed target dose constraints (D<sub>95%</sub> &#x2265; 95% and D<sub>2%</sub>&#x2264; 105%) but achieved clinically acceptable doses to organs at risk (OARs). IMPT plans adhered to all target and OAR dose constraints. FLASH-PT plans showed a reduction in both target homogeneity (p <italic>&lt;</italic> 0.001) and dose conformity (non-significant) compared to IMPT.</p>
</sec>
<sec>
<title>Conclusion</title>
<p>Without accounting for a sparing effect, IMPT plans were superior in target coverage, dose conformity, target homogeneity, and OAR sparing compared to FLASH-PT. Further research is warranted in treatment planning optimisation and beam delivery for clinical implementation of Bragg peak FLASH-PT.</p>
</sec>
</abstract>
<kwd-group>
<kwd>FLASH radiotherapy</kwd>
<kwd>proton therapy</kwd>
<kwd>feasibility studies</kwd>
<kwd>intensity-modulated proton therapy</kwd>
<kwd>Bragg peak</kwd>
</kwd-group>
<contract-sponsor id="cn001">Medical Research Council<named-content content-type="fundref-id">10.13039/501100000265</named-content>
</contract-sponsor>
<counts>
<fig-count count="7"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="38"/>
<page-count count="10"/>
<word-count count="5236"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Radiation Oncology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>FLASH radiotherapy (FLASH-RT) is an emerging technique aiming to deliver ultra-high dose rate (&#x2265; 40 Gy/s) radiation to treat patients. Pre-clinically, the use of FLASH-RT has been shown to induce a normal tissue sparing effect whilst maintaining the anti-tumour effectiveness of conventional radiotherapy techniques, known as the &#x201c;FLASH effect&#x201d; (<xref ref-type="bibr" rid="B1">1</xref>&#x2013;<xref ref-type="bibr" rid="B14">14</xref>). Compared to conventional radiotherapy, the use of ultra-high dose rates for FLASH-RT allows for a substantial reduction in overall treatment times. FLASH-RT may therefore be able to mitigate uncertainties associated with intra-fractional motion, whilst also potentially increasing the availability and accessibility of radiotherapy treatments. These aspects of FLASH-RT could benefit treatments of targets where organ motion is a prevalent issue, normally requiring the use of large margins, gating, or tracking, such as the lungs, bowels, and liver (<xref ref-type="bibr" rid="B15">15</xref>). Treatments of regions where critical structures are at risk of receiving increased doses due to their proximity to the target, such as the brain and spinal cord, could also benefit from the normal tissue sparing of FLASH-RT.</p>
<p>Dose rates can be defined in different ways; two common definitions include the average dose rate (ADR) which accounts for the dead-time between each spot delivered, and the dose-averaged dose rate (DADR) which only considers the instantaneous dose rate delivered (<xref ref-type="bibr" rid="B16">16</xref>). The &#x2018;dose rate coverage&#x2019; or &#x2018;FLASH coverage&#x2019; are terms used to quantify the volume of an irradiated structure receiving dose at ultra-high dose rates. It is defined as the ratio of the organ at risk (OAR) volume receiving doses at ultra-high dose rates to the OAR volume receiving a defined threshold dose. Several studies report to achieve a FLASH dose rate coverage of approximately 80% to OARs when evaluating the ADR, however, this is highly dependent on how the study has defined &#x2018;FLASH dose rate coverage&#x2019; as dose and dose rate thresholds vary (<xref ref-type="bibr" rid="B17">17</xref>&#x2013;<xref ref-type="bibr" rid="B22">22</xref>)</p>
<p>Due to limited depth-penetration of clinically used electron beams, and vast power requirements and heating effects of producing ultra-high dose rate photons, protons are one of the most promising candidates for treatments of deep-seated targets. Heavier ion-beams and very-high energy electrons may also be contenders; however, ultra-high dose rate protons are readily available given few modifications to clinical delivery systems (<xref ref-type="bibr" rid="B23">23</xref>, <xref ref-type="bibr" rid="B24">24</xref>). FLASH proton therapy (FLASH-PT) has therefore become a rapidly progressing branch within FLASH-RT research. Following the successful completion of the FAST-01 clinical trial which demonstrated the safe delivery of transmission beam FLASH-PT to patients, the FAST-02 trial is currently investigating the toxicities and pain relief associated with the treatment technique (<xref ref-type="bibr" rid="B25">25</xref>&#x2013;<xref ref-type="bibr" rid="B27">27</xref>).</p>
<p>Despite transmission beam FLASH-PT showing promising results, the technique does not maximise the benefits associated with the use of the Bragg peak in standard proton therapy. Bragg peak proton therapy allows for improved dose conformity, with the use of fewer beams, and reduced doses to OARs when compared to transmission beam proton therapy (<xref ref-type="bibr" rid="B20">20</xref>, <xref ref-type="bibr" rid="B28">28</xref>). The use of fewer beams allows for an increase in dose rates delivered by each beam, making Bragg peak FLASH-PT more favourable than transmission beam FLASH-PT. However, without accounting for a sparing effect, previous Bragg peak and transmission beam FLASH-PT treatment planning studies have shown inferiority to conventional proton therapy plans (<xref ref-type="bibr" rid="B20">20</xref>, <xref ref-type="bibr" rid="B29">29</xref>). This is linked to both the beam delivery and treatment planning process. As the incorporation of dose rates into the plan optimisation has not been needed prior to FLASH-RT, the software capabilities and limitations are not well understood.</p>
<p>With regards to beam delivery, conventional energy-modulation techniques for standard Bragg peak proton therapy are not fast enough for ultra-high dose rates to be achieved. Instead, Bragg peak FLASH-PT can make use of patient-specific static energy degrades. Previous studies have adopted the use of range shifters or ridge filters, with the inclusion of collimators, along the path of a mono-energetic proton beam (<xref ref-type="bibr" rid="B20">20</xref>&#x2013;<xref ref-type="bibr" rid="B22">22</xref>, <xref ref-type="bibr" rid="B30">30</xref>). However, these systems have been associated with enlarged spot sizes due to multiple Coulomb scattering, and an increase in spot-to-spot distance to reduce overlap in the ridge filters, when compared to conventional proton therapy beam delivery (<xref ref-type="bibr" rid="B20">20</xref>&#x2013;<xref ref-type="bibr" rid="B22">22</xref>, <xref ref-type="bibr" rid="B29">29</xref>). How these increases impact the FLASH-PT treatments have not been well established; alternative beam setups and treatment planning optimisations could reduce their importance and thereby the deviations of FLASH-PT treatment planning from that of standard proton therapy. In addition to the difference in spot sizes, the time scale of gantry rotation for standard proton therapy machines are drastically longer than the treatment delivery time for FLASH-PT. Multi-beam fractions can therefore not be utilised for the ultra-high dose rate technique. To bypass this limitation, and to simplify the treatment setup due to the use of patient-and field-specific static energy degraders, it may be beneficial to deliver Bragg peak FLASH-PT using one beam per fraction.</p>
<p>The rapid progression of FLASH-PT, from pre-clinical studies to clinical trials, necessitates a full understanding of the capabilities, limitations, and safety of the treatment modality before it is implemented into clinical practice. Few studies focusing on treatment plans for FLASH-PT delivered in one beam per fraction have been conducted for Bragg peak proton beams, and no guidelines regarding the treatment planning process for this technique exists. The aim of this study was to investigate the feasibility and constraints of Bragg peak FLASH-PT treatment planning, using a novel beam setup and one beam per fraction, and determine if these plans are comparable to clinical intensity-modulated proton therapy (IMPT) plans.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<p>Ten patient cases were planned with IMPT and FLASH-PT and consisted of bone (n = 3), brain (n = 4), and lung (n = 3) targets. Target sizes ranged between 1.42-324.09 cm<sup>3</sup>, with a median size of 16.18 cm<sup>3</sup>. Details of the target types, defined OARs, target sizes, prescribed doses, and number of fractions for the individual patient cases are shown in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Table outlining the target type, defined organs at risk, target size, prescribed dose (Gy), and number of fractions used in the treatment plans for each patient case.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Patient Case</th>
<th valign="top" align="left">Target Type</th>
<th valign="top" align="left">Defined Organs at Risk</th>
<th valign="top" align="left">Target Size (cm<sup>3</sup>)</th>
<th valign="top" align="left">Prescribed Dose (Gy)</th>
<th valign="top" align="left">Number of Fractions</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">1</td>
<td valign="top" align="left">Bone</td>
<td valign="top" align="left">Body</td>
<td valign="top" align="left">62.20</td>
<td valign="top" align="left">8</td>
<td valign="top" align="left">1</td>
</tr>
<tr>
<td valign="top" align="left">2</td>
<td valign="top" align="left">Bone</td>
<td valign="top" align="left">Body</td>
<td valign="top" align="left">324.09</td>
<td valign="top" align="left">8</td>
<td valign="top" align="left">1</td>
</tr>
<tr>
<td valign="top" align="left">3</td>
<td valign="top" align="left">Bone</td>
<td valign="top" align="left">Left Lung, Right Lung, Body</td>
<td valign="top" align="left">207.07</td>
<td valign="top" align="left">8</td>
<td valign="top" align="left">1</td>
</tr>
<tr>
<td valign="top" align="left">4</td>
<td valign="top" align="left">Brain</td>
<td valign="top" align="left">Body</td>
<td valign="top" align="left">19.71</td>
<td valign="top" align="left">20</td>
<td valign="top" align="left">2</td>
</tr>
<tr>
<td valign="top" align="left">5</td>
<td valign="top" align="left">Brain</td>
<td valign="top" align="left">Brain, Body</td>
<td valign="top" align="left">1.42</td>
<td valign="top" align="left">30</td>
<td valign="top" align="left">3</td>
</tr>
<tr>
<td valign="top" align="left">6</td>
<td valign="top" align="left">Brain</td>
<td valign="top" align="left">Brain, Body</td>
<td valign="top" align="left">2.26</td>
<td valign="top" align="left">30</td>
<td valign="top" align="left">3</td>
</tr>
<tr>
<td valign="top" align="left">7</td>
<td valign="top" align="left">Brain</td>
<td valign="top" align="left">Brain, Body</td>
<td valign="top" align="left">12.65</td>
<td valign="top" align="left">34</td>
<td valign="top" align="left">2</td>
</tr>
<tr>
<td valign="top" align="left">8</td>
<td valign="top" align="left">Lung</td>
<td valign="top" align="left">Left Lung, Thoracic Wall, Body</td>
<td valign="top" align="left">11.98</td>
<td valign="top" align="left">45</td>
<td valign="top" align="left">3,4,5</td>
</tr>
<tr>
<td valign="top" align="left">9</td>
<td valign="top" align="left">Lung</td>
<td valign="top" align="left">Left Lung, Bronchial Tree, Body</td>
<td valign="top" align="left">10.66</td>
<td valign="top" align="left">45</td>
<td valign="top" align="left">3,4,5</td>
</tr>
<tr>
<td valign="top" align="left">10</td>
<td valign="top" align="left">Lung</td>
<td valign="top" align="left">Right Lung, Thoracic Wall, Great Vessel, Heart, Body</td>
<td valign="top" align="left">21.61</td>
<td valign="top" align="left">45</td>
<td valign="top" align="left">3,4,5</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Treatment plans were produced using a MATLAB-based simulation software, utilising an open-source Monte Carlo engine known as MCsquare (many-core Monte Carlo) (<xref ref-type="bibr" rid="B31">31</xref>, <xref ref-type="bibr" rid="B32">32</xref>). An open-source robust IMPT research treatment planning system (TPS), openMIROpt, produced by UCLouvain and accessed through OpenReggui was used to produce the IMPT plans (<xref ref-type="bibr" rid="B33">33</xref>). FLASH-PT plans were generated using the same open-source software, along with a ConformalFLASH library (<ext-link ext-link-type="uri" xlink:href="http://www.openFLASH.software">www.openFLASH.software</ext-link>), developed by Ion Beam Applications SA (IBA), allowing for dose rate calculation and optimisation (<xref ref-type="bibr" rid="B34">34</xref>). The treatment plans were simulated for use on the IBA Proteus<sup>&#xae;</sup>Plus (Louvain-la-Neuve, Belgium) machine. The maximum field-size currently available from IBA for this beam setup is 8&#xd7;8&#xd7;8 cm<sup>3</sup>. The IMPT spots were placed in a hexagonal pattern for optimal dose coverage, and a constant relative biological effectiveness (RBE) of 1.1 was assumed. Robust planning is not yet available for ConformalFLASH, and so planning target volume-based treatment planning was adopted for all treatment plans to allow for a fair comparison. Additionally, a couch angle of 0<sup>&#xb0;</sup> was used as the research TPS does not allow for uses of other angles. Proton beams were planned using a cyclotron with a clinically used beam current of 300 nA for IMPT, and the maximum available beam current of 500 nA for FLASH-PT.</p>
<p>FLASH-PT beams were simulated using monoenergetic spot-scanned protons at an energy of 230 MeV. The ultra-high dose rate protons traverse through a conformal energy modulator (CEM), a range shifter, and an aperture, before reaching the target volume. The CEM is a patient- and field-specific &#x2018;hedgehog&#x2019;, which modulates the protons to the 3D shape of the target; it is made from a 3D printed plastic-like material called &#x2018;Tusk&#x2019;. The range shifter is used to set the distal range of the delivered dose to the target, and the aperture to reduce dose contributions from scattered particles and to sharpen the penumbra.</p>
<p>Standard IMPT treatment plans are optimised using either single-field optimisation (SFO) or multi-field optimisation (MFO). Plans are optimised using SFO when a single beam is used to cover the entire target and each beam is optimised independently from the others. MFO allows for use of multiple beams to cover the entire target, and the dose contribution of each beam is considered in the dose calculation and optimisation of the other beams. This improves target coverage whilst minimising dose to surrounding structures, making MFO the standard optimisation technique used in treatment planning (<xref ref-type="bibr" rid="B29">29</xref>). As Bragg peak FLASH-PT requires dose delivery to be made using one beam per fraction, due to ultra-high dose rate requirements and to improve target homogeneity, all FLASH-PT plans in this study were optimised using SFO. For treatments of simpler structures, such as the vertebrae cases in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>, a single-field fraction was used. However, for more advanced structures, such as the brain and lung targets, the use of a single field would lead to increased doses to OARs and worsened dose conformity. In these cases, multiple single-field fractions were used. Different fraction numbers were investigated for the three lung cases, as seen in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>. Initially, the idea was to generate plans in line with standard fractionations schemes for this target type by using 3 and 5 fractions. However, with an increase in number of fractions there is a decrease in intra-fraction dose, and thus dose rates, to OARs. To determine if there is a middle ground between reducing doses to OARs and increasing the dose rates delivered, the use of 4 fractions was also investigated. In addition to the standard MFO IMPT plans being generated, SFO IMPT plans were also produced to allow for fair comparisons of the SFO FLASH-PT plans to IMPT. This will show any discrepancies between the plans which are solely due to the use of ultra-high dose rates and those which may be due to the use of SFO instead of MFO. The sum of the SFO plans, for all fractions, will be used to compare the FLASH-PT and IMPT SFO plans to those of IMPT MFO.</p>
<p>A percentile dose rate definition, similar to that of the ADR, is adopted by the ConformalFLASH library; it is defined as the ratio of 98% of dose delivered to the time taken to deliver that dose (<xref ref-type="bibr" rid="B34">34</xref>). A minimum dose rate constraint of &#x2265; 40 Gy/s is implemented for each OAR volume receiving a dose &#x2265; 2 Gy. Dose rate coverage in this study is therefore defined as the volume of OAR which adheres to this constraint. As pre-clinical studies have found varying values for the dose modifying factor (DMF) for the FLASH effect, a DMF of 1 is used in all treatment plans (<xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B10">10</xref>, <xref ref-type="bibr" rid="B11">11</xref>). A &#x201c;FLASH effect&#x201d; is therefore not assumed in the FLASH-PT treatment plans. Each plan was normalised such that D<sub>50%</sub> = D<italic>
<sub>p</sub>
</italic>, where D<italic>
<sub>p</sub>
</italic> is the prescribed dose. Dose constraints used for the target structures were D<sub>95%</sub> &#x2265; 95% and D<sub>2%</sub> &#x2264; 105% of the relative target volumes.</p>
<p>Dose distribution maps, dose volume histograms (DVHs), and dose statistics were used to qualitatively compare the dose distributions of IMPT MFO, IMPT SFO, and FLASH-PT plans. Quantitative comparisons are made using DVH parameters and dose statistics; this allows for conformity and homogeneity indices to be determined, along with statistical comparisons of the DVH parameters using the two-sided Wilcoxon Rank Sum Test, Spearman&#x2019;s Rank Correlation Coefficient (<italic>&#x3c1;</italic>
<sub>SRC</sub>), and Pearson Correlation Coefficient (<italic>&#x3c1;</italic>
<sub>PCC</sub>), with significance defined for p&#xa0;<italic>&lt;</italic>&#xa0;0.05. Different samples involving the mean dose, V<sub>95%</sub>, V<sub>100%</sub>, V<sub>105%</sub>, V<sub>107%</sub>, D<sub>2%</sub>, D<sub>5%</sub>, D<sub>95%</sub>, and D<sub>98%</sub> were tested. Clinical acceptability for doses delivered to OARs were defined using the European Particle Therapy Network (EPTN) consensus and the Quantitative Analysis of Normal Tissue Effects in the Clinic (QUANTEC) review (<xref ref-type="bibr" rid="B35">35</xref>, <xref ref-type="bibr" rid="B36">36</xref>).</p>
<p>The radiation conformity index (RCI) was evaluated using the ratio of the target volume to the total volume receiving 95% of the prescribed dose (RCI = V<italic>
<sub>target</sub>
</italic>/V<sub>95%</sub>) (<xref ref-type="bibr" rid="B37">37</xref>). The homogeneity index (HI) was calculated using HI =(D<sub>2%</sub>-D<sub>98%</sub>)/D<sub>50%</sub> as defined in the ICRU Report 83 (<xref ref-type="bibr" rid="B38">38</xref>). All treatment planning data was included in these comparisons.</p>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<p>The field size currently available from IBA imposed a size constraint on the targets. This is due to limitation in the space between the CEM and range shifter in the FLASH-PT adapted beam nozzle. Treatment plans could therefore not be simulated for targets with width (inline or crossline) above 70&#xa0;mm. Patient cases 1, 2, 3, and 4 all exceeded this size limit and were therefore excluded from the study. A minimum angle criterion of &#x2265; 40&#xb0; between each beam was chosen to reduce beam overlap and hotspots for the SFO plans. However, the space limitation also affected which beam angles that could be used, likely due to the target sizes along with the margins for spot placements exceeding 70&#xa0;mm. Three FLASH-PT plans could therefore not adhere to the minimum angle criterion: Case 8 (5 fractions) and Case 10 (3 and 5 fractions).</p>
<p>IMPT MFO achieves the clinical dose constraints imposed, with D<sub>95%</sub> <italic>&gt;</italic> 95% and D<sub>2%</sub> <italic>&lt;</italic> 105%, as seen in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>. IMPT SFO has two cases where these dose limits are not achieved: Cases 5 and 6 which have the two smallest target volumes, as seen in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>. FLASH-PT, however, only achieves similar mean doses to the other two IMPT techniques, with D<sub>95%</sub> and D<sub>2%</sub> being significantly different to both IMPT MFO and IMPT SFO, as shown in <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>. Since IMPT MFO and IMPT SFO are comparable, the discrepancy seen for FLASH-PT is not only due to the use of SFO but also due to the limitations imposed for the use of ultra-high dose rate protons. In addition to this, FLASH-PT is associated with a significant increase in dose to OARs, as seen in <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>. Nevertheless, the doses are still within clinically acceptable limits, as defined by the EPTN consensus and the QUANTEC review (<xref ref-type="bibr" rid="B35">35</xref>, <xref ref-type="bibr" rid="B36">36</xref>). An increase in the dose at the distal edge of the targets were observed for all FLASH-PT plans compared to IMPT, with examples shown in <xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>. This is likely due to range straggling of the 230 MeV protons following interactions with the energy degraders used in this beam setup; and would also contribute to the increase in dose to OARs observed for FLASH-PT. Dose volume histograms (DVHs) for the different cases, fractionation schemes, and treatment techniques are shown in the <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Appendix</bold>
</xref>.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Plot summarising target dose metrics. D<sub>95%</sub> (square), mean dose (circle), and D<sub>2%</sub> (diamond) are normalised with respect to the prescribed dose for all patient cases and for each treatment technique. The dashed lines indicate 95%, 100%, and 105% of the prescribed dose.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1369065-g001.tif"/>
</fig>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Table of median values of: mean dose (%), D<sub>98%</sub> (%), D<sub>95%</sub> (%), D<sub>5%</sub> (%), D<sub>2%</sub> (%), V<sub>95%</sub> (%), V<sub>100%</sub> (%), V<sub>105%</sub> (%), V<sub>107%</sub> (%), radiation conformity index, and homogeneity index for IMPT MFO, IMPT SFO, and FLASH-PT.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="bottom" align="left">Dose Parameters</th>
<th valign="bottom" align="left">IMPT MFO</th>
<th valign="bottom" align="left">IMPT SFO</th>
<th valign="bottom" align="left">FLASH-PT</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="bottom" align="left">Mean Dose (%)</td>
<td valign="bottom" align="left">99.86</td>
<td valign="bottom" align="left">99.79</td>
<td valign="bottom" align="left">99.98</td>
</tr>
<tr>
<td valign="bottom" align="left">D<sub>98%</sub> (%)</td>
<td valign="bottom" align="left">94.40</td>
<td valign="bottom" align="left">94.37</td>
<td valign="bottom" align="left">91.64 <sup>*</sup>
</td>
</tr>
<tr>
<td valign="bottom" align="left">D<sub>95%</sub> (%)</td>
<td valign="bottom" align="left">95.79</td>
<td valign="bottom" align="left">95.54</td>
<td valign="bottom" align="left">93.06 <sup>**</sup>
</td>
</tr>
<tr>
<td valign="bottom" align="left">D<sub>5%</sub> (%)</td>
<td valign="bottom" align="left">103.69</td>
<td valign="bottom" align="left">103.30</td>
<td valign="bottom" align="left">106.76 <sup>**</sup>
</td>
</tr>
<tr>
<td valign="bottom" align="left">D<sub>2%</sub> (%)</td>
<td valign="bottom" align="left">104.47</td>
<td valign="bottom" align="left">103.98</td>
<td valign="bottom" align="left">107.87 <sup>***</sup>
</td>
</tr>
<tr>
<td valign="bottom" align="left">V<sub>95%</sub> (%)</td>
<td valign="bottom" align="left">96.87</td>
<td valign="bottom" align="left">96.76</td>
<td valign="bottom" align="left">87.26 <sup>**</sup>
</td>
</tr>
<tr>
<td valign="bottom" align="left">V<sub>100%</sub> (%)</td>
<td valign="bottom" align="left">49.99</td>
<td valign="bottom" align="left">50.46</td>
<td valign="bottom" align="left">50.24</td>
</tr>
<tr>
<td valign="bottom" align="left">V<sub>105%</sub> (%)</td>
<td valign="bottom" align="left">0.89</td>
<td valign="bottom" align="left">0.41</td>
<td valign="bottom" align="left">10.95 <sup>***</sup>
</td>
</tr>
<tr>
<td valign="bottom" align="left">V<sub>107%</sub> (%)</td>
<td valign="bottom" align="left">0.02</td>
<td valign="bottom" align="left">0.01</td>
<td valign="bottom" align="left">4.20 <sup>**</sup>
</td>
</tr>
<tr>
<td valign="bottom" align="left">Radiation Conformity Index (RCI)</td>
<td valign="bottom" align="left">0.93</td>
<td valign="bottom" align="left">0.87</td>
<td valign="bottom" align="left">0.80</td>
</tr>
<tr>
<td valign="bottom" align="left">Homogeneity Index (HI)</td>
<td valign="bottom" align="left">0.10</td>
<td valign="bottom" align="left">0.09</td>
<td valign="bottom" align="left">0.16 <sup>***</sup>
</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>The Wilcoxon Rank Sum Test was performed using the dose parameters for each patient case and for each treatment technique with the median value being used as a representative value. FLASH-PT plans are significantly inferior in the majority of dose metrics compared to both IMPT techniques. <sup>&#x2217;&#x2217;&#x2217;</sup>p &lt; 0.001, <sup>&#x2217;&#x2217;</sup>p &lt; 0.01, <sup>&#x2217;</sup>p &lt; 0.05.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Bar charts representing the mean doses delivered to OARs for IMPT MFO, IMPT SFO and FLASH-PT, normalised with respect to the prescribed dose for each patient case. The considered OAR structures are the total brain (top) and the lung side in which the target is located (bottom). For the lung cases, the doses to the lung are shown for the 3-,4-, and 5-fraction plans.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1369065-g002.tif"/>
</fig>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Dose maps for an example brain (Case 7) and lung case (Case 10, 3 fractions). Comparisons of IMPT MFO, IMPT SFO, and FLASH-PT are made in the axial, sagittal, and coronal directions. Doses are normalised with respect to the prescribed dose.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1369065-g003.tif"/>
</fig>
<p>Spearman&#x2019;s Correlation Coefficient showed a very strong (p&#xa0;<italic>&lt;</italic>&#xa0;0.01) negative correlation (<italic>&#x3c1;<sub>SRC</sub>
</italic>) = -0.80) between the RCI and the target size in the IMPT SFO case. The use of SFO therefore appears to reduce conformity, however, this result was not found for the Bragg peak FLASH-PT plans, likely due to the inclusion of the aperture in the beam setup. Pearson&#x2019;s Correlation Coefficient showed a significant (p <italic>&lt;</italic> 0.05) positive correlation (<italic>&#x3c1;<sub>PCC</sub>
</italic>) = 0.63) between the HI and the target sizes for the FLASH-PT plans. An increase in target size is therefore related to a linear increase in HI, and thus worsened target homogeneity. This relationship was not found for IMPT MFO or IMPT SFO.</p>
<p>FLASH-PT plans appeared to be dependent on the target shape, position, smoothness of distal edge of the target, and distance from the nozzle of the beam to the target. Therefore, treatment plans generated using the same treatment planning parameters (spot spacings, layer spacings, and margins for spot placements) varied greatly with respect to the dose distributions and subsequently the plan quality achieved. To circumvent this, treatment planning parameters suitable for each patient case and each treatment fraction were used and obtained through manual parameter optimisation. Different treatment planning parameters were therefore required for each treatment fraction, with significantly larger spot spacings, layer spacings, and margins for spot placement needed to achieve similar mean doses and conformity to that of standard IMPT. Spot and layer spacings were at least twice as large for FLASH-PT plans compared to those used for IMPT. An increase in the spacings lead to a decrease in the number of spots placed on the target, making it more difficult for the TPS to generate and optimise conformal and homogeneous treatment plans (<xref ref-type="fig" rid="f4">
<bold>Figures&#xa0;4A, C</bold>
</xref>). Additionally, larger margins for spot placements improved the target coverage of the FLASH-PT plans, albeit at the expense of increasing the dose to nearby OARs. Details on the spot spacings, layer spacings, and margins of spot placements for each patient case and treatment fraction are shown in <xref ref-type="supplementary-material" rid="SM1">
<bold>Tables S1&#x2013;S3</bold>
</xref> in the <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Appendix</bold>
</xref>.</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Left, boxplots of the <bold>(A)</bold> radiation conformity index (RCI) and <bold>(C)</bold> homogeneity index (HI) evaluated for each patient case and each treatment technique. Right, bar charts of the <bold>(B)</bold> RCI and <bold>(D)</bold> HI for 3,4, and 5 fractions shown for Case 10. Each point (x) shows the exact RCI and HI for each patient case. The line within the box represents the median value, the bottom and top of the box represents the 25<italic>
<sup>th</sup>
</italic> and 75<italic>
<sup>th</sup>
</italic> percentiles, the whiskers show the range of the data points (excluding outliers), and outliers are shown using individual (o) markers.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1369065-g004.tif"/>
</fig>
<p>FLASH-PT plans were less conformal than both IMPT plans (<xref ref-type="fig" rid="f3">
<bold>Figures&#xa0;3</bold>
</xref>, <xref ref-type="fig" rid="f4">
<bold>4A</bold>
</xref>), although not significantly. Dose conformity for Bragg peak FLASH-PT can be improved by increasing the resolution at which the CEM is calculated, however, at the expense of drastically increasing the simulation time using this TPS. IMPT MFO and IMPT SFO plans produced comparable HI; this indicates that the inferior target homogeneity of FLASH-PT, seen in <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref> and <xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4C</bold>
</xref>, is mainly attributed to the limitations imposed for use of ultra-high dose rate protons. More specifically, the worsened target homogeneity is likely due to the increase in spot sizes and subsequent decrease in number of spots than can be placed, in addition to the impact of target size. Furthermore, as there are multiple objectives in the TPS, the optimiser appears to struggle to fulfil all criteria with fewer spot numbers than what is conventionally used for IMPT.</p>
<p>Case 10, having the largest target size and the largest number of OARs (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>), was selected as an example case to investigate the capabilities of the treatment planning optimiser and the subsequent treatment planning quality. The large target size, and the number of OARs, for this case led to inferior FLASH-PT treatment plans and so it is shown as a &#x201c;worse-case&#x201d; scenario in comparison to the rest of the samples. As shown in <xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4B</bold>
</xref>, the FLASH-PT RCI is comparable to both IMPT techniques, for all fractionation schemes. <xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4D</bold>
</xref>, however, shows a drastic reduction in target homogeneity for FLASH-PT when compared to IMPT (MFO and SFO). The 3-fraction FLASH-PT plan appeared to be the most inferior in terms of homogeneity; this is expected as reducing the number of beams will reduce the target homogeneity. This was also observed for the FLASH-PT plans generated for the other lung cases using 3 fractions. However, comparing DVHs shown in the <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Appendix (Figures S2&#x2013;S4)</bold>
</xref>, increasing the number of fractions does not generally improve target coverage or OAR doses.</p>
<p>For the same example case, variations in dose rates appeared to be dependent on the OARs as well as the number of fractions, and thus beams, used to deliver the treatment (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5</bold>
</xref>, <xref ref-type="fig" rid="f6">
<bold>6</bold>
</xref>). Comparing the brain and lung case shown in <xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>, the use of two fractions to treat a target with few density differences along the path of the beam (brain) generates higher dose rates to the surrounding normal tissues compared to the use of three fractions treating a target in a location where more density differences are present (lung). A further decrease in dose rates was also observed when increasing the number of fractions used for the lung cases.</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Dose rate maps for an example brain (Case 7) and lung case (Case 10, 3 fractions), showing the percentile dose rate in Gy/s. Comparisons of the dose rates achieved in each fraction are made in the axial, sagittal, and coronal directions.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1369065-g005.tif"/>
</fig>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Dose rate volume histograms of the Right Lung, Thoracic Wall, Great Vessel, and Heart structures for Case 10, separated for the 3 fractions planned. For simplicity, as no OAR receives FLASH dose rates to relative volumes &#x2265; 30% in this study, only relative volumes in the range of 0-30% are included in these plots. Solid lines represent the dose rate, different colours indicate the different fractions planned, and dashed lines represent the volume of the OAR which receives doses &#x2265; 2 Gy (below line). In the case of the Right Lung, the optimisation criterion was not fully achieved, as there are volumes receiving doses &#x2265; 2 Gy (dashed lines) that do not receive the dose at dose rates &#x2265; 40 Gy/s (solid lines); fraction 2 nearly fulfils it. Similar, for the Thoracic Wall, only fraction 2 fully achieves the dose rate optimisation criterion, while fractions 1 and 3 do not. The criterion is fully achieved for fractions 2 and 3 for the Great Vessel, and for all fractions for the Heart.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1369065-g006.tif"/>
</fig>
<p>Focusing more on the &#x201c;worst-case&#x201d; scenario (Case 10), the minimum dose rate constraint was not fully achieved for the Right Lung (for any fraction); this means that there were volumes receiving &#x2265; 2 Gy (dashed lines) that did not receive these doses at dose rates &#x2265; 40 Gy/s (solid lines), as shown in <xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>. This is likely due to low-dose, and therefore dose rate, regions at the edges of the beam, as seen in <xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref> and is a consequence of strict OAR dose constraints for this structure. Similarly for the Thoracic Wall, only fraction 2 fully achieved the dose rate optimisation criterion, while fractions 1 and 3 did not. However, the criterion was fully achieved for fractions 2 and 3 for the Great Vessel and for all fractions for the Heart. Larger OARs, such as the Right Lung and Thoracic Wall will have more low-dose rate regions due to the beam-path chosen and beam-spread compared to structures such as the Heart and Great Vessel, which due to their proximity to the target receive most of the dose at much higher dose rates.</p>
<p>Dose rate maps and dose rate volume histograms (DRVHs) give a good representation of the intra-fraction dose rates achieved for a patient case, however, to allow for comparisons between target types and fraction numbers, dose rate coverages were evaluated using DRVH dose rate metrics and the corresponding DVH dose metrics. Median values of dose rate coverages for the brain plans and 3-, 4-, and 5- fraction lung plans are 87%, 71%, 69%, and 56%, respectively. For the different lung fractionation schemes investigated, the 5-fraction plans had significantly (p <italic>&lt;</italic> 0.05) lower dose rates than the 3- and 4-fraction plans. This is likely due to the reduced intra-fraction dose making it harder for the optimiser to fulfil the ultra-high dose rate objective, as it paradoxically must increase the dose to achieve FLASH dose rates whilst also adhering to a maximum OAR dose constraint. Increasing the number of fractions lead to an increase in the spread of dose rate coverages achieved to OARs, as seen in <xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>. This was expected as the dose per fraction decreased.</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Boxplot showing the dose rate coverage of the OARs defined for the brain and lung cases for each fractionation scheme investigated. Each point shows the exact dose rate coverage for each OAR in each patient case, the line across the box represents the median value, the bottom and top of the box represent the 25<italic>
<sup>th</sup>
</italic> and 75<italic>
<sup>th</sup>
</italic> percentiles, the whiskers show the range of the data point (excluding outliers).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1369065-g007.tif"/>
</fig>
<p>Significantly higher dose rate coverages were achieved for the brain cases compared to the lung cases for all fractionation schemes investigated. This could be due to the increase in density differences along the path of the beam for the lung cases compared to the brain cases.</p>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>Due to the novelty of Bragg peak FLASH-PT, few feasibility studies focusing on the clinical implementation, and specifically the treatment planning process, have been conducted. The aim of this study was to determine this feasibility, and outline any constraints, of producing Bragg Peak FLASH-PT treatment plans delivered in one beam per fraction, and whether the generated plans were comparable to those for a clinically used technique such as IMPT. Without accounting for the sparing effect attributed to ultra-high dose rates, IMPT plans had better target coverage, dose conformity, target homogeneity, and OAR sparing compared to Bragg peak FLASH-PT using this beam setup and its implementation into the MIROpt TPS.</p>
<p>An interplay between the Bragg peak FLASH-PT beam setup, target location, fractionation scheme, and treatment planning optimisation is evident. Density differences along the path of the beam leads to variations in the treatment plans; target areas with few density differences, such as the brain, showed less variations whilst also achieving better target coverage, reduced doses to OARs, and improved dose rate coverage compared to the lung cases. In this study, brain cases appear to be more favourable than lung when using Bragg peak FLASH-PT, although the sample size is not large enough to ensure statistical significance.</p>
<p>The FLASH coverage for the brain and lung (3 fractions) cases, with OAR volumes receiving doses &#x2265; 2 Gy at dose rates &#x2265; 40 Gy/s, is comparable to the coverages of approximately 80% reported in previous treatment planning studies. The achieved FLASH coverage appears to be independent of patient case and target type when comparing brain and lung targets. In terms of fractionated treatment delivery, increasing the number of fractions for the lung cases did not result in any major differences in the target dose distributions or doses to OARs. A reduction in the intra-fraction dose rate coverage was found, as expected due to the intra-fraction dose being reduced. Fewer fractions will be more favourable if the FLASH sparing effect is to be fully exploited, however, ensuring clinically acceptable target coverage is still a priority.</p>
<p>The implementation of the CEM from the Conformal FLASH library to the MIROpt research TPS, the sample size, and the imposed target size constraint, are clear limitations in this study. However, in addition to being open-source and thereby very accessible, the use of the Conformal FLASH library in conjunction with MIROpt provided valuable information with regards to the treatment planning process of Bragg peak FLASH-PT. This includes the evaluation of fractionated beam delivery, dose comparisons between target types, the effect of target size on the planning quality, and the capability of achieving FLASH dose rates to the majority of the irradiated OAR volumes. Nevertheless, implementing the same beam setup into a more clinically used and validated system, such as RayStation<sup>&#xae;</sup> (RaySearch AB, Stockholm, Sweden), will allow for more efficient optimisation with additional degrees of freedom. This should improve the Bragg peak FLASH-PT plan quality. In addition to this, further studies should include a greater sample size, potentially focus on target types of varying sizes, and explore the use of hypo-fractionated regimens for this novel technique.</p>
</sec>
<sec id="s5" sec-type="conclusions">
<label>5</label>
<title>Conclusion</title>
<p>The aim of this study was to determine the feasibility and constraint of producing IMPT equivalent Bragg Peak FLASH-PT treatment plans using an open-source TPS. Without accounting for the FLASH sparing effect attributed to delivery at ultra-high dose rates, IMPT plans were superior in target coverage, dose conformity, target homogeneity, and OAR sparing compared to Bragg peak FLASH-PT using this beam setup. Establishing whether the limitations found are software or hardware related will be key in further developing, improving, and demonstrating the feasibility of clinically implementing Bragg peak FLASH-PT.</p>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s7" sec-type="ethics-statement">
<title>Ethics statement</title>
<p>EPN Lund 2013/742" allowed us the use of patient CT scans in radiotherapy research and development.</p>
</sec>
<sec id="s8" sec-type="author-contributions">
<title>Author contributions</title>
<p>NL: Writing &#x2013; original draft. KP: Writing &#x2013; review &amp; editing, Supervision, Funding acquisition. IF: Writing &#x2013; review &amp; editing, Resources.</p>
</sec>
</body>
<back>
<sec id="s9" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. NL&#x2019;s position is funded by the Medical Research Council (MRC) UK. The support by Cancer Research UK&#x2014;RadNet (C6078/A28736) is gratefully acknowledged. We thank the MRC for current financial support through a programme grant (MR/X006611/1). We also thank the National Cancer Institute/NIH/DHHS for their support towards this work (1P01CA257904).</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>Many thanks to Rudi Labarbe and Lucian Hotoiu for their help and advice related to the use of the Conformal FLASH library. The authors would also like to acknowledge the use of the University of Oxford Advanced Research Computing (ARC) facility in carrying out this work.</p>
</ack>
<sec id="s10" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s11" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s12" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fonc.2024.1369065/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fonc.2024.1369065/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="DataSheet_1.pdf" id="SM1" mimetype="application/pdf"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Montay-Gruel</surname> <given-names>P</given-names>
</name>
<name>
<surname>Petersson</surname> <given-names>K</given-names>
</name>
<name>
<surname>Jaccard</surname> <given-names>M</given-names>
</name>
<name>
<surname>Boivin</surname> <given-names>G</given-names>
</name>
<name>
<surname>Germond</surname> <given-names>JF</given-names>
</name>
<name>
<surname>Petit</surname> <given-names>B</given-names>
</name>
<etal/>
</person-group>. <article-title>Irradiation in a flash: Unique sparing of memory in mice after whole brain irradiation with dose rates above 100gy/s</article-title>. <source>Radiotherapy Oncol</source>. (<year>2017</year>) <volume>124</volume>:<page-range>365&#x2013;9</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.radonc.2017.05.003</pub-id>
</citation>
</ref>
<ref id="B2">
<label>2</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Montay-Gruel</surname> <given-names>P</given-names>
</name>
<name>
<surname>Acharya</surname> <given-names>MM</given-names>
</name>
<name>
<surname>Petersson</surname> <given-names>K</given-names>
</name>
<name>
<surname>Alikhani</surname> <given-names>L</given-names>
</name>
<name>
<surname>Allen</surname> <given-names>CYBD</given-names>
</name>
<name>
<surname>Ollivier</surname> <given-names>J</given-names>
</name>
<etal/>
</person-group>. <article-title>Long-term neurocognitive benefits of flash radiotherapy driven by reduced reactive oxygen species</article-title>. <source>Proc Natl Acad Sci USA</source>. (<year>2019</year>) <volume>116</volume>:<page-range>10943&#x2013;51</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.1901777116</pub-id>
</citation>
</ref>
<ref id="B3">
<label>3</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dokic</surname> <given-names>I</given-names>
</name>
<name>
<surname>Meister</surname> <given-names>S</given-names>
</name>
<name>
<surname>Bojcevski</surname> <given-names>J</given-names>
</name>
<name>
<surname>Tessonnier</surname> <given-names>T</given-names>
</name>
<name>
<surname>Walsh</surname> <given-names>D</given-names>
</name>
<name>
<surname>Knoll</surname> <given-names>M</given-names>
</name>
<etal/>
</person-group>. <article-title>Neuroprotective effects of ultra-high dose rate flash bragg peak proton irradiation</article-title>. <source>Int J Radiat Oncol Biol Phys</source>. (<year>2022</year>) <volume>113</volume>:<page-range>614&#x2013;23</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijrobp.2022.02.020</pub-id>
</citation>
</ref>
<ref id="B4">
<label>4</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alaghband</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Cheeks</surname> <given-names>SN</given-names>
</name>
<name>
<surname>Allen</surname> <given-names>BD</given-names>
</name>
<name>
<surname>Montay-Gruel</surname> <given-names>P</given-names>
</name>
<name>
<surname>Doan</surname> <given-names>NL</given-names>
</name>
<name>
<surname>Petit</surname> <given-names>B</given-names>
</name>
<etal/>
</person-group>. <article-title>Neuroprotection of radiosensitive juvenile mice by ultra-high dose rate flash irradiation</article-title>. <source>Cancers</source>. (<year>2020</year>) <volume>12</volume>:<fpage>1671</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/cancers12061671</pub-id>
</citation>
</ref>
<ref id="B5">
<label>5</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Diffenderfer</surname> <given-names>ES</given-names>
</name>
<name>
<surname>Verginadis</surname> <given-names>II</given-names>
</name>
<name>
<surname>Kim</surname> <given-names>MM</given-names>
</name>
<name>
<surname>Shoniyozov</surname> <given-names>K</given-names>
</name>
<name>
<surname>Velalopoulou</surname> <given-names>A</given-names>
</name>
<name>
<surname>Goia</surname> <given-names>D</given-names>
</name>
<etal/>
</person-group>. <article-title>Design, implementation, and <italic>in vivo</italic> validation of a novel proton flash radiation therapy system</article-title>. <source>Int J Radiat Oncol Biol Phys</source>. (<year>2020</year>) <volume>106</volume>:<page-range>440&#x2013;8</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijrobp.2019.10.049</pub-id>
</citation>
</ref>
<ref id="B6">
<label>6</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Adrian</surname> <given-names>G</given-names>
</name>
<name>
<surname>Konradsson</surname> <given-names>E</given-names>
</name>
<name>
<surname>Beyer</surname> <given-names>S</given-names>
</name>
<name>
<surname>Wittrup</surname> <given-names>A</given-names>
</name>
<name>
<surname>Butterworth</surname> <given-names>KT</given-names>
</name>
<name>
<surname>McMahon</surname> <given-names>S</given-names>
</name>
<etal/>
</person-group>. <article-title>Cancer cells can exhibit a sparing flash effect at low doses under normoxic <italic>in vitro</italic>-conditions</article-title>. <source>Front Oncol</source>. (<year>2021</year>) <volume>11</volume>:<elocation-id>686142</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fonc.2021.686142</pub-id>
</citation>
</ref>
<ref id="B7">
<label>7</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Velalopoulou</surname> <given-names>A</given-names>
</name>
<name>
<surname>Karagounis</surname> <given-names>IV</given-names>
</name>
<name>
<surname>Cramer</surname> <given-names>GM</given-names>
</name>
<name>
<surname>Kim</surname> <given-names>MM</given-names>
</name>
<name>
<surname>Skoufos</surname> <given-names>G</given-names>
</name>
<name>
<surname>Goia</surname> <given-names>D</given-names>
</name>
<etal/>
</person-group>. <article-title>Flash proton radiotherapy spares normal epithelial and mesenchymal tissues while preserving sarcoma response</article-title>. <source>Cancer Res</source>. (<year>2021</year>) <volume>81</volume>:<page-range>4808&#x2013;21</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1158/0008-5472.CAN-21-1500</pub-id>
</citation>
</ref>
<ref id="B8">
<label>8</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ruan</surname> <given-names>JL</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>C</given-names>
</name>
<name>
<surname>Wouters</surname> <given-names>S</given-names>
</name>
<name>
<surname>Tullis</surname> <given-names>I</given-names>
</name>
<name>
<surname>Verslegers</surname> <given-names>M</given-names>
</name>
<name>
<surname>Mysara</surname> <given-names>M</given-names>
</name>
<etal/>
</person-group>. <article-title>Irradiation at ultra-high (flash) dose rates reduces acute normal tissue toxicity in the mouse gastrointestinal system</article-title>. <source>Int J Radiat Oncology Biology Phys</source>. (<year>2021</year>) <volume>111</volume>:<page-range>1250&#x2013;61</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijrobp.2021.08.004</pub-id>
</citation>
</ref>
<ref id="B9">
<label>9</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Levy</surname> <given-names>K</given-names>
</name>
<name>
<surname>Natarajan</surname> <given-names>S</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J</given-names>
</name>
<name>
<surname>Chow</surname> <given-names>S</given-names>
</name>
<name>
<surname>Eggold</surname> <given-names>JT</given-names>
</name>
<name>
<surname>Loo</surname> <given-names>PE</given-names>
</name>
<etal/>
</person-group>. <article-title>Abdominal flash irradiation reduces radiation-induced gastrointestinal toxicity for the treatment of ovarian cancer in mice</article-title>. <source>Sci Rep</source>. (<year>2020</year>) <volume>10</volume>:<fpage>21600</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-020-78017-7</pub-id>
</citation>
</ref>
<ref id="B10">
<label>10</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>S&#xf8;rensen</surname> <given-names>BS</given-names>
</name>
<name>
<surname>Sitarz</surname> <given-names>MK</given-names>
</name>
<name>
<surname>Ankj&#xe6;rgaard</surname> <given-names>C</given-names>
</name>
<name>
<surname>Johansen</surname> <given-names>J</given-names>
</name>
<name>
<surname>Andersen</surname> <given-names>CE</given-names>
</name>
<name>
<surname>Kanouta</surname> <given-names>E</given-names>
</name>
<etal/>
</person-group>. <article-title>
<italic>In vivo</italic> validation and tissue sparing factor for acute damage of pencil beam scanning proton flash</article-title>. <source>Radiotherapy Oncol J Eur Soc Ther Radiol Oncol</source>. (<year>2022</year>) <volume>167</volume>:<page-range>109&#x2013;15</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.radonc.2021.12.022</pub-id>
</citation>
</ref>
<ref id="B11">
<label>11</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>S&#xf8;rensen</surname> <given-names>BS</given-names>
</name>
<name>
<surname>Sitarz</surname> <given-names>MK</given-names>
</name>
<name>
<surname>Ankj&#xe6;rgaard</surname> <given-names>C</given-names>
</name>
<name>
<surname>Johansen</surname> <given-names>JG</given-names>
</name>
<name>
<surname>Andersen</surname> <given-names>CE</given-names>
</name>
<name>
<surname>Kanouta</surname> <given-names>E</given-names>
</name>
<etal/>
</person-group>. <article-title>Pencil beam scanning proton flash maintains tumor control while normal tissue damage is reduced in a mouse model</article-title>. <source>Radiotherapy Oncol J Eur Soc Ther Radiol Oncol</source>. (<year>2022</year>) <volume>175</volume>:<page-range>178&#x2013;84</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.radonc.2022.05.014</pub-id>
</citation>
</ref>
<ref id="B12">
<label>12</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cao</surname> <given-names>X</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>R</given-names>
</name>
<name>
<surname>Esipova</surname> <given-names>TV</given-names>
</name>
<name>
<surname>Allu</surname> <given-names>SR</given-names>
</name>
<name>
<surname>Ashraf</surname> <given-names>R</given-names>
</name>
<name>
<surname>Rahman</surname> <given-names>M</given-names>
</name>
<etal/>
</person-group>. <article-title>Quantification of oxygen depletion during flash irradiation <italic>in vitro</italic> and <italic>in vivo</italic>
</article-title>. <source>Int J Radiat Oncol Biol Phys</source>. (<year>2021</year>) <volume>111</volume>:<page-range>240&#x2013;8</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijrobp.2021.03.056</pub-id>
</citation>
</ref>
<ref id="B13">
<label>13</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname> <given-names>MM</given-names>
</name>
<name>
<surname>Verginadis</surname> <given-names>II</given-names>
</name>
<name>
<surname>Goia</surname> <given-names>D</given-names>
</name>
<name>
<surname>Haertter</surname> <given-names>A</given-names>
</name>
<name>
<surname>Shoniyozov</surname> <given-names>K</given-names>
</name>
<name>
<surname>Zou</surname> <given-names>W</given-names>
</name>
<etal/>
</person-group>. <article-title>Comparison of flash proton entrance and the spread-out bragg peak dose regions in the sparing of mouse intestinal crypts and in a pancreatic tumor model</article-title>. <source>Cancers</source>. (<year>2021</year>) <volume>13</volume>:<fpage>4244</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/cancers13164244</pub-id>
</citation>
</ref>
<ref id="B14">
<label>14</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Iturri</surname> <given-names>L</given-names>
</name>
<name>
<surname>Bertho</surname> <given-names>A</given-names>
</name>
<name>
<surname>Lamirault</surname> <given-names>C</given-names>
</name>
<name>
<surname>Juchaux</surname> <given-names>M</given-names>
</name>
<name>
<surname>Gilbert</surname> <given-names>C</given-names>
</name>
<name>
<surname>Espenon</surname> <given-names>J</given-names>
</name>
<etal/>
</person-group>. <article-title>Proton flash radiation therapy and immune infiltration: Evaluation in an orthotopic glioma rat model</article-title>. <source>Int J Radiat Oncol Biol Phys</source>. (<year>2022</year>) <volume>116</volume>(<issue>3</issue>):<page-range>655&#x2013;65</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijrobp.2022.12.018</pub-id>
</citation>
</ref>
<ref id="B15">
<label>15</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname> <given-names>S</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>H</given-names>
</name>
<name>
<surname>Choi</surname> <given-names>JI</given-names>
</name>
<name>
<surname>Press</surname> <given-names>RH</given-names>
</name>
<name>
<surname>Lazarev</surname> <given-names>S</given-names>
</name>
<name>
<surname>Kabarriti</surname> <given-names>R</given-names>
</name>
<etal/>
</person-group>. <article-title>Flash radiotherapy using single-energy proton PBS transmission beams for hypofractionation liver cancer: Dose and dose rate quantification</article-title>. <source>Front Oncol</source>. (<year>2022</year>) <volume>11</volume>:<elocation-id>813063</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fonc.2021.813063</pub-id>
</citation>
</ref>
<ref id="B16">
<label>16</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>van de Water</surname> <given-names>S</given-names>
</name>
<name>
<surname>Safai</surname> <given-names>S</given-names>
</name>
<name>
<surname>Schippers</surname> <given-names>JM</given-names>
</name>
<name>
<surname>Weber</surname> <given-names>DC</given-names>
</name>
<name>
<surname>Lomax</surname> <given-names>AJ</given-names>
</name>
</person-group>. <article-title>Towards flash proton therapy: The impact of treatment planning and machine characteristics on achievable dose rates</article-title>. <source>Acta Oncol</source>. (<year>2019</year>) <volume>58</volume>:<page-range>1463&#x2013;9</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/0284186X.2019.1627416</pub-id>
</citation>
</ref>
<ref id="B17">
<label>17</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname> <given-names>S</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>H</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>S</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>C</given-names>
</name>
<name>
<surname>Xiong</surname> <given-names>W</given-names>
</name>
<name>
<surname>Zhai</surname> <given-names>H</given-names>
</name>
<etal/>
</person-group>. <article-title>Dose rate and dose robustness for proton transmission flash-rt treatment in lung cancer</article-title>. <source>Front Oncol</source>. (<year>2022</year>) <volume>12</volume>:<elocation-id>970602</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fonc.2022.970602</pub-id>
</citation>
</ref>
<ref id="B18">
<label>18</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kang</surname> <given-names>M</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>S</given-names>
</name>
<name>
<surname>Choi</surname> <given-names>JI</given-names>
</name>
<name>
<surname>Simone</surname> <given-names>CB</given-names>
<suffix>II</suffix>
</name>
<name>
<surname>Lin</surname> <given-names>H</given-names>
</name>
</person-group>. <article-title>Quantitative assessment of 3d dose rate for proton pencil beam scanning flash radiotherapy and its application for lung hypofractionation treatment planning</article-title>. <source>Cancers</source>. (<year>2021</year>) <volume>13</volume>:<fpage>3549</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/cancers13143549</pub-id>
</citation>
</ref>
<ref id="B19">
<label>19</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname> <given-names>S</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>H</given-names>
</name>
<name>
<surname>Choi</surname> <given-names>JI</given-names>
</name>
<name>
<surname>Simone</surname> <given-names>CB</given-names>
<suffix>II</suffix>
</name>
<name>
<surname>Kang</surname> <given-names>M</given-names>
</name>
</person-group>. <article-title>A novel proton pencil beam scanning FLASH RT delivery method enables optimal oar sparing and ultra-high dose rate delivery: A comprehensive dosimetry study for lung tumors</article-title>. <source>Cancers</source>. (<year>2021</year>) <volume>13</volume>:<fpage>5790</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/cancers13225790</pub-id>
</citation>
</ref>
<ref id="B20">
<label>20</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname> <given-names>S</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>H</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>C</given-names>
</name>
<name>
<surname>Xiong</surname> <given-names>W</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>CC</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>S</given-names>
</name>
<etal/>
</person-group>. <article-title>Use of single-energy proton pencil beam scanning bragg peak for intensity-modulated proton therapy flash treatment planning in liverhypofractionated radiation therapy</article-title>. <source>Med Phys</source>. (<year>2022</year>) <volume>49</volume>:<page-range>6560&#x2013;74</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/mp.15894</pub-id>
</citation>
</ref>
<ref id="B21">
<label>21</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>G</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>W</given-names>
</name>
<name>
<surname>Peng</surname> <given-names>H</given-names>
</name>
</person-group>. <article-title>Design of static and dynamic ridge filters for flash&#x2013;impt: A simulation study</article-title>. <source>Med Phys</source>. (<year>2022</year>) <volume>49</volume>:<page-range>5387&#x2013;99</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/mp.15717</pub-id>
</citation>
</ref>
<ref id="B22">
<label>22</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname> <given-names>C</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>X</given-names>
</name>
<name>
<surname>Chang</surname> <given-names>C</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>R</given-names>
</name>
<name>
<surname>Bohannon</surname> <given-names>D</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>L</given-names>
</name>
<etal/>
</person-group>. <article-title>Feasibility study of hybrid inverse planning with transmission beams and single-energy spread-out bragg peaks for proton flash radiotherapy</article-title>. <source>Med Phys</source>. (<year>2023</year>) <volume>50</volume>(<issue>6</issue>):<page-range>3687&#x2013;700</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/mp.16370</pub-id>
</citation>
</ref>
<ref id="B23">
<label>23</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tessonnier</surname> <given-names>T</given-names>
</name>
<name>
<surname>Mein</surname> <given-names>S</given-names>
</name>
<name>
<surname>Walsh</surname> <given-names>DWM</given-names>
</name>
<name>
<surname>Schuhmacher</surname> <given-names>N</given-names>
</name>
<name>
<surname>Liew</surname> <given-names>H</given-names>
</name>
<name>
<surname>Cee</surname> <given-names>R</given-names>
</name>
<etal/>
</person-group>. <article-title>Flash dose rate helium ion beams: First <italic>in vitro</italic> investigations</article-title>. <source>Int J Radiat Oncol Biol Phys</source>. (<year>2021</year>) <volume>111</volume>:<elocation-id>12345</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fnins.2013.12345</pub-id>
</citation>
</ref>
<ref id="B24">
<label>24</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Whitmore</surname> <given-names>L</given-names>
</name>
<name>
<surname>Mackay</surname> <given-names>RI</given-names>
</name>
<name>
<surname>van Herk</surname> <given-names>M</given-names>
</name>
<name>
<surname>Jones</surname> <given-names>JK</given-names>
</name>
<name>
<surname>Jones</surname> <given-names>RM</given-names>
</name>
</person-group>. <article-title>Focused vhee (very high energy electron) beams and dose delivery for radiotherapy applications</article-title>. <source>Sci Rep</source>. (<year>2021</year>) <volume>11</volume>:<fpage>14013</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-021-93276-8</pub-id>
</citation>
</ref>
<ref id="B25">
<label>25</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mascia</surname> <given-names>AE</given-names>
</name>
<name>
<surname>Daugherty</surname> <given-names>EC</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>E</given-names>
</name>
<name>
<surname>Xiao</surname> <given-names>Z</given-names>
</name>
<name>
<surname>Sertorio</surname> <given-names>M</given-names>
</name>
<etal/>
</person-group>. <article-title>Proton flash radiotherapy for the treatment of symptomatic bone metastases: The fast-01 nonrandomized trial</article-title>. <source>JAMA Oncol</source>. (<year>2023</year>) <volume>9</volume>(<issue>1</issue>):<page-range>62&#x2013;9</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1001/jamaoncol.2022.5843</pub-id>
</citation>
</ref>
<ref id="B26">
<label>26</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<collab>FAST-01 trial</collab>
</person-group>. <article-title>Feasibility study of flash radiotherapy for the treatment of symptomatic bone metastases (fast-01)</article-title>. (<year>2020</year>). ClinicalTrials.gov Identifier: NCT04592887.</citation>
</ref>
<ref id="B27">
<label>27</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<collab>FAST-02 trial</collab>
</person-group>. <article-title>Flash radiotherapy for the treatment of symptomatic bone metastases in the thorax (fast-02)</article-title>. (<year>2022</year>). ClinicalTrials.gov Identifier: NCT05524064.</citation>
</ref>
<ref id="B28">
<label>28</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Krieger</surname> <given-names>M</given-names>
</name>
<name>
<surname>van de Water</surname> <given-names>S</given-names>
</name>
<name>
<surname>Folkerts</surname> <given-names>MM</given-names>
</name>
<name>
<surname>Mazal</surname> <given-names>A</given-names>
</name>
<name>
<surname>Fabiano</surname> <given-names>S</given-names>
</name>
<name>
<surname>Bizzocchi</surname> <given-names>N</given-names>
</name>
<etal/>
</person-group>. <article-title>A quantitative flash effectiveness model to reveal potentials and pitfalls of high dose rate proton therapy</article-title>. <source>Med Phys</source>. (<year>2022</year>) <volume>49</volume>:<page-range>2026&#x2013;38</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/mp.15459</pub-id>
</citation>
</ref>
<ref id="B29">
<label>29</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kang</surname> <given-names>M</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>S</given-names>
</name>
<name>
<surname>Choi</surname> <given-names>JI</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>H</given-names>
</name>
<name>
<surname>Simone</surname> <given-names>CB</given-names>
<suffix>II</suffix>
</name>
</person-group>. <article-title>A universal range shifter and range compensator can enable proton pencil beam scanning single-energy bragg peak FLASH-RT treatment using current commercially available proton systems</article-title>. <source>Int J Radiat Oncol Biol Phys</source>. (<year>2022</year>) <volume>113</volume>:<page-range>203&#x2013;13</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijrobp.2022.01.009</pub-id>
</citation>
</ref>
<ref id="B30">
<label>30</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ramesh</surname> <given-names>P</given-names>
</name>
<name>
<surname>Gu</surname> <given-names>W</given-names>
</name>
<name>
<surname>Ruan</surname> <given-names>D</given-names>
</name>
<name>
<surname>Sheng</surname> <given-names>K</given-names>
</name>
</person-group>. <article-title>Dose and dose rate objectives in bragg peak and shoot-through beam orientation optimization for flash proton therapy</article-title>. <source>Med Phys</source>. (<year>2022</year>) <volume>49</volume>:<page-range>7826&#x2013;37</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/mp.16009</pub-id>
</citation>
</ref>
<ref id="B31">
<label>31</label>
<citation citation-type="other">
<person-group person-group-type="author">
<collab>MATLAB</collab>
</person-group>. (<year>2019</year>). <source>Version: 9.7.0 (R2019b)</source>. <publisher-loc>Natick, Massachusetts</publisher-loc>: <publisher-name>The MathWorks Inc.</publisher-name> Available at: https://uk.mathworks.com/products/matlab.html?s_tid=hp_products_matlab.</citation>
</ref>
<ref id="B32">
<label>32</label>
<citation citation-type="other">
<person-group person-group-type="author">
<collab>Mcsquare</collab>
</person-group>. <publisher-name>Universit&#xe9; catholique de Louvain</publisher-name>, <publisher-loc>Belgium</publisher-loc> (<year>2022</year>). Available at: <uri xlink:href="http://www.openmcsquare.org/">http://www.openmcsquare.org/</uri> (Accessed <access-date>24/02/2023</access-date>).</citation>
</ref>
<ref id="B33">
<label>33</label>
<citation citation-type="other">
<person-group person-group-type="author">
<collab>Openreggui</collab>
</person-group>. (<year>2022</year>). Available at: <uri xlink:href="https://openreggui.org/">https://openreggui.org/</uri> (Accessed <access-date>24/02/2023</access-date>).</citation>
</ref>
<ref id="B34">
<label>34</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<collab>Ion Beam Applications SA (IBA SA)</collab>
</person-group>. <article-title>ConformalFLASH</article-title>. (<year>2022</year>). Available at: <uri xlink:href="http://www.openflash.software">http://www.openflash.software</uri> (Accessed <access-date>24/02/2023</access-date>).</citation>
</ref>
<ref id="B35">
<label>35</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roeck</surname> <given-names>LD</given-names>
</name>
<name>
<surname>van der Weide</surname> <given-names>HL</given-names>
</name>
<name>
<surname>Eekers</surname> <given-names>DB</given-names>
</name>
<name>
<surname>Kramer</surname> <given-names>MC</given-names>
</name>
<name>
<surname>Alapetite</surname> <given-names>C</given-names>
</name>
<name>
<surname>Blomstrand</surname> <given-names>M</given-names>
</name>
<etal/>
</person-group>. <article-title>The European particle therapy network (EPTN) consensus on the follow-up of adult patients with brain and skull base tumours treated with photon or proton irradiation</article-title>. <source>Radiotherapy Oncol</source>. (<year>2022</year>) <volume>168</volume>:<page-range>241&#x2013;9</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.radonc.2022.01.018</pub-id>
</citation>
</ref>
<ref id="B36">
<label>36</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Marks</surname> <given-names>LB</given-names>
</name>
<name>
<surname>Yorke</surname> <given-names>ED</given-names>
</name>
</person-group>. <article-title>Use of normal tissue complication probability models in the clinic</article-title>. <source>Int J Radiat Oncology,Biology,Physics</source>. (<year>2010</year>) <volume>76</volume>:<page-range>S10&#x2013;9</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijrobp.2009.07.1754</pub-id>
</citation>
</ref>
<ref id="B37">
<label>37</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kn&#xf6;&#xf6;s</surname> <given-names>T</given-names>
</name>
<name>
<surname>Kristensen</surname> <given-names>I</given-names>
</name>
<name>
<surname>Nilsson</surname> <given-names>P</given-names>
</name>
</person-group>. <article-title>Volumetric and dosimetric evaluation of radiation treatment plans: radiation conformity index</article-title>. <source>Int J Radiat Oncol Biol Physics</source>. (<year>1998</year>) <volume>42</volume>:<page-range>1169&#x2013;76</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0360-3016(98)00239-9</pub-id>
</citation>
</ref>
<ref id="B38">
<label>38</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<collab>ICRU 83 report</collab>
</person-group>. <article-title>The international commission on radiation units and measurements</article-title>. <source>J ICRU</source>. (<year>2010</year>) <volume>10</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/jicru/ndq001</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>