Elliptic flow and R AA of D mesons at FAIR comparing the UrQMD hybrid model and the coarse-graining approach.

We present a study of the elliptic flow and R AA of D and D ¯ mesons in Au+Au collisions at FAIR energies. We propagate the charm quarks and the D mesons following a previously applied Langevin dynamics. The evolution of the background medium is modeled in two different ways: (I) we use the UrQMD hydrodynamics + Boltzmann transport hybrid approach including a phase transition to QGP and (II) with the coarse-graining approach employing also an equation of state with QGP. The latter approach has previously been used to describe di-lepton data at various energies very successfully. This comparison allows us to explore the effects of partial thermalization and viscous effects on the charm propagation. We explore the centrality dependencies of the collisions, the variation of the decoupling temperature and various hadronization parameters. We find that the initial partonic phase is responsible for the creation of most of the D / D ¯ mesons elliptic flow and that the subsequent hadronic interactions seem to play only a minor role. This indicates that D / D ¯ mesons elliptic flow is a smoking gun for a partonic phase at FAIR energies. However, the results suggest that the magnitude and the details of the elliptic flow strongly depend on the dynamics of the medium and on the hadronization procedure, which is related to the medium properties as well. Therefore, even at FAIR energies the charm quark might constitute a very useful tool to probe the quark-gluon plasma and investigate its physics.

Rapidity gap survival in enhanced Pomeron scheme.

We apply the phenomenological Reggeon field theory framework to investigate rapidity gap survival (RGS) probability for diffractive dijet production in proton-proton collisions. In particular, we study in some detail rapidity gap suppression due to elastic rescatterings of intermediate partons in the underlying parton cascades, described by enhanced (Pomeron-Pomeron interaction) diagrams. We demonstrate that such contributions play a subdominant role, compared to the usual, so-called "eikonal", rapidity gap suppression due to elastic rescatterings of constituent partons of the colliding protons. On the other hand, the overall RGS factor proves to be sensitive to color fluctuations in the proton. Hence, experimental data on diffractive dijet production can be used to constrain the respective model approaches.

Lateral variations of radiobiological properties of therapeutic fields of 1H, 4He, 12C and 16O ions studied with Geant4 and microdosimetric kinetic model.

As known, in cancer therapy with ion beams the relative biological effectiveness (RBE) of ions changes in the course of their propagation in tissues. Such changes are caused not only by increasing the linear energy transfer (LET) of beam particles with the penetration depth towards the Bragg peak, but also by nuclear reactions induced by beam nuclei leading to the production of various secondary particles. Although the changes of RBE along the beam axis have been studied quite well, much less attention has been paid to the evolution of RBE in the transverse direction, perpendicular to the beam axis. In order to fill this gap, we simulated radiation fields of 1H, 4He, 12C and 16O nuclei of 20 mm in diameter by means of a Geant4-based Monte Carlo model for heavy-ion therapy connected with the modified microdosimetric kinetic model to describe the response of normal ([Formula: see text] Gy) and early-responding ([Formula: see text] Gy) tissues. Depth and radial distributions of saturation-corrected dose-mean lineal energy, RBE and RBE-weighted dose are investigated for passive beam shaping and active beam scanning. The field of 4He has a small lateral spread as compared with 1H field, and it is characterised by a modest lateral variation of RBE suggesting the use of fixed RBE values across the field transverse cross section at each depth. Reduced uncertainties of RBE on the boundary of a 4He treatment field can be advantageous in a specific case of an organ at risk located in lateral proximity to the target volume. It is found that the lateral distributions of RBE calculated for 12C and 16O fields demonstrate fast variations in the radial direction due to changes of dose and composition of secondary fragments in the field penumbra. Nevertheless, the radiation fields of all four projectiles at radii larger than 20 mm can be characterized by a common RBE value defined by tissue radiosensitivity. These findings can help, in particular, in accessing the transverse homogeneity of radiation fields of ions used in studies in vitro.

