Towards a Field Theoretical Stochastic Model for Description of Tumour Growth

We develop a field theory-inspired stochastic model for description of tumour growth based on an analogy with an SI epidemic model, where the susceptible individuals (S) would represent the healthy cells and the infected ones (I), the cancer cells. From this model, we obtain a curve describing the tumour volume as a function of time, which can be compared to available experimental data.


Introduction
In a recent letter [1], we have shown how the standard field theoretical language based on creation and annihilation operators (building blocks of the second quantization method [2]) may be used for a straightforward derivation of closed master equations [3] to describe the population dynamics of multivariate stochastic epidemic models. This was mainly motivated by the observation that, as remarked in [4], for the kinds of model studied in population biology and epidemiology, a field theoretical description is notationally neater and more manageable than standard methods, in often replacing sets of equations with single equations with the same content. Indeed, a single hamiltonian function sums up the dynamics compactly, even when births and deaths allow the population size to change, and it may be easily written down from a verbal description of the transitions presented in these models.
In the present work, we employ the very same methodology established in [1] to develop a field-theory inspired stochastic model for description of tumour growth based on an analogy with an SI epidemic model [5], where the susceptible individuals (S) would represent the healthy cells and the infected ones (I), the cancer cells. From this model, we were able to obtain a curve to describe the time evolution of the tumour volume, which is then compared to available experimental data for liver cancer. Our main motivation comes from the pioneering ideas about cancer as a phase transition presented in [6], more specifically, from the following observations: (1) progression of cancer must also involve a population-level shift due to a competition between two co-existing phenotypes: normal and cancerous ; (2) given enough time and resources, cancer cells will usually outcompete healthy cells in the organ or tissue where they coexist, in the competition for space and resources.
The rest of this work is organized as follows. In Section 2, we introduce the basic aspects of our model which allow us to obtain differential equations to describe the time evolution of the mean number of individuals in the interacting populations we are dealing with, i.e., normal (healthy) and cancer cells. The analytical solutions for these differential equations, besides the specific curve which describes the time evolution of the tumour volume and its comparison to experimental data are presented in Section 3. Finally, in Section 4 we present our concluding remarks.

Building the Model
We will start by considering interacting populations, whose total sizes are allowed to change, composed of two types of individuals: the normal (healthy) Since the configurations are given entirely in terms of occupation numbers ( ) , n m , which calls for a representation in terms of second-quantized bosonic operators [7], we will introduce creation and annihilation operators for the normal cells, respectively, † h and h , and for the cancer cells, namely, † c and c . These operators must satisfy the following commutation relations † † , As usual in the second quantization framework, we say that † h and † c "create", respectively, normal and cancer cells when applied over the reference (vacuum) state 0, 0 . This allows us to build our space from basis vectors of the This vacuum state has the following properties: 0, 0 0, 0 0 h c = = (from which "annihilation" operators) and 0, 0 0, 0 1 = (inner product).
Following the above definitions, we also have † † , 1, ; , The time evolution of our system will then be generated by a linear operator  (called hamiltonian) which may be constructed directly from the transition rates present in our model according to The notational advantage of this field theoretical description is made clear at this point if we observe that, from the above definitions, the equation which represents the flux of probability between states at rates defined by our model (the so-called master equation or forward Kolmogorov differential equation [5]) takes the very compact form of an imaginary-time Schrödinger-type linear equation, namely d , , . dt We can get, after some algebra, a more common representation for the master equation by substituting the expressions for the hamiltonian, Equation (7), and the vector state, Equation (3), into Equation (8) .
and, in general .
After multiplying Equation (9) Finally, by differentiating the above equation with respect to φ and evaluating the result at 0 φ θ = = we get the following differential equation for On the other hand, if we differentiate Equation (13) with respect to θ and evaluate at 0 φ θ = = we get the following differential equation for

Comparison to Liver Cancer Data. Fitting the Parameters in the Analytical Solution for
is the initial volume of the normal tissue, The above expression is compared to experimental data for liver cancer (average tumour volume) in Figure 1. The corresponding data have been obtained through the analysis of computed tomography (CT) scans for a set of 34 patients available in [8].

Concluding Remarks
We would like to finish this work by making a few comments about Equation (18). As far as we know, this is the first time that an expression describing the time evolution of the volume of tumours 1 is obtained from basic assumptions about cancer as a phase transition and, most certainly, this is the very first model for description of tumour growth built by using standard field theoretical language commonly found in models describing fundamental interactions of elementary particles. In a future work we are going to present a qualitative 1 A comprehensive review of the most common mathematical models for description of tumour growth may be found in [10]. . This curve was obtained by fitting the experimental data available in [8] (square points in the figure), using the software OriginPro 8 [9].
analysis of the behaviour of solutions for the system of first order linear differential equations composed by (14) and (15) in the phase plane, which we believe will shed more light on how our model is indeed connected to the idea of cancer as a phase transition.
Last but not least, other possible extensions for the present model should consider the inclusion of other kinds of dependence in the birth/death rates of cells, as temperature and/or concentration gradients of toxic carcinogens, for example.