Rydberg-Blockade-Based Parity Quantum Optimization.

We present a scalable architecture for solving higher-order constrained binary optimization (HCBO) problems on current neutral-atom hardware operating in the Rydberg blockade regime. In particular, we formulate the recently developed parity encoding of arbitrary connected HCBO problems as a maximum-weight independent set (MWIS) problem on disk graphs, that are directly encodable on such devices. Our architecture builds from small MWIS modules in a problem-independent way, crucial for practical scalability.

Universal Parity Quantum Computing.

We propose a universal gate set for quantum computing with all-to-all connectivity and intrinsic robustness to bit-flip errors based on parity encoding. We show that logical controlled phase gate and R_{z} rotations can be implemented in parity encoding with single-qubit operations. Together with logical R_{x} rotations, implemented via nearest-neighbor controlled-NOT gates and an R_{x} rotation, these form a universal gate set. As the controlled phase gate requires only single-qubit rotations, the proposed scheme has advantages for several cornerstone quantum algorithms, e.g., the quantum Fourier transform. We present a method to switch between different encoding variants via partial on-the-fly encoding and decoding.

Quantum Optimization via Four-Body Rydberg Gates.

A large ongoing research effort focuses on obtaining a quantum advantage in the solution of combinatorial optimization problems on near-term quantum devices. A particularly promising platform implementing quantum optimization algorithms are arrays of trapped neutral atoms, laser coupled to highly excited Rydberg states. However, encoding combinatorial optimization problems in atomic arrays is challenging due to limited interqubit connectivity of the native finite-range interactions. Here, we present a four-body Rydberg parity gate, enabling a direct and straightforward implementation of the parity architecture, a scalable architecture for encoding arbitrarily connected interaction graphs. Our gate relies on adiabatic laser pulses and is fully programmable by adjusting two hold times during operation. We numerically demonstrate implementations of the quantum approximate optimization algorithm (QAOA) for small-scale test problems. Variational optimization steps can be implemented with a constant number of system manipulations, paving the way for experimental investigations of QAOA beyond the reach of numerical simulations.

Perspectives of quantum annealing: methods and implementations.

Quantum annealing is a computing paradigm that has the ambitious goal of efficiently solving large-scale combinatorial optimization problems of practical importance. However, many challenges have yet to be overcome before this goal can be reached. This perspectives article first gives a brief introduction to the concept of quantum annealing, and then highlights new pathways that may clear the way towards feasible and large scale quantum annealing. Moreover, since this field of research is to a strong degree driven by a synergy between experiment and theory, we discuss both in this work. An important focus in this article is on future perspectives, which complements other review articles, and which we hope will motivate further research.

The influence of errors in small field dosimetry on the dosimetric accuracy of treatment plans.

Background: Dosimetric effects of inaccuracies of output factors (OFs) implemented in treatment planning systems (TPSs) were investigated.Materials and methods: Modified beam models (MBM) for which the OFs of small fields (down to 1 × 1 cm2) were increased by up to 12% compared to the original beam models (OBM) were created for two TPSs. These beam models were used to recalculate treatment plans of different complexity. Treatment plans using stereotactic 3D-conformal (s3D-CRT) for brain metastasis as well as VMAT plans for head and neck and prostate cancer patients were generated. Dose distributions calculated with the MBM and the OBM were compared to measured dose distributions acquired using film dosimetry and a 2D-detector-array. For the s3D-CRT plans the calculated and measured dose at the isocenter was evaluated. For VMAT, gamma pass rates (GPRs) were calculated using global gamma index with 3%/3 mm, 2%/3 mm, 1%/3 mm and 2%/2 mm with a 20% threshold. Contribution of small fields to the total fluence was expressed as the ratio (F) of fluence trough leaf openings smaller than 2 cm to the total fluence.Results: Using film dosimetry for the s3D-CRT plans, the average of the ratio of calculated dose to measured dose at the isocenter was 1.01 and 1.06 for the OBM and MBM model, respectively. A significantly lower GPR of the MBM compared to the OBM was only found for the localized prostate cases (F = 12.4%) measured with the 2D-detector-array and an acceptance criterion of 1%/3 mm.Conclusion: The effects of uncertainties in small field OFs implemented in TPSs are most pronounced for s3D-CRT cases and can be clearly identified using patient specific quality assurance. For VMAT these effects mainly remain undetected using standard patient specific quality assurance. Using tighter acceptance criteria combined with an analysis of the fluence generated by small fields can help identifying inaccuracies of OFs implemented in TPSs.