Distributions of deposited energy and ionization clusters around ion tracks studied with Geant4 toolkit.

The Geant4-based Monte Carlo model for Heavy-Ion Therapy (MCHIT) was extended to study the patterns of energy deposition at sub-micrometer distance from individual ion tracks. Dose distributions for low-energy (1)H, (4)He, (12)C and (16)O ions measured in several experiments are well described by the model in a broad range of radial distances, from 0.5 to 3000 nm. Despite the fact that such distributions are characterized by long tails, a dominant fraction of deposited energy (∼80%) is confined within a radius of about 10 nm. The probability distributions of clustered ionization events in nanoscale volumes of water traversed by (1)H, (2)H, (4)He, (6)Li, (7)Li, and (12)C ions are also calculated. A good agreement of calculated ionization cluster-size distributions with the corresponding experimental data suggests that the extended MCHIT can be used to characterize stochastic processes of energy deposition to sensitive cellular structures.

Spatio-Temporal Dynamics of Hypoxia during Radiotherapy.

Tumour hypoxia plays a pivotal role in cancer therapy for most therapeutic approaches from radiotherapy to immunotherapy. The detailed and accurate knowledge of the oxygen distribution in a tumour is necessary in order to determine the right treatment strategy. Still, due to the limited spatial and temporal resolution of imaging methods as well as lacking fundamental understanding of internal oxygenation dynamics in tumours, the precise oxygen distribution map is rarely available for treatment planing. We employ an agent-based in silico tumour spheroid model in order to study the complex, localized and fast oxygen dynamics in tumour micro-regions which are induced by radiotherapy. A lattice-free, 3D, agent-based approach for cell representation is coupled with a high-resolution diffusion solver that includes a tissue density-dependent diffusion coefficient. This allows us to assess the space- and time-resolved reoxygenation response of a small subvolume of tumour tissue in response to radiotherapy. In response to irradiation the tumour nodule exhibits characteristic reoxygenation and re-depletion dynamics which we resolve with high spatio-temporal resolution. The reoxygenation follows specific timings, which should be respected in treatment in order to maximise the use of the oxygen enhancement effects. Oxygen dynamics within the tumour create windows of opportunity for the use of adjuvant chemotherapeutica and hypoxia-activated drugs. Overall, we show that by using modelling it is possible to follow the oxygenation dynamics beyond common resolution limits and predict beneficial strategies for therapy and in vitro verification. Models of cell cycle and oxygen dynamics in tumours should in the future be combined with imaging techniques, to allow for a systematic experimental study of possible improved schedules and to ultimately extend the reach of oxygenation monitoring available in clinical treatment.

Comparative study of dose distributions and cell survival fractions for 1H, 4He, 12C and 16O beams using Geant4 and Microdosimetric Kinetic model.

Depth and radial dose profiles for therapeutic (1)H, (4)He, (12)C and (16)O beams are calculated using the Geant4-based Monte Carlo model for Heavy-Ion Therapy (MCHIT). (4)He and (16)O ions are presented as alternative options to (1)H and (12)C broadly used for ion-beam cancer therapy. Biological dose profiles and survival fractions of cells are estimated using the modified Microdosimetric Kinetic model. Depth distributions of cell survival of healthy tissues, assuming 10% and 50% survival of tumor cells, are calculated for 6 cm SOBPs at two tumor depths and for different tissues radiosensitivities. It is found that the optimal ion choice depends on (i) depth of the tumor, (ii) dose levels and (iii) the contrast of radiosensitivities of tumor and surrounding healthy tissues. Our results indicate that (12)C and (16)O ions are more appropriate to spare healthy tissues in the case of a more radioresistant tumor at moderate depths. On the other hand, a sensitive tumor surrounded by more resistant tissues can be better treated with (1)H and (4)He ions. In general, (4)He beam is found to be a good candidate for therapy. It better spares healthy tissues in all considered cases compared to (1)H. Besides, the dose conformation is improved for deep-seated tumors compared to (1)H, and the damage to surrounding healthy tissues is reduced compared to heavier ions due to the lower impact of nuclear fragmentation. No definite advantages of (16)O with respect to (12)C ions are found in this study.

