On the Inability of Markov Models to Capture Criticality in Human Mobility

We examine the non-Markovian nature of human mobility by exposing the inability of Markov models to capture criticality in human mobility. In particular, the assumed Markovian nature of mobility was used to establish a theoretical upper bound on the predictability of human mobility (expressed as a minimum error probability limit), based on temporally correlated entropy. Since its inception, this bound has been widely used and empirically validated using Markov chains. We show that recurrent-neural architectures can achieve significantly higher predictability, surpassing this widely used upper bound. In order to explain this anomaly, we shed light on several underlying assumptions in previous research works that has resulted in this bias. By evaluating the mobility predictability on real-world datasets, we show that human mobility exhibits scale-invariant long-range correlations, bearing similarity to a power-law decay. This is in contrast to the initial assumption that human mobility follows an exponential decay. This assumption of exponential decay coupled with Lempel-Ziv compression in computing Fano's inequality has led to an inaccurate estimation of the predictability upper bound. We show that this approach inflates the entropy, consequently lowering the upper bound on human mobility predictability. We finally highlight that this approach tends to overlook long-range correlations in human mobility. This explains why recurrent-neural architectures that are designed to handle long-range structural correlations surpass the previously computed upper bound on mobility predictability

Applications of Monte Carlo methods to special radiotherapeutic techniques

Monte Carlo (MC) methods are considered one of the most powerful and precise approaches to study and solve medical physics issues. They, indeed, can be applied in all the situations where to use deterministic algorithms is infeasible or impossible. Surprising improvements in computer technology have promoted a wide diffusion of this technique, giving rise to the born of several Monte Carlo codes, such as the GEANT4 toolkit. In this paper we show some of the applications we developed using GEANT4. In particular, the simulation of two different radiotherapy techniques, such as proton/ion therapy and stereotactic radiosurgery will be discussed. In the first case we show the main features of our last public version of the GEANT4 Hadrontherapy program, also discussing the issues related to the nuclear fragmentation. In the second case, we show the procedures followed for the simulation of a Gamma Knife device, in order to validate the Treatment Planning System (TPS) used for the dose computation

Alteration of p53 and Bax/ Bcl-2 ratio by fotemustine and proton irradiation

Deregulation of apoptosis commonly occurs in melanoma cells and could be a reason for resistance. The effectiveness of different treatments depends on their ability to activate this process. In this study the effects of combined treatments with fotemustine (FM) and proton irradiation on the regulators of apoptosis were analyzed. Sub-confluent HTB140 human melanoma cells were treated with FM (100, 250 µM) 24 h prior to irradiation (12, 16 Gy). Cells were irradiated in the middle of the therapeutic 62 MeV proton spread out Bragg peak. Flow cytometric analysis of apoptosis and the Western blot analysis of apoptotic regulators were performed 6 or 48 h after treatments. Percent of apoptotic nuclei increased after applied treatments, reaching the level of 4 to 41 %. Induction of apoptosis was associated with p53 and Bax up regulation and Bcl-2 down regulation. The obtained results imply that analyzed treatments induce apoptosis through the activation of the mitochondrial apoptotic pathway, with better pro-apoptotic effects achieved by combined treatments

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

Study of the 12C + 12C reaction at 62 A MeV for hadrontherapy applications

The largest uncertainty on the physical dose deposition in hadrontherapy is due to ion’s nuclear interaction through the traversed material. Today the simulation codes are not able to reproduce the fragmentation process with the required precision. To improve the knowledge of 12C fragmentation at intermediate energies we have measured at the Laboratori Nazionali del Sud in Catania production cross sections, energy spectra and angular distributions of fragments produced in 12C fragmentation on thin 12C target, at 62A MeV

Charged particle's flux measurement from PMMA irradiated by 80 MeV/u carbon ion beam

Hadrontherapy is an emerging technique in cancer therapy that uses beams of charged particles. To meet the improved capability of hadrontherapy in matching the dose release with the cancer position, new dose monitoring techniques need to be developed and introduced into clinical use. The measurement of the fluxes of the secondary particles produced by the hadron beam is of fundamental importance in the design of any dose monitoring device and is eagerly needed to tune Monte Carlo simulations. We report the measurements done with charged secondary particles produced from the interaction of a 80 MeV/u fully stripped carbon ion beam at the INFN Laboratori Nazionali del Sud, Catania, with a Poly-methyl methacrylate target. Charged secondary particles, produced at 90$\degree$ with respect to the beam axis, have been tracked with a drift chamber, while their energy and time of flight has been measured by means of a LYSO scintillator. Secondary protons have been identified exploiting the energy and time of flight information, and their emission region has been reconstructed backtracking from the drift chamber to the target. Moreover a position scan of the target indicates that the reconstructed emission region follows the movement of the expected Bragg peak position. Exploting the reconstruction of the emission region, an accuracy on the Bragg peak determination in the submillimeter range has been obtained. The measured differential production rate for protons produced with $E^{\rm Prod}_{\rm kin} >$ 83 MeV and emitted at 90$\degree$ with respect to the beam line is: $dN_{\rm P}/(dN_{\rm C}d\Omega)(E^{\rm Prod}_{\rm kin} > 83 {\rm ~MeV}, \theta=90\degree)= (2.69\pm 0.08_{\rm stat} \pm 0.12_{\rm sys})\times 10^{-4} sr^{-1}$.Comment: 13 pages, 9 figure

First full-beam PET acquisitions in proton therapy with a modular dual-head dedicated system

During particle therapy irradiation, positron emitters with half-lives ranging from 2 to 20 min are generated from nuclear processes. The half-lives are such that it is possible either to detect the positron signal in the treatment room using an in-beam positron emission tomography (PET) system, right after the irradiation, or to quickly transfer the patient to a close PET/CT scanner. Since the activity distribution is spatially correlated with the dose, it is possible to use PET imaging as an indirect method to assure the quality of the dose delivery. In this work, we present a new dedicated PET system able to operate in-beam. The PET apparatus consists in two 10 cm × 10 cm detector heads. Each detector is composed of four scintillating matrices of 23 × 23 LYSO crystals. The crystal size is 1.9 mm × 1.9 mm × 16 mm. Each scintillation matrix is read out independently with a modularized acquisition system. The distance between the two opposing detector heads was set to 20 cm. The system has very low dead time per detector area and a 3 ns coincidence window, which is capable to sustain high single count rates and to keep the random counts relatively low. This allows a new full-beam monitoring modality that includes data acquisition also while the beam is on. The PET system was tested during the irradiation at the CATANA (INFN, Catania, Italy) cyclotron-based proton therapy facility. Four acquisitions with different doses and dose rates were analysed. In all cases the random to total coincidences ratio was equal or less than 25%. For each measurement we estimated the accuracy and precision of the activity range on a set of voxel lines within an irradiated PMMA phantom. Results show that the inclusion of data acquired during the irradiation, referred to as beam-on data, improves both the precision and accuracy of the range measurement with respect to data acquired only after irradiation. Beam-on data alone are enough to give precisions better than 1 mm when at least 5 Gy are delivered

Proton range monitoring with in-beam PET: Monte Carlo activity predictions and comparison with cyclotron data

n/