IAEA methodology for on-site end-to-end IMRT/VMAT audits: an international pilot study.

Background: The IAEA has developed and tested an on-site, end-to-end IMRT/VMAT dosimetry audit methodology for head and neck cases using an anthropomorphic phantom. The audit methodology is described, and the results of the international pilot testing are presented.Material and methods: The audit utilizes a specially designed, commercially available anthropomorphic phantom capable of accommodating a small volume ion chamber (IC) in four locations (three in planning target volumes (PTVs) and one in an organ at risk (OAR)) and a Gafchromic film in a coronal plane for the absorbed dose to water and two-dimensional dose distribution measurements, respectively. The audit consists of a pre-visit and on-site phases. The pre-visit phase is carried out remotely and includes a treatment planning task and a set of computational exercises. The on-site phase aims at comparing the treatment planning system (TPS) calculations with measurements in the anthropomorphic phantom following an end-to-end approach. Two main aspects were tested in the pilot study: feasibility of the planning constraints and the accuracy of IC and film results in comparison with TPS calculations. Treatment plan quality was scored from 0 to 100.Results: Forty-two treatment plans were submitted by 14 institutions from 10 countries, with 79% of them having a plan quality score over 90. Seventeen sets of IC measurement results were collected, and the average measured to calculated dose ratio was 0.988 ± 0.016 for PTVs and 1.020 ± 0.029 for OAR. For 13 film measurement results, the average gamma passing rate was 94.1% using criteria of 3%/3 mm, 20% threshold and global gamma.Conclusions: The audit methodology was proved to be feasible and ready to be adopted by national dosimetry audit networks for local implementation.

Testing the methodology for a dosimetric end-to-end audit of IMRT/VMAT: results of IAEA multicentre and national studies.

Introduction: Within an International Atomic Energy Agency (IAEA) co-ordinated research project (CRP), a remote end-to-end dosimetric quality audit for intensity modulated radiation therapy (IMRT)/ volumetric arc therapy (VMAT) was developed to verify the radiotherapy chain including imaging, treatment planning and dose delivery. The methodology as well as the results obtained in a multicentre pilot study and national trial runs conducted in close cooperation with dosimetry audit networks (DANs) of IAEA Member States are presented.Material and methods: A solid polystyrene phantom containing a dosimetry insert with an irregular solid water planning target volume (PTV) and organ at risk (OAR) was designed for this audit. The insert can be preloaded with radiochromic film and four thermoluminescent dosimeters (TLDs). For the audit, radiotherapy centres were asked to scan the phantom, contour the structures, create an IMRT/VMAT treatment plan and irradiate the phantom. The dose prescription was to deliver 4 Gy to the PTV in two fractions and to limit the OAR dose to a maximum of 2.8 Gy. The TLD measured doses and film measured dose distributions were compared with the TPS calculations.Results: Sixteen hospitals from 13 countries and 64 hospitals from 6 countries participated in the multicenter pilot study and in the national runs, respectively. The TLD results for the PTV were all within ±5% acceptance limit for the multicentre pilot study, whereas for national runs, 17 participants failed to meet this criterion. All measured doses in the OAR were below the treatment planning constraint. The film analysis identified seven plans in national runs below the 90% passing rate gamma criteria.Conclusion: The results proved that the methodology of the IMRT/VMAT dosimetric end-to-end audit was feasible for its intended purpose, i.e., the phantom design and materials were suitable; the phantom was easy to use and it was robust enough for shipment. Most importantly the audit methodology was capable of identifying suboptimal IMRT/VMAT delivery.

Glass Transitions in Monodisperse Cluster-Forming Ensembles: Vortex Matter in Type-1.5 Superconductors.

At low enough temperatures and high densities, the equilibrium configuration of an ensemble of ultrasoft particles is a self-assembled, ordered, cluster crystal. In the present Letter, we explore the out-of-equilibrium dynamics for a two-dimensional realization, which is relevant to superconducting materials with multiscale intervortex forces. We find that, for small temperatures following a quench, the suppression of the thermally activated particle hopping hinders the ordering. This results in a glass transition for a monodispersed ensemble, for which we derive a microscopic explanation in terms of an "effective polydispersity" induced by multiscale interactions. This demonstrates that a vortex glass can form in clean systems of thin films of "type-1.5" superconductors. An additional setup to study this physics can be layered superconducting systems, where the shape of the effective vortex-vortex interactions can be engineered.

Spatial Patterns in Rydberg Excitations from Logarithmic Pair Interactions.

