The STATFLUX code: a statistical method for calculation of flow and set of parameters, based on the Multiple-Compartment Biokinetical Model

The code STATFLUX, implementing a new and simple statistical procedure for the calculation of transfer coefficients in radionuclide transport to animals and plants, is proposed. The method is based on the general multiple-compartment model, which uses a system of linear equations involving geometrical volume considerations. Flow parameters were estimated by employing two different least-squares procedures: Derivative and Gauss-Marquardt methods, with the available experimental data of radionuclide concentrations as the input functions of time. The solution of the inverse problem, which relates a given set of flow parameter with the time evolution of concentration functions, is achieved via a Monte Carlo Simulation procedure.Program summaryTitle of program: STATFLUXCatalogue identifier: ADYS_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADYS_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputer for which the program is designed and others on which it has been tested: Micro-computer with Intel Pentium III, 3.0 GHzInstallation: Laboratory of Linear Accelerator, Department of Experimental Physics, University of São Paulo, BrazilOperating system: Windows 2000 and Windows XPProgramming language used: Fortran-77 as implemented in Microsoft Fortran 4.0. NOTE: Microsoft Fortran includes non-standard features which are used in this program. Standard Fortran compilers such as, g77, f77, ifort and NAG95, are not able to compile the code and therefore it has not been possible for the CPC Program Library to test the program.Memory, required to execute with typical data: 8 Mbytes of RAM memory and 100 MB of Hard disk memoryNo. of bits in a word: 16No. of lines in distributed program, including test data, etc.: 6912No. of bytes in distributed Program, including test data, etc.: 229 541Distribution format: tar.gzNature of the physical problem: the investigation of transport mechanisms for radioactive substances, through environmental pathways, is very important for radiological protection of populations. One such pathway, associated with the food chain, is the grass-animal-man sequence. The distribution of trace elements in humans and laboratory animals has been intensively studied over the past 60 years [R.C. Pendlenton, C.W. Mays, R.D. Lloyd, A.L. Brooks, Differential accumulation of iodine-131 from local fallout in people and milk, Health Phys. 9 (1963) 1253-1262]. In addition, investigations on the incidence of cancer in humans, and a possible causal relationship to radioactive fallout, have been undertaken [E.S. Weiss, M.L. Rallison, W.T. London, W.T. Carlyle Thompson, Thyroid nodularity in southwestern Utah school children exposed to fallout radiation, Amer. J. Public Health 61 (1971) 241-249; M.L. Rallison, B.M. Dobyns, F.R. Keating, J.E. Rall, F.H. Tyler, Thyroid diseases in children, Amer. J. Med. 56 (1974) 457-463; J.L. Lyon, M.R. Klauber, J.W. Gardner, K.S. Udall, Childhood leukemia associated with fallout from nuclear testing, N. Engl. J. Med. 300 (1979) 397-402]. From the pathways of entry of radionuclides in the human (or animal) body, ingestion is the most important because it is closely related to life-long alimentary (or dietary) habits. Those radionuclides which are able to enter the living cells by either metabolic or other processes give rise to localized doses which can be very high. The evaluation of these internally localized doses is of paramount importance for the assessment of radiobiological risks and radiological protection. The time behavior of trace concentration in organs is the principal input for prediction of internal doses after acute or chronic exposure. The General Multiple-Compartment Model (GMCM) is the powerful and more accepted method for biokinetical studies, which allows the calculation of concentration of trace elements in organs as a function of time, when the flow parameters of the model are known. However, few biokinetics data exist in the literature, and the determination of flow and transfer parameters by statistical fitting for each system is an open problem.Restriction on the complexity of the problem: This version of the code works with the constant volume approximation, which is valid for many situations where the biological half-live of a trace is lower than the volume rise time. Another restriction is related to the central flux model. The model considered in the code assumes that exist one central compartment (e.g., blood), that connect the flow with all compartments, and the flow between other compartments is not included.Typical running time: Depends on the choice for calculations. Using the Derivative Method the time is very short (a few minutes) for any number of compartments considered. When the Gauss-Marquardt iterative method is used the calculation time can be approximately 5-6 hours when similar to 15 compartments are considered. (C) 2006 Elsevier B.V. All rights reserved.

Estudo comparativo entre o tratamento radioterápico com elétrons e betaterapia, após cirurgia de queloides

Introduction: The surgery of keloids as an individual therapy with high rates of recurrence needs adjuvant therapies as complement. Recent keloids are rich in fibroblasts, which are highly radiosensitive and what explains the good results of percutaneous radiotherapy after surgery. Objective: To compare the beta-therapy results with electron-beam technique in newly operated keloids. Material and Methods: a prospective, comparative and randomized study was outlined including patients in the immediate postoperative period of surgical excision of keloids. Divided in G1 and G2 groups respectively, they received treatment with Sr90 (beta-therapy)plates or with electrons from the linear accelerator. The monitoring lasted 10 years. The results were performed using the information from patients, photograph parameters, observation and measurement of injuries, according to the criteria: Unchanged, Regular, Good and Excellent. Results: There were 26 patients, 13 in each group. In G1, 54% presented regular and unchanged improvement criteria and 46% had good or excellent criteria. In G2, the results were respectively 23% and 77%. Conclusion: The irradiation with electron-beam is better than beta-teraphy for the treatment of keloids surgery, due to the distribution in the tissue. No radio-induced tumors were observed.

Avaliação da resposta à quimiorradioterapia neoadjuvante em pacientes com adenocarcinomas retais

Pós-graduação em Bases Gerais da Cirurgia - FMB

Tratamento de quelóides usando radioterapia pós-operatória com elétrons: estudo comparativo de dois esquemas

Pós-graduação em Bases Gerais da Cirurgia - FMB

Radioterapia com elétrons

Radiotherapy is a field of medical physics, which has been going through a scientific and technological advancement, it is extreme important the professionals that work in this area continuous their study to improve the quality of service provided. For it, they should know the physical principles related as well their employment in radiotherapy. The electrons have been using in therapy of superficial tumors, because they show rapid decrease of dose in depth, they don’t expose the deeper tissues to radiation. The electron therapy has replaced the surface therapy with x-ray and mold brachytherapy, and nowadays, it represents 10 to 15% of radiotherapy treatments. This study brings the definition of physical parameters used in the dosimetry of electron beam, describes the tests for quality control of linear accelerator with a electrons beam, as well the activities realized in the Radiotherapy Division, of Hospital de Clinicas, in the Universidade Estadual of Campinas, in the Clinac 2100C machine

Uma avaliação do uso do NCRP 151 na realidade brasileira

Medical Physics is an interdisciplinary field that applies concepts and laws of physics in medical practices. Currently, one of its main applications is the use of ionizing radiation in the treatment of oncological diseases. Due to its wide use and highly dangerous, many of radioprotection procedures should be adopted with the objective of protecting human beings from harmful effects of radiation. Thus, you can better enjoy the benefits that the practice can offer. The methodology proposed by the National Council on Radiation Protection 151 (NCRP 151), relates technical information necessary to Structural Shielding Design and Evaluation for Megavoltage X- and Gamma- Ray Radiotherapy Facilities. However, many parameters used to calculate the shield are based on estimates only, and it is an international standard that may not be adequate to the Brazilian reality. Thus, the central idea of this study is the collection of data from the routine of the Radiotherapy Service of the Real e Benemérita Associação Portuguesa de Beneficência, in particular equipment cobalt therapy Theratron 780 (Atomic Energy of Canada Ltd.) and the linear accelerator Varian Clinac 2100C for measurement of workload, number of patients, fields, and dose factors to determine the best use of barrier protection. Furthermore, this work features a profile of radiotherapy treatments carried out closer to the Brazilian reality

Sistema de cálculo computacional da unidade monitora para feixes de elétrons em radioterapia

