Research article

A phase-field model for non-small cell lung cancer under the effects of immunotherapy


  • Formulating mathematical models that estimate tumor growth under therapy is vital for improving patient-specific treatment plans. In this context, we present our recent work on simulating non-small-scale cell lung cancer (NSCLC) in a simple, deterministic setting for two different patients receiving an immunotherapeutic treatment. At its core, our model consists of a Cahn-Hilliard-based phase-field model describing the evolution of proliferative and necrotic tumor cells. These are coupled to a simplified nutrient model that drives the growth of the proliferative cells and their decay into necrotic cells. The applied immunotherapy decreases the proliferative cell concentration. Here, we model the immunotherapeutic agent concentration in the entire lung over time by an ordinary differential equation (ODE). Finally, reaction terms provide a coupling between all these equations. By assuming spherical, symmetric tumor growth and constant nutrient inflow, we simplify this full 3D cancer simulation model to a reduced 1D model. We can then resort to patient data gathered from computed tomography (CT) scans over several years to calibrate our model. Our model covers the case in which the immunotherapy is successful and limits the tumor size, as well as the case predicting a sudden relapse, leading to exponential tumor growth. Finally, we move from the reduced model back to the full 3D cancer simulation in the lung tissue. Thereby, we demonstrate the predictive benefits that a more detailed patient-specific simulation including spatial information as a possible generalization within our framework could yield in the future.

    Citation: Andreas Wagner, Pirmin Schlicke, Marvin Fritz, Christina Kuttler, J. Tinsley Oden, Christian Schumann, Barbara Wohlmuth. A phase-field model for non-small cell lung cancer under the effects of immunotherapy[J]. Mathematical Biosciences and Engineering, 2023, 20(10): 18670-18694. doi: 10.3934/mbe.2023828

    Related Papers:

    [1] YongKyung Oh, JiIn Kwak, Sungil Kim . Time delay estimation of traffic congestion propagation due to accidents based on statistical causality. Electronic Research Archive, 2023, 31(2): 691-707. doi: 10.3934/era.2023034
    [2] Ming Wei, Shaopeng Zhang, Bo Sun . Airport passenger flow, urban development and nearby airport capacity dynamic correlation: 2006-2019 time-series data analysis for Tianjin city, China. Electronic Research Archive, 2022, 30(12): 4447-4468. doi: 10.3934/era.2022226
    [3] Yuxia Liu, Qi Zhang, Wei Xiao, Tianguang Chu . Characteristic period analysis of the Chinese stock market using successive one-sided HP filter. Electronic Research Archive, 2023, 31(10): 6120-6133. doi: 10.3934/era.2023311
    [4] Jiping Xing, Xiaohong Jiang, Yu Yuan, Wei Liu . Incorporating mobile phone data-based travel mobility analysis of metro ridership in aboveground and underground layers. Electronic Research Archive, 2024, 32(7): 4472-4494. doi: 10.3934/era.2024202
    [5] Moutian Liu, Lixia Duan . In-phase and anti-phase spikes synchronization within mixed Bursters of the pre-Bözinger complex. Electronic Research Archive, 2022, 30(3): 961-977. doi: 10.3934/era.2022050
    [6] Zhihui Wang, Yanying Yang, Lixia Duan . Dynamic mechanism of epileptic seizures induced by excitatory pyramidal neuronal population. Electronic Research Archive, 2023, 31(8): 4427-4442. doi: 10.3934/era.2023226
    [7] Li Li, Zhiguo Zhao . Inhibitory autapse with time delay induces mixed-mode oscillations related to unstable dynamical behaviors near subcritical Hopf bifurcation. Electronic Research Archive, 2022, 30(5): 1898-1917. doi: 10.3934/era.2022096
    [8] Jianmin Hou, Quansheng Liu, Hongwei Yang, Lixin Wang, Yuanhong Bi . Stability and bifurcation analyses of p53 gene regulatory network with time delay. Electronic Research Archive, 2022, 30(3): 850-873. doi: 10.3934/era.2022045
    [9] Jieqiong Yang, Panzhu Luo . Study on the spatial correlation network structure of agricultural carbon emission efficiency in China. Electronic Research Archive, 2023, 31(12): 7256-7283. doi: 10.3934/era.2023368
    [10] Ariel Leslie, Jianzhong Su . Modeling and simulation of a network of neurons regarding Glucose Transporter Deficiency induced epileptic seizures. Electronic Research Archive, 2022, 30(5): 1813-1835. doi: 10.3934/era.2022092
  • Formulating mathematical models that estimate tumor growth under therapy is vital for improving patient-specific treatment plans. In this context, we present our recent work on simulating non-small-scale cell lung cancer (NSCLC) in a simple, deterministic setting for two different patients receiving an immunotherapeutic treatment. At its core, our model consists of a Cahn-Hilliard-based phase-field model describing the evolution of proliferative and necrotic tumor cells. These are coupled to a simplified nutrient model that drives the growth of the proliferative cells and their decay into necrotic cells. The applied immunotherapy decreases the proliferative cell concentration. Here, we model the immunotherapeutic agent concentration in the entire lung over time by an ordinary differential equation (ODE). Finally, reaction terms provide a coupling between all these equations. By assuming spherical, symmetric tumor growth and constant nutrient inflow, we simplify this full 3D cancer simulation model to a reduced 1D model. We can then resort to patient data gathered from computed tomography (CT) scans over several years to calibrate our model. Our model covers the case in which the immunotherapy is successful and limits the tumor size, as well as the case predicting a sudden relapse, leading to exponential tumor growth. Finally, we move from the reduced model back to the full 3D cancer simulation in the lung tissue. Thereby, we demonstrate the predictive benefits that a more detailed patient-specific simulation including spatial information as a possible generalization within our framework could yield in the future.



    In quantitative finance, volatility is a crucial factor in pricing derivatives, dynamic hedging, and portfolio management within the financial market. For FX (foreign exchange) options, prices are quoted in volatility. Given its significance in evaluating financial derivatives, volatility has been a focal point in both academia and practice.

    For many years, academics have been studying how to model volatility. One of the most famous volatility models is the Black-Scholes model [1]. They described the dynamic of the asset price in terms of geometric Brownian model (GBM) in which volatility of the underlying asset is assumed to be a constant and the stock price has a log-normal distribution. Because of this, the Black-Scholes model has the advantage of being tractable, but has the disadvantage of not reflecting the time-varying volatility in the real market. In fact, in the real market, there is high probability that excess skewness and leptokurtosis occur as properties of the risk-neutral probability distribution unlike the log-normal distribution. Since then, to overcome these shortcomings of the Black-Scholes model, two major types of volatility models have been developed: Stochastic volatility (SV) models and local volatility models.

    In SV models, taking account of many extraordinary volatility behaviors in the financial market after 1987 Financial Crash, the existence of a nonflat implied volatility has become remarkable. Therefore, participants interested in financial transactions have taken notice of the models that can predict the movement of financial assets, noting that the volatility of an underlying asset follows a stochastic process, and then the (pure) SV model has been proposed for describing and reflecting real situations in financial markets. In fact, the Heston model (cf. Heston [2]) and the fast mean-reverting SV model given by Fouque et al. [3] have become representative SV models, which are designed to capture the phenomenon of the mean-reversion of volatility in the real market. In addition, the Hull and White model [4] adopted the instantaneous variance process as a geometric Brownian motion. The mean-reverting Ornstein-Uhlenbeck process was also used in several models as a stochastic process (see Scott [5]; Chesney and Scott [6]; Schöbel and Zhu [7]). In addition, Heston model [2] assumed that the volatility is a Cox-Ingersoll-Ross (CIR) process and the Heston model is the one of most popular stochastic models because of the tractability.

    The local volatility models were developed by Dupire [8] and Derman and Kani [9] for the continuous and discrete cases, respectively, which are also called non-parametric local volatility models. In this model, the volatility depends on the both asset price and time, addressing the importance of the correlation between the change in the underlying asset price and the randomness of volatility to price options. In addition, Cox and Ross [10] proposed the constant elasticity of variance (CEV) model as the parametric local volatility models. Especially, as seen in Tian et al. [11], the CEV model has been known to generate the U-shaped implied volatility contrary to the Black-Scholes model. However, in the CEV model mentioned by Cox and Ross [10], the movement in volatility and underlying asset price has perfect correlation, either positively or negatively, relying on the elasticity parameter. However, in the empirical study, there is definite correlation at all times, showing volatility is time-varying as shown in Ghysels et al. [12].

    Recently, considering that the local volatility models or the SV models fail to capture the empirical evidence demonstrating that the implied volatility of equity options exhibits smile and skew curves simultaneously. A model that combines these two volatility models, called a hybrid model, has been proposed by many researchers, concentrating that the mixed models were designed to have the advantages of these two type of volatility models. First, one hybrid model was developed based on the non-parametric local volatility model and the stochastic model (see Van der Stoep et al. [13]; Tian et al. [11]; Cui et al. [14]), namely stochastic local volatility models. Second, another hybrid model that combines the CEV model and the SV model has also been developed (e.g., Andersen and Piterbarg [15]; Lord et al. [16]; Choi et al. [17]; Cui et al. [18]).

    As volatility has a direct impact on pricing and hedging performance, several volatility derivatives have been introduced. In particular, these products enable investors to manage the volatility risk, that is, they are used to hedge volatility risk. Among the volatility derivatives, variance swap is the most famous [19,20]. The payoff of variance swap is calculated from the difference between the realized variance and a strike price of the contract. Starting with the studies of Carr and Madan [19] and Demeterfi et al. [21], who evaluated variance swap by replicating the swap payoff with a portfolio of European options, various studies have been introduced. Zhu and Lian [22] priced the variance swap under the Heston model. Zheng and Kwok [23] developed a probabilistic approach to evaluate the variance swap under an SV model with jumps in the underlying asset price. A closed-form pricing formula for the discrete-time variance swap under general affine GARCH (generalized autoregressive conditional heteroskedasticity) type models was introduced by Badescu et al. [20]. Issaka [24] evaluate the variance swap by assuming that the underlying assets follow the multiple SV models. Xi and Wong [25] priced the variance swap in a rough Heston model [26] and investigated the sensitivity to roughness on the variance swap.

    In addition, as another volatility derivative, there is a timer option. The payoff of the timer option is similar to that of the vanilla option, but it has a random maturity contrary to the vanilla option, whose payoff is determined by only the fixed maturity embedded in the option. The timer option expires when the accumulated realized variance of the underlying asset reaches a predetermined level. In particular, timer options were first studied by Neuberger [27] as the "mileage option", another name for timer options, in academia and launched by the Société Générale Corporate and Investment Banking (SG CIB) in 2007. Since April 2007, timer options have been traded on the market (see [28]). After the trade began, several studies have been conducted on the evaluation of time options. Bick [29] considered the continuous version of timer options, and Li [30] implemented the Monte Carlo simulation to price timer options. Saunders [31] investigated the pricing of timer options under the fast mean-reverting SV and provided the closed-form approximation formulas. Liang et al. [32] obtained closed-form pricing formulas for timer options under the 3/2 and Heston models by using the path-integral techniques. Bernard and Cui [33] developed a one-dimensional problem for pricing of a timer call option under a general SV model. In addition, they showed the empirical results in the Hull-White and Heston models. Li [34] studied perpetual timer options under general SV models. Thus, they obtained pricing formulas for the Heston model and the 3/2 model. Ma et al. [35] derived a fast approximate analytic method to evaluate timer options under the Vasicek interest rate model. Li [36] applied the Heston model to price the timer options. In addition, as mentioned in Zhang et al. [37], the perpetual timer option was examined based on the Hull-White SV model. [38] developed concise pricing formulas for timer options using a probabilistic method to explicitly derive probability densities in stochastic volatility models, including Heston and 3/2. Recently, Wang et al. [39] employed multiscale volatility models to price timer options, stressing that the factor of the volatility is assumed to follow both fast mean-reverting and slowly varying factors. In particular, [40] introduced the return timer option, a financial derivative with a random expiry triggered by the first occurrence of a return exceeding upper or lower barriers.

    In this study, we investigate the timer option prices under the CEV model. In previous studies, the pricing of the time options has been only evaluated using various stochastic volatility models. To the best of our knowledge, there has been no research about the timer option pricing based upon CEV model in which the volatility process is considered to be the SV given by Heston [2] or Fouque et al. [3]. In fact, since the model dynamics of the underlying asset for the classical timer option is a form of the SV model, including the fast-mean reverting process given by Fouque et al. [3] or CIR process described by Heston [2], the pricing procedure of the timer options based upon CEV model is very similar to that of European options under the stochastic and local volatility (SVCEV) model introduced by Choi et al. [17], whose model is driven by the fast-mean reverting volatility process and constant elasticity of variance process. According to Choi et al. [17], the SVCEV model shows some improvements over the traditional CEV model. The SVCEV model has been widely used to evaluate various contingent claims. For example, Bock et al. [41] developed the pricing formula of a European barrier option under the SVCEV model. In addition, Kim et al. [42] applied the SVCEV model to evaluate the real options. Furthermore, several extended models based on SVCEV have been developed. For instance, Choi et al. [43] introduced a multiscale hybrid model consisting of fast and slow factors and evaluated an equity-linked annuity under the multiscale hybrid model. The previous SVCEV model [17] derived an approximation solution using up to the first order correction. Recently, an extended model has been developed using up to the second order correction term [44].

    In this article, by using asymptotic analysis on the partial differential equations (PDEs) for timer options, we derive the approximation formulas of timer options with the CEV model. Later, we call the timer options under the CEV model TO-CEVs. To demonstrate that our approximated solution for the timer options has been obtained accurately and effectively, we compare the option value with a Monte-Carlo price. Moreover, in the numerical implications, we analyze the price sensitivity of the SV for the given elasticity value on the timer options. Above all else, the most remarkable result of our paper is that the standard European option under CEV may or may not be overvalued compared with timer option with CEV, depending on the elasticity value, which is contrary to the results of Sawyer [28] emphasizing that the price of the European vanilla call options can be quite overvalued compared to the price of timer call options in the region of in-the-money (ITM).

    The remainder of this paper is organized as follows: In Section 2, we model stochastic differential equations (SDEs) to price a timer option under CEV models. Next, in Section 3, we derive the approximation solution of the option value taking advantage of asymptotic analysis. In Section 4, we examine the pricing accuracy of the approximated formula for the timer options through a Monte-Carlo simulation and analyze the parameter sensitivity on the first approximated option price. Finally, Section 5 provides concluding remarks.

    In this section, we provide the stochastic dynamics for the price of a timer call option based on hybrid stochastic and local volatility. According to Choi et al. [17], we assume that the underlying asset price model, Xt, and its instantaneous variance process, Yt, satisfy the following stochastic differential equation (SDE) model under a real probability measure P:

    dXt=μXtdt+f(Yt)Xθ/2tdWXt, (2.1)
    dYt=α(mYt)dt+βdWYt, (2.2)

    where μ is the expected return rate of the underlying asset process, Xt, f() is a smooth function satisfying 0<c1fc2< for some constants c1 and c2, θ is an elasticity parameter, α and β are positive constants, and m is the long-run mean level of Y. In addition, WXt and WYt are standard Brownian motions such that dWX,WYt=ρdt, where |ρ|1.

    Herein, based on the fact that the volatility of the underlying asset has been demonstrated to be mean-reverting through data-driven study [45], we construct the model of the volatility as a fast mean-reverting Ornstein-Uhlenbeck process, as depicted in SDE (2.2). Since the ergodic process Yt is explicitly presented by Yt=m+(Y0m)eαt+βt0eα(ts)dWYt, the process Yt follows the normal distribution. i.e., YN(m+(Y0m)eαt,β22α(1e2αt)). Interestingly, as t goes to , its invariant distribution is given by N(m,β22α), which is independent of the initial volatility level, Y0. Nevertheless, empirical studies grounded in financial market data (refer to Choi et al. [17] and Choi et al. [44]) indicate that a (pure) stochastic volatility process (in the case of θ=2) may not accurately capture the geometric properties of the implied volatility, specifically, the implied-volatility-smile and implied-volatility-skew phenomena. Therefore, to overcome the aforementioned limitations observed in the SV model, we consider the timer options based upon the CEV diffusion, which is very similar to that of European options under stochastic and local volatility (in brief, SVCEV), proposed by Choi et al. [17].

    The most distinctive feature of the SVCEV model proposed by Choi et al. [17] lies in the significant impact of the elasticity parameter θ. This economic parameter, when set to θ=2, represents the classical Black-Scholes market, whereas for θ>2, it depicts the phenomenon noticed in commodity markets (especially, coal, copper, or gold [46]) where an increase in commodity prices leads to a concurrent rise in volatility, known as the inverse leverage effect. On the other hand, when θ<2, it captures the leverage effect commonly observed in stock markets, where a fall in stock prices results in an increase in the underlying asset's volatility. Therefore, the SVCEV model of volatility for the underlying asset, as adopted from Choi et al. [17], can be considered as a suitable and realistic volatility model, efficiently capturing the random structures across various financial markets.

    Then, by the Girsanov theorem in Øksendal [47], under the risk-neutral measure Q, the stochastic dynamics in (2.1) and (2.2) are transformed into

    dXt=rXtdt+f(Yt)Xθ/2td¯WXt, (2.3)
    dYt={1ϵ(mYt)ν2ϵΛ(Yt)}dt+ν2ϵd¯WYt, (2.4)

    where r is risk-free interest rate, ν is the standard deviation of the invariant distribution of Y, and Λ(Yt) is the market price of volatility risk. In addition, ¯WX and ¯WY are transformed Brownian motions such that ¯WX,¯WYt=ρdt.

    We now demonstrate the timer call option and derive its PDE. Let us define the accumulated variance process at time t as

    It=t0f2(Yu)Xθ2udu.

    We note that if θ=2, then it reduces to the form of the accumulated variance process found in standard timer options. Next, if we denote a pre-determined variance budget as B, depending on the option buyer's choice, then the stopping timer or random maturity, τB, is given by

    τB=inf{t>0,It=B}. (2.5)

    In other words, the expiration of the timer option coincides with the initial moment at which the accumulated variance budget, It, reaches the predetermined variance budget level, B.

    Then, the payoff of a timer call option is max(XτBK,0), where K is the strike price. Therefore, we obtain the price of the timer call option at time t as follows:

    P(tτB,x,y,I)=EQ[er(τBtτB)max(XτBK,0)|XtτB=x,YtτB=y,ItτB=I]=EQ[erτBImax(XτBIK,0)|X0=x,Y0=y], (2.6)

    where EQ is the expectation with respect to risk-neutral measure Q.

    Now, by using the well-known Feynman-Kac formula and employing the well-known fact, as described in Li [36], that the price of timer options does not depend on the time variable t, we can transform this conditional expectation (2.6) into the following partial differential equation:

    (1ϵL0+1ϵL1+L2)P(x,y,I)=0, (2.7)
    P(x,y,B)=max(xK,0), (2.8)

    where

    L0=(my)y+ν22y2,L1=ν2Λ(y)y+ν2ρf(y)x2xy,L2=r(xx)+xθ2f2(y)I+12f2(y)xθ2x2.

    In this section, we aim to build an approximate solution utilizing the asymptotic analysis proposed by Fouque et al. [45]. Referring to Fouque et al. [45], for small parameter ϵ, we expand the original solution P(x,y,I) with respect to ϵ in the following form:

    P(x,y,I)=P0(x,y,I)+ϵP1(x,y,I)+ϵP2(x,y,I)+, (3.1)

    where P0,P1, are functions of (x,y,I) satisfying P0(x,B)=(xK)+ and Pn(x,B)=0 for all n1. Now, if we substitute the expansion in (3.1) into the singular perturbed PDE presented in (2.7), we obtain

    1ϵL0P0+1ϵ(L0P1+L1P0)+(L0P2+L1P1+L2P0)+ϵ(L0P3+L1P2+L2P1)=O(ϵ). (3.2)

    Now, we present two useful and fundamental results, known as the growth condition and the centering condition, for asymptotic analysis.

    Lemma 3.1. (Growth condition) As y, if the leading-order price P0(x,y,I) and the correction term P1(x,y,I) do not grow as much as P0yey2/2 and P1yey2/2, respectively, then P0 and P1 are independent of the unobservable variable y.

    Proof. The proof is similar to that in Theorem 3.1 of Ha et al. [48] or Lemma 1 and Theorem 4.1 in the work of Choi et al. [17]. From (3.2), we have the following PDEs:

    L0P0=0,L0P1+L1P0=0. (3.3)

    Next, the ODE L0P0=0 in (3.3) yields

    P0(x,I)=c1(x,I)y0e(mu)22ν2du+c2(x,I)

    for some functions c1(x,I) and c2(x,I). Then, by using the assumption, as y, if the leading-order price P0 does not grow as much as P0yey2/2, then c1(x,I)=0. Therefore, P0 is independent of the unobservable y. Similarly, if we apply this assumption to L0P1+L1P0=0, the correction term P1 is independent of y.

    Lemma 3.2. (Centering condition) If the ODE L0G(y)+G=0, a Poisson equation for the function G(y), has a unique solution, then the following centering condition must be satisfied:

    GRG(y)Φ(y)dy=0,

    where Φ(y)=12πνe(ym)22ν2 is the invariant density distribution of the ergodic process Y.

    Proof. Refer to Section 5.2 in Fouque et al. [45].

    Proposition 3.1. Under the growth condition presented in Lemma 3.1, the leading-order price P0(x,I) satisfies the following PDE:

    LTOCEV(ˆσ)P0(x,I)=0, (3.4)
    P0(x,B)=(xK)+, (3.5)

    where LTOCEV is given by

    LTOCEV=ˆσ2xθ2I+ˆσ22xθ2x2+r(xx),ˆσ=f2.

    Proof. According to Lemma 3.2, since L0P2+L2P0=0 is a Poisson equation of P2 with respect to y and P0 is independent of the variable y, the following centering condition must be hold:

    L2P0=L2P0=0,

    where

    L2=ˆσ2xθ2I+ˆσ22xθ2x2+r(xx)LTOCEV(ˆσ)

    and the effective volatility ˆσ is given by

    ˆσ=f2.

    Therefore, we can obtain the PDE (3.4) with (3.5).

    Now, we provide the analytic form solution to the PDE presented in Proposition 3.1.

    Theorem 3.1. The leading-order price, P0(x,I), is given by

    P0(x,I)=erBIˆσ2xθ2x+˜K(˜xs)12(2θ)e(˜x+s)B12θ(2˜xs)ds+erBIˆσ2sθ2K+˜K(s˜x)12(2θ)e(˜x+s)B12θ(2˜xs)ds,˜x=2xe(2θ)BIˆσ2xθ2(2θ)2χ,χ=ˆσ2(2θ)r(er(2θ)Bˆσ2sθ2er(2θ)Iˆσ2sθ2),˜K=2K2θ(2θ)2χ, (3.6)

    where the modified Bessel function of the first kind of order q, denoted as Bq(x), is given by

    Bq(x)=n=0(x2)2n+qr!Γ(n+1+q).

    Proof. If we define the new state variable ξ as ξ=Iˆσ2xθ2, then P0I=1ˆσ2xθ2P0ξ holds. Therefore, the PDE problem in (3.4) and (3.5) can be transformed into

    P0ξ+ˆσ22xθ2P0x2+r(xP0xP0)=0, (3.7)
    P0(x,Bˆσ2xθ2)=max(xK,0). (3.8)

    Surprisingly, this PDE is related to the pricing formula of a European call option under the CEV diffusion, except that the current time t and expiry date T are now replaced by the initial accumulated variance budget ξ and variance budget B, respectively. Therefore, referring to Lipton [49], the solution of PDE (3.7) with boundary condition (3.8) is given by

    P0(x,I)=er(Bˆσ2xθ2ξ)x+˜K(˜xs)12(2θ)e(˜x+s)B12θ(2˜xs)ds+er(Bˆσ2xθ2ξ)K+˜K(s˜x)12(2θ)e(˜x+s)B12θ(2˜xs)ds,

    where

    ˜x=2xe(2θ)(Bˆσ2xθ2ξ)(2θ)2χ,χ=ˆσ2(2θ)r(er(2θ)Bˆσ2xθ2er(2θ)ξ)and˜K=2K2θ(2θ)2χ.

    Finally, if we replace ξ with I and use the change of variables ξ=Iˆσ2xθ2, then we can obtain the leading-order price P0 as in (3.6).

    Next, we continue our asymptotic approaches to obtain an explicit expression for correction term P1, which plays a key role in this research and will reflect the impact of stochastic volatility on timer prices.

    Proposition 3.2. Under the growth condition given in Lemma 3.1, the first-order correction term P1(x,I) satisfies the following non-homogeneous PDE:

    LTOCEV(ˆσ)P1(x,I)=H(x,I;θ), (3.9)
    P1(x,B)=0, (3.10)

    where the non-homogeneous term H(x,I;θ) is given by

    H(x,I)=ν2Λϕxθ2P0I+νρ2fϕxx(xθ2P0I)+12{ν2Λϕxθ2P0x2+νρ2fϕxx(xθ2P0x2)}.

    Proof. In (3.2), the fourth term (or ϵ-order term) yields

    L0P3+L1P2+L2P1=0. (3.11)

    Now, by using the centering condition for the Poisson equation P3, we can obtain L1P2+L2P1=0. Then, since L2P0=0, we deduce

    L2P0=L2P0L2P0=(f2(y)ˆσ2)xθ2P0I+12(f2(y)ˆσ2)xθ2P0x2.

    In addition, from the relation L0P2+L2P0=0, P2 satisfies

    P2=L10L2P0=(ϕ(y)+C)xθ2P0I12(ϕ(y)+C)xθ2P0x2,forCR,

    where the function ϕ(y), at most polynomially growing (refer to Lemma 3.1 in Fouque et al. [3]), is solution to the following Poisson equation:

    L0ϕ(y)=f2(y)ˆσ2.

    Then, L1P2+L2P1=0 becomes

    L2P1=L1P2=L1ϕ(y)xθ2P0I+12L1ϕ(y)xθ2P0x2=ν2Λϕxθ2P0I+νρ2fϕxx(xθ2P0I)+12{ν2Λϕxθ2P0x2+νρ2fϕxx(xθ2P0x2)}H(x,I;θ).

    Therefore, combining the results L2P1=H(x,I;θ) and the boundary condition P1(x,B)=0, we complete the proof.

    Next, we present the first-order correction term P1(x,I), which is the solution to the PDE given in Proposition 3.2.

    Theorem 3.2. The first-order correction term P1(x,I) is given by

    P1(x,I)=RBˆσ2xθ2Iˆσ2xθ2er(sIˆσ2xθ2)H(x,s;θ)D(x,s;x,I)dsdx, (3.12)

    where the transition probability density function, D(x,s;x,I), is described by

    D(x,s;x,I)=(2θ)k12θ(u1u1θ2)12θeu1u2,k=2rˆσ2(2θ)(er(2θ)(sI)1),u1=kx2θer(2θ)(sI),u2=(x)2θk.

    Proof. As shown in Theorem 3.1, if we define the new state variable ξ as ξ=Iˆσ2xθ2, then P1I=1ˆσ2xθ2P1ξ holds. Therefore, the PDE problem in (3.9) and (3.10) can be transformed into

    P1ξ+ˆσ22xθ2P1x2+r(xP1xP1)=0, (3.13)
    P1(x,Bˆσ2xθ2)=0. (3.14)

    Now, using the Feynman-Kac formula and referring to Choi et al. [17], the solution to the PDE (3.13) with the source term, H(x,I;θ), as defined in Proposition 3.2, and with a zero boundary condition (3.14), is given by

    P1(x,ξ)=EQ[Bˆσ2xθ2ξer(sξ)H(x,s)ds|Xt=x].

    Finally, if we replace ξ with I, then we can obtain the analytic-form solution P1(x,I) as described in (3.12).

    This section investigates the error of the approximation of the timer call option price P(x,y,I) given by (2.6) under the CEV model and examines the price sensitivities on the TO-CEV with regard to the model parameters.

    If we combine the results of leading-order term P0(x,I) presented in Theorem 3.1 and the correction term P1(x,I) described in Theorem 3.2, then we can obtain a first-order approximation to the TO-CEV price P(x,y,I), which is denoted by ˜Pϵ, given by

    P(x,y,I)P0(x,I)+ϵP1(x,I)˜Pϵ(x,I). (4.1)

    Then, one may want to examine the accuracy of the above approximation ˜Pϵ(x,I). In this regard, by referring to Fouque and Lorig [50], if we assume that the payoff function and its derivatives are continuously differentiable and bounded, utilizing the asymptotic expansions methods, then the accuracy of ˜Pϵ(x,I) is given by

    |P(x,y,I)˜Pϵ(x,I)|O(ϵ). (4.2)

    However, for a non-smooth payoff, that is, the payoff function is not continuously differentiable, it is required to introduce a regularization argument because of the singularity at x=K. In this research, we numerically investigate the accuracy of the analytic form solution to the PDE (2.7) via Monte-Carlo simulations (refer to Table 1 in Section 4.2.1), instead of providing the mathematical proof for the accuracy.

    Table 1.  For given θ=1.8,1.9,2.1 and θ=2.2, solutions ˜Pϵ and errors are given against ϵ, where "R.Error" is the relative percentage error. Note that the baseline parameters are given by X0=100,T=1.0,r=0.03,K=100,B=0.3,v=0.05,Λ(y)ψ(y)=0.0141,f(y)ψ(y)=0.0023,ρ=0.1,m=0.5,u=1.0,f(y)=ey and f(y)2=0.0272.
    θ=1.8 θ=1.9
    Simul. ϵ 0.01 0.001 0.0001 0.01 0.001 0.0001
    ˜Pϵ 17.0555 17.5684 17.7306 21.3202 21.8599 22.0306
    10000 31.1461 18.5245 17.0744 25.8019 21.0699 22.6008
    R.Error(%) 45.24 5.16 3.84 17.36 3.74 2.52
    50000 31.8303 18.3064 18.1076 25.7835 21.1684 21.7066
    R.Error(%) 46.41 4.03 2.08 17.31 3.26 1.49
    100000 31.9046 18.6246 17.7975 25.8621 21.3048 22.184
    R.Error(%) 46.54 5.67 0.37 17.56 2.61 0.69
    θ=2.1 θ=2.2
    Simul. ϵ 0.01 0.001 0.0001 0.01 0.001 0.0001
    ˜Pϵ 31.0636 31.9086 32.1757 35.8661 36.5744 36.7984
    10000 37.2880 31.5311 33.0563 44.0128 36.9723 37.2356
    R.Error(%) 16.69 1.19 2.66 18.51 1.07 1.17
    50000 37.5411 32.9560 32.8798 44.2231 36.8926 37.1461
    R.Error(%) 17.25 3.17 2.14 18.89 0.86 0.93
    100000 37.9496 32.5929 32.4217 44.7175 37.0560 36.6976
    R.Error(%) 18.14 2.09 0.75 19.79 1.29 0.27

     | Show Table
    DownLoad: CSV

    In this section, we provide numerical implications by using the asymptotic approximations (4.1) of the TO-CEV prices. In addition, we compare the option price calculated by the closed-form with that obtained from the Monte-Carlo method. The parameters we have chosen in the numerical experiments are described by X0=100,T=1.0,r=0.03,K=100,B=0.3,v=0.05,Λ(y)ψ(y)=0.0141,f(y)ψ(y)=0.0023,ρ=0.1,m=0.5,u=1.0,f(y)=ey and f(y)2=0.0272, referring to Ha et al. [48].

    We show the accuracy of the model formula by comparing the price of the timer option, presented by the model formula and the option value through Monte Carlo (MC) simulation. Table 1 demonstrates a comparison of the results. In the table, PMC represents the Monte Carlo price and ˜Pϵ is the approximate option price given by (4.2). In addition, we calculate the relative percentage error given by 100×|PMC˜Pϵ(x,I)PMC| for the numerical implications. According to Table 1, by performing simulations 10000, 50000, or 100000 times, we obtain stable MC simulation results. Furthermore, we confirm that the value of relative error approaches zero as ϵ goes to zero. The average CPU times (in seconds) required for the execution of PMC are 27.5863, 4894.1, and 34745 for simulations conducted 10000, 50000, and 100000 times, respectively. In contrast, the CPU time for ˜Pϵ is 0.6485. Consequently, our options' formula verifies the excellence of the model not only with regards to the computational accuracy but also in terms of the computational efficiency.

    Figure 1 displays the changes in the correction term price P1 with respect to the accumulated variance under the variance budget for a given θ value. As seen in both figures, as the value of the accumulated variance is less than the variance budget (B), the impact of the correction price is significant, but when the v is greater than the B, the values of the correction term are zero and have no influence on the option price. In fact, in the case of timer options, as the accumulated variance reaches the variance budget, the options get exercised. Therefore, when the aggregated variance is larger than the budget level, the effect of the stochastic volatility (SV) must be zero, and then the correction term price remains zero. Moreover, it can be seen that the greater the variance budget level, the lower the price impact of the correction term is for the accumulated variance. It implies that as the time-to-maturity of the timer options is higher (as the budget levels get bigger), the price sensitivities become smaller, which ultimately decreases the influence of the stochastic volatility (SV) on the timer option. What's more remarkable is that, as shown in Figure 1, the more the elasticity value increases, the more the price of the correction term is sensitive to the variance budget. It means that the effect of the SV on the timer option tends to rise as the elasticity parameter gets larger, and in particular, it becomes more prominent as the expiration is longer for the larger elasticity value.

    Figure 1.  Value of the correction term Pϵ1 in terms of the accumulated variance v for the given θ. Note: X0=100,T=1.0,r=0.03,K=100,Λ(y)ψ(y)=0.0141, f(y)ψ(y)=0.0023, ρ=0.1,ϵ=0.001,m=log(0.1),u=1.0,f(y)=exp(y), f(y)2=0.0272.

    Figure 2 shows the changes in the correction term P1 of the timer option with constant elasticity of variance (CEV) with regard to the underlying asset for each strike price and for a given elasticity parameter. As seen in Figure 2(a)(d), the correction prices have a tendency to decrease sharply as the underlying asset is near the strike price (K) regardless of the value of the elasticity θ. It implies that the SV has a more significant impact on the time option if the value of the underlying asset is adjacent to the strike value (K) for the elasticity parameter θ, especially, exhibiting a hump phenomenon for the graph of the correction term for θ>2.

    Figure 2.  Value of the correction term Pϵ1 with respect to the underlying asset x for the given θ. Note: B=0.03,v=0.01. The values of the remaining parameters are the same as in Figure 1.

    Figure 3 presents the comparison of the price of the timer option (TO) and the price of the timer CEV options (TO-CEVs). Taking account of the fact that the classical timer option has the elasticity θ=2, note that the price of TO-CEVs tends to decrease as the elasticity value gets bigger, and the sensitivities of the option value get larger when the elasticity decreases.

    Figure 3.  Comparison of the values of timer CEV option (TO-CEV) with respect to the underlying asset x for a given θ. Note: B=0.3,v=0.05. The values of the remaining parameters are the same as in Figure 1.

    Figure 4 exhibits the changes of the TO-CEV prices against the underlying asset for a given elasticity value, comparing them with those of the European CEV options (EO-CEVs). As shown in Figure 4(a)(f), one can observe that the price gap between the TO-CEVs and EO-CEVs depends on the value of the elasticity. We can see that, for θ>2 or θ=2, definitely, the prices of EO-CEVs tend to be bigger than that of the TO-CEVs, and as the elasticity parameter increases, the phenomenon is more pronounced, demonstrating that our findings that TO-CEVs are noticeably underestimated compared with EO-CEVs are consistent with those of Sawyer [28]. However, in the case of θ<2, the price difference between TO-CEVs and EO-CEVs becomes smaller as the elasticity value grows down-after all, the TO-CEVs tend to be overpriced compared to the EO-CEVs in terms of underlying asset price, contrary to the results of Sawyer [28]. In particular, as seen in Figure 4(c), the EO-CEVs and the TO-CEVs have almost the same values in the region of deep ITM. In fact, referring to the empirical evidence verifying that the elasticity value calibrated from the historical data analysis of the volatility of the S&P 500 index is between 1.9 and 2, given by Choi et al. [17] and Figure 4(c), (d), it can be observed that the price of the EO-CEVs tend to get larger than that of the TO-CEVs in the real financial situation, especially, in the area of ITM. It implies that, even in the case of θ<2, if we take account of the elasticities obtained in the real situation, the value of TO-CEVs becomes undervalued compared to that of EO-CEVs in the region of ITM, similar to the findings described by Sawyer [28].

    Figure 4.  Value of Vanilla and timer option prices by CEV models with respect to the underlying asset x for the given θ. Note: B=0.022,v=0.02. The values of the remaining parameters are the same as in Figure 1.

    In this paper, we investigated the valuing problem of timer options under the constant elasticity of variance (CEV). In the real financial market, timer options are one of the financial derivatives that take account of the level of volatility so that the options are expired as the realized variance arrives at the variance budget. The purpose of this research is to analyze the price sensitivities between the timer options with the CEV (TO-CEVs) and the standard European CEV options (EO-CEVs).

    First, under the partial differential equations (PDEs) obtained from the underlying asset model and Feynman-Kac formula, we derive the approximated analytic solutions for the options by making use of the method of asymptotic analysis given by Fouque et al. [3]. Second, we implemented the numerical experiments to verify the pricing accuracy of our analytic-form formulas for the timer CEV options, compared to Monte Carlo prices. Third, we observed the pricing sensitivities of the options in terms of some model parameters, emphasizing the impact of the SV on the timer options in terms of the elasticity parameter. Fourth, we found that the timer CEV options tend to be underpriced compared with the standard European CEV options for θ>2 or the elasticity values that reflect the real market situation, in line with the results of [28]. Finally, our research can be more widely applied to several challengeable volatility managing problems or risk management strategies, including the problem for the timer options under diverse underlying asset models.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The work of S. Y. Choi was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2021R1F1A1046138) and the Gachon University research fund of 2023 (2023039990001). The research by J. H. Yoon was supported by the the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2022R1A5A1033624 and No. 2023R1A2C1006600).

    The authors have no conflicts of interest to declare.



    [1] R. C. Rockne, J. G. Scott, Introduction to mathematical oncology, JCO Clin. Cancer Inform., 3 (2019), 1–4. https://doi.org/10.1200/CCI.19.00010 doi: 10.1200/CCI.19.00010
    [2] R. A. Weinberg, The Biology of Cancer, W.W. Norton & Company, (2006). https://doi.org/10.1201/9780203852569
    [3] T. A. Graham, A. Sottoriva, Measuring cancer evolution from the genome, J. Pathol., 241 (2017), 183–191. https://doi.org/10.1002/path.4821 doi: 10.1002/path.4821
    [4] D. Hanahan, R. A. Weinberg, Hallmarks of cancer: The next generation, Cell, 144 (2011), 646–674. https://doi.org/10.1016/j.cell.2011.02.013 doi: 10.1016/j.cell.2011.02.013
    [5] R. D. Schreiber, L. J. Old, M. J. Smyth, Cancer immunoediting: Integrating immunity's roles in cancer suppression and promotion, Science, 331 (2011), 1565–1570. https://doi.org/10.1126/science.1203486 doi: 10.1126/science.1203486
    [6] Y. Zhang, Z. Zhang, The history and advances in cancer immunotherapy: Understanding the characteristics of tumor-infiltrating immune cells and their therapeutic implications, Cell. Mol. Immunol., 17 (2020), 807–821. https://doi.org/10.1038/s41423-020-0488-6 doi: 10.1038/s41423-020-0488-6
    [7] A. Rounds, J. Kolesar, Nivolumab for second-line treatment of metastatic squamous non-small-cell lung cancer, Am. J. Health-Syst. Pharm., 72 (2015), 1851–1855. https://doi.org/10.2146/ajhp150235 doi: 10.2146/ajhp150235
    [8] G. M. Keating, Nivolumab: A review in advanced nonsquamous non-small cell lung cancer, Drugs, 76 (2016), 969–978. https://doi.org/10.1007/s40265-016-0589-9 doi: 10.1007/s40265-016-0589-9
    [9] Y. Iwai, J. Hamanishi, K. Chamoto, T. Honjo, Cancer immunotherapies targeting the PD-1 signaling pathway, J. Biomed. Sci., 24 (2017), 26. https://doi.org/10.1186/s12929-017-0329-9 doi: 10.1186/s12929-017-0329-9
    [10] N. Ghaffari Laleh, C. M. L. Loeffler, J. Grajek, K. Staňková, A. T. Pearson, H. S. Muti, et al., Classical mathematical models for prediction of response to chemotherapy and immunotherapy, PLOS Comput. Biol., 18 (2022), 1–18. https://doi.org/10.1371/journal.pcbi.1009822 doi: 10.1371/journal.pcbi.1009822
    [11] I. Ezhov, K. Scibilia, K. Franitza, F. Steinbauer, S. Shit, L. Zimmer, et al., Learn-Morph-Infer: A new way of solving the inverse problem for brain tumor modeling, Med. Image Anal., 83 (2023), 102672. https://doi.org/10.1016/j.media.2022.102672 doi: 10.1016/j.media.2022.102672
    [12] A. K. Laird, Dynamics of tumour growth: Comparison of growth rates and extrapolation of growth curve to one cell, Br. J. Cancer, 19 (1965), 278–291. https://doi.org/10.1038/bjc.1965.32 doi: 10.1038/bjc.1965.32
    [13] L. Norton, A Gompertzian model of human breast cancer growth, Cancer Res., 48 (1988), 7067–7071.
    [14] S. Benzekry, C. Lamont, A. Beheshti, A. Tracz, J. M. L. Ebos, L. Hlatky, et al., Classical mathematical models for description and prediction of experimental tumor growth, PLoS Comput. Biol., 10 (2014), e1003800. https://doi.org/10.1371/journal.pcbi.1003800 doi: 10.1371/journal.pcbi.1003800
    [15] M. Bilous, C. Serdjebi, A. Boyer, P. Tomasini, C. Pouypoudat, D. Barbolosi, et al., Quantitative mathematical modeling of clinical brain metastasis dynamics in non-small cell lung cancer, Sci. Rep., 9 (2019), 13018. https://doi.org/10.1038/s41598-019-49407-3 doi: 10.1038/s41598-019-49407-3
    [16] P. Schlicke, C. Kuttler, C. Schumann, How mathematical modeling could contribute to the quantification of metastatic tumor burden under therapy: Insights in immunotherapeutic treatment of non-small cell lung cancer, Theor. Biol. Med. Model., 18 (2021), 1–15. https://doi.org/10.1186/s12976-021-00142-1 doi: 10.1186/s12976-021-00142-1
    [17] S. Benzekry, C. Sentis, C. Coze, L. Tessonnier, N. André, Development and validation of a prediction model of overall survival in high-risk neuroblastoma using mechanistic modeling of metastasis, JCO Clin. Cancer Inf., 5 (2021), 81–90. https://doi.org/10.1200/CCI.20.00092 doi: 10.1200/CCI.20.00092
    [18] S. Benzekry, P. Schlicke, P. Tomasini, E. Simon, Mechanistic modeling of brain metastases in NSCLC provides computational markers for personalized prediction of outcome, medRxiv preprint, 2023. https://doi.org/10.1101/2023.01.10.23284189
    [19] F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, A. Jemal, Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries, CA. Cancer J. Clin., 68 (2018), 394–424. https://doi.org/10.3322/caac.21492 doi: 10.3322/caac.21492
    [20] C. Zappa, S. A. Mousa, Non-small cell lung cancer: Current treatment and future advances, Transl. Lung Cancer Res., 5 (2016). https://doi.org/10.21037/tlcr.2016.06.07
    [21] W. D. Travis, E. Brambilla, A. G. Nicholson, Y. Yatabe, J. H. Austin, M. B. Beasley, et al., The 2015 world health organization classification of lung tumors: Impact of genetic, clinical and radiologic advances since the 2004 classification, J. Thorac. Oncol., 10 (2015), 1243–1260. https://doi.org/10.1097/JTO.0000000000000630 doi: 10.1097/JTO.0000000000000630
    [22] J. S. Lowengrub, H. B. Frieboes, F. Jin, Y. L. Chuang, X. Li, P. Macklin, et al., Nonlinear modelling of cancer: Bridging the gap between cells and tumours, Nonlinearity, 23 (2009). https://doi.org/10.1088/0951-7715/23/1/r01
    [23] V. Cristini, J. Lowengrub, Multiscale Modeling of Cancer: An Integrated Experimental and Mathematical Modeling Approach, Cambridge University Press, (2010). https://doi.org/10.1017/cbo9780511781452
    [24] O. Clatz, M. Sermesant, P. Y. Bondiau, H. Delingette, S. K. Warfield, G. Malandain, et al., Realistic simulation of the 3-D growth of brain tumors in MR images coupling diffusion with biomechanical deformation, IEEE Trans. Med. Imaging, 24 (2005), 1334–1346. https://doi.org/10.1109/tmi.2005.857217 doi: 10.1109/tmi.2005.857217
    [25] S. Subramanian, A. Gholami, G. Biros, Simulation of glioblastoma growth using a 3D multispecies tumor model with mass effect, J. Math. Biol., 79 (2019), 941–967. https://doi.org/10.1007/s00285-019-01383-y doi: 10.1007/s00285-019-01383-y
    [26] H. J. Bowers, E. E. Fannin, A. Thomas, J. A. Weis, Characterization of multicellular breast tumor spheroids using image data-driven biophysical mathematical modeling, Sci. Rep., 10 (2020), 1–12. https://doi.org/10.1038/s41598-020-68324-4 doi: 10.1038/s41598-020-68324-4
    [27] P. Friedl, D. Gilmour, Collective cell migration in morphogenesis, regeneration and cancer, Nat. Rev. Mol. Cell Biol., 10 (2009), 445–457. https://doi.org/10.1038/nrm2720 doi: 10.1038/nrm2720
    [28] H. Garcke, K. F. Lam, E. Sitka, V. Styles, A Cahn–Hilliard–Darcy model for tumour growth with chemotaxis and active transport, Math. Models Method Appl. Sci., 26 (2016), 1095–1148. https://doi.org/10.1142/s0218202516500263 doi: 10.1142/s0218202516500263
    [29] S. Frigeri, M. Grasselli, E. Rocca, On a diffuse interface model of tumour growth, Eur. J. Appl. Math., 26 (2015), 215–243. https://doi.org/10.1017/s0956792514000436 doi: 10.1017/s0956792514000436
    [30] H. G. Lee, Y. Kim, J. Kim, Mathematical model and its fast numerical method for the tumor growth, Math. Biosci. Eng., 12 (2015), 1173–1187. https://doi.org/10.3934/mbe.2015.12.1173 doi: 10.3934/mbe.2015.12.1173
    [31] M. Ebenbeck, H. Garcke, Analysis of a Cahn–Hilliard–Brinkman model for tumour growth with chemotaxis, J. Differ. Equations, 266 (2019), 5998–6036. https://doi.org/10.1016/j.jde.2018.10.045 doi: 10.1016/j.jde.2018.10.045
    [32] M. Ebenbeck, H. Garcke, On a Cahn–Hilliard–Brinkman model for tumor growth and its singular limits, SIAM J. Math. Anal., 51 (2019), 1868–1912. https://doi.org/10.1137/18m1228104 doi: 10.1137/18m1228104
    [33] M. Fritz, E. Lima, J. T. Oden, B. Wohlmuth, On the unsteady Darcy–Forchheimer–Brinkman equation in local and nonlocal tumor growth models, Math. Models Method Appl. Sci., 29 (2019), 1691–1731. https://doi.org/10.1142/S0218202519500325 doi: 10.1142/S0218202519500325
    [34] K. F. Lam, H. Wu, Thermodynamically consistent Navier–Stokes–Cahn–Hilliard models with mass transfer and chemotaxis, Eur. J. Appl. Math., 29 (2018), 595–644. https://doi.org/10.1017/s0956792517000298 doi: 10.1017/s0956792517000298
    [35] G. Lorenzo, A. M. Jarrett, C. T. Meyer, V. Quaranta, D. R. Tyson, T. E. Yankeelov, Identifying mechanisms driving the early response of triple negative breast cancer patients to neoadjuvant chemotherapy using a mechanistic model integrating in vitro and in vivo imaging data, arXiv preprint, (2022), arXiv: 2212.04270. https://doi.org/10.48550/arXiv.2212.04270
    [36] H. Garcke, K. F. Lam, R. Nürnberg, E. Sitka, A multiphase Cahn–Hilliard–Darcy model for tumour growth with necrosis, Math. Models Method Appl. Sci., 28 (2018), 525–577. https://doi.org/10.1142/s0218202518500148 doi: 10.1142/s0218202518500148
    [37] J. T. Oden, A. Hawkins, S. Prudhomme, General diffuse-interface theories and an approach to predictive tumor growth modeling, Math. Models Method Appl. Sci., 20 (2010), 477–517. https://doi.org/10.1142/s0218202510004313 doi: 10.1142/s0218202510004313
    [38] S. M. Wise, J. S. Lowengrub, H. B. Frieboes, V. Cristini, Three-dimensional multispecies nonlinear tumor growth – I: Model and numerical method, J. Theor. Biol., 253 (2008), 524–543. https://doi.org/10.1016/j.jtbi.2008.03.027 doi: 10.1016/j.jtbi.2008.03.027
    [39] V. Cristini, X. Li, J. S. Lowengrub, S. M. Wise, Nonlinear simulations of solid tumor growth using a mixture model: Invasion and branching, J. Math. Biol., 58 (2009), 723–763. https://doi.org/10.1007/s00285-008-0215-x doi: 10.1007/s00285-008-0215-x
    [40] H. B. Frieboes, F. Jin, Y. L. Chuang, S. M. Wise, J. S. Lowengrub, V. Cristini, Three-dimensional multispecies nonlinear tumor growth – II: Tumor invasion and angiogenesis, J. Theor. Biol., 264 (2010), 1254–1278. https://doi.org/10.1016/j.jtbi.2010.02.036 doi: 10.1016/j.jtbi.2010.02.036
    [41] E. Lima, J. T. Oden, R. Almeida, A hybrid ten-species phase-field model of tumor growth, Math. Models Method Appl. Sci., 24 (2014), 2569–2599. https://doi.org/10.1142/s0218202514500304 doi: 10.1142/s0218202514500304
    [42] A. Hawkins-Daarud, K. G. van der Zee, J. T. Oden, Numerical simulation of a thermodynamically consistent four-species tumor growth model, Int. J. Numer. Meth. Biol., 28 (2012), 3–24. https://doi.org/10.1002/cnm.1467 doi: 10.1002/cnm.1467
    [43] M. Fritz, P. K. Jha, T. Köppl, J. T. Oden, A. Wagner, B. Wohlmuth, Modeling and simulation of vascular tumors embedded in evolving capillary networks, Comput. Methods Appl. Mech. Eng., 384 (2021), 113975. https://doi.org/10.1016/j.cma.2021.113975 doi: 10.1016/j.cma.2021.113975
    [44] M. Fritz, P. K. Jha, T. Köppl, J. T. Oden, B. Wohlmuth, Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network, Nonlinear Anal. Real World Appl., 61 (2021), 103331. https://doi.org/10.1016/j.nonrwa.2021.103331 doi: 10.1016/j.nonrwa.2021.103331
    [45] G. Lorenzo, M. A. Scott, K. Tew, T. J. Hughes, Y. J. Zhang, L. Liu, et al., Tissue-scale, personalized modeling and simulation of prostate cancer growth, Proc. Natl. Acad. Sci., 113 (2016), E7663–E7671. https://doi.org/10.1073/pnas.1615791113 doi: 10.1073/pnas.1615791113
    [46] G. Song, T. Tian, X. Zhang, A mathematical model of cell-mediated immune response to tumor, Math. Biosci. Eng., 18 (2021), 373–385. https://doi.org/10.3934/mbe.2021020 doi: 10.3934/mbe.2021020
    [47] D. Kirschner, A. Tsygvintsev, On the global dynamics of a model for tumor immunotherapy, Math. Biosci. Eng., 6 (2009), 573–583. https://doi.org/10.3934/mbe.2009.6.573 doi: 10.3934/mbe.2009.6.573
    [48] K. R. Fister, J. H. Donnelly, Immunotherapy: An optimal control theory approach, Math. Biosci. Eng., 2 (2005), 499–510. https://doi.org/10.3934/mbe.2005.2.499 doi: 10.3934/mbe.2005.2.499
    [49] A. Soboleva, A. Kaznatcheev, R. Cavill, K. Schneider, K. Stankova, Polymorphic gompertzian model of cancer validated with in vitro and in vivo data, bioRxiv preprint, 2023. https://doi.org/10.1101/2023.04.19.537467
    [50] G. G. Powathil, D. J. Adamson, M. A. Chaplain, Towards predicting the response of a solid tumour to chemotherapy and radiotherapy treatments: Clinical insights from a computational model, PLoS Comput. Biol., 9 (2013), 1–14. https://doi.org/10.1371/journal.pcbi.1003120 doi: 10.1371/journal.pcbi.1003120
    [51] C. Wu, D. A. Hormuth, G. Lorenzo, A. M. Jarrett, F. Pineda, F. M. Howard, et al., Towards patient-specific optimization of neoadjuvant treatment protocols for breast cancer based on image-guided fluid dynamics, IEEE Trans. Biomed. Eng., 69 (2022), 3334–3344. https://doi.org/10.1109/tbme.2022.3168402 doi: 10.1109/tbme.2022.3168402
    [52] A. M. Jarrett, D. A. Hormuth, S. L. Barnes, X. Feng, W. Huang, T. E. Yankeelov, Incorporating drug delivery into an imaging-driven, mechanics-coupled reaction diffusion model for predicting the response of breast cancer to neoadjuvant chemotherapy: theory and preliminary clinical results, Phys. Med. Biol., 63 (2018), 105015. https://doi.org/10.1088/1361-6560/aac040 doi: 10.1088/1361-6560/aac040
    [53] R. C. Rockne, A. D. Trister, J. Jacobs, A. J. Hawkins-Daarud, M. L. Neal, K. Hendrickson, et al., A patient-specific computational model of hypoxia-modulated radiation resistance in glioblastoma using 18F-FMISO-PET, J. R. Soc. Interface, 12 (2015). https://doi.org/10.1098/rsif.2014.1174
    [54] P. Colli, H. Gomez, G. Lorenzo, G. Marinoschi, A. Reali, E. Rocca, Optimal control of cytotoxic and antiangiogenic therapies on prostate cancer growth, Math. Models Method Appl. Sci., 31 (2021), 1419–1468. https://doi.org/10.1142/s0218202521500299 doi: 10.1142/s0218202521500299
    [55] M. Fritz, C. Kuttler, M. L. Rajendran, L. Scarabosio, B. Wohlmuth, On a subdiffusive tumour growth model with fractional time derivative, IMA J. Appl. Math., 86 (2021), 688–729. https://doi.org/10.1093/imamat/hxab009 doi: 10.1093/imamat/hxab009
    [56] S. A. Quezada, K. S. Peggs, Exploiting CTLA-4, PD-1 and PD-L1 to reactivate the host immune response against cancer, Br. J. Cancer, 108 (2013), 1560–1565. https://doi.org/10.1038/bjc.2013.117 doi: 10.1038/bjc.2013.117
    [57] A. Ribas, Tumor immunotherapy directed at PD-1, N. Engl. J. Med., 366 (2012), 2517–2519. https://doi.org/10.1056/NEJMe1205943 doi: 10.1056/NEJMe1205943
    [58] D. M. Pardoll, The blockade of immune checkpoints in cancer immunotherapy, Nat. Rev. Cancer, 12 (2012), 252–264. https://doi.org/10.1038/nrc3239 doi: 10.1038/nrc3239
    [59] E. N. Rozali, S. V. Hato, B. W. Robinson, R. A. Lake, W. J. Lesterhuis, Programmed death ligand 2 in cancer-induced immune suppression, Clin. Dev. Immunol., 2012 (2012), 1–8. https://doi.org/10.1155/2012/656340 doi: 10.1155/2012/656340
    [60] S. P. Patel, R. Kurzrock, PD-L1 expression as a predictive biomarker in cancer immunotherapy, Mol. Cancer Ther., 14 (2015), 847–856. https://doi.org/10.1158/1535-7163.MCT-14-0983 doi: 10.1158/1535-7163.MCT-14-0983
    [61] Y. Viossat, R. Noble, A theoretical analysis of tumour containment, Nat. Ecol. Evol., 5 (2021), 826–835. https://doi.org/10.1038/s41559-021-01428-w doi: 10.1038/s41559-021-01428-w
    [62] T. Hillen, K. J. Painter, M. Winkler, Convergence of a cancer invasion model to a logistic chemotaxis model, Math. Models Method Appl. Sci., 23 (2013), 165–198.
    [63] B. Gompertz, On the nature of the function expressive of the law of human mortality, and on the new mode of determining the value of life contingencies, Philos. Trans. R. Soc., 115 (1825), 513–585. https://doi.org/10.1098/rstl.1825.0026 doi: 10.1098/rstl.1825.0026
    [64] K. Erbertseder, J. Reichold, B. Flemisch, P. Jenny, R. Helmig, A coupled discrete/continuum model for describing cancer-therapeutic transport in the lung, PLOS ONE, 7 (2012), 1–17. https://doi.org/10.1371/journal.pone.0031966 doi: 10.1371/journal.pone.0031966
    [65] H. Garcke, K. F. Lam, Well-posedness of a Cahn–Hilliard system modelling tumour growth with chemotaxis and active transport, Eur. J. Appl. Math., 28 (2017), 284–316. https://doi.org/10.1017/S0956792516000292 doi: 10.1017/S0956792516000292
    [66] H. Garcke, K. F. Lam, Analysis of a Cahn–Hilliard system with non-zero Dirichlet conditions modeling tumor growth with chemotaxis, Discrete Contin. Dyn. Syst. Ser. A, 37 (2017), 4277–4308. https://doi.org/10.3934/dcds.2017183 doi: 10.3934/dcds.2017183
    [67] P. Colli, H. Gomez, G. Lorenzo, G. Marinoschi, A. Reali, E. Rocca, Optimal control of cytotoxic and antiangiogenic therapies on prostate cancer growth, Math. Models Method Appl. Sci., 31 (2021), 1419–1468. https://doi.org/10.1142/s0218202521500299 doi: 10.1142/s0218202521500299
    [68] D. J. Eyre, Unconditionally gradient stable time marching the Cahn–Hilliard equation, MRS Online Proc. Lib., 529 (1998), 39–46. https://doi.org/10.1557/proc-529-39 doi: 10.1557/proc-529-39
    [69] S. C. Brenner, A. E. Diegel, L. Y. Sung, A robust solver for a mixed finite element method for the Cahn–Hilliard equation, J. Sci. Comput., 77 (2018), 1234–1249. https://doi.org/10.1007/s10915-018-0753-3 doi: 10.1007/s10915-018-0753-3
    [70] S. Balay, S. Abhyankar, M. F. Adams, S. Benson, J. Brown, P. Brune, et al., PETSc/TAO} Users Manual, Technical Report ANL-21/39 - Revision 3.18, Argonne National Laboratory, 2022.
    [71] A. Logg, K. A. Mardal, G. Wells, Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book, Springer Science & Business Media, (2012). https://doi.org/10.1007/978-3-642-23099-8
    [72] R. Kikinis, S. D. Pieper, K. G. Vosburgh, 3D Slicer: A platform for subject-specific image analysis, visualization, and clinical support, in Intraoperative Imaging and Image-guided Therapy, Springer, (2013), 277–289. https://doi.org/10.1007/978-1-4614-7657-3_19
    [73] Blender, Accessed: 2022-12-02. Available from: https://www.blender.org/.
    [74] C. Geuzaine, J. F. Remacle, Gmsh: A 3-D finite element mesh generator with built-in pre-and post-processing facilities, Int. J. Numer. Meth. Eng., 79 (2009), 1309–1331. https://doi.org/10.1002/nme.2579 doi: 10.1002/nme.2579
    [75] The Vascular Modeling Toolkit website, Accessed: 2022-12-02. www.vmtk.org.
    [76] L. Antiga, M. Piccinelli, L. Botti, B. Ene-Iordache, A. Remuzzi, D. A. Steinman, An image-based modeling framework for patient-specific computational hemodynamics, Med. Biol. Eng. Comput., 46 (2008), 1097–1112. https://doi.org/10.1007/s11517-008-0420-1 doi: 10.1007/s11517-008-0420-1
    [77] E. Eisenhauer, P. Therasse, J. Bogaerts, L. Schwartz, D. Sargent, R. Ford, et al., New response evaluation criteria in solid tumours: Revised RECIST guideline (version 1.1), Eur. J. Cancer, 45 (2009), 228–247. https://doi.org/10.1016/j.ejca.2008.10.026 doi: 10.1016/j.ejca.2008.10.026
    [78] L. Hanin, J. Rose, Suppression of metastasis by primary tumor and acceleration of metastasis following primary tumor resection: A natural law, Bull. Math. Biol., 80 (2018), 519–539. https://doi.org/10.1007/s11538-017-0388-9 doi: 10.1007/s11538-017-0388-9
  • This article has been cited by:

    1. Ping Wang, Changgui Gu, Huijie Yang, Haiying Wang, Identify causality by multi-scale structural complexity, 2024, 633, 03784371, 129398, 10.1016/j.physa.2023.129398
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2701) PDF downloads(143) Cited by(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog