{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,3]],"date-time":"2026-02-03T19:05:44Z","timestamp":1770145544368,"version":"3.49.0"},"reference-count":40,"publisher":"IOP Publishing","issue":"21","license":[{"start":{"date-parts":[[2020,11,13]],"date-time":"2020-11-13T00:00:00Z","timestamp":1605225600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/publishingsupport.iopscience.iop.org\/iop-standard\/v1"},{"start":{"date-parts":[[2020,11,13]],"date-time":"2020-11-13T00:00:00Z","timestamp":1605225600000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/iopscience.iop.org\/info\/page\/text-and-data-mining"}],"funder":[{"name":"Collaborative Health Research Projects","award":["523394-18"],"award-info":[{"award-number":["523394-18"]}]},{"DOI":"10.13039\/501100000038","name":"Natural Sciences and Engineering Research Council of Canada","doi-asserted-by":"crossref","award":["241018"],"award-info":[{"award-number":["241018"]}],"id":[{"id":"10.13039\/501100000038","id-type":"DOI","asserted-by":"crossref"}]}],"content-domain":{"domain":["iopscience.iop.org"],"crossmark-restriction":false},"short-container-title":["Phys. Med. Biol."],"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>\n                    The purpose of this work was to develop an efficient quadratic mixed integer programming algorithm for high dose rate (HDR) brachytherapy treatment planning problems and integrate the algorithm into an open-source Monte Carlo based treatment planning software, RapidBrachyMCTPS. The mixed-integer algorithm yields a globally optimum solution to the dose volume histogram (DVH) based problem and, unlike other methods, is not susceptible to local minimum trapping. A hybrid linear-quadratic penalty model coupled to a mixed integer programming model was used to optimize treatment plans for 10 prostate cancer patients. Dose distributions for each dwell position were calculated with RapidBrachyMCTPS with type A uncertainties less than 0.2% in voxels within the planning target volume (PTV). The optimization process was divided into two parts. First, the data was preprocessed, in which the problem size was reduced by eliminating voxels that had negligible impact on the solution (e.g. far from the dwell position). Second, the best combination of dwell times to obtain a plan with the highest score was found. The dwell positions and dose volume constraints were used as input to a commercial mixed integer optimizer (Gurobi Optimization, Inc.). A penalty-based criterion was adopted for the scoring. The voxel-reduction technique successfully reduced the problem size by an average of 91%, without loss of quality. The preprocessing of the optimization process required on average 4 s and solving for the global maximum required on average 33 s. The total optimization time averaged 37 s, which is a substantial improvement over the \u223c15 min optimization time reported in published literature. The plan quality was evaluated by evaluating dose volume metrics, including PTV\n                    <jats:italic>D<\/jats:italic>\n                    <jats:sub>90<\/jats:sub>\n                    , rectum and bladder\n                    <jats:italic>D<\/jats:italic>\n                    <jats:sub>\n                      1\n                      <jats:italic>cc<\/jats:italic>\n                    <\/jats:sub>\n                    and urethra\n                    <jats:italic>D<\/jats:italic>\n                    <jats:sub>\n                      0.1\n                      <jats:italic>cc<\/jats:italic>\n                    <\/jats:sub>\n                    . In conclusion, fast mixed integer optimization is an order of magnitude faster than current mixed-integer approaches for solving HDR brachytherapy treatment planning problems with DVH based metrics.\n                  <\/jats:p>","DOI":"10.1088\/1361-6560\/aba317","type":"journal-article","created":{"date-parts":[[2020,7,6]],"date-time":"2020-07-06T18:30:13Z","timestamp":1594060213000},"page":"215005","update-policy":"https:\/\/doi.org\/10.1088\/crossmark-policy","source":"Crossref","is-referenced-by-count":14,"title":["Fast mixed integer optimization (FMIO) for high dose rate brachytherapy"],"prefix":"10.1088","volume":"65","author":[{"given":"Majd","family":"Antaki","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Christopher","family":"L Deufel","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Shirin A","family":"Enger","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"266","published-online":{"date-parts":[[2020,11,13]]},"reference":[{"key":"pmbaba317bib1","doi-asserted-by":"publisher","first-page":"250","DOI":"10.1016\/S0168-9002(03)01368-8","type":"journal-article","article-title":"GEANT4-a simulation toolkit","volume":"506","author":"Agostinelli","year":"2003","journal-title":"Nucl. Instrum. Methods. Phys. Res. A"},{"key":"pmbaba317bib2","doi-asserted-by":"publisher","first-page":"270","DOI":"10.1109\/TNS.2006.869826","type":"journal-article","article-title":"Geant4 developments and applications","volume":"53","author":"Allison","year":"2006","journal-title":"IEEE Trans. Nucl. Sci."},{"key":"pmbaba317bib3","doi-asserted-by":"publisher","first-page":"4012","DOI":"10.1118\/1.2349685","type":"journal-article","article-title":"Optimization of HDR brachytherapy dose distributions using linear programming with penalty costs","volume":"33","author":"Alterovitz","year":"2006","journal-title":"Med. Phys."},{"key":"pmbaba317bib4","doi-asserted-by":"publisher","first-page":"283","DOI":"10.1007\/978-3-642-03474-9_81","type":"journal-article","article-title":"Influence of modulation restriction in inverse optimization with HIPO of prostate implants on plan quality: Analysis using dosimetric and radiobiological indices","author":"Baltas","year":"2009","journal-title":"World Congress on Medical Physics and Biomedical Engineering (Berlin: Springer)"},{"key":"pmbaba317bib5","doi-asserted-by":"publisher","first-page":"6208","DOI":"10.1118\/1.4747264","type":"journal-article","article-title":"Report of the Task Group 186 on model-based dose calculation methods in brachytherapy beyond the TG-43 formalism: Current status and recommendations for clinical implementation","volume":"39","author":"Beaulieu","year":"2012","journal-title":"Med. Phys."},{"key":"pmbaba317bib6","author":"Beli\u00ebn","year":"2009","type":"journal-article"},{"key":"pmbaba317bib7","doi-asserted-by":"publisher","first-page":"955","DOI":"10.1016\/j.ijrobp.2007.07.2321","type":"journal-article","article-title":"Inverse planning approach for 3-D MRI-based pulse-dose rate intracavitary brachytherapy in cervix cancer","volume":"69","author":"Chajon","year":"2007","journal-title":"Int. J. Radiat. Oncol. Biol. Phys."},{"key":"pmbaba317bib8","author":"Cullen","year":"1997","type":"journal-article"},{"key":"pmbaba317bib9","doi-asserted-by":"publisher","first-page":"10","DOI":"10.1016\/j.brachy.2004.11.002","type":"journal-article","article-title":"Ct-guided multi-catheter insertion technique for partial breast brachytherapy: reliable target coverage and dose homogeneity","volume":"4","author":"Cuttino","year":"2005","journal-title":"Brachytherapy"},{"key":"pmbaba317bib10","doi-asserted-by":"publisher","first-page":"2200","DOI":"10.1118\/1.598418","type":"journal-article","article-title":"Monte Carlo-aided dosimetry of a new high dose-rate brachytherapy source","volume":"25","author":"Daskalov","year":"1998","journal-title":"Med. Phys."},{"key":"pmbaba317bib11","doi-asserted-by":"publisher","first-page":"525","DOI":"10.1088\/0031-9155\/59\/3\/525","type":"journal-article","article-title":"Quality assurance for high dose rate brachytherapy treatment planning optimization: using a simple optimization to verify a complex optimization","volume":"59","author":"Deufel","year":"2014","journal-title":"Phys. Med. Biol."},{"key":"pmbaba317bib12","doi-asserted-by":"publisher","first-page":"71","DOI":"10.1016\/0360-3016(91)90168-4","type":"journal-article","article-title":"Dose-volume histograms","volume":"21","author":"Drzymala","year":"1991","journal-title":"Int. J. Radiat. Oncol. Biol. Phys."},{"key":"pmbaba317bib13","doi-asserted-by":"publisher","first-page":"297","DOI":"10.1088\/0031-9155\/46\/2\/303","type":"journal-article","article-title":"An iterative sequential mixed-integer approach to automated prostate brachytherapy treatment plan optimization","volume":"46","author":"D\u2019Souza","year":"2001","journal-title":"Phys. Med. Biol."},{"key":"pmbaba317bib14","doi-asserted-by":"publisher","first-page":"311","DOI":"10.1007\/s00291-007-0106-3","type":"journal-article","article-title":"Interactive decision support in radiation therapy treatment planning","volume":"30","author":"Ehrgott","year":"2008","journal-title":"OR Spectrum"},{"key":"pmbaba317bib15","doi-asserted-by":"publisher","DOI":"10.1088\/1361-6560\/aad97a","type":"journal-article","article-title":"RapidBrachyMCTPS: a Monte Carlo-based treatment planning system for brachytherapy applications","volume":"63","author":"Famulari","year":"2018","journal-title":"Phys. Med. Biol."},{"key":"pmbaba317bib16","first-page":"278","type":"journal-article","article-title":"Mixed integer programming optimization models for brachytherapy treatment planning","author":"Gallagher","year":"1997","journal-title":"Proc. AMIA Annu. Fall Symp."},{"key":"pmbaba317bib17","article-title":"RapidBrachyMCTPS 2.0: A comprehensive and flexible Monte Carlo-based treatment planning system for brachytherapy applications","author":"Glickman","year":"2020","type":"preprint"},{"key":"pmbaba317bib18","doi-asserted-by":"publisher","first-page":"1041","DOI":"10.1088\/0031-9155\/58\/4\/1041","type":"journal-article","article-title":"Mixed integer programming improves comprehensibility and plan quality in inverse optimization of prostate HDR brachytherapy","volume":"58","author":"Gorissen","year":"2013","journal-title":"Phys. Med. Biol."},{"key":"pmbaba317bib19","author":"Holm","year":"2013","journal-title":"Dissertation Department of Mathematics, Link\u00f6ping University","type":"journal-article"},{"key":"pmbaba317bib20","doi-asserted-by":"publisher","first-page":"1021","DOI":"10.1118\/1.3676179","type":"journal-article","article-title":"Impact of using linear optimization models in dose planning for HDR brachytherapy","volume":"39","author":"Holm","year":"2012","journal-title":"Med. Phys."},{"key":"pmbaba317bib21","doi-asserted-by":"publisher","first-page":"820","DOI":"10.1016\/j.ijrobp.2008.02.009","type":"journal-article","article-title":"Anatomy-based inverse planning simulated annealing optimization in high-dose-rate prostate brachytherapy: significant dosimetric advantage over other optimization techniques","volume":"72","author":"Jacob","year":"2008","journal-title":"Int. J. Radiat. Oncol. Biol. Phys."},{"key":"pmbaba317bib22","doi-asserted-by":"publisher","first-page":"612","DOI":"10.1007\/978-3-642-03474-9_172","type":"journal-article","article-title":"Optimization of catheter position and dwell time in prostate HDR brachytherapy using HIPO and linear programming","author":"Karabis","year":"2009","journal-title":"World Congress on Medical Physics and Biomedical Engineering September 7\u201312, 2009, Munich, Germany"},{"key":"pmbaba317bib23","doi-asserted-by":"publisher","first-page":"957","DOI":"10.1118\/1.597857","type":"journal-article","article-title":"A comparison of mixed integer programming and fast simulated annealing for optimizing beam weights in radiation therapy","volume":"23","author":"Langer","year":"1996","journal-title":"Med. Phys."},{"key":"pmbaba317bib24","doi-asserted-by":"publisher","first-page":"145","DOI":"10.1088\/0031-9155\/44\/1\/012","type":"journal-article","article-title":"Treatment planning for brachytherapy: an integer programming model, two computational approaches and experiments with permanent prostate implant planning","volume":"44","author":"Lee","year":"1999","journal-title":"Phys. Med. Biol."},{"key":"pmbaba317bib25","doi-asserted-by":"publisher","first-page":"147","DOI":"10.1023\/A:1022986523863","type":"journal-article","article-title":"Mixed integer programming approaches to treatment planning for brachytherapy\u2013application to permanent prostate implants","volume":"119","author":"Lee","year":"2003","journal-title":"Ann. Oper. Res."},{"key":"pmbaba317bib26","doi-asserted-by":"publisher","first-page":"773","DOI":"10.1118\/1.1368127","type":"journal-article","article-title":"Inverse planning anatomy-based dose optimization for HDR-brachytherapy of the prostate using fast simulated annealing algorithm and dedicated objective function","volume":"28","author":"Lessard","year":"2001","journal-title":"Med. Phys."},{"key":"pmbaba317bib27","doi-asserted-by":"publisher","first-page":"394","DOI":"10.1016\/j.radonc.2018.01.007","type":"journal-article","article-title":"The influence of tissue composition uncertainty on dose distributions in brachytherapy","volume":"126","author":"Mann-Krzisnik","year":"2018","journal-title":"Radiother. Oncol."},{"key":"pmbaba317bib28","doi-asserted-by":"publisher","first-page":"167","DOI":"10.1016\/j.radonc.2006.12.004","type":"journal-article","article-title":"Quality of interstitial PDR-brachytherapy-implants of head-and-neck-cancers: Predictive factors for local control and late toxicity?","volume":"82","author":"Melzner","year":"2007","journal-title":"Radiother. Oncol."},{"key":"pmbaba317bib29","doi-asserted-by":"publisher","first-page":"23","DOI":"10.1023\/A:1021386030224","type":"journal-article","article-title":"MIP models and BB strategies in brachytherapy treatment optimization","volume":"25","author":"Meyer","year":"2003","journal-title":"J Glob Optim."},{"key":"pmbaba317bib30","doi-asserted-by":"publisher","first-page":"2556","DOI":"10.1002\/mp.13533","type":"journal-article","article-title":"An extended dose\u2013volume model in high dose-rate brachytherapy\u2013Using mean-tail-dose to reduce tumor underdosage","volume":"46","author":"Mor\u00e9n","year":"2019","journal-title":"Med. Phys."},{"key":"pmbaba317bib31","type":"journal-article","article-title":"Tables and graphs of electron-interaction cross sections from 10 eV to 100 GeV derived from the LLNL evaluated data library (EEDL) Z\u2009=\u20091\u2009\u2212\u2009100","author":"Perkins","year":"1991","journal-title":"Technical Report No. UCRL-50400 vol 31 21\u20134 (Livermore, CA: Lawrence Livermore National Laboratory)"},{"key":"pmbaba317bib32","doi-asserted-by":"publisher","first-page":"71","DOI":"10.5114\/jcb.2010.14405","type":"journal-article","article-title":"Dose volume uniformity index: a simple tool for treatment plan evaluation in brachytherapy","volume":"2","author":"Prabhakar","year":"2010","journal-title":"J. Contemp. Brachytherapy"},{"key":"pmbaba317bib33","doi-asserted-by":"publisher","first-page":"e1\u2013e5","DOI":"10.1002\/mp.12702","type":"journal-article","article-title":"RECORDS: improved Reporting of montE CarlO RaDiation transport studies: report of the AAPM Research Committee Task Group 268","volume":"45","author":"Sechopoulos","year":"2018","journal-title":"Med. Phys."},{"key":"pmbaba317bib34","author":"Tuli","year":"2001","journal-title":"Technical Report No. BNL\u2013NCS-51655-01\/02-REV (Upton, NY: Brookhaven National Laboratory)","type":"journal-article"},{"key":"pmbaba317bib35","doi-asserted-by":"publisher","first-page":"170","DOI":"10.1016\/j.ijrobp.2004.09.021","type":"journal-article","article-title":"High-dose regions versus likelihood of cure after prostate brachytherapy","volume":"62","author":"Wallner","year":"2005","journal-title":"Int. J. Radiat. Oncol. Biol. Phys."},{"key":"pmbaba317bib36","type":"journal-article","article-title":"DOSXYZnrc Users Manual","author":"Walters","year":"2005","journal-title":"Ionizing Radiation Standards NRCC Report PIRS-794revB (Ottawa: National Research Council of Canada)"},{"key":"pmbaba317bib37","doi-asserted-by":"publisher","first-page":"489","DOI":"10.1016\/j.ijrobp.2005.06.028","type":"journal-article","article-title":"Accelerated partial breast irradiation: an analysis of variables associated with late toxicity and long-term cosmetic outcome after high-dose-rate interstitial brachytherapy","volume":"64","author":"Wazer","year":"2006","journal-title":"Int. J. Radiat. Oncol. Biol. Phys."},{"key":"pmbaba317bib38","doi-asserted-by":"publisher","first-page":"567","DOI":"10.1118\/1.596069","type":"journal-article","article-title":"Monte Carlo evaluation of kerma at a point for photon transport problems","volume":"14","author":"Williamson","year":"1987","journal-title":"Med. Phys."},{"key":"pmbaba317bib39","doi-asserted-by":"publisher","first-page":"28","DOI":"10.1016\/j.ijrobp.2007.09.011","type":"journal-article","article-title":"Postimplantation analysis enables improvement of dose\u2013volume histograms and reduction of toxicity for permanent seed implantation","volume":"71","author":"Wust","year":"2008","journal-title":"Int. J. Radiat. Oncol. Biol. Phys."},{"key":"pmbaba317bib40","doi-asserted-by":"publisher","first-page":"20","DOI":"10.1016\/j.brachy.2011.09.008","type":"journal-article","article-title":"American Brachytherapy Society consensus guidelines for high-dose-rate prostate brachytherapy","volume":"11","author":"Yamada","year":"2012","journal-title":"Brachytherapy"}],"container-title":["Physics in Medicine &amp; Biology"],"original-title":[],"link":[{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317","content-type":"text\/html","content-version":"am","intended-application":"text-mining"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317\/pdf","content-type":"application\/pdf","content-version":"am","intended-application":"text-mining"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317\/pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317\/pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317\/pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317\/pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317\/pdf","content-type":"application\/pdf","content-version":"am","intended-application":"similarity-checking"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317\/pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317\/pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"},{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317","content-type":"text\/html","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,12,13]],"date-time":"2025-12-13T14:43:03Z","timestamp":1765636983000},"score":1,"resource":{"primary":{"URL":"https:\/\/iopscience.iop.org\/article\/10.1088\/1361-6560\/aba317"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2020,11,13]]},"references-count":40,"journal-issue":{"issue":"21","published-online":{"date-parts":[[2020,11,13]]}},"URL":"https:\/\/doi.org\/10.1088\/1361-6560\/aba317","relation":{},"ISSN":["0031-9155","1361-6560"],"issn-type":[{"value":"0031-9155","type":"print"},{"value":"1361-6560","type":"electronic"}],"subject":[],"published":{"date-parts":[[2020,11,13]]},"assertion":[{"value":"Fast mixed integer optimization (FMIO) for high dose rate brachytherapy","name":"article_title","label":"Article Title"},{"value":"Physics in Medicine & Biology","name":"journal_title","label":"Journal Title"},{"value":"paper","name":"article_type","label":"Article Type"},{"value":"\u00a9 2020 Institute of Physics and Engineering in Medicine. All rights, including for text and data mining, AI training, and similar technologies, are reserved.","name":"copyright_information","label":"Copyright Information"},{"value":"2020-04-17","name":"date_received","label":"Date Received","group":{"name":"publication_dates","label":"Publication dates"}},{"value":"2020-07-06","name":"date_accepted","label":"Date Accepted","group":{"name":"publication_dates","label":"Publication dates"}},{"value":"2020-11-13","name":"date_epub","label":"Online publication date","group":{"name":"publication_dates","label":"Publication dates"}}]}}