In radiation theraphy with electron beam, the electrons are produced in linear accelerators, and energy the most used have between 4MeV and 20MeV. Generally, the treatments are done for superficial injuries, because the low penetration of these particles. In this work a system for calculation of monitor units (U.M.) for cases of treatments with electron beam was developed. The Excel program of Microsoft was used and is easily found in the operational system of the personal microcomputers. In the Excel has been inserted the pertinent data of the linear accelerator of Varian, model 2100C, used in the Service of radiation theraphy of the Hospital of the Clinics of the College of Medicine of the UNESP of Botucatu. For some values of the physical parameters, such as: factors field and factors calibration, not supplied in the tests of acceptance of the machine, still proceeded calculations from interpolation and extrapolation. The mathematical formulas for automatic search of these and others factors used in the calculations of the determination of the U.M had been developed in agreement available routines in Excel. For this the functions had been used the function IF (that it imposes search condition) and the PROCH (that looks a value in a column from determined line), beyond the basic functions of addition, multiplication and division. It is intended to optimize the routine of the Services of radiation theraphy that perform through eletrontheraphy procedures, speeding the calculations and minimizing the occurrence of errors and uncertainties deriving of the maken a mistake manipulation of the parameters gotten in tables of data of electron beams

Measurement of the elastic electron-proton cross section and separation of the electric and magnetic form factor in the Q 2 range from 0.004 to 1 (GeV/c) 2

The electromagnetic form factors of the proton are fundamental quantities sensitive to the distribution of charge and magnetization inside the proton. Precise knowledge of the form factors, in particular of the charge and magnetization radii provide strong tests for theory in the non-perturbative regime of QCD. However, the existing data at Q^2 below 1 (GeV/c)^2 are not precise enough for a hard test of theoretical predictions.rnrnFor a more precise determination of the form factors, within this work more than 1400 cross sections of the reaction H(e,e′)p were measured at the Mainz Microtron MAMI using the 3-spectrometer-facility of the A1-collaboration. The data were taken in three periods in the years 2006 and 2007 using beam energies of 180, 315, 450, 585, 720 and 855 MeV. They cover the Q^2 region from 0.004 to 1 (GeV/c)^2 with counting rate uncertainties below 0.2% for most of the data points. The relative luminosity of the measurements was determined using one of the spectrometers as a luminosity monitor. The overlapping acceptances of the measurements maximize the internal redundancy of the data and allow, together with several additions to the standard experimental setup, for tight control of systematic uncertainties.rnTo account for the radiative processes, an event generator was developed and implemented in the simulation package of the analysis software which works without peaking approximation by explicitly calculating the Bethe-Heitler and Born Feynman diagrams for each event.rnTo separate the form factors and to determine the radii, the data were analyzed by fitting a wide selection of form factor models directly to the measured cross sections. These fits also determined the absolute normalization of the different data subsets. The validity of this method was tested with extensive simulations. The results were compared to an extraction via the standard Rosenbluth technique.rnrnThe dip structure in G_E that was seen in the analysis of the previous world data shows up in a modified form. When compared to the standard-dipole form factor as a smooth curve, the extracted G_E exhibits a strong change of the slope around 0.1 (GeV/c)^2, and in the magnetic form factor a dip around 0.2 (GeV/c)^2 is found. This may be taken as indications for a pion cloud. For higher Q^2, the fits yield larger values for G_M than previous measurements, in agreement with form factor ratios from recent precise polarized measurements in the Q2 region up to 0.6 (GeV/c)^2.rnrnThe charge and magnetic rms radii are determined as rn⟨r_e⟩=0.879 ± 0.005(stat.) ± 0.004(syst.) ± 0.002(model) ± 0.004(group) fm,rn⟨r_m⟩=0.777 ± 0.013(stat.) ± 0.009(syst.) ± 0.005(model) ± 0.002(group) fm.rnThis charge radius is significantly larger than theoretical predictions and than the radius of the standard dipole. However, it is in agreement with earlier results measured at the Mainz linear accelerator and with determinations from Hydrogen Lamb shift measurements. The extracted magnetic radius is smaller than previous determinations and than the standard-dipole value.