The collective excitations in ensembles of dissipative, laser driven ultracold atoms exhibit crystal-like patterns, a many-body effect of the Rydberg blockade mechanism. These crystalline structures are revealed in an experiment from a postselection of configurations with fixed numbers of excitations. Here, we show that these subensembles can be well represented by ensembles of effective particles that interact via logarithmic pair potentials. This allows one to study the emergent patterns with a small number of effective particles to determine the phases of Rydberg crystals and to systematically study contributions from N-body terms.

Role of quantum fluctuations in the hexatic phase of cold polar molecules.

Two-dimensional crystals melt via an intermediate hexatic phase, which is characterized by an anomalous scaling of spatial and orientational correlation functions and the absence of an attraction between dislocations. We propose a protocol to study the effect of quantum fluctuations on the nature of this phase with a model system of strongly correlated ultracold polar molecules. Dislocations can be located in experiment from local energy differences which induce internal stark shifts in the molecules. We present a criterion to identify the hexatic phase from the statistics of the end points of topological defect strings and find a hexatic phase, which is dominated by quantum fluctuations, between the crystal and superfluid phases.

Equilibrium and non-equilibrium cluster phases in colloids with competing interactions.

The phase behavior of colloids that interact via competing interactions - short-range attraction and long-range repulsion - is studied by computer simulation. In particular, for a fixed strength and range of repulsion, the effect of the strength of an attractive interaction (ε) on the phase behavior is investigated at various colloid densities (ρ). A thermodynamically stable equilibrium colloidal cluster phase, consisting of compact crystalline clusters, is found below the fluid-solid coexistence line in the ε-ρ parameter space. The mean cluster size is found to linearly increase with the colloid density. At large ε and low densities, and at small ε and high densities, a non-equilibrium cluster phase, consisting of elongated Bernal spiral-like clusters, is observed. Although gelation can be induced either by increasing ε at constant density or vice versa, the gelation mechanism is different in either route. While in the ρ route gelation occurs via a glass transition of compact clusters, gelation in the ε route is characterized by percolation of elongated clusters. This study both provides the location of equilibrium and non-equilibrium cluster phases with respect to the fluid-solid coexistence, and reveals the dependencies of the gelation mechanism on the preparation route.

From classical to quantum glasses with ultracold polar molecules.

We study the dynamics of a bilayer system of ultracold polar molecules, which exhibits classical and quantum glassy behavior, characterized by long tails in the relaxation time and dynamical heterogeneity. In the proposed setup, quantum fluctuations are of the order of thermal fluctuations and the degree of frustration can be tuned by the interlayer distance. We discuss the possible observation of a glassy anomalous diffusion and dynamical heterogeneity in experiment using internal degrees of freedom of the molecules in combination with optical detection.

Cluster Luttinger liquids of Rydberg-dressed atoms in optical lattices.

We investigate the zero-temperature phases of bosonic and fermionic gases confined to one dimension and interacting via a class of finite-range soft-shoulder potentials (i.e., soft-core potentials with an additional hard-core onsite interaction). Using a combination of analytical and numerical methods, we demonstrate the stabilization of critical quantum liquids with qualitatively new features with respect to the Tomonaga-Luttinger liquid paradigm. These features result from frustration and cluster formation in the corresponding classical ground state. Characteristic signatures of these liquids are accessible in state-of-the-art experimental setups with Rydberg-dressed ground-state atoms trapped in optical lattices.

Scalable quantum processors empowered by the Fermi scattering of Rydberg electrons.

Quantum computing promises exponential speed-up compared to its classical counterpart. While the neutral atom processors are the pioneering platform in terms of scalability, the dipolar Rydberg gates impose the main bottlenecks on the scaling of these devices. This article presents an alternative scheme for neutral atom quantum processing, based on the Fermi scattering of a Rydberg electron from ground-state atoms in spin-dependent lattice geometries. Instead of relying on Rydberg pair-potentials, the interaction is controlled by engineering the electron cloud of a sole Rydberg atom. The present scheme addresses the scaling obstacles in Rydberg processors by exponentially suppressing the population of short-lived states and by operating in ultra-dense atomic lattices. The restoring forces in molecule type Rydberg-Fermi potential preserve the trapping over a long interaction period. Furthermore, the proposed scheme mitigates different competing infidelity criteria, eliminates unwanted cross-talks, and significantly suppresses the operation depth in running complicated quantum algorithms.

Scalable set of reversible parity gates for integer factorization.

Classical microprocessors operate on irreversible gates, that, when combined with AND, half-adder and full-adder operations, execute complex tasks such as multiplication of integers. We introduce parity versions of all components of a multiplication circuit. The parity gates are reversible quantum gates based on the recently introduced parity transformation and build on ground-space encoding of the corresponding gate logic. Using a quantum optimization heuristic, e.g., an adiabatic quantum computing protocol, allows one to quantum mechanically reverse the process of multiplication and thus factor integers, which has applications in cryptography. Our parity approach builds on nearest-neighbor constraints equipped with local fields, able to encode the logic of a binary multiplication circuit in a modular and scalable way.

Evaluation of a novel CBCT conversion method implemented in a treatment planning system.

BACKGROUND: To evaluate a novel CBCT conversion algorithm for dose calculation implemented in a research version of a treatment planning system (TPS). METHODS: The algorithm was implemented in a research version of RayStation (v. 11B-DTK, RaySearch, Stockholm, Sweden). CBCTs acquired for each ten head and neck (HN), gynecology (GYN) and lung cancer (LNG) patients were collected and converted using the new algorithm (CBCTc). A bulk density overriding technique implemented in the same version of the TPS was used for comparison (CBCTb). A deformed CT (dCT) was created by using deformable image registration of the planning CT (pCT) to the CBCT to reduce anatomical changes. All treatment plans were recalculated on the pCT, dCT, CBCTc and the CBCTb. The resulting dose distributions were analyzed using the MICE toolkit (NONPIMedical AB Sweden, Umeå) with local gamma analysis, with 1% dose difference and 1 mm distance to agreement criteria. A Wilcoxon paired rank sum test was applied to test the differences in gamma pass rates (GPRs). A p value smaller than 0.05 considered statistically significant. RESULTS: The GPRs for the CBCTb method were systematically lower compared to the CBCTc method. Using the 10% dose threshold and the dCT as reference the median GPRs were for the CBCTc method were 100% and 99.8% for the HN and GYN cases, respectively. Compared to that the GPRs of the CBCTb method were lower with values of 99.8% and 98.0%, for the HN and GYN cases, respectively. The GPRs of the LNG cases were 99.9% and 97.5% for the CBCTc and CBCTb method, respectively. These differences were statistically significant. The main differences between the dose calculated on the CBCTs and the pCTs were found in regions near air/tissue interfaces, which are also subject to anatomical variations. CONCLUSION: The dose distribution calculated using the new CBCTc method showed excellent agreement with the dose calculated using dCT and pCT and was superior to the CBCTb method. The main reasons for deviations of the calculated dose distribution were caused by anatomical variations between the pCT and the corrected CBCT.

Impact of log file source and data frequency on accuracy of log file-based patient specific quality assurance.

Performing phantom measurements for patient-specific quality assurance (PSQA) adds a significant amount of time to the adaptive radiotherapy procedure. Log file based PSQA can be used to increase the efficiency of this process. This study compared the dosimetric accuracy of high-frequency linear accelerator (Linac) log files and low-frequency log data stored in the oncology information system (OIS). Thirty patients were included, that were recently treated in the head and neck (HN), brain, and prostate region with volumetric modulated arc therapy (VMAT) and an additional ten patients treated using stereotactic body radiation therapy (SBRT) with 3D-conformal radiotherapy (3D-CRT) technique. Log data containing a single fraction were used to calculate the dose distributions. The dosimetric differences between Linac log files and OIS logs were evaluated with a gamma analysis with 2%/2 mm criterion and dose threshold of 30%. The original treatment plan was used as a reference. Moreover, DVH parameters of D98%, D50%, and D2% of the planning-target volume (PTV) and dose to several organs at risk (OARs) were reported. Significant differences in dose distributions between the two log types and the original dose were observed for PTV D98% and D2% (r < 0.001) for HN cases, PTV D98% (r = 0.005) for brain cases, and PTV D50% (r = 0.015) for prostate cases. No significant differences were found between the two log types with respect to D50%. The root mean square (RMS) error of the leaf positions of the OIS log was approximately twice the RMS error of the Linac log file for VMAT plans, but identical for 3D-CRT plans. The relationship between the gamma pass rate and the RMS error showed a moderate correlation for the Linac log files (r = -0.58, p < 0.001) and strong correlation for OIS logs (r = -0.71, p < 0.001). Furthermore, all doses calculated using Linac log files and OIS log data had a GPR >90% for an RMS error < 3.3 mm. Based on these findings, a tolerance limit of RMS error of 3.3 mm for considering OIS log based PSQA was established. Nevertheless, the OIS log data quality should be improved to achieve adequate PSQA.

Generating deliverable DICOM RT treatment plans for prostate VMAT by predicting MLC motion sequences with an encoder-decoder network.

BACKGROUND: Deep learning-based auto-planning is an active research field; however, for some tasks a treatment planning system (TPS) is still required. PURPOSE: To introduce a deep learning-based model generating deliverable DICOM RT treatment plans that can be directly irradiated by a linear accelerator (LINAC). The model was based on an encoder-decoder network and can predict multileaf collimator (MLC) motion sequences for prostate VMAT radiotherapy. METHODS: A total of 619 treatment plans from 460 patients treated for prostate cancer with single-arc VMAT were included in this study. An encoder-decoder network was trained using 465 clinical treatment plans and validated on 77 plans. The performance was analyzed on a separate test set of 77 treatment plans. Separate L1 losses were computed for the leaf and jaw positions as well as the monitor units, with the leaf loss being weighted by a factor of 100 before being added to the other losses. The generated treatment plans were recalculated in a treatment planning system and the dose-volume metrics and gamma passing rates were compared to the original dose. RESULTS: All generated treatment plans showed good agreement with the original data, with an average gamma passing rate (3%/3 mm) of 91.9 ± 7.1%. However, the coverage of the PTVs. was slightly lower for the generated plans (D98%  = 92.9 ± 2.6%) in comparison to the original plans (D98%  = 95.7 ± 2.2%). There was no significant difference in mean dose to the bladder between the predicted and original plan (Dmean of 28.0 ± 13.5 vs. 28.1 ± 13.3% of prescribed dose) or rectum (Dmean of 42.3 ± 7.4 vs. 42.6 ± 7.5%). The maximum dose to bladder was only slightly higher in the predicted plans (D2% of 100.7 ± 5.3 vs. 99.8 ± 4.0%) and for the rectum it was even lower (D2% of 100.5 ± 3.7 vs. 100.1 ± 4.3). CONCLUSIONS: The deep learning-based model could predict MLC motion sequences in prostate VMAT plans, eliminating the need for sequencing inside a TPS, thus revolutionizing autonomous treatment planning workflows. This research completes the loop in deep learning-based treatment planning processes, enabling more efficient workflows for real-time or online adaptive radiotherapy.

Review and recommendations on deformable image registration uncertainties for radiotherapy applications.

Deformable image registration (DIR) is a versatile tool used in many applications in radiotherapy (RT). DIR algorithms have been implemented in many commercial treatment planning systems providing accessible and easy-to-use solutions. However, the geometric uncertainty of DIR can be large and difficult to quantify, resulting in barriers to clinical practice. Currently, there is no agreement in the RT community on how to quantify these uncertainties and determine thresholds that distinguish a good DIR result from a poor one. This review summarises the current literature on sources of DIR uncertainties and their impact on RT applications. Recommendations are provided on how to handle these uncertainties for patient-specific use, commissioning, and research. Recommendations are also provided for developers and vendors to help users to understand DIR uncertainties and make the application of DIR in RT safer and more reliable.

Clinical Implementation and Evaluation of Auto-Segmentation Tools for Multi-Site Contouring in Radiotherapy.

Background and purpose: Tools for auto-segmentation in radiotherapy are widely available, but guidelines for clinical implementation are missing. The goal was to develop a workflow for performance evaluation of three commercial auto-segmentation tools to select one candidate for clinical implementation. Materials and Methods: One hundred patients with six treatment sites (brain, head-and-neck, thorax, abdomen, and pelvis) were included. Three sets of AI-based contours for organs-at-risk (OAR) generated by three software tools and manually drawn expert contours were blindly rated for contouring accuracy. The dice similarity coefficient (DSC), the Hausdorff distance, and a dose/volume evaluation based on the recalculation of the original treatment plan were assessed. Statistically significant differences were tested using the Kruskal-Wallis test and the post-hoc Dunn Test with Bonferroni correction. Results: The mean DSC scores compared to expert contours for all OARs combined were 0.80 ± 0.10, 0.75 ± 0.10, and 0.74 ± 0.11 for the three software tools. Physicians' rating identified equivalent or superior performance of some AI-based contours in head (eye, lens, optic nerve, brain, chiasm), thorax (e.g., heart and lungs), and pelvis and abdomen (e.g., kidney, femoral head) compared to manual contours. For some OARs, the AI models provided results requiring only minor corrections. Bowel-bag and stomach were not fit for direct use. During the interdisciplinary discussion, the physicians' rating was considered the most relevant. Conclusion: A comprehensive method for evaluation and clinical implementation of commercially available auto-segmentation software was developed. The in-depth analysis yielded clear instructions for clinical use within the radiotherapy department.