In hospitals, the nurse schedules are sensitive to disruptions such as sick leave since the absence of nurses with the right skills can have a severe impact on healthcare quality. We consider nurse rerostering from a strategic perspective in the case of having multiple skills and varying staffing demand. Our aim is to design a decision support tool to be used on a strategic level to analyse how the distribution of skills affects the sick-leave robustness of schedules. Schedules are evaluated based on generated scenarios with sick leave and possible rerostering of tasks. We have performed a case study at a ward in a Swedish hospital, and show how our tool can be used to indicate which of the available skills can be removed, with only a small impact on the understaffing in the nominal scenario and in scenarios with sick leave.

BackgroundIntensity-modulated brachytherapy (IMBT) is an emerging technology for cancer treatment, in which radiation sources are shielded to shape the dose distribution. The rotatable shields provide an additional degree of freedom, but also introduce an additional, directional, type of uncertainty, compared to conventional high-dose-rate brachytherapy (HDR BT). PurposeWe propose and evaluate a robust optimization approach to mitigate the effects of rotational uncertainty in the shields with respect to planning criteria. MethodsA previously suggested prototype for platinum-shielded prostate Yb-169-based dynamic IMBT is considered. We study a retrospective patient data set (anatomical contours and catheter placement) from two clinics, consisting of six patients that had previously undergone conventional Ir-192 HDR BT treatment. The Monte Carlo-based treatment planning software RapidBrachyMCTPS is used for dose calculations. In our computational experiments, we investigate systematic rotational shield errors of +/- 10 degrees and +/- 20 degrees, and the same systematic error is applied to all dwell positions in each scenario. This gives us three scenarios, one nominal and two with errors. The robust optimization approach finds a compromise between the average and worst-case scenario outcomes. ResultsWe compare dose plans obtained from standard models and their robust counterparts. With dwell times obtained from a linear penalty model (LPM), for 10 degrees errors, the dose to urethra (D0.1cc) and rectum (D0.1cc and D1cc) increase with up to 5% and 7%, respectively, in the worst-case scenario, while with the robust counterpart, the corresponding increases were 3% and 3%. For all patients and all evaluated criteria, the worst-case scenario outcome with the robust approach had lower deviation compared to the standard model, without compromising target coverage. We also evaluated shield errors up to 20 degrees and while the deviations increased to a large extent with the standard models, the robust models were capable of handling even such large errors. ConclusionsWe conclude that robust optimization can be used to mitigate the effects from rotational uncertainty and to ensure the treatment plan quality of IMBT.

BackgroundSpatial properties of a dose distribution, such as volumes of contiguous hot spots, are of clinical importance in treatment planning for high dose-rate brachytherapy (HDR BT). We have in an earlier study developed an optimization model that reduces the prevalence of contiguous hot spots by modifying a tentative treatment plan. PurposeThe aim of this study is to incorporate the correction of hot spots in a standard inverse planning workflow and to validate the integrated model in a clinical treatment planning system. The spatial function is included in the objective function for the inverse planning, as opposed to in the previous study where it was applied as a separate post-processing step. Our aim is to demonstrate that fine-adjustments of dose distributions, which are often performed manually in todays clinical practice, can be automated. MethodsA spatial optimization function was introduced in the treatment planning system RayStation (RaySearch Laboratories AB, Stockholm, Sweden) via a research interface. A series of 10 consecutive prostate patients treated with HDR BT was retrospectively replanned with and without the spatial function. ResultsOptimization with the spatial function decreased the volume of the largest contiguous hot spot by on average 31%, compared to if the function was not included. The volume receiving at least 200% of the prescription dose decreased by on average 11%. Target coverage, measured as the fractions of the clinical target volume (CTV) and the planning target volume (PTV) receiving at least the prescription dose, was virtually unchanged (less than a percent change for both metrics). Organs-at-risk received comparable or slightly decreased doses if the spatial function was included in the optimization model. ConclusionsOptimization of spatial properties such as the volume of contiguous hot spots can be integrated in a standard inverse planning workflow for brachytherapy, and need not be conducted as a separate post-processing step.