Laser driven proton acceleration and beam shaping

In the race to obtain protons with higher energies, using more compact systems at the same time, laser-driven plasma accelerators are becoming an interesting possibility. But for now, only beams with extremely broad energy spectra and high divergence have been produced. The driving line of this PhD thesis was the study and design of a compact system to extract a high quality beam out of the initial bunch of protons produced by the interaction of a laser pulse with a thin solid target, using experimentally reliable technologies in order to be able to test such a system as soon as possible. In this thesis, different transport lines are analyzed. The first is based on a high field pulsed solenoid, some collimators and, for perfect filtering and post-acceleration, a high field high frequency compact linear accelerator, originally designed to accelerate a 30 MeV beam extracted from a cyclotron. The second one is based on a quadruplet of permanent magnetic quadrupoles: thanks to its greater simplicity and reliability, it has great interest for experiments, but the effectiveness is lower than the one based on the solenoid; in fact, the final beam intensity drops by an order of magnitude. An additional sensible decrease in intensity is verified in the third case, where the energy selection is achieved using a chicane, because of its very low efficiency for off-axis protons. The proposed schemes have all been analyzed with 3D simulations and all the significant results are presented. Future experimental work based on the outcome of this thesis can be planned and is being discussed now.

Developments for transactinide chemistry experiments behind the gas-filled separator TASCA

Topic of this thesis is the development of experiments behind the gas-filled separator TASCA(TransActinide Separator and Chemistry Apparatus) to study the chemical properties of the transactinide elements.rnIn the first part of the thesis, the electrodepositions of short-lived isotopes of ruthenium and osmium on gold electrodes were studied as model experiments for hassium. From literature it is known that the deposition potential of single atoms differs significantly from the potential predicted by the Nernst equation. This shift of the potential depends on the adsorption enthalpy of therndeposited element on the electrode material. If the adsorption on the electrode-material is favoured over the adsorption on a surface made of the same element as the deposited atom, the electrode potential is shifted to higher potentials. This phenomenon is called underpotential deposition.rnPossibilities to automatize an electro chemistry experiment behind the gas-filled separator were explored for later studies with transactinide elements.rnThe second part of this thesis is about the in-situ synthesis of transition-metal-carbonyl complexes with nuclear reaction products. Fission products of uranium-235 and californium-249 were produced at the TRIGA Mainz reactor and thermalized in a carbon-monoxide containing atmosphere. The formed volatile metal-carbonyl complexes could be transported in a gas-stream.rnFurthermore, short-lived isotopes of tungsten, rhenium, osmium, and iridium were synthesised at the linear accelerator UNILAC at GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt. The recoiling fusion products were separated from the primary beam and the transfer products in the gas-filled separator TASCA. The fusion products were stopped in the focal plane of TASCA in a recoil transfer chamber. This chamber contained a carbon-monoxide – helium gas mixture. The formed metal-carbonyl complexes could be transported in a gas stream to various experimental setups. All synthesised carbonyl complexes were identified by nuclear decay spectroscopy. Some complexes were studied with isothermal chromatography or thermochromatography methods. The chromatograms were compared with Monte Carlo Simulations to determine the adsorption enthalpyrnon silicon dioxide and on gold. These simulations based on existing codes, that were modified for the different geometries of the chromatography channels. All observed adsorption enthalpies (on silcon oxide as well as on gold) are typical for physisorption. Additionally, the thermalstability of some of the carbonyl complexes was studied. This showed that at temperatures above 200 °C therncomplexes start to decompose.rnIt was demonstrated that carbonyl-complex chemistry is a suitable method to study rutherfordium, dubnium, seaborgium, bohrium, hassium, and meitnerium. Until now, only very simple, thermally stable compounds have been synthesized in the gas-phase chemistry of the transactindes. With the synthesis of transactinide-carbonyl complexes a new compound class would be discovered. Transactinide chemistry would reach the border between inorganic and metallorganic chemistry.rnFurthermore, the in-situ synthesised carbonyl complexes would allow nuclear spectroscopy studies under low background conditions making use of chemically prepared samples.

