Simple beam models for Monte Carlo photon beam dose calculations in radiotherapy

Monte Carlo (code GEANT) produced 6 and 15 MV phase space (PS) data were used to define several simple photon beam models. For creating the PS data the energy of starting electrons hitting the target was tuned to get correct depth dose data compared to measurements. The modeling process used the full PS information within the geometrical boundaries of the beam including all scattered radiation of the accelerator head. Scattered radiation outside the boundaries was neglected. Photons and electrons were assumed to be radiated from point sources. Four different models were investigated which involved different ways to determine the energies and locations of beam particles in the output plane. Depth dose curves, profiles, and relative output factors were calculated with these models for six field sizes from 5x5 to 40x40cm2 and compared to measurements. Model 1 uses a photon energy spectrum independent of location in the PS plane and a constant photon fluence in this plane. Model 2 takes into account the spatial particle fluence distribution in the PS plane. A constant fluence is used again in model 3, but the photon energy spectrum depends upon the off axis position. Model 4, finally uses the spatial particle fluence distribution and off axis dependent photon energy spectra in the PS plane. Depth dose curves and profiles for field sizes up to 10x10cm2 were not model sensitive. Good agreement between measured and calculated depth dose curves and profiles for all field sizes was reached for model 4. However, increasing deviations were found for increasing field sizes for models 1-3. Large deviations resulted for the profiles of models 2 and 3. This is due to the fact that these models overestimate and underestimate the energy fluence at large off axis distances. Relative output factors consistent with measurements resulted only for model 4.

Practicality of time-optimal two-qubit Hamiltonian simulation

What is the time-optimal way of using a set of control Hamiltonians to obtain a desired interaction? Vidal, Hammerer, and Cirac [Phys. Rev. Lett. 88, 237902 (2002)] have obtained a set of powerful results characterizing the time-optimal simulation of a two-qubit quantum gate using a fixed interaction Hamiltonian and fast local control over the individual qubits. How practically useful are these results? We prove that there are two-qubit Hamiltonians such that time-optimal simulation requires infinitely many steps of evolution, each infinitesimally small, and thus is physically impractical. A procedure is given to determine which two-qubit Hamiltonians have this property, and we show that almost all Hamiltonians do. Finally, we determine some bounds on the penalty that must be paid in the simulation time if the number of steps is fixed at a finite number, and show that the cost in simulation time is not too great.

Spin squeezing as a measure of entanglement in a two-qubit system

We show that the two definitions of spin squeezing extensively used in the literature [M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993) and D.J. Wineland , Phys. Rev. A 50, 67 (1994)] give different predictions of entanglement in the two-atom Dicke system. We analyze differences between the definitions and show that the spin squeezing parameter of Kitagawa and Ueda is a better measure of entanglement than the commonly used spectroscopic spin squeezing parameter. We illustrate this relation by examining different examples of a driven two-atom Dicke system in which spin squeezing and entanglement arise dynamically. We give an explanation of the source of the difference using the negativity criterion for entanglement.

Fast nonadiabatic two-qubit gates for the Kane quantum computer

In this paper, we apply the canonical decomposition of two-qubit unitaries to find pulse schemes to control the proposed Kane quantum computer. We explicitly find pulse sequences for the controlled-NOT, swap, square root of swap, and controlled Z rotations. We analyze the speed and fidelity of these gates, both of which compare favorably to existing schemes. The pulse sequences presented in this paper are theoretically faster, with higher fidelity, and simpler. Any two-qubit gate may be easily found and implemented using similar pulse sequences. Numerical simulation is used to verify the accuracy of each pulse scheme.

Numerical analysis of capacities for two-qubit unitary operations

We present numerical results on the capacities of two-qubit unitary operations for performing communication and creating entanglement. The capacities for communication considered are based upon the increase in Holevo information of an ensemble. Our results indicate that the capacity may be accurately estimated using ensemble sizes and ancilla dimensions of 4. In addition, the calculated values of these capacities were close to, and in some cases equal to, the similarly defined entangling capacities; this result indicates connections between these capacities.

Monte Carlo simulations for dosimetric verification in photon and electron beam radiotherapy

Dissertação para obtenção do Grau de Doutor em Engenharia Biomédica

Brazilian experience with two conditioning regimens in patients with multiple sclerosis: BEAM/horse ATG and CY/rabbit ATG

