The FLUKA code for application of Monte Carlo methods to promote high precision ion beam therapy

Monte Carlo (MC) methods are increasingly being utilized to support several aspects of commissioning and clinical operation of ion beam therapy facilities. In this contribution two emerging areas of MC applications are outlined. The value of MC modeling to promote accurate treatment planning is addressed via examples of application of the FLUKA code to proton and carbon ion therapy at the Heidelberg Ion Beam Therapy Center in Heidelberg, Germany, and at the Proton Therapy Center of Massachusetts General Hospital (MGH) Boston, USA. These include generation of basic data for input into the treatment planning system (TPS) and validation of the TPS analytical pencil-beam dose computations. Moreover, we review the implementation of PET/CT (Positron-Emission-Tomography / Computed- Tomography) imaging for in-vivo verification of proton therapy at MGH. Here, MC is used to calculate irradiation-induced positron-emitter production in tissue for comparison with the +-activity measurement in order to infer indirect information on the actual dose delivery

The role of nuclear reactions in Monte Carlo calculations of absorbed and biological effective dose distributions in hadron therapy

Monte Carlo codes are rapidly spreading among hadron therapy community due to their sophisticated nuclear/electromagnetic models which allow an improved description of the complex mixed radiation field produced by nuclear reactions in therapeutic irradiation. In this contribution results obtained with the Monte Carlo code FLUKA are presented focusing on the production of secondary fragments in carbon ion interaction with water and on CT-based calculations of absorbed and biological effective dose for typical clinical situations. The results of the simulations are compared with the available experimental data and with the predictions of the GSI analytical treatment planning code TRiP

Large-Eddy Simulation of a Gas Turbine Model Combustor

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/106433/1/AIAA2013-172.pd

Study of the time and space distribution of beta+ emitters from 80 MeV/u carbon ion beam irradiation on PMMA

Proton and carbon ion therapy is an emerging technique used for the treatment of solid cancers. The monitoring of the dose delivered during such treatments and the on-line knowledge of the Bragg peak position is still a matter of research. A possible technique exploits the collinear 511\ \kilo\electronvolt photons produced by positrons annihilation from $\beta^+$ emitters created by the beam. This paper reports rate measurements of the 511\ \kilo\electronvolt photons emitted after the interactions of a 80\ \mega\electronvolt / u fully stripped carbon ion beam at the Laboratori Nazionali del Sud (LNS) of INFN, with a Poly-methyl methacrylate target. The time evolution of the $\beta^+$ rate was parametrized and the dominance of $^{11}C$ emitters over the other species ($^{13}N$, $^{15}O$, $^{14}O$) was observed, measuring the fraction of carbon ions activating $\beta^+$ emitters $A_0=(10.3\pm0.7)\cdot10^{-3}$. The average depth in the PMMA of the positron annihilation from $\beta^+$ emitters was also measured, D_{\beta^+}=5.3\pm1.1\ \milli\meter, to be compared to the expected Bragg peak depth D_{Bragg}=11.0\pm 0.5\ \milli\meter obtained from simulations

Update On the Status of the FLUKA Monte Carlo Transport Code*

The FLUKA Monte Carlo transport code is a well-known simulation tool in High Energy Physics. FLUKA is a dynamic tool in the sense that it is being continually updated and improved by the authors. We review the progress achieved since the last CHEP Conference on the physics models, some technical improvements to the code and some recent applications. From the point of view of the physics, improvements have been made with the extension of PEANUT to higher energies for p, n, pi, pbar/nbar and for nbars down to the lowest energies, the addition of the online capability to evolve radioactive products and get subsequent dose rates, upgrading of the treatment of EM interactions with the elimination of the need to separately prepare preprocessed files. A new coherent photon scattering model, an updated treatment of the photo-electric effect, an improved pair production model, new photon cross sections from the LLNL Cullen database have been implemented. In the field of nucleus-- nucleus interactions the electromagnetic dissociation of heavy ions has been added along with the extension of the interaction models for some nuclide pairs to energies below 100 MeV/A using the BME approach, as well as the development of an improved QMD model for intermediate energies. Both DPMJET 2.53 and 3 remain available along with rQMD 2.4 for heavy ion interactions above 100 MeV/A. Technical improvements include the ability to use parentheses in setting up the combinatorial geometry, the introduction of pre-processor directives in the input stream. a new random number generator with full 64 bit randomness, new routines for mathematical special functions (adapted from SLATEC). Finally, work is progressing on the deployment of a user-friendly GUI input interface as well as a CAD-like geometry creation and visualization tool. On the application front, FLUKA has been used to extensively evaluate the potential space radiation effects on astronauts for future deep space missions, the activation dose for beam target areas, dose calculations for radiation therapy as well as being adapted for use in the simulation of events in the ALICE detector at the LHC

Decay constants, semi-leptonic and non-leptonic decays in a Bethe-Salpeter Model

We evaluate the decay constants for the B and $D$ mesons and the form factors for the semileptonic decays of the B meson to $D$ and $D^*$ mesons in a Bethe-Salpeter model. From data we extract $V_{cb}=0.039 \pm 0.002$ from ${\bar B} \to D^* l {\bar{\nu}}$ and $V_{cb}=0.037 \pm 0.004$ from ${\bar B} \to D l {\bar{\nu}}$ decays. The form factors are then used to obtain non-leptonic decay partial widths for $B\to D \pi (K)$ and $B \to D D (D_s)$ in the factorization approximation.Comment: 15 Pages, 3 Postscript figures (available also from [email protected]

Excitable media in open and closed chaotic flows

We investigate the response of an excitable medium to a localized perturbation in the presence of a two-dimensional smooth chaotic flow. Two distinct types of flows are numerically considered: open and closed. For both of them three distinct regimes are found, depending on the relative strengths of the stirring and the rate of the excitable reaction. In order to clarify and understand the role of the many competing mechanisms present, simplified models of the process are introduced. They are one-dimensional baker-map models for the flow and a one-dimensional approximation for the transverse profile of the filaments.Comment: 14 pages, 16 figure

Relativistic Calculation of the Meson Spectrum: a Fully Covariant Treatment Versus Standard Treatments

A large number of treatments of the meson spectrum have been tried that consider mesons as quark - anti quark bound states. Recently, we used relativistic quantum "constraint" mechanics to introduce a fully covariant treatment defined by two coupled Dirac equations. For field-theoretic interactions, this procedure functions as a "quantum mechanical transform of Bethe-Salpeter equation". Here, we test its spectral fits against those provided by an assortment of models: Wisconsin model, Iowa State model, Brayshaw model, and the popular semi-relativistic treatment of Godfrey and Isgur. We find that the fit provided by the two-body Dirac model for the entire meson spectrum competes with the best fits to partial spectra provided by the others and does so with the smallest number of interaction functions without additional cutoff parameters necessary to make other approaches numerically tractable. We discuss the distinguishing features of our model that may account for the relative overall success of its fits. Note especially that in our approach for QCD, the resulting pion mass and associated Goldstone behavior depend sensitively on the preservation of relativistic couplings that are crucial for its success when solved nonperturbatively for the analogous two-body bound-states of QED.Comment: 75 pages, 6 figures, revised content