Development and in vitro testing of liposomal gadolinium-formulations for neutron capture therapy of glioblastoma multiforme

For the improvement of current neutron capture therapy, several liposomal formulations of neutron capture agent gadolinium were developed and tested in a glioma cell model. Formulations were analyzed regarding physicochemical and biological parameters, such as size, zeta potential, uptake into cancer cells and performance under neutron irradiation. The neutron and photon dose derived from intracellular as well as extracellular Gd was calculated via Monte Carlo simulations and set in correlation with the reduction of cell survival after irradiation. To investigate the suitability of Gd as a radiosensitizer for photon radiation, cells were also irradiated with synchrotron radiation in addition to clinically used photons generated by linear accelerator.rnIrradiation with neutrons led to significantly lower survival for Gd-liposome-treated F98 and LN229 cells, compared to irradiated control cells and cells treated with non-liposomal Gd-DTPA. Correlation between Gd-content and -dose and respective cell survival displayed proportional relationship for most of the applied formulations. Photon irradiation experiments showed the proof-of-principle for the radiosensitizer approach, although the photon spectra currently used have to be optimized for higher efficiency of the radiosensitizer. In conclusion, the newly developed Gd-liposomes show great potential for the improvement of radiation treatment options for highly malignant glioblastoma.rn

Comparison of monte carlo collimator transport methods for photon treatment planning in radiotherapy

The aim of this work was a Monte Carlo (MC) based investigation of the impact of different radiation transport methods in collimators of a linear accelerator on photon beam characteristics, dose distributions, and efficiency. Thereby it is investigated if it is possible to use different simplifications in the radiation transport for some clinical situations in order to save calculation time.

[Cosmetic results after breast conserving carcinoma treatment in patients with intramammarian seromas]

INTRODUCTION: There were 373 patients irradiated after breast conserving carcinoma treatment. A planning computed tomography revealed in 97 of these patients seromas and tissue defects exceeding 2cm in diameter. The cosmetic results in those patients and the impact of seromas herein had to be evaluated. PATIENTS AND METHODS: Mean age was 59 years. A quadrant resection was performed in 17.5 percent of the patients, a segmental resection in 27.8 percent and a tumour excision in 54.6 percent. Radiation therapy was applied with the Linear accelerator and 6 MeV photons up to a total dose in the residual breast of 50 Gy followed by a boost dose to the former tumour bed of 10 Gy. A distinct evaluation and documentation of therapy related side effects and the resulting cosmesis was done in 51 patients. RESULTS: In all the examined seroma patients there were moderate acute skin reactions grade 1 to 3. As late effects in 82.3 percent scar indurations were noticed. At the skin 51 percent showed enhanced pigmentation, 68.6 percent atrophia and only 11 percent teleangiectasia. Subcutaneous fibrosis occurred in 56,9 percent of the patients, 78.4 percent of the women had breast asymmetries. In 41.2 percent there were a lymphedema subcutaneously, in 72.5 percent impaired sensibility. The overall cosmetic result documented with a 5 point score was "very good" (score 1) in 19.6 percent and "good" (score 2) in 33.3 percent, 37.3 percent were "satisfactory" (score 3) and 9.8 percent "bad" (score 4) respectively. No "very bad" results (score 5) were seen. CONCLUSIONS: The cosmetic results in the examined group of seroma and hematoma patients were inferior to those reported in the literature. We conclude that postoperative seroma and hematoma have an adverse effect on the resulting cosmesis and that their frequency and extent have to be reduced in future by the responsible surgeons.

Adjuvant therapy after resection of brain metastases. Frameless image-guided LINAC-based radiosurgery and stereotactic hypofractionated radiotherapy