Studies have shown that autologous hematopoietic SCT (HSCT) can be used as an intensive immunosuppressive therapy to treat refractory patients and to prevent the progression of multiple sclerosis (MS). This is a prospective multicentric Brazilian MS trial comparing two conditioning regimens: BEAM/horse ATG and CY/rabbit ATG. Most (80.4%) of the 41 subjects in the study had the secondary progressive MS subtype and the mean age was 42 years. The baseline EDSS score in 58.5% of the subjects was 6.5 and 78% had a score of 6.0 or higher, respectively. The complication rate during the intra-transplantation period was 56% for all patients: 71.4% of the patients in the BEAM/hATG group and 40% in the CY/rATG group (P = 0.04). Three subjects (7.5%) died of cardiac toxicity, sepsis and alveolar hemorrhage, all of them in the BEAM/ATG group. EFS was 58.54% for a ll patients: 47% in the BEAM/hATG group and 70% in the CY/rATG group (P = 0.288). In conclusion, the CY/rATG regimen seems to be associated with similar outcome results, but presented less toxicity when compared with the BEAM/hATG regimen. Long-term follow-up would be required to fully assess the differences in therapeutic effectiveness between the two regimens. Bone Marrow Transplantation (2010) 45, 239-248; doi:10.1038/bmt.2009.127; published online 6 July 2009

Investigation of Two Bridge Alternatives for Low Volume Roads Volume 2 of 2 Concept 2: Beam In Slab Bridge, HR-382, 1997

Recent reports have indicated that 23.5% of the nation's highway bridges are structurally deficient and 17.7% are functionally obsolete. A significant number of these bridges are on the Iowa secondary road system where over 86% of the rural bridge management responsibilities are assigned to the counties. Some of the bridges can be strengthened or otherwise rehabilitated, but many more are in need of immediate replacement. In a recent investigation (HR-365 "Evaluation of Bridge Replacement Alternatives for the County Bridge System") several types of replacement bridges that are currently being used on low volume roads were identified. It was also determined that a large number of counties (69%) have the ability and are interested in utilizing their own forces to design and construct short span bridges. After reviewing the results from HR-365, the research team developed one "new" bridge replacement concept and a modification of a replacement system currently being used. Both of these bridge replacement alternatives were investigated in this study, the results of which are presented in two volumes. This volume (Volume 2) presents the results of Concept 2 - Modification of the Beam-in-Slab Bridge. Concept 1 - Steel Beam Precast Units is presented in Volume 1. Concept 2 involves various laboratory tests of the Beam-in-Slab bridge (BISB) currently being used by Benton County and several other Iowa counties. In this investigation, the behavior and strength of the BISB were determined; a new method of obtaining composite action between the steel beams and concrete was also tested. Since the Concept 2 bridge is primarily intended for use on low-volume roads, the system can be constructed with new or used beams. In the experimental part of the investigation, there were three types of laboratory tests: push-out tests, service and ultimate load tests of models of the BISB, and composite beam tests utilizing the newly developed shear connection. In addition to the laboratory tests, there was a field test in which an existing BISB was service load tested. An equation was developed for predicting the strength of the shear connection investigated; in addition, a finite element model for analyzing the BISB was also developed. Push-out tests were completed to determine the strength of the recently developed shear connector. A total of 36 specimens were tested, with variables such as hole diameter, hole spacing, presence of reinforcement, etc. being investigated. In the model tests of the BISB, two and four beam specimens [L=9,140 mm (30 ft)] were service load tested for behavior and load distribution data. Upon completion of these tests, both specimens were loaded to failure. In the composite beam tests, four beams, one with standard shear studs and three using the shear connection developed, were tested. Upon completion of the service load tests, all four beams were loaded to failure. The strength and behavior of the beams with the new shear connection were found to be essentially the same as that of the specimen with standard shear studs.

Stability of the Helical TomoTherapy Hi·Art II detector for treatment beam irradiations.

The Hi·Art II Helical TomoTherapy (HT) unit is equipped with a built-in onboard MVCT detector used for patient imaging and beam monitoring. Our aim was to study the detector stability for treatment beam measurements. We studied the MVCT detector response with the 6 MV photon beam over time, throughout short-term (during an irradiation) and long-term (two times 50 days) periods. Our results show a coefficient of variation ≤ 1% for detector chambers inside the beam (excluding beam gradients) for short- and long-term response of the MVCT detector. Larger variations were observed in beam gradients and an influence of the X-ray target where degradation was found. The results assume that an 'air scan' procedure is performed daily to recalibrate the detector with the imaging beam. On short term, the detector response stability is comparable to other devices. Long-term measure- ments during two 50-day periods show a good reproducibility. 