PURPOSE: The aim was to evaluate a postprocessing optimization algorithms ability to improve the spatial properties of a clinical treatment plan while preserving the target coverage and the dose to the organs at risk. The goal was to obtain a more homogenous treatment plan, minimizing the need for manual adjustments after inverse treatment planning. MATERIALS AND METHODS: The study included 25 previously treated prostate cancer pa-tients. The treatment plans were evaluated on dose-volume histogram parameters established clin-ical and quantitative measures of the high dose volumes. The volumes of the four largest hot spots were compared and complemented with a human observer study with visual grading by eight oncologists. Statistical analysis was done using ordinal logistic regression. Weighted kappa and Fleiss kappa were used to evaluate intra-and interobserver reliability. RESULTS: The quantitative analysis showed that there was no change in planning target volume (PTV) coverage and dose to the rectum. There were significant improvements for the adjusted treatment plan in: V150% and V200% for PTV, dose to urethra, conformal index, and dose nonhomogeneity ratio. The three largest hot spots for the adjusted treatment plan were significantly smaller compared to the clinical treatment plan. The observers preferred the adjusted treatment plan in 132 cases and the clinical in 83 cases. The observers preferred the adjusted treatment plan on homogeneity and organs at risk but preferred the clinical plan on PTV coverage. CONCLUSIONS: Quantitative analysis showed that the postadjustment optimization tool could improve the spatial properties of the treatment plans while maintaining the target coverage. (c) 2022 The Authors. Published by Elsevier Inc. on behalf of American Brachytherapy Society. This is an open access article under the CC BY license ( http://creativecommons.org/licenses/by/4.0/ )

Brachytherapy is a mature treatment modality. The literature is abundant in terms of review articles and comprehensive books on the latest established as well as evolving clinical practices. The intent of this article is to part ways and look beyond the current state-of-the-art and review emerging technologies that are noteworthy and perhaps may drive the future innovations in the field. There are plenty of candidate topics that deserve a deeper look, of course, but with practical limits in this communicative platform, we explore four topics that perhaps is worthwhile to review in detail at this time. First, intensity modulated brachytherapy (IMBT) is reviewed. The IMBT takes advantage of anisotropic radiation profile generated through intelligent high-density shielding designs incorporated onto sources and applicators such to achieve high quality plans. Second, emerging applications of 3D printing (i.e. additive manufacturing) in brachytherapy are reviewed. With the advent of 3D printing, interest in this technology in brachytherapy has been immense and translation swift due to their potential to tailor applicators and treatments customizable to each individual patient. This is followed by, in third, innovations in treatment planning concerning catheter placement and dwell times where new modelling approaches, solution algorithms, and technological advances are reviewed. And, fourth and lastly, applications of a new machine learning technique, called deep learning, which has the potential to improve and automate all aspects of brachytherapy workflow, are reviewed. We do not expect that all ideas and innovations reviewed in this article will ultimately reach clinic but, nonetheless, this review provides a decent glimpse of what is to come. It would be exciting to monitor as IMBT, 3D printing, novel optimization algorithms, and deep learning technologies evolve over time and translate into pilot testing and sensibly phased clinical trials, and ultimately make a difference for cancer patients. Todays fancy is tomorrows reality. The future is bright for brachytherapy.

Treatment planning in high dose‐rate brachytherapy has traditionally been conducted with manual forward planning, but inverse planning is today increasingly used in clinical practice. There is a large variety of proposed optimization models and algorithms to model and solve the treatment planning problem. Two major parts of inverse treatment planning for which mathematical optimization can be used are the decisions about catheter placement and dwell time distributions. Both these problems as well as integrated approaches are included in this review. The proposed models include linear penalty models, dose–volume models, mean‐tail dose models, quadratic penalty models, radiobiological models, and multiobjective models. The aim of this survey is twofold: (i) to give a broad overview over mathematical optimization models used for treatment planning of brachytherapy and (ii) to provide mathematical analyses and comparisons between models. New technologies for brachytherapy treatments and methods for treatment planning are also discussed. Of particular interest for future research is a thorough comparison between optimization models and algorithms on the same dataset, and clinical validation of proposed optimization approaches with respect to patient outcome.

Cancer is a widespread class of diseases that each year affects millions of people. It is mostly treated with chemotherapy, surgery, radiation therapy, or combinations thereof. High doserate (HDR) brachytherapy (BT) is one modality of radiation therapy, which is used to treat for example prostate cancer and gynecologic cancer. In BT, catheters (i.e., hollow needles) or applicators are used to place a single, small, but highly radioactive source of ionizing radiation close to or within a tumour, at dwell positions. An emerging technique for HDR BT treatment is intensity modulated brachytherapy (IMBT), in which static or dynamic shields are used to further shape the dose distribution, by hindering the radiation in certain directions. 

The topic of this thesis is the application of mathematical optimization to model and solve the treatment planning problem. The treatment planning includes decisions on catheter placement, that is, how many catheters to use and where to place them, as well as decisions for dwell times. Our focus is on the latter decisions. The primary treatment goals are to give the tumour a sufficiently high radiation dose while limiting the dose to the surrounding healthy organs, to avoid severe side effects. Because these aims are typically in conflict, optimization models of the treatment planning problem are inherently multiobjective. Compared to manual treatment planning, there are several advantages of using mathematical optimization for treatment planning. First, the optimization of treatment plans requires less time, compared to the time-consuming manual planning. Secondly, treatment plan quality can be improved by using optimization models and algorithms. Finally, with the use of sophisticated optimization models and algorithms the requirements of experience and skill level for the planners are lower. The use of optimization for treatment planning of IMBT is especially important because the degrees of freedom are too many for manual planning. 

The contributions of this thesis include the study of properties of treatment planning models, suggestions for extensions and improvements of proposed models, and the development of new optimization models that take clinically relevant, but uncustomary aspects, into account in the treatment planning. A common theme is the modelling of constraints on dosimetric indices, each of which is a restriction on the portion of a volume that receives at least a specified dose, or on the lowest dose that is received by a portion of a volume. Modelling dosimetric indices explicitly yields mixed-integer programs which are computationally demanding to solve. We have therefore investigated approximations of dosimetric indices, for example using smooth non-linear functions or convex functions. Contributions of this thesis are also a literature review of proposed treatment planning models for HDR BT, including mathematical analyses and comparisons of models, and a study of treatment planning for IMBT, which shows how robust optimization can be used to mitigate the risks from rotational errors in the shield placement. 

Cancer är en grupp av sjukdomar som varje år drabbar miljontals människor. De vanligaste behandlingsformerna är cellgifter, kirurgi, strålbehandling eller en kombination av dessa. I denna avhandling studeras högdosrat brachyterapi (HDR BT), vilket är en form av strålbehandling som till exempel används vid behandling av prostatacancer och gynekologisk cancer. Vid brachyterapibehandling används ihåliga nålar eller applikatorer för att placera en millimeterstor strålkälla antingen inuti eller intill en tumör. I varje nål finns det ett antal så kallade dröjpositioner där strålkällan kan stanna en viss tid för att bestråla den omkringliggande vävnaden, i alla riktningar. Genom att välja lämpliga tider för dröjpositionerna kan dosfördelningen formas efter patientens anatomi. Utöver HDR BT studeras också den nya tekniken intensitetsmodulerad brachyterapi (IMBT) vilket är en variation på HDR BT där skärmning används för att minska strålningen i vissa riktningar vilket gör det möjligt att forma dosfördelningen bättre. 

Planeringen av en behandling med HDR BT omfattar hur många nålar som ska användas, var de ska placeras samt hur länge strålkällan ska stanna i de olika dröjpositionerna. För HDR BT kan dessa vara flera hundra stycken medan det för IMBT snarare handlar om tusentals möjliga kombinationer av dröjpositioner och inställningar av skärmarna. Planeringen resulterar i en dosplan som beskriver hur hög stråldos som tumören och intilliggande frisk vävnad och riskorgan utsätts för. Dosplaneringen kan formuleras som ett matematiskt optimeringsproblem vilket är ämnet för avhandlingen. De övergripande målsättningarna för behandlingen är att ge en tillräckligt hög stråldos till tumören, för att döda alla cancerceller, samt att undvika att bestråla riskorgan eftersom det kan ge allvarliga biverkningar. Då alla målsättningarna inte samtidigt kan uppnås fullt ut så fås optimeringsproblem där flera målsättningar behöver prioriteras mot varandra. Utöver att dosplanen uppfyller kliniska behandlingsriktlinjer så är också tidsaspekten av planeringen viktig eftersom det är vanligt att den görs medan patienten är bedövad eller sövd. 

Vid utvärdering av en dosplan används dos-volymmått. För en tumör anger ett dosvolymmått hur stor andel av tumören som får en stråldos som är högre än en specificerad nivå. Dos-volymmått utgör en viktig del av målen för dosplaner som tas upp i kliniska behandlingsriktlinjer och ett exempel på ett sådant mål vid behandling av prostatacancer är att 95% av prostatans volym ska få en stråldos som är minst den föreskrivna dosen. Dos-volymmått utläses ur de kliniskt betydelsefulla dos-volym histogrammen som för varje stråldosnivå anger motsvarande volym som erhåller den dosen. 

En fördel med att använda matematisk optimering för dosplanering är att det kan spara tid jämfört med manuell planering. Med väl utvecklade modeller så finns det också möjlighet att skapa bättre dosplaner, till exempel genom att riskorganen nås av en lägre dos men med bibehållen dos till tumören. Vidare så finns det även fördelar med en process som inte är lika personberoende och som inte kräver erfarenhet i lika stor utsträckning som manuell dosplanering i dagsläget gör. Vid IMBT är det dessutom så många frihetsgrader att manuell planering i stort sett blir omöjligt. 

I avhandlingen ligger fokus på hur dos-volymmått kan användas och modelleras explicit i optimeringsmodeller, så kallade dos-volymmodeller. Detta omfattar såväl analys av egenskaper hos befintliga modeller, utvidgningar av tidigare använda modeller samt utveckling av nya optimeringsmodeller. Eftersom dos-volymmodeller modelleras som heltalsproblem, vilka är beräkningskrävande att lösa, så är det också viktigt att utveckla algoritmer som kan lösa dem tillräckligt snabbt för klinisk användning. Ett annat mål för modellutvecklingen är att kunna ta hänsyn till fler kriterier som är kliniskt relevanta men som inte ingår i dos-volymmodeller. En sådan kategori av mått är hur dosen är fördelad rumsligt, exempelvis att volymen av sammanhängande områden som får en alldeles för hög dos ska vara liten. Sådana områden går dock inte att undvika helt eftersom det är typiskt för dosplaner för brachyterapi att stråldosen fördelar sig ojämnt, med väldigt höga doser till små volymer precis intill strålkällorna. Vidare studeras hur små fel i inställningarna av skärmningen i IMBT påverkar dosplanens kvalitet och de olika utvärderingsmått som används kliniskt. Robust optimering har använts för att säkerställa att en dosplan tas fram som är robust sett till dessa möjliga fel i hur skärmningen är placerad. 

Slutligen ges en omfattande översikt över optimeringsmodeller för dosplanering av HDR BT och speciellt hur optimeringsmodellerna hanterar de motstridiga målsättningarna.  

High dose-rate brachytherapy is a modality of radiation therapy used for cancer treatment, in which the radiation source is placed within the body. The treatment goal is to give a high enough dose to the tumour while sparing nearby healthy tissue and organs (organs-at-risk). The most common criteria for evaluating dose distributions are dosimetric indices. For the tumour, such an index is the portion of the volume that receives at least a specified dose level (e.g. the prescription dose), while for organs-at-risk it is instead the portion of the volume that receives at most a specified dose level. Dosimetric indices are aggregate criteria and do not consider spatial properties of the dose distribution. Further, there are neither any established evaluation criteria for characterizing spatial properties, nor have such properties been studied in the context of mathematical optimization of brachytherapy. Spatial properties are however of clinical relevance and therefore dose plans are sometimes adjusted manually to improve them. We propose an optimization model for reducing the prevalence of contiguous volumes with a too high dose (hot spots) or a too low dose (cold spots) in a tentative dose plan. This model is independent of the process of constructing the tentative plan. We conduct computational experiments with tentative plans obtained both from optimization models and from clinical practice. The objective function considers pairs of dose points and each pair is given a distance-based penalty if the dose is either too high or too low at both dose points. Constraints are included to retain dosimetric indices at acceptable levels. Our model is designed to automate the manual adjustment step in the planning process. In the automatic adjustment step large-scale optimization models are solved. We show reductions of the volumes of the largest hot and cold spots, and the computing times are feasible in clinical practice.

Purpose High dose-rate brachytherapy is a method of radiotherapy for cancer treatment in which the radiation source is placed within the body. In addition to give a high enough dose to a tumor, it is also important to spare nearby healthy organs [organs at risk (OAR)]. Dose plans are commonly evaluated using the so-called dosimetric indices; for the tumor, the portion of the structure that receives a sufficiently high dose is calculated, while for OAR it is instead the portion of the structure that receives a sufficiently low dose that is of interest. Models that include dosimetric indices are referred to as dose-volume models (DVMs) and have received much interest recently. Such models do not take the dose to the coldest (least irradiated) volume of the tumor into account, which is a distinct weakness since research indicates that the treatment effect can be largely impaired by tumor underdosage even to small volumes. Therefore, our aim is to extend a DVM to also consider the dose to the coldest volume. Methods An improved DVM for dose planning is proposed. In addition to optimizing with respect to dosimetric indices, this model also takes mean dose to the coldest volume of the tumor into account. Results Our extended model has been evaluated against a standard DVM in ten prostate geometries. Our results show that the dose to the coldest volume could be increased, while also computing times for the dose planning were improved. Conclusion While the proposed model yields dose plans similar to other models in most aspects, it fulfils its purpose of increasing the dose to cold tumor volumes. An additional benefit is shorter solution times, and especially for clinically relevant times (of minutes) we show major improvements in tumour dosimetric indices.