BACKGROUND: Tumor bed stereotactic radiosurgery (SRS) after resection of brain metastases is a new strategy to delay or avoid whole-brain irradiation (WBRT) and its associated toxicities. This retrospective study analyzes results of frameless image-guided linear accelerator (LINAC)-based SRS and stereotactic hypofractionated radiotherapy (SHRT) as adjuvant treatment without WBRT. MATERIALS AND METHODS: Between March 2009 and February 2012, 44 resection cavities in 42 patients were treated with SRS (23 cavities) or SHRT (21 cavities). All treatments were delivered using a stereotactic LINAC. All cavities were expanded by ≥ 2 mm in all directions to create the clinical target volume (CTV). RESULTS: The median planning target volume (PTV) for SRS was 11.1 cm(3). The median dose prescribed to the PTV margin for SRS was 17 Gy. Median PTV for SHRT was 22.3 cm(3). The fractionation schemes applied were: 4 fractions of 6 Gy (5 patients), 6 fractions of 4 Gy (6 patients) and 10 fractions of 4 Gy (10 patients). Median follow-up was 9.6 months. Local control (LC) rates after 6 and 12 months were 91 and 77 %, respectively. No statistically significant differences in LC rates between SRS and SHRT treatments were observed. Distant brain control (DBC) rates at 6 and 12 months were 61 and 33 %, respectively. Overall survival (OS) at 6 and 12 months was 87 and 63.5 %, respectively, with a median OS of 15.9 months. One patient treated by SRS showed symptoms of radionecrosis, which was confirmed histologically. CONCLUSION: Frameless image-guided LINAC-based adjuvant SRS and SHRT are effective and well tolerated local treatment strategies after resection of brain metastases in patients with oligometastatic disease.

Monte Carlo based beam model using a photon MLC for modulated electron radiotherapy

PURPOSE Modulated electron radiotherapy (MERT) promises sparing of organs at risk for certain tumor sites. Any implementation of MERT treatment planning requires an accurate beam model. The aim of this work is the development of a beam model which reconstructs electron fields shaped using the Millennium photon multileaf collimator (MLC) (Varian Medical Systems, Inc., Palo Alto, CA) for a Varian linear accelerator (linac). METHODS This beam model is divided into an analytical part (two photon and two electron sources) and a Monte Carlo (MC) transport through the MLC. For dose calculation purposes the beam model has been coupled with a macro MC dose calculation algorithm. The commissioning process requires a set of measurements and precalculated MC input. The beam model has been commissioned at a source to surface distance of 70 cm for a Clinac 23EX (Varian Medical Systems, Inc., Palo Alto, CA) and a TrueBeam linac (Varian Medical Systems, Inc., Palo Alto, CA). For validation purposes, measured and calculated depth dose curves and dose profiles are compared for four different MLC shaped electron fields and all available energies. Furthermore, a measured two-dimensional dose distribution for patched segments consisting of three 18 MeV segments, three 12 MeV segments, and a 9 MeV segment is compared with corresponding dose calculations. Finally, measured and calculated two-dimensional dose distributions are compared for a circular segment encompassed with a C-shaped segment. RESULTS For 15 × 34, 5 × 5, and 2 × 2 cm(2) fields differences between water phantom measurements and calculations using the beam model coupled with the macro MC dose calculation algorithm are generally within 2% of the maximal dose value or 2 mm distance to agreement (DTA) for all electron beam energies. For a more complex MLC pattern, differences between measurements and calculations are generally within 3% of the maximal dose value or 3 mm DTA for all electron beam energies. For the two-dimensional dose comparisons, the differences between calculations and measurements are generally within 2% of the maximal dose value or 2 mm DTA. CONCLUSIONS The results of the dose comparisons suggest that the developed beam model is suitable to accurately reconstruct photon MLC shaped electron beams for a Clinac 23EX and a TrueBeam linac. Hence, in future work the beam model will be utilized to investigate the possibilities of MERT using the photon MLC to shape electron beams.