Understanding the Influence of Photon Energy on 6MV Non Reference Dosimetry and CT Dosimetry using TLD and OSLD

Measurement of the absorbed dose from ionizing radiation in medical applications is an essential component to providing safe and reproducible patient care. There are a wide variety of tools available for measuring radiation dose; this work focuses on the characterization of two common, solid-state dosimeters in medical applications: thermoluminescent dosimeters (TLD) and optically stimulated luminescent dosimeters (OSLD). There were two main objectives to this work. The first objective was to evaluate the energy dependence of TLD and OSLD for non-reference measurement conditions in a radiotherapy environment. The second objective was to fully characterize the OSLD nanoDot in a CT environment, and to provide validated calibration procedures for CT dose measurement using OSLD. Current protocols for dose measurement using TLD and OSLD generally assume a constant photon energy spectrum within a nominal beam energy regardless of measurement location, tissue composition, or changes in beam parameters. Variations in the energy spectrum of therapeutic photon beams may impact the response of TLD and OSLD and could thereby result in an incorrect measure of dose unless these differences are accounted for. In this work, we used a Monte Carlo based model to simulate variations in the photon energy spectra of a Varian 6MV beam; then evaluated the impact of the perturbations in energy spectra on the response of both TLD and OSLD using Burlin Cavity Theory. Energy response correction factors were determined for a range of conditions and compared to measured correction factors with good agreement. When using OSLD for dose measurement in a diagnostic imaging environment, photon energy spectra are often referenced to a therapy-energy or orthovoltage photon beam – commonly 250kVp, Co-60, or even 6MV, where the spectra are substantially different. Appropriate calibration techniques specifically for the OSLD nanoDot in a CT environment have not been presented in the literature; furthermore the dependence of the energy response of the calibration energy has not been emphasized. The results of this work include detailed calibration procedures for CT dosimetry using OSLD, and a full characterization of this dosimetry system in a low-dose, low-energy setting.

Conditional production of superpositions of coherent states with inefficient photon detection

It is shown that a linear superposition of two macroscopically distinguishable optical coherent states can be generated using a single photon source and simple all-optical operations. Weak squeezing on a single photon, beam mixing with an auxiliary coherent state, and photon detecting with imperfect threshold detectors are enough to generate a coherent state superposition in a free propagating optical field with a large coherent amplitude (alpha>2) and high fidelity (F>0.99). In contrast to all previous schemes to generate such a state, our scheme does not need photon number resolving measurements nor Kerr-type nonlinear interactions. Furthermore, it is robust to detection inefficiency and exhibits some resilience to photon production inefficiency.

NMR analog of Bell's inequalities violation test

In this paper, we present an analog of Bell's inequalities violation test for N qubits to be performed in a nuclear magnetic resonance (NMR) quantum computer. This can be used to simulate or predict the results for different Bell's inequality tests, with distinct configurations and a larger number of qubits. To demonstrate our scheme, we implemented a simulation of the violation of the Clauser, Horne, Shimony and Holt (CHSH) inequality using a two-qubit NMR system and compared the results to those of a photon experiment. The experimental results are well described by the quantum mechanics theory and a local realistic hidden variables model (LRHVM) that was specifically developed for NMR. That is why we refer to this experiment as a simulation of Bell's inequality violation. Our result shows explicitly how the two theories can be compatible with each other due to the detection loophole. In the last part of this work, we discuss the possibility of testing some fundamental features of quantum mechanics using NMR with highly polarized spins, where a strong discrepancy between quantum mechanics and hidden variables models can be expected.

Calculation of correction factors for ionization chamber measurements with small fields in low-density media.

The quantity of interest for high-energy photon beam therapy recommended by most dosimetric protocols is the absorbed dose to water. Thus, ionization chambers are calibrated in absorbed dose to water, which is the same quantity as what is calculated by most treatment planning systems (TPS). However, when measurements are performed in a low-density medium, the presence of the ionization chamber generates a perturbation at the level of the secondary particle range. Therefore, the measured quantity is close to the absorbed dose to a volume of water equivalent to the chamber volume. This quantity is not equivalent to the dose calculated by a TPS, which is the absorbed dose to an infinitesimally small volume of water. This phenomenon can lead to an overestimation of the absorbed dose measured with an ionization chamber of up to 40% in extreme cases. In this paper, we propose a method to calculate correction factors based on the Monte Carlo simulations. These correction factors are obtained by the ratio of the absorbed dose to water in a low-density medium □D(w,Q,V1)(low) averaged over a scoring volume V₁ for a geometry where V₁ is filled with the low-density medium and the absorbed dose to water □D(w,QV2)(low) averaged over a volume V₂ for a geometry where V₂ is filled with water. In the Monte Carlo simulations, □D(w,QV2)(low) is obtained by replacing the volume of the ionization chamber by an equivalent volume of water, according to the definition of the absorbed dose to water. The method is validated in two different configurations which allowed us to study the behavior of this correction factor as a function of depth in phantom, photon beam energy, phantom density and field size.