In silico analysis of cell cycle synchronisation effects in radiotherapy of tumour spheroids.

Tumour cells show a varying susceptibility to radiation damage as a function of the current cell cycle phase. While this sensitivity is averaged out in an unperturbed tumour due to unsynchronised cell cycle progression, external stimuli such as radiation or drug doses can induce a resynchronisation of the cell cycle and consequently induce a collective development of radiosensitivity in tumours. Although this effect has been regularly described in experiments it is currently not exploited in clinical practice and thus a large potential for optimisation is missed. We present an agent-based model for three-dimensional tumour spheroid growth which has been combined with an irradiation damage and kinetics model. We predict the dynamic response of the overall tumour radiosensitivity to delivered radiation doses and describe corresponding time windows of increased or decreased radiation sensitivity. The degree of cell cycle resynchronisation in response to radiation delivery was identified as a main determinant of the transient periods of low and high radiosensitivity enhancement. A range of selected clinical fractionation schemes is examined and new triggered schedules are tested which aim to maximise the effect of the radiation-induced sensitivity enhancement. We find that the cell cycle resynchronisation can yield a strong increase in therapy effectiveness, if employed correctly. While the individual timing of sensitive periods will depend on the exact cell and radiation types, enhancement is a universal effect which is present in every tumour and accordingly should be the target of experimental investigation. Experimental observables which can be assessed non-invasively and with high spatio-temporal resolution have to be connected to the radiosensitivity enhancement in order to allow for a possible tumour-specific design of highly efficient treatment schedules based on induced cell cycle synchronisation.

Hadron formation in relativistic nuclear collisions and the QCD phase diagram.

We study particle production in ultrarelativistic nuclear collisions at CERN SPS and LHC energies and the conditions of chemical freeze-out. We determine the effect of the inelastic reactions between hadrons occurring after hadronization and before chemical freeze-out employing the ultrarelativistic quantum molecular dynamics hybrid model. The differences between the initial and the final hadronic multiplicities after the rescattering stage resemble the pattern of data deviation from the statistical equilibrium calculations. By taking these differences into account in the statistical model analysis of the data, we are able to reconstruct the original hadrochemical equilibrium points in the (T, µ(B)) plane which significantly differ from chemical freeze-out ones and closely follow the parton-hadron phase boundary recently predicted by lattice quantum chromodynamics.

Nonthermal p/π ratio at LHC as a consequence of hadronic final state interactions.

Recent LHC data on Pb+Pb reactions at sqrt[s](NN) = 2.7 TeV suggests that the p/π is incompatible with thermal models. We explore several hadron ratios (K/π, p/π, Λ/π, Ξ/π) within a hydrodynamic model with a hadronic after burner, namely the ultrarelativistic quantum molecular dynamics model 3.3, and show that the deviations can be understood as a final state effect. We propose the p/π as an observable sensitive on whether final state interactions take place or not. The measured values of the hadron ratios do then allow us to gauge the transition energy density from hydrodynamics to the Boltzmann description. We find that the data can be explained with transition energy densities of 840 ± 150 MeV/fm(3).

Nanolesions induced by heavy ions in human tissues: Experimental and theoretical studies.

The biological effects of energetic heavy ions are attracting increasing interest for their applications in cancer therapy and protection against space radiation. The cascade of events leading to cell death or late effects starts from stochastic energy deposition on the nanometer scale and the corresponding lesions in biological molecules, primarily DNA. We have developed experimental techniques to visualize DNA nanolesions induced by heavy ions. Nanolesions appear in cells as "streaks" which can be visualized by using different DNA repair markers. We have studied the kinetics of repair of these "streaks" also with respect to the chromatin conformation. Initial steps in the modeling of the energy deposition patterns at the micrometer and nanometer scale were made with MCHIT and TRAX models, respectively.