Radiothérapie asservie à la respiration en combinaison avec l'utilisation d'un faisceau sans filtre égalisateur

La radiothérapie stéréotaxique corporelle (SBRT) est une technique couramment employée pour le traitement de tumeurs aux poumons lorsque la chirurgie n’est pas possible ou refusée par le patient. Une complication de l’utilisation de cette méthode provient du mouvement de la tumeur causé par la respiration. Dans ce contexte, la radiothérapie asservie à la respiration (RGRT) peut être bénéfique. Toutefois, la RGRT augmente le temps de traitement en raison de la plus petite proportion de temps pour laquelle le faisceau est actif. En utilisant un faisceau de photons sans filtre égalisateur (FFF), ce problème peut être compensé par le débit de dose plus élevé d’un faisceau FFF. Ce mémoire traite de la faisabilité d’employer la technique de RGRT en combinaison avec l’utilisation un faisceau FFF sur un accélérateur Synergy S (Elekta, Stockholm, Suède) avec une ceinture pneumatique, le Bellows Belt (Philips, Amsterdam, Pays-Bas), comme dispositif de suivi du signal respiratoire. Un Synergy S a été modifié afin de pouvoir livrer un faisceau 6 MV FFF. Des mesures de profils de dose et de rendements en profondeur ont été acquises en cuve à eau pour différentes tailles de champs. Ces mesures ont été utilisées pour créer un modèle du faisceau 6 MV FFF dans le système de planification de traitement Pinnacle3 de Philips. Les mesures ont été comparées au modèle à l’aide de l’analyse gamma avec un critère de 2%, 2 mm. Par la suite, cinq plans SBRT avec thérapie en arc par modulation volumétrique (VMAT) ont été créés avec le modèle 6 MV du Synergy S, avec et sans filtre. Une comparaison des paramètres dosimétriques a été réalisée entre les plans avec et sans filtre pour évaluer la qualité des plans FFF. Les résultats révèlent qu’il est possible de créer des plans SBRT VMAT avec le faisceau 6 MV FFF du Synergy S qui sont cliniquement acceptables (les crières du Radiation Therapy Oncology Group 0618 sont respectés). Aussi, une interface physique de RGRT a été mise au point pour remplir deux fonctions : lire le signal numérique de la ceinture pneumatique Bellows Belt et envoyer une commande d’irradiation binaire au linac. L’activation/désactivation du faisceau du linac se fait par l’entremise d’un relais électromécanique. L’interface comprend un circuit électronique imprimé fait maison qui fonctionne en tandem avec un Raspberry Pi. Un logiciel de RGRT a été développé pour opérer sur le Raspberry Pi. Celui-ci affiche le signal numérique du Bellows Belt et donne l’option de choisir les limites supérieure et inférieure de la fenêtre d’irradiation, de sorte que lorsque le signal de la ceinture se trouve entre ces limites, le faisceau est actif, et inversement lorsque le signal est hors de ces limites. Le logiciel envoie donc une commande d’irradiation au linac de manière automatique en fonction de l’amplitude du signal respiratoire. Finalement, la comparaison entre la livraison d’un traitement standard sans RGRT avec filtre par rapport à un autre plan standard sans RGRT sans filtre démontre que le temps de traitement en mode FFF est réduit en moyenne de 54.1% pour un arc. De la même manière, la comparaison entre la livraison d’un traitement standard sans RGRT avec filtre par rapport à un plan de RGRT (fenêtre d’irradiation de 75%) sans filtre montre que le temps de traitement de RGRT en mode FFF est réduit en moyenne de 27.3% par arc. Toutefois, il n’a pas été possible de livrer des traitements de RGRT avec une fenêtre de moins de 75%. Le linac ne supporte pas une fréquence d’arrêts élevée.

Extracting Information from Qubit-Environment Correlations

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)