
Long-term stability belongs to the main problems of dye-sensitized solar cells (DSSCs), impeding their practical application. Especially the usually fluid electrolyte tends to evaporation, thus drying the cells if they are not perfectly sealed. While gelling the electrolyte with different polymers often reduces the efficiency, using a glycerol-based electrolyte was already shown to result in similar or even improved efficiency. At the same time, drying of the cells was significantly reduced. Here we report on improving glycerol-based electrolytes further by varying the iodine-triiodide ratio and the overall concentration in the electrolyte. Long-term tests with unsealed glass-based DSSCs were performed over more than 1 year, showing that most of the cells increased efficiency during this time, opposite to cells with a commercial solvent-based iodine-triiodide electrolyte which completely dried after 2–3 months.
Citation: Katrin Gossen, Marius Dotter, Bennet Brockhagen, Jan Lukas Storck, Andrea Ehrmann. Long-term investigation of unsealed DSSCs with glycerol-based electrolytes of different compositions[J]. AIMS Materials Science, 2022, 9(2): 283-296. doi: 10.3934/matersci.2022017
[1] | Nada Almuallem, Dumitru Trucu, Raluca Eftimie . Oncolytic viral therapies and the delicate balance between virus-macrophage-tumour interactions: A mathematical approach. Mathematical Biosciences and Engineering, 2021, 18(1): 764-799. doi: 10.3934/mbe.2021041 |
[2] | Philip Gerlee, Alexander R. A. Anderson . Diffusion-limited tumour growth: Simulations and analysis. Mathematical Biosciences and Engineering, 2010, 7(2): 385-400. doi: 10.3934/mbe.2010.7.385 |
[3] | Bethan Morris, Lee Curtin, Andrea Hawkins-Daarud, Matthew E. Hubbard, Ruman Rahman, Stuart J. Smith, Dorothee Auer, Nhan L. Tran, Leland S. Hu, Jennifer M. Eschbacher, Kris A. Smith, Ashley Stokes, Kristin R. Swanson, Markus R. Owen . Identifying the spatial and temporal dynamics of molecularly-distinct glioblastoma sub-populations. Mathematical Biosciences and Engineering, 2020, 17(5): 4905-4941. doi: 10.3934/mbe.2020267 |
[4] | Qingfeng Tang, Guohong Zhang . Stability and Hopf bifurcations in a competitive tumour-immune system with intrinsic recruitment delay and chemotherapy. Mathematical Biosciences and Engineering, 2021, 18(3): 1941-1965. doi: 10.3934/mbe.2021101 |
[5] | Farinaz Forouzannia, Sivabal Sivaloganathan, Mohammad Kohandel . A mathematical study of the impact of cell plasticity on tumour control probability. Mathematical Biosciences and Engineering, 2020, 17(5): 5250-5266. doi: 10.3934/mbe.2020284 |
[6] | Abdulhamed Alsisi, Raluca Eftimie, Dumitru Trucu . Non-local multiscale approach for the impact of go or grow hypothesis on tumour-viruses interactions. Mathematical Biosciences and Engineering, 2021, 18(5): 5252-5284. doi: 10.3934/mbe.2021267 |
[7] | Shihong Yin, Qifang Luo, Yanlian Du, Yongquan Zhou . DTSMA: Dominant Swarm with Adaptive T-distribution Mutation-based Slime Mould Algorithm. Mathematical Biosciences and Engineering, 2022, 19(3): 2240-2285. doi: 10.3934/mbe.2022105 |
[8] | Abdulhamed Alsisi, Raluca Eftimie, Dumitru Trucu . Nonlocal multiscale modelling of tumour-oncolytic viruses interactions within a heterogeneous fibrous/non-fibrous extracellular matrix. Mathematical Biosciences and Engineering, 2022, 19(6): 6157-6185. doi: 10.3934/mbe.2022288 |
[9] | J. Leonel Rocha, Sandra M. Aleixo . An extension of Gompertzian growth dynamics: Weibull and Fréchet models. Mathematical Biosciences and Engineering, 2013, 10(2): 379-398. doi: 10.3934/mbe.2013.10.379 |
[10] | Nara Bobko, Jorge P. Zubelli . A singularly perturbed HIV model with treatment and antigenic variation. Mathematical Biosciences and Engineering, 2015, 12(1): 1-21. doi: 10.3934/mbe.2015.12.1 |
Long-term stability belongs to the main problems of dye-sensitized solar cells (DSSCs), impeding their practical application. Especially the usually fluid electrolyte tends to evaporation, thus drying the cells if they are not perfectly sealed. While gelling the electrolyte with different polymers often reduces the efficiency, using a glycerol-based electrolyte was already shown to result in similar or even improved efficiency. At the same time, drying of the cells was significantly reduced. Here we report on improving glycerol-based electrolytes further by varying the iodine-triiodide ratio and the overall concentration in the electrolyte. Long-term tests with unsealed glass-based DSSCs were performed over more than 1 year, showing that most of the cells increased efficiency during this time, opposite to cells with a commercial solvent-based iodine-triiodide electrolyte which completely dried after 2–3 months.
The beginning of a primary solid tumour is the result of a single normal cell that is transformed as a result of mutations in certain key genes. Cells can mutate spontaneously, or mutations can be environmentally induced. Mutations occur during cell division, and most of the time the immune system can recognise mutated cells and eliminate them. When the immune system fails to eliminate cells with mutations in genes that control cell proliferation, cells become cancerous. It is known that all cancer cell lines have at least one mutation, with most cancer cell lines having more than one mutation; e.g., in [1] 137 oncogenic mutations were identified in 14 out of 24 known cancer genes in 60 human cancer cell lines. Moreover, cancer cells are genetically unstable and cells inside the solid tumours keep mutating leading to very heterogeneous tumour masses.
The mechanisms behind the mutation pressure are still not fully understood. Experimental studies have shown that some changes in the extracellular matrix (ECM) can correlate with sustained cell proliferative signalling and an increased risk of developing cancer [2]. Other studies have shown that culturing cells for long times in stiff hydrogels can lead to the subclonal selection of genomic aberrations in cells [3], thus suggesting that the ECM properties could impact the mutation status of cells in solid tumours. Other studies suggested that the maintenance of cells at high density in the absence of proliferation leads to an increase in mutagenesis following cell division [4]. Therefore, there seem to be different mechanisms that can trigger and influence the mutation rate of cells.
Mutated cells not only incur sustained proliferation (see Figure 1), but can also exhibit migrational and invasion properties [5], eventually leading to cancer metastasis. The invasion of surrounding tissue is the result of ECM degradation and remodelling by the cancer cells (which can secrete various proteolytic enzymes, such as matrix metalloproteinases (MMPs)) as well as other cells in the microenvironment. The last decades have seen the development of numerous theoretical studies based on mathematical models, which investigate computationally the biological mechanisms behind cancer cell invasion into the tissue. Most of these models are single-scale models; see [6,7,8,9,10,11] and references therein. More recently, multi-scale mathematical models have started to be developed to consider also the multi-scale aspects of various biological processes occurring during cancer invasion [12,13,14,15,16].
The main issue faced by all these single-scale and especially multi-scale mathematical models for cancer growth and invasion is parameter estimation. The last few years have seen the publication of various mathematical studies that try to estimate numerically different model parameters using inverse problem formulations [17,18,19,20,21]. These studies focus on identifying the magnitude of isotropic [21] and anisotropic tumour diffusion [18], the magnitude of tumour growth rate [19,21], the strength/location of tumour-induced tissue deformation [19], the position of the blood vessels that act as a source for the oxygen concentration that influences tumour growth [17], the location of the source of tumours [20].
To our knowledge, no studies have tried to estimate the mutation laws for cells inside heterogeneous tumours. To address this aspect, in this current study we consider a simplified problem with only two cancer cell populations (one mutating into the second one), that can exhibit both random and directed haptotactic movement. Moreover, we assume that the cells can degrade and remodel the surrounding ECM density. Furthermore, the type of interactions among cancer cells and between cells and ECM is not always very clear: some simpler mathematical models for cancer growth and invasion consider local cell-cell and cell-ECM interactions [8,9], while other more complex models consider non-local cell-cell and cell-ECM interactions [11,16,22,23]. In this study, we consider both type of models, local and non-local, and estimate the mutation laws for both cases. Finally, we assume here the knowledge of additional information in terms of both exact and noisy measurements of the tumour constituent density at some later time in the tumour evolution. We test our inversion approach on several cancer proliferation laws that are usually used in cancer modelling: logistic and Gompertz proliferation.
This paper is structured as follows. In Section 2 we describe the local mathematical model for the dynamics of the two cancer cell sub-populations and the extracellular matrix (ECM). In Sections 3 and 3.1 we formulate the inverse problem for this local model under the assumptions of logistic cell proliferation and mutation depending only on cancer cell density. In Section 3.4 we formulate the inverse problem for this local model under the assumption of logistic cell proliferation and mutation depending on both cell and ECM density. In Section 4 we describe a non-local model for two cancer cell populations and their interactions with the ECM, while in Section A we present the numerical approach for this forward nonlocal model. In Section 4.1 we formulate the inverse problem for this nonlocal model for cancer invasion. We summarise and discuss our results in Section 5.
In this study we consider two populations of tumour cells: a primary cell population c1(x,t) that can mutate into a more aggressive secondary cell population c2(x,t). These cell populations exercise a spatial redistribution via random movement (with diffusion coefficients D1 and D2) and directed movement towards extracellular matrix (ECM) gradients v(x,t) (with haptotactic coefficients η1 and η2). Moreover, the two populations undergo logistic growth (at the same rate μc), up to a carrying capacity (Kc). We also assume that the ECM undergoes degradation (at a rate ρ) and remodelling (at a rate μv). The above assumptions are described by the equations:
∂c1∂t=D1Δc1⏟diffusion−η1∇⋅(c1∇v)⏟haptotacticmovement+μcc1(1−c1+c2+vKc)⏟logisticproliferation−ω(t)⏟mutationswitchQ(⋅,⋅)⏟unknownmutation, | (2.1a) |
∂c2∂t=D2Δc2⏟diffusion−η2∇⋅(c2∇v)⏟haptotacticmovement+μcc2(1−c1+c2+vKc)⏟logisticproliferation+ω(t)⏟mutationswitchQ(⋅,⋅)⏟unknownmutation, | (2.1b) |
∂v∂t=−ρ(c1+c2)v⏟degradation+μv(Kc−v−c1−c2)+⏟remodellingterm. | (2.1c) |
where (Kc−v−c1−c2)+:=max{(Kc−v−c1−c2),0}. Finally, the unknown term Q(⋅,⋅) represents the mutation law of cell subpopulation c1 into cell subpopulation c2, which is assumed to be mediated by a time-dependant mutation enhancement ω(t) that is known a priori and is taken here of the form
ω(t):=(1+tanh(t−t1,2ts))2, |
where t1,2 is the time at which mutations from c1 to c2 start occurring, and ts>0 is a time-steepness coefficient for this mutation law.
The mutation law Q(⋅,⋅) is considerd here unknown due to either unknown dependance on the primary cell population c1, or unknown dependance on the ECM v, or unknown dependance on both primary tumour cell population and ECM. In this study we investigate three assumptions related to this mutation term, namely: (i) mutation depends linearly on the density of primary tumour cells as experimental studies have shown that tumour hypoxia, generated by high tumour cell density, is linearly correlated with an increase in genomic changes toward more aggressive tumours [24]; (ii) mutation depends linearly on the density of primary tumour cells, and nonlinearly on the ECM density [3]; (iii) mutation law is very general and depends autonomously on both the primary tumour and ECM, as suggested by various experimental studies regarding the role of the tumour microenvironment in cancer cells mutation process [3,24,25]. Thus mathematically, these cases correspond to three inverse problems that seek to identify the unknown mutation law Q(⋅,⋅) in the following three situations:
(i) Mutation depends linearly on the density of primary tumour cell sub-population c1 but does not depend at all on ECM, and so this is given by the unknown term
Q(c1,v):=˜Q1(c1)=δ0c1, | (2.2) |
with δ0 representing the unknown mutation rate.
(ii) Mutation depends in a known linear manner on c1 and in an unknown nonlinear way on the density of ECM. The unknown dependence on v is denoted mathematically by the unknown function ˜Q2(v), and so the entire mutation law is therefore of the form
Q(c1,v):=δ0c1˜Q2(v), | (2.3) |
with δ0>0 here being considered known. A usual choice for ˜Q2(v) is of the form [26,27]:
˜Q2(v):={exp(−1κ2−(1−v(x,t))2)exp(−1κ2)if1−κ<v(x,t)<1,0,otherwise, | (2.4) |
where κ>0 is a certain level of ECM beyond which mutations can occur.
(iii) Mutation is given by an unknown general nonlinear law Q(c1,v) that is autonomous in c1 and v, which will be reconstructed from the data available at a specific later time in the tumour evolution.
To complete the description of the model, we mention that the initial conditions for the cancer-ECM dynamics described by model (2.1) are:
c1(x,0):=c1,0(x),c2(x,0):=c2,0(x)andv(x,0):=v0(x),forx∈Ω. | (2.5) |
Here, c1,0(⋅), c2,0(⋅) and v0(⋅) give the initial distributions of the primary cell subpopulation, mutated cell subpopulation and ECM, respectively. Furthermore, we assume that the cells and the ECM components do not leave the tissue region Ω, and therefore we consider zero Neumann boundary conditions:
∂c1∂n|∂Ω=0,∂c2∂n|∂Ω=0and∂v∂n|∂Ω=0, | (2.6) |
where n(ξ) is the usual normal direction at any given tissue boundary point ξ∈∂Ω.
Throughout the following sections, we refer to the tumour dynamics model (2.1) together with the initial and boundary conditions (2.5) and (2.6) as the "forward model".
We start with the forward model defined by the tumour dynamics model (2.1) in the presence of the initial and boundary conditions (2.5) and (2.6). Our goal is to reconstruct the unknown cancer cells mutation law Q(⋅) from additional information given by measurements of the cancer cells and ECM densities taken at some later time tf:=T>0 in the tumour evolution. These measurements are therefore given in the form of the following functions on Ω, which are considered to be known:
c∗1(⋅):Ω→Rforthecancersubpopulationc1, | (3.1a) |
c∗2(⋅):Ω→Rforthecancersubpopulationc2, | (3.1b) |
v∗(⋅):Ω→RfortheECMdensity. | (3.1c) |
In the following, we will explore the reconstruction of the unknown cancer mutation law Q(⋅) when the known measurements c∗1(x), c∗2(x) and v∗(x) will be given both as exact (accurate) data and as noisy data, ∀x∈Ω.
In this section we outline in a unitary manner the forward solver involved in the retrieval of the mutation laws corresponding to Cases (i) and (ii) that require the retrieval of ˜Q1(c1) and ˜Q2(v), respectively. To that end, for r=1,2, denoting by er either the primary tumour or the ECM, i.e., er∈{c1,v}, enables us to proceed with addressing simultaneously both cases by simply referring to the retrieval of the term compactly denoted as ˜Qr(er) that is specified by
er:={c1,ifr=1,v,ifr=2,andsubsequently˜Qr(er):={˜Q1(c1),ifr=1,˜Q2(v),ifr=2. | (3.2) |
We start by considering a uniform discretisation GΩ:={(xi,yj)}i,j=1…N of step size Δx=Δy>0 for a square maximal tissue region Ω⊂R2 where the tumour exercises its dynamics. At any given time t∈[0,tf] the discretisations of cancer cells densities c1(⋅,t) and c2(⋅,t) as well as the density of ECM v(⋅,t) are therefore given by the N×N matrices ˜c1(t):={˜c1,i,j(t)}i,j=1…N, ˜c2(t):={˜c2,i,j(t)}i,j=1…N and ˜v(t):={˜vi,j(t)}i,j=1…N, with ˜c1,i,j(t):=c1((xi,yj),t), ˜c2,i,j(t):=c2((xi,yj),t) and ˜vi,j(t):=v((xi,yj),t), ∀i,j=1…N. Correspondingly, in the following ˜er(t) will denote either ˜c1(t) or ˜v(t), i.e., ˜er(t)∈{˜c1(t),˜v(t)}, as required by model (3.2).
Throughout this study we assume that we have a priori knowledge that the cumulated ECM and cancer densities do not exceed the tissue carrying capacity Kc. Under this assumption, the unknown mutation law can be written in terms of an unknown (for the moment) function mc∗1,c∗2,v∗:[0,Kc]→[0,∞). Moreover, this unknown function mc∗1,c∗2,v∗ will be appropriately identified within a suitable family of functions M1 such that the corresponding solution for the tumour model (2.1) will match the measurements given in model (3.1). Thus, denoting by ˜Qc∗1,c∗2,v∗r(⋅) the unknown mutation term for which the corresponding solution of model (2.1) matches measurement model (3.1), at each (xi,yj) we can write this as
˜Qc∗1,c∗2,v∗r(˜emc∗1,c∗2,v∗r,i,j(t)):=M1i,j(˜emc∗1,c∗2,v∗r(t),mc∗1,c∗2,v∗), |
where M1(⋅,⋅):={M1i,j(⋅,⋅)}i,j=1…N, with M1(⋅,⋅):RN×N×M1→RN×N representing a "trial mutation operator" that will be specified below alongside the family of functions M1. Furthermore, ˜emc∗1,c∗2,v∗r(t):={˜emc∗1,c∗2,v∗r,i,j(t)}i,j=1..N represents the solution for the density of either the primary cell population (if r=1) or the ECM (if r=2) that is obtained for model (2.1) when, instead of the unknown term ˜Qr(er), in the mutation law we use the trial mutation term M1(⋅,mc∗1,c∗2,v∗).
Next, we consider an uniform discretisation for the domain [0,Kc] that is given by an equally spaced grid G1M:={ηl}l=1…M of step size Δη>0. On this discretised domain, the unknown function mc∗1,c∗2,v∗ is identified through a suitable approximation within the following M−dimensional space of functions associated with G1M, namely
M1:={m:[0,Kc]→R|m|El=∑p=0,1m(ηl+p)ϕl+p,∀El∈G1,tilesM}withthefamilyofintervalsG1,tilesM:={El:=[ηl,ηl+1]|l=1…M−1},and∀El∈G1,tilesM,{ϕl+p}p=0,1describetheusuallinearshapefunctionsonEl. | (3.3) |
Thus, for any candidate function m∈M1, the corresponding trial mutation operator M1 has each of its components Mi,j, ∀i,j=1…N, given by
M1i,j(˜emr(t),m):=m|El(˜emr,i,j(t)),withindexlbeingindependentofitschoicewithintheassociatedsetofindicesΛi,j,namely:Λi,j:={l∈{1,…,M−1}|∃El∈G1,tilesMsuchthat˜emr,i,j(t)∈El}. | (3.4) |
Here, as per model (3.2), ˜emr(t):={˜emri,j(t)}i,j=1…N represents either the solutions for the density of primary cancer cell population, {˜cm1,i,j(t)}i,j=1…N, or for the density of ECM, {˜vmi,j(t)}i,j=1…N, which alongside the density of mutated cell population, {˜cm1,i,j(t)}i,j=1…N, are obtained with model (2.1) when this uses within the mutation law the trial mutation term M1(⋅,m):={M1i,j(⋅,m)}i,j=1…N given in model (3.4) instead of the unknown term ˜Qr(⋅). Finally, the trial mutation form for the full mutation law given in models (2.2) and (2.3) for Cases (i) and (ii), respectively, is denoted by ¯Mr(˜cm1(t),˜vm(t),m):={¯Mr,i,j(˜cm1(t),˜vm(t),m)}i.j=1…N and is given by:
¯Mr,i,j(˜cm1(t),˜vm(t),m):={M1i,j(˜cm1(t),m),ifr=1,δ0˜cm1,i,j(t)M1i,j(˜vm(t),m),ifr=2. | (3.5) |
Therefore, in space-discretised form, model (3.4) that uses M1(⋅,m) can therefore be written as
∂∂t[˜cm1˜cm2˜vm]=[H1(˜cm1,cm2,˜vm,m)H2(˜cm1,cm2,˜vm,m)H3(˜cm1,cm2,˜vm)], | (3.6) |
Here, H1(⋅,⋅,⋅,⋅)={H1i,j(⋅,⋅,⋅,⋅)}i,j=1…N represents the spatial discretisation corresponding to Eq (2.1a). Each of its components H1i,j(⋅,⋅,⋅,⋅), ∀i,j=1…N, are given by
H1i,j(˜cm1(t),˜cm2(t),˜vm(t),m):=D1(Δx)2(˜cm1,i−1,j(t)+˜cm1,i+1,j(t)+˜cm1,i,j−1(t)+˜cm1,i,j+1(t)−4˜cm1,i,j(t))−η12(Δx)2((˜cm1,i,j(t)+˜cm1,i+1,j(t))(˜vmi+1,j(t)−˜vmi,j(t))−(˜cm1,i,j(t)+˜cm1,i−1,j(t))(˜vmi,j(t)−˜vmi−1,j(t))+(˜cm1,i,j(t)+˜cm1,i,j+1(t))(˜vmi,j+1(t)−˜vmi,j(t))−(˜cm1,i,j(t)+˜cm1,i,j−1(t))(˜vmi,j(t)−˜vmi,j−1(t)))+μc˜cm1,i,j(t)(Kc−˜cm1,i,j(t)−˜cm2,i,j(t)−˜vmi,j(t))−ω(t)¯Mr,i,j(˜cm1(t),˜vm(t),m). | (3.7) |
Similarly, H2(⋅,⋅,⋅,⋅)={H2i,j(⋅,⋅,⋅,⋅)}i,j=1…N represents the spatial discretisation corresponding to Eq (2.1d), and each of its components H2i,j(⋅,⋅,⋅,⋅), ∀i,j=1…N, are given by
H2i,j(˜cm2(t),˜cm2(t),˜vm(t),m):=D2(Δx)2(˜cm2,i−1,j(t)+˜cm2,i+1,j(t)+˜cm2,i,j−1(t)+˜cm2,i,j+1(t)−4˜cm2,i,j(t))−η22(Δx)2((˜cm2,i,j(t)+˜cm2,i+1,j(t))(˜vmi+1,j(t)−˜vmi,j(t))−(˜cm2,i,j(t)+˜cm2,i−1,j(t))(˜vmi,j(t)−˜vmi−1,j(t))+(˜cm2,i,j(t)+˜cm2,i,j+1(t))(˜vmi,j+1(t)−˜vmi,j(t))−(˜cm2,i,j(t)+˜cm2,i,j−1(t))(˜vmi,j(t)−˜vmi,j−1(t)))+μc˜cm2,i,j(t)(Kc−˜cm1,i,j(t)−˜cm2,i,j(t)−˜vmi,j(t))+ω(t)¯Mr,i,j(˜cm1(t),˜vm(t),m). | (3.8) |
Finally, H3(⋅,⋅,⋅)={H3i,j(⋅,⋅,⋅)}i,j=1…N represents the discretisation of the ECM Eq (2.1g), and each of its components H3i,j(⋅,⋅), ∀i,j=1…N, are given by
H3i,j(˜cm1(t),˜cm2(t),˜vm(t)):=−ρ(˜cm1,i,j(t)+˜cm2,i,j(t))˜vmi,j(t)+μv(Kc−˜cm1,i,j(t)−˜cm2,i,j(t)−˜vmi,j(t))+. | (3.9) |
Consider now a time discretisation {tn}n=0…L with time step Δt:=T/L. For each n∈{0,…,L}, a simple Euler time-marching scheme can be written for model (3.6) via the associated operator
Km:RN×N×RN×N×RN×N→RN×N×RN×N×RN×NgivenbyKm([˜cm,n1˜cm,n2˜vm,n]):=[˜cm,n1˜cm,n2˜vm,n]+Δt[H1(˜cm,n1,˜cm,n2,˜vm,n,m)H2(˜cm,n1,˜cm,n2,˜vm,n,m)H3(˜cm,n1,˜cm,n2,˜vm,n)], | (3.10) |
where ˜cm,n1:=˜cm1(tn), ˜cm,n2:=˜cm2(tn) and ˜vm,n:=˜vm(tn). The right-hand-side operators are given by
H1(˜cm,n1,˜cm,n2,˜vm,n,m):=H1(˜cm1(tn),˜cm2(tn),˜vm(tn),m),H2(˜cm,n1,˜cm,n2,˜vm,n,m):=H2(˜cm1(tn),˜cm2(tn),˜vm(tn),m),H3(˜cm,n1,˜cm,n2,˜vm,n):=H3(˜cm1(tn),˜cm2(tn),˜vm(tn)). |
This allows us to formulate the "forward operator" K between the family of functions M1 (where we search for the appropriate cancer cells mutation function mc∗1,c∗2,v∗) and the space where the discretised measurements model (3.1) are recorded. Hence, the "forward operator" K is defined as
K:S→RN×N×RN×N×RN×NgivenbyK(m):=Km∘Km∘⋯⋯∘Km⏟Ltimes([˜c1,0˜c2,0˜v0]) | (3.11) |
where ˜c1,0:={c1,0(xi,yj)}i,j=1,…,N, ˜c2,0:={c2,0(xi,yj)}i,j=1,…,N and ˜v0:={v0(xi,yj)}i,j=1,…,N are the discretised initial conditions (2.5) for the governing tumour forward model (2.1). Hence, for each m∈M1, the forward operator K gives the spatio-temporal progression of the initial condition [˜c1,0,˜c2,0,˜v0]T under the invasion model (2.1), which is obtained when the cell mutation law at each instance of time t>0 involves the corresponding trial mutation operator M(⋅,m) instead of the unknown mutation terms ˜Qr(⋅).
From models (3.10) and (3.11) we have that our forward operator K is given as a finite composition of affine functions of the form
M1∋m⟼Km∈ℓ2(ℓ2(E×E×E);ℓ2(E×E×E)). | (3.12) |
Here, ℓ2(ℓ2(E×E×E);ℓ2(E×E×E)) is the usual finite-dimensional Bochner space of square integrable vector-value functions [28] with respect to the counting measure (see [29]) that are defined on ℓ2(E×E×E) and take values in ℓ2(E×E×E), and E:={Ei,j}i,j=1…N represents the standard basis of elementary matrices associated with the grid GΩ. As a direct consequence, we immediately obtain that this operator is both continuous and compact, from where we obtain that K is also closed and sequentially bounded [28]. Therefore, K satisfies the hypotheses assumed in [30] that ensure convergence for the nonlinear Tikhonov regularisation strategy given by the functionals {Jα}α>0,
Jα:M1→R,∀α>0,definedbyJα(m):=‖K(m)−[˜c∗1˜c∗2˜v∗]‖22+α‖m‖22,∀m∈M1. | (3.13) |
The minimisation of these functionals enable us to identify mc∗1,c∗2,v∗ as the limit α→0 of the points of minimum mα of Jα (these points correspond to the smallest discrepancy between the data measurements and the solution of our system that uses mα as a mutation law). The two norms involved in model (3.13) represent the usual Euclidean norms on the corresponding finite dimensional spaces. Indeed, while the first is the standard Euclidean norm on RN×N×RN×N×RN×N, the second is also the Euclidean norm induced on the M−dimensional space of functions M1 via the standard isomorphism that we establish between M1 and RM by which each m∈M1 is uniquely represented through its nodal values {m(ηl)}l=1…M with respect to the linear basis functions {ˉϕl}l=1…M associated to GM [31]:
sincem=∑l=1…Mm(ηl)ˉϕl,wethereforemaketheidentification:m≡{m(ηl)}l=1…M. | (3.14) |
Finally, in model (3.13), ˜c∗1, ˜c∗2 and ˜v∗ represent the discretised measurements of the densities of cancer cells and ECM given in Eqs (3.1a)–(3.1c), i.e., ˜c∗1:={c∗1(xi,yj)}i,j=1,…,N, ˜c∗2:={c∗2(xi,yj)}i,j=1,…,N and ˜v∗:={v∗(xi,yj)}i,j=1,…,N. We assume that these data measurements are either exact or are corrupted by a certain noise level δ≥0. Thus, maintaining for simplicity the measurements notation unchanged, these measurements are given by
˜c∗1(x)=˜c∗1exact(x)+δγc1(x), | (3.15a) |
˜c∗2(x)=˜c∗2exact(x)+δγc2(x), | (3.15b) |
˜v∗(x)=˜v∗exact(x)+δγv(x), | (3.15c) |
where, ∀x∈Ω, we have that ˜c∗1exact(x), ˜c∗2exact(x) and ˜v∗exact(x) describe the exact data, and γc1(x), γc2(x) and γv(x) are signal-independent noise generated from a Gaussian normal distribution with mean zero and standard deviations σc1, σc2 and σv, respectively, given by
{σc1:=1λ(Ω)∫Ω˜c∗1exact(x)dx,σc2:=1λ(Ω)∫Ω˜c∗2exact(x)dx,σv:=1λ(Ω)∫Ω˜v∗exact(x)dx, | (3.16) |
with λ(⋅) being the usual Lebesgue measure. In the numerical results below, we generate the random variables γc1(x), γc2(x) and γv(x) via MATLAB function normrnd by taking {γc1(xi,yj)}i,j=1…N:=normrnd(0,σc1,N×N), {γc2(xi,yj)}i,j=1…N:=normrnd(0,σc2,N×N) and {γv(xi,yj)}i,j=1…N:=normrnd(0,σv,N×N).
We explore now the inversion approach that we formulated so far in the context of the forward model (2.1) by proceeding with the reconstruction of the unknown terms involved in the mutation laws in Cases (i) and (ii), namely ˜Q1(c1) and ˜Q2(v), respectively.
Initial Conditions The initial conditions (2.5) that we consider in the computations for the forward model (2.1) are as follows:
c1,0(x):=0.5(exp(−‖x−(2,2)‖220.03)−exp(−9.407)), | (3.17a) |
c2,0(x):=0, | (3.17b) |
v0(x):=0.5+0.3⋅sin(4π⋅‖x‖2),∀x∈Ω. | (3.17c) |
Here, we assume that c2,0(x)=0 because this second cell population will arise after a period of time following mutations of the first cell population c1.
To identify the cancer cells mutation law, we consider both exact and noisy measurement data model (3.15) as additional information for the forward model (2.1) in the presence of initial conditions (3.17) and boundary conditions (2.6). Specifically, in each of the two cases, we assume that the exact data (namely ˜c∗1,exact(x), ˜c∗2,exact(x) and ˜v∗exact(x)) that appear in model (3.15) are given by the solutions densities for primary tumour cells population, ˉc1(x,t), mutated cells population, ˉc2(x,t), and ECM, and ˉv(x,t), evaluated at the final time tf>0, i.e.,
˜c∗1,exact(x):=ˉc1(x,tf),˜c∗2,exact(x):=ˉc2(x,tf)and˜v∗exact(x):=ˉv(x,tf),∀x∈Ω, | (3.18) |
which are obtained from the forward model (2.1) as follows:
● for Case (i): we assume that the mutation law is of the form given in model (2.2) but when parameter δ0>0 is considered known and has the value given in parameter Table 1;
Parameter | Value | Description | Reference |
D1 | 0.00675 | diffusion of primary tumour | [11] |
D2 | 0.00675 | diffusion of secondary tumour | [11] |
η1 | 2.85×10−2 | haptotaxis to ECM from c1 | [37] |
η2 | 2.85×10−2 | haptotaxis to ECM from c2 | [37] |
μc | 0.25 | proliferation of tumour cells c | [38] |
Kc | 1 | tissue carrying capacity | [16] |
ρ | 2 | ECM degradation coefficient | [39] |
μv | 0.40 | ECM remodelling coefficient | [16] |
Scc | (0.5000.3) | cell-cell adhesion function | [11] |
Scv | (0.1000.5) | cell matrix adhesion function | [11] |
t1,2 | 10 | time initiation for mutations | [38,40] |
ts | 3 | time-steepness coefficient | [38] |
δ0 | 0.3 | mutation from primary tumour | [38,40] |
Δx | 0.03125 | discretisation step size for GΩ | [16] |
Δt | 10−3 | time step size | [16] |
Δη | 0.0625 | mesh size used for G1M and G2M | Estimated |
● for Case (ii): we assume that the mutation law is of the form given by model (2.3) with the known term ˜Q2(v) specified in model (2.4).
For each regularisation parameter α>0 considered here, the minimisation process for Jα is initiated with m0=I×10−3, (where I represents the M vector of ones), and for the actual minimisation we employed here the nonlinear minimisation MATLAB function lsqnonlin. Finally, since there are no data to test the trial mutation operators beyond the maximal accessible region Ac and Av defined by the minimum and maximum values of the solution, i.e.,
Ac:=[ˉcmin1,ˉcmax1],with:ˉcmin1:=min(x,t)∈Ω×[0,T]c1(x,t),ˉcmax1:=max(x,t)∈Ω×[0,T]c1(x,t), | (3.19) |
Av:=[ˉvmin,ˉvmax],with:ˉvmin:=min(x,t)∈Ω×[0,T]v(x,t),ˉvmax:=max(x,t)∈Ω×[0,T]v(x,t), | (3.20) |
the reconstructions in this section will be attempted only for the restriction of the sought mutation laws to Ac and Av. We need to emphasise that Ac is used to reconstruct ˜Q1(c1), while Av is used to reconstruct ˜Q2(v).
An acceptable numerical reconstruction of the mutation law mc∗1,c∗2,v∗, i.e.,
mc∗1,c∗2,v∗:=mα∗, | (3.21) |
is obtained for the choice of the regularisation parameter α∗, which throughout this work is selected based on a standard discrepancy principle argument [32].
Figure 2 shows the reconstruction of the cancer cell mutation law ˜Q1(c1) for model (2.1) in the presence of the measurements given by models (3.15) and (3.18) that are considered here both exact and affected by a level of noise δ∈{1%,3%}. The first row of this figure shows from left to right the true mutation law restricted at the maximal accessible region Ac (where the reconstruction is being attempted) and the reconstruction of the mutation law equation ˜Q1(c1) given by model (2.2) on Ac with no noise, respectively. The second row of the figure show from left to right the reconstruction of the mutation law on Ac with 1%, and 3% of noise in the measured data, respectively.
Figure 3 shows the reconstruction of unknown nonlinear cancer cell mutation term ˜Q2(v) (that is part of cell mutation law introduced in model (2.3)) for model (2.1) in the presence of the measurements given by models (3.15) and (3.18) that are considered here both exact and affected by a level of noise δ∈{1%,3%}. Again, the first row shows from left to right the true mutation law restricted to Av where the reconstruction is attempted, and the reconstruction of the mutation law model (2.3) on Av with no noise, respectively. The second row of the figure show from left to right the mutation reconstruction on Av with 1%, and 3% of noise in the measured data, respectively.
From Figures 2 and 3 we observe that when the measurement data is not affected by noise, we obtain good mutation laws reconstructions in both cases, i.e., (i) mutation depends on primary tumour only, and (ii) mutation depend on the ECM only. However, as expected, as soon as the level of noise in the measurements increases, the reconstruction gradually looses accuracy in both cases (i.e., Case (i) shown in Figure 2 lower panels, and Case (ii) shown in Figure 3 lower panels). This loss in accuracy increases with the increase in the c1 density (horizontal axis in Figure 2) and in the v density (horizontal axis in Figure 3).
In this section, we study the inverse problem of identifying the general mutation law in case (iii), namely Q(c1,v), which is unknown and depends nonlinearly in an autonomous manner on both primary tumour cell density c1 and on ECM density v.
Maintaining here the same spatio-temporal discretisation notations introduced in Section 3.1 for the cancer cell populations and ECM, we seek to identify the unknown nonlinear mutation law Q(c1,v) such that we match the measurements model (3.1). To achieve this, similar to the previous two cases addressed in Sections 3.1-3.3, within an suitable space of functions M2 that will be detailed below, we proceed to identify an appropriate a mutation approximating function defined on the 2-dimensional region where the pair (c1,v) ranges during its spatio-temporal evolution, namely ¯mc∗1,c∗2,v∗:[0,Kc]×[0,Kc]→[0,∞). The function ¯mc∗1,c∗2,v∗ will be identified such that this will directly determine an acceptable "mutation law candidate", denoted here by Qc∗1,c∗2,v∗(⋅,⋅), which will enable a solution in forward model (2.1) that matches the measurements model (3.1). To select this mutation law candidate Qc∗1,c∗2,v∗(⋅,⋅), adopting a similar approach to the one in Cases (i) and (ii), we involve again an appropriately constructed mutation trial operator M2(⋅,⋅,⋅):={M2i,j(⋅,⋅,⋅)}i,j=1…N, M2(⋅,⋅,⋅):RN×N×RN×N×M2→RN×N that will be specified in a moment and will enable us the express the acceptable mutation law candidate as
Qc∗1,c∗2,v∗(˜c¯mc∗1,c∗2,v∗1,i,j(t),˜v¯mc∗1,c∗2,v∗i,j(t)):=M2i,j(˜c¯mc∗1,c∗2,v∗1(t),˜v¯mc∗1,c∗2,v∗(t),¯mc∗1,c∗2,v∗) | (3.22) |
where ˜c¯mc∗1,c∗2,v∗1(t):={˜c¯mc∗1,c∗2,v∗1(t)}i,j=1…N and ˜v¯mc∗1,c∗2,v∗(t):={˜v¯mc∗1,c∗2,v∗(t)}i,j=1…N represents the solution for the density of the primary cell population and of the ECM, respectively, which are obtained for model (2.1) when, instead of the unknown term Q(⋅,⋅), in the mutation law we use the trial mutation term M2(⋅,⋅,mc∗1,c∗2,v∗)
Denoting by G2M:={(ηl,ζk)}l,k=1…M the equally spaced grid given by the uniform discretisation of the [0,Kc]×[0,Kc] with step size Δη=Δζ>0, the space of functions where we seek to identify the function ¯mc∗1,c∗2,v∗ is an M×M- dimensional space of potential mutation law shape candidates, which is given by
M2:={¯m:[0,Kc]×[0,Kc]→R|¯m|El,k=∑p,q=0,1¯m(ηl+p,ζk+q)ϕl+p,k+q,∀El,k∈G2,tilesM}withG2,tilesM:={El,k:=[ηl,ηl+1]×[ζk,ζk+1]|l,k=1…M−1},and∀El,k∈G2,tilesM,{ϕl+p,k+q}p,q=0,1aretheusualbilinearshapefunctionsonEl,k. | (3.23) |
Therefore, for any ¯m∈M2, the trial proliferation operator M2 is given by
M2i,j(˜c¯m1(t),˜v¯m(t),¯m):=¯m|El,k(˜c¯m1(t),˜v¯m(t))with(l,k)∈Λi,j:={l,k∈{1,…,M−1}|∃El′,k′∈GtilesMsuchthat(˜c¯m1,i,j(t),˜v¯mi,j(t))∈El′,k′},andnotingalsoherethat(l,k)isindependentofitschoicewithinΛi,j. | (3.24) |
Here, ˜c¯m1(t):={˜c¯m1,i,j(t)}i,j=1…N, ˜c¯m2(t):={˜c¯m2,i,j(t)}i,j=1…N and ˜v¯m(t):={˜v¯mi,j(t)}i,j=1…N represent the solutions at the grid points and time t>0 for the cancer cells and ECM densities obtained with model (2.1) when this uses M2i,j(⋅,⋅,¯m) as mutation law given in model (3.24). Therefore, in space-discretised form, model (2.1) can be recasted also in this case as
∂∂t[˜c¯m1˜c¯m2˜v¯m]=[H1(˜c¯m1,˜c¯m2,˜v¯m,¯m)H2(˜c¯m1,˜c¯m2,˜v¯m,¯m)H3(˜c¯m1,˜c¯m2,˜v¯m)], | (3.25) |
Also in in this case (i.e., case (iii)) we have that H1(⋅,⋅,⋅,⋅) and H2(⋅,⋅,⋅,⋅) are correspondingly defined through Eqs (3.7) and (3.8) when the trial mutation form for the full mutation law (which in Cases (i) and (ii) was given by operator ¯Mr) is given here by M2i,j(˜c¯m1(t),˜v¯m(t),¯m).
Finally, adopting again the same time discretisation as in Section 3.1 and using the Euler method, a time marching step for model (3.25) is given by following operator
ˉKˉm:RN×N×RN×N×RN×N→RN×N×RN×N×RN×NgivenbyˉKˉm([˜cm1,n1˜cm1,n2˜vm1,n]):=[˜cm1,n1˜cm1,n2˜vm1,n]+Δt[H1(˜cˉm,n1,˜cˉm,n2,˜vˉm,n,ˉm)H2(˜cˉm,n1,˜cˉm,n2,˜vˉm,n,ˉm)H3(˜cˉm,n1,˜cˉm,n2,˜vˉm,n)], | (3.26) |
where, for any n∈{0,…,L}, we have ˜c¯m,n1:=˜c¯m1(tn), ˜c¯m,n2:=˜c¯m2(tn), and ˜v¯m,n:=˜v¯m(tn), while
H1(˜c¯m,n1,˜c¯m,n2,˜v¯m,n,¯m):=H1(˜c¯m1(tn),c¯m2(tn),˜v¯m(tn),¯m),H2(˜c¯m,n1,˜c¯m,n2,˜v¯m,n,¯m):=H2(˜c¯m1(tn),c¯m2(tn),˜v¯m(tn),¯m),H3(˜c¯m,n1,˜c¯m,n2,˜v¯m,n):=H3(˜c¯m1(tn),˜c¯m2(tn),˜v¯m(tn)). |
Therefore, this allows us formulate "forward operator" denoted here by ¯K defined by
¯K:M2→RN×N×RN×N×RN×Ngivenby¯K(¯m):=¯K¯m∘¯K¯m∘⋯⋯∘¯K¯m⏟Ltimes([˜c1,0˜c2,0˜v0]) | (3.27) |
where ˜c1,0, ˜c2,0 and ˜v0 are the discretised initial conditions for model (2.1) (as introduced in model (3.11)).
Similar to the previous two inverse problem cases, the forward operator ¯K will enable us to identify the ¯mc∗1,c∗2,v∗∈M2 such that the solution of the resulting model matches measurements model (3.1). Thus, similar to Section 3.2, we observe again that ¯K can be written down as
M2∋¯m⟼¯K¯m∈ℓ2(ℓ2(E×E×E);ℓ2(E×E×E)) | (3.28) |
Here, ℓ2(ℓ2(E×E×E);ℓ2(E×E×E)) is again the finite-dimensional Bochner space of square integrable vector-value functions defined on ℓ2(E×E×E) with values in the same space. As in Section 3.2, also here we have that the mappings defined in model (3.28) are continuous and compact. Therefore, we obtain that ¯K is also closed sequentially bounded, and as a consequence we have that also in this case we satisfy the inverse problems hypotheses adopted [30], and thus, the convergence of the subsequent nonlinear Tikhonov regularisation is ensured. The Tikhonov functionals {¯Jα}α>0 will have essentially the same form as in the first two cases, except that the space where these are defined is different (i.e., in this case we have M2 of dimension M×M rather that M1 of dimension M that we had in Cases (i) and (ii)), namely:
¯Jα:M2→R,∀α>0,definedby¯Jα(¯m):=‖¯K(¯m)−[˜c∗1˜c∗2˜v∗]‖22+α‖¯m‖22,∀¯m∈M2. | (3.29) |
Computationally, we address here the reconstruction of the general mutation law Q(c1,v) that appears in Case (iii) for model (2.1) in the presence of zero flux boundary conditions and the initial conditions prescribed in Eqs (3.17a)–(3.17c). Furthermore, the inverse problem is addressed in the presence of both exact and noisy measurements of the form detailed in Eqs (3.15) and (3.16), with the exact measurements given by
˜c∗1,exact(x):=ˉc1(x,tf),˜c∗2,exact(x):=ˉc2(x,tf),and˜v∗exact(x):=ˉv(x,tf),∀x∈Ω, | (3.30) |
where ˉc1(x,t),ˉc2(x,t) and ˉv(x,t) represent the solution densities for primary cell population, mutated cell population and ECM, respectively, which are obtained when model (2.1) uses as mutation law the expression given by model (2.3) with the known term ˜Q2(v) specified in model (2.4).
Although the dimensionality is different with respect to Cases (i) and (ii), the inversion method in this case follows the same steps and numerical minimisation procedure and steps that were described in Section 3.3 (reason for which we skip here details already outlined there). Indeed, for any α∈{10−i|i=1,…12}, the minimisation of the Tikhonov functional Jα starts with an initial guess ¯m0=I×10−3, (where I represents the M×M matrix of ones), and leads to the numerical identification of the associated point of minimum mα that correspond to the smalles mismatch between the associated solution (that is obtained when model (2.1) uses mα as mutation law) and the measurements. Furthermore, since no data is available beyond the maximal region maximal accessible regions Acv given by
Acv:=Ac×Av=[ˉcmin1,ˉcmax1]×[ˉvmin,ˉvmax],withAcandAvdescribedinEqs(3.19)and(3.20). |
the reconstruction of the mutation law is explored only on Acv. Finally, using again an discrepancy principle-based argument to choose the regularisation parameter α∗∈{10−i|i=1,…12}, the reconstructed mutation law will be given by the corresponding M×M matrix ¯mc∗1,c∗2,v∗:=¯mα∗, which in turn will determine the precise shape of the reconstructed mutation law Qc∗1,c∗2,v∗(˜c¯mc∗1,c∗2,v∗1,i,j(t),˜v¯mc∗1,c∗2,v∗i,j(t)) that is defined as per model (3.22). Figure 4 shows the reconstruction of the most general cancer cell mutation law Q(c1,v), starting from the measurements given by models (3.15) and (3.18) for exact data as well as data impacted by noise δ∈{1%,3%}. For comparison, the first row of this figure shows the true mutation law restricted at the maximal accessible region Acv where the reconstruction is being attempted. The second row of the figure show from left to right the reconstruction of the mutation law on Acv with no noise (left), with 1% noise (center), and 3% noise (right) in the measured data. From this figure, we observe that we obtain a good reconstruction of the mutation law when the measurement data are not affected by noise. However, as expected, as soon as the level of noise increases in the measurements, the reconstruction gradually looses accuracy.
Throughout the previous sections we assumed that both cancer cell populations proliferate logistically: μcc1,2(1c1+c2+vKc). However, the sigmoid shape of tumour growth that is given by the logistic term can be obtained also with other proliferation rules, such as the Gompertz rule [33,34,35]: μcc1,2log(Kcc1+c2+v). This raises the question as to what happens when different cancer cell populations use different proliferation laws.
In this subsection we reconstruct the general mutation law Q(c1,v,t) in case (iii) when we assume that the primary c1 population proliferation proliferates logistically, while the second (mutated) cell population c2 proliferates according to Gompertz law. In Figure 5 we present the numerical reconstruction results. We observe that the results are similar to those in Figure 4; this could be explained by the fact that the mutation starts in the primary tumour which proliferate logistically.
Since various mathematical studies have assumed nonlocal cell-cell and cell-ECM interactions to explain the invasion of cancer cells [11,16,22,23], in the following we generalise the model (2.1) by replacing the local haptotaxis towards local ECM gradients with a nonlocal haptotaxis flux generated by these cell-cell and cell-ECM interactions.
As before, we consider a primary cancer cell subpopulation c1(x,t) and a secondary mutated cancer cell subpopulation c2(x,t). These cancer cell populations interact with each other as well as with the extracellular matrix (ECM), v(x,t), which they degrade and remodel. For compact notation, we consider the combined vector of primary cancer cells c1, mutated cancer cells c2 and extracellular matrix v defined as
u(x,t):=[c1(x,t),c2(x,t),v(x,t)]T. |
We use this vector to describe in a compact manner the flux term generated by the nonlocal cell-cell and cell-ECM interactions (A1,2(x,t,u⋅,t)). Therefore the coupled tumour dynamics in this case is given by:
∂c1∂t=D1Δc1⏟randommotility−∇⋅[c1A1(x,t,u(⋅,t))]⏟adhesion+μcc1(1−c1+c2+vKc)⏟logisticproliferation−ω(t)⏟mutationswitchremodQ(⋅,⋅)⏟unknownmutation, | (4.1a) |
∂c2∂t=D2Δc2⏟randommotility−∇⋅[c2A2(x,t,u(⋅,t))]⏟adhesion+μcc2(1−c1+c2+vKc)⏟logisticproliferation+ω(t)⏟mutationswitchremodQ(⋅,⋅)⏟unknownmutation, | (4.1b) |
∂v∂t=−ρ(c1+c2)v⏟degradation+μv(Kc−c1−c2−v)+⏟ECMremodelling. | (4.1c) |
As before, D1,2 are the diffusion coefficients, μc is the cancer cell proliferation rate, ρ is the ECM degradation rate, and μv is the ECM remodelling rate. The flux term Ap(x,t,u(⋅,t)), p=1,2, has been proposed in [11,23] to describe the directed movement of cells due to cell-cell and cell-matrix adhesion:
Ap(x,t,u(⋅,t)):=1R∫B((0,0),R)n(y)⋅K(‖y‖2)⋅gp(u(x+y,t),t)χΩ(x+y) dy. | (4.2) |
By using these notations the approximation for DpΔcp−∇⋅[cpAp(x,t,u(⋅,t))] in (4.1) is as follows:
It is assumed that the interactions between a cell and its neighbouring cells as well as the components of the ECM occur inside a sensing region B((0,0),R)⊂R, where R>0 is the sensing radius, as illustrated in Figure 6. In the above equation, χΩ(⋅) represents the characteristic function of Ω. Further, n(y) is the unit radial vector given by
n(y):={y/‖y‖2if y∈B(0,R)∖{(0,0)},(0,0)otherwise. | (4.3) |
Furthermore, gp(u(x+y,t),t) represents the p−th component of the adhesion function g(u(x+y,t),t) that is given by
g(u,t)=[Sccc+Scvv]⋅(Kc−c1−c2−v)+,withScc=(Sc1c1Sc1c2Sc2c2Sc2c2)andScv=(Sc1v00Sc2v). | (4.4) |
Here {Scicj}i=1,2 are the non-negative cell-cell adhesion strengths of the adhesion bonds established between the primary and mutated cancer cell populations, while {Sciv}i=1,2 are the non-negative cell-matrix adhesion stands for the strengths for the adhesive interactions between each of the two cancer subpopulation and the ECM. Finally, no overcrowding of the cell population and ECM over the tumour domain is ensured here through the term (Kc−c1−c2−v)+:=max(Kc−c1−c2−v,0).
The nonlocal dynamics model (4.1) is assumed to take place in the presence of the same initial and boundary condition as those assumed in for the local model and given by models (2.5) and (2.6) in Section 2. The numerical approximation of the forward cancer invasion model (4.1) is presented in Appendix Section A. This includes the discretisation of the main spatial operators as well as the off grid numerical approach for the nonlocal adhesion flux terms Ap(x,t,u(⋅,t)) whose schematic is given by the right side of Figure 6. Finally, the general mutation law Q(c1,v) is assumed to be unknown and its identification will be our main focus in this section.
Building on the inverse problems approach developed for the local cancer invasion model in Section 3.4 and applying this to the nonlocal model (4.1), we proceed now to address the reconstruction of the unknown mutation law Q(c1,v) within the nonlocal cancer invasion model (4.1) from exact and noisy measured data that are considered here again to be of the form prescribed in models (3.15) and (3.16). Here, the exact measurement data that are given by
˜c∗1,exact(x):=ˉc1(x,tf),˜c∗2,exact(x):=ˉc2(x,tf),and˜v∗exact(x):=ˉv(x,tf),∀x∈Ω, | (4.5) |
where ˉcp(x,t) p=1,2 and ˉv(x,t) are the corresponding solution densities of the primary and mutated cancer cell populations as well as that of ECM that are obtained for model (4.1) obtained when this uses as mutation law given by model (2.3) with the known term ˜Q2(v) detailed in model (2.4). For the numerical reconstruction we consider here the same initial and boundary conditions as in Section 3.3. Figure 7 shows the reconstruction of the cancer cell mutation law Q(c1,v) for model (4.1) in the presence of the measurements given by model (4.5) that are considered here both exact and affected by a level of noise δ∈{1%,3%}. As for the previous figures, the first row shows the true mutation law restricted at the maximal accessible region Acv where the reconstruction is being attempted. The second row of the figure shows, from left to right, the reconstruction of the mutation law on Acv with no noise, with 1%, and with 3% of noise in the measured data, respectively. From this figure, we observe that we obtain a good mutation law reconstruction when the measurement data are not affected by noise. However, as expected, as soon as the level of noise increases in the measurements, the reconstruction gradually looses accuracy.
As in Section 3.4.3, here we investigate the reconstruction of the general mutation law Q(c1,v) when cancer cell growth is described by different rules: logistic proliferation for the c1 cells, and Gompertz proliferation for the second (mutated) cancer cell population c2. Again, we do not see any significant differences between the reconstruction of mutation law in this case (see Figure 8) and in the previous case where both cell populations grow logistically (see Figure 7). We believe that is because in both cases, the population that starts mutating (i.e., the c1 population) grows logistically.
In this work we explored a new inverse problem that addresses the reconstruction of the cancer cells mutation law in two heterotypic cancer invasion models: a model with a local cancer cell haptotaxis flux towards ECM (see model (2.1)), and a model with a nonlocal haptotaxis flux generated by cell-cell and cell-ECM adhesion forces (see model (4.1)). For the reconstruction of this mutation law through an inverse problem Tikhonov regularisation-based approach, we used a numerically-generated spatial tumour snapshot data assumed to be acquired at a later stage in the tumour evolution (in practice, the data can be provided through a medical imaging scan).
This inverse problem approach was implemented computationally via a mixed finite differences - finite element numerical scheme. Specifically, on one hand, we used a Crank-Nicholson-type finite difference scheme for the discretisation of the forward models that arises in each of the two tumour invasion dynamics considered here (i.e., local cancer cell invasion, and nonlocal cancer cell invasion). On the other hand, we developed a finite element approach involving a bilinear shape functions on a square mesh for the discretisation of mutation law candidates recruited from a proposed space of functions S, as well as for their evaluation on a maximal accessible regions where the mutation law reconstruction was performed. Finally, these two parts were appropriately assembled in an optimisation solver that sought to reconstruct the cancer cell mutation law by minimising over the space S the emerging Tikhonov functionals that were formulated in each of the two cases considered.
Finally, this inversion approach was explored and tested on the reconstruction of different cancer cell mutation laws used in local and nonlocal cancer modelling: (i) reconstructing the cancer cell mutation assuming a linear dependence only on the c1 cell population; (ii) reconstructing the cancer cell mutation assuming a linear dependence on the c1 population and a nonlinear dependence on the ECM density, and (iii) reconstructing a very general mutation law, while assuming no prior knowledge about the mutation. To explore the robustness of the proposed inverse problem approach we used "in silico" constructed data that mathematically could be expressed in the general form given by models (3.15) and (3.16), but that were specified individually for each of the cases considered for the local and non-local models (2.1) and (4.1), respectively. Indeed, while the noise was always added the same way (as per model (3.16)), the exact measurements considered here were maps given by solutions at the final time for the densities of the cancer cell populations (i.e, both primary and mutated populations) and the ECM, namely (ˉc1(⋅,tf),ˉc2(⋅,tf),ˉv(⋅,tf)), which were obtained by equipping the local and nonlocal forward models with a selection of known mutation laws (given either by model (2.2) or by models (2.3) and (2.4)) considered as appropriate for each of the cases explored in the context of both local and nonlocal models for tumour dynamics. We note that the results in Figures 4, 5, 7, and 8 look correspondingly similar because we reconstruct the same mutation law Q(c, v), while for the second (mutated) cancer cell population, c2, we use logistic proliferation for Figures 4 and 7 and Gompertz proliferation law for Figures 5 and 8. While there are small numerical differences, these are hard to be distinguished among the figures without zooming in. Furthermore, this similarity of the results is not significantly affected when switching between the Logistic and Gompertz proliferation laws for the second cancer cells population. An explanation for this similarity is that in addition to that fact that the mutation starts in the primary tumour which proliferate logistically, the logistic and Gompertz laws considered for c2 are part of the same class of size growth laws [33,34,35], with similar profiles on the parts of their domain where these are evaluated by the cancer models (2.1) and (4.1). All numerical reconstruction results for the local and nonlocal models showed that for exact measurements we obtained a good reconstruction of the mutation law, while for increasingly noisy measurements the reconstruction gradually deteriorates. Finally, these in silico reconstruction tests recommend the inversion framework developed in this work also for the context of real measurements data that come in the form of imaging data (such as MRI) with levels of noise that may not be a priori known. Nevertheless, the accuracy of the reconstruction in the real data case is expected to deteriorate for significant levels of noise.
The first author would like to thank Saudi Arabian Cultural Bureau in the United Kingdom (UKSACB) on behalf of Taif University and the University College of Al-Khurmah in Saudi Arabia for supporting and sponsoring his PhD studies at the University of Dundee.
The authors declare no conflict of interest.
In this section, we briefly discuss the numerical methods used to solve the forward models (4.1a)–(4.1c). To discretise the system in space we use the method of lines approach. For the time-evolution of the system, we use a predictor-corrector scheme introduced in [16], using the Euler method as the predictor and a trapezoidal-type rule as the corrector. In the reaction-diffusion models (4.1a) and (4.1b), the terms DpΔcp−∇⋅[cpAp(x,t,u(⋅,t))], p=1,2, will be approximated through a second-order mid-point rule [36] as detailed below. In brief, ∀p∈{1,2},n=¯0,L,i,j=¯1,N, these midpoint approximations are given by:
{cnp,i,j+12:=cnp,i,j+cnp,i,j+12,cnp,i,j−12:=cnp,i,j+cnp,i,j−12,cnp,i+12,j:=cnp,i,j+cnp,i+1,j2,cnp,i−12,j:=cnp,i,j+cnp,i−1,j2.and {Anp,i,j+12:=Anp,i,j+Anp,i,j+12,Anp,i,j−12:=Anp,i,j+Anp,i,j−12,Anp,i+12,j:=Anp,i,j+Anp,i+1,j2,Anp,i−12,j:=Anp,i,j+Anp,i−1,j2. | (A.1) |
Further, the central differences are given by
{[cp,y]ni,j+12:=cnp,i,j+1−cnp,i,jΔy,[cp,y]ni,j−12:=cnp,i,j−cnp,i,j−1Δy,[cp,x]ni+12,j:=cnp,i+1,j−cnp,i,jΔx,[cp,x]ni−12,j:=cnp,i,j−cnp,i−1,jΔx. | (A.2) |
DpΔcp−∇⋅[cpAp(x,t,u(⋅,t))]=≃Dp[cp,x]ni+12,j−Dp[cp,x]ni−12,j−cnp,i+12,j⋅Anp,i+12,j+cnp,i−12,j⋅Anp,i−12,jΔx+Dp[cp,y]ni,j+12−Dp[cp,y]ni,j−12−cnp,i,j+12⋅Anp,i,j+12+cnp,i,j−12⋅Anp,i,j−12Δy. | (A.3) |
Next, we shift our attention to the numerical approach of the adhesive flux Ap (that explores the effects of cell-cell and cell-matrix adhesion of cancer cells subpopulations). We carry out these computations off-grid by decomposing the sensing region B(x,R) in s∈N∗ annuli centred at x (with the inner most central circle being of numerically negligible size), and for each annuli k∈{1,…,s} (counted from the inner most to the outer most), we consider a radial decomposition of this in 2h+(k−1) uniformly distributed radial sectors (with h∈N fixed), as illustrated in the right side of Figure 6 (for h=2). This leads to a collection of sectors {Sν}ν=1,Ns, where
Ns:=s∑k=12h+(k−1). |
Then for each annulus sector Sν, we evaluate the total primary cancer cell population c1, the total mutated cancer cell population c2 and the total ECM mass distributed on Sν that are given by
ωSν,c1:=1λ(Sν)∫Sνc1(x,t)dx,ωSν,c2:=1λ(Sν)∫Sνc2(x,t)dx,andωSν,v:=1λ(Sν)∫Sνv(x,t)dx. | (A.4) |
Further, on each sector Sν, we consider the off-grid barycentres of each annulus sector by bSν and the values of each of the three densities at these barycentres, namely c1(bSν,t) c2(bSν,t) and v(bSν,t), are obtained via interpolation with bilinear shape functions on the grid rectangles {yibSν}i=1,4 that contain bSν. Therefore, the approximation of the adhesion integral at each instance of time t>0 is given by
Ap(x,t,u(⋅,t))=s2m∑ν=1λ(Sν)Rn(bSν)⋅K(bSν)gp(˜u(bSν,t)), | (A.5) |
where
˜u(bSν,t):=[ωSν,c1(t),ωSν,c2(t),ωSν,v(t)]T, |
and gp(˜u(bSν,t)) is the p−th component of
g(˜u(bSν,t))=[Scc[ωSν,c1(t),ωSν,c2(t)]T+ωSν,v(t)Scv]⋅(Kc−c(bSν,t)−v(bSν,t))+. |
For all the cancer cells mutation laws reconstructions considered in this work, we use the non-dimensional parameter set specified in Table 1.
[1] |
O'Reagan B, Grätzel M (1991) A low cost, high efficiency solar cell based on dye sensitized colloidal TiO2 films. Nature 353: 737-740. https://doi.org/10.1038/353737a0 doi: 10.1038/353737a0
![]() |
[2] |
Gong J, Sumathy K, Qiao Q, et al. (2017) Review on dye-sensitized solar cells (DSSCs): Advanced techniques and research trends. Renewable Sustainable Energy Rev 68: 234-246. https://doi.org/10.1016/j.rser.2016.09.097 doi: 10.1016/j.rser.2016.09.097
![]() |
[3] |
Ehrmann A, Blachowicz T (2019) Recent coating materials for textile-based solar cells. AIMS Mater Sci 6: 234-251. https://doi.org/10.3934/matersci.2019.2.234 doi: 10.3934/matersci.2019.2.234
![]() |
[4] | Ehrmann A, Blachowicz T (2019) Comment on 'Dye-sensitized solar cells using Aloe vera and cladode of cactus extracts as natural sensitizers' [Chem. Phys. Lett. 679 (2017) 97-101]. Chem Phys Lett 714: 227-229. https://doi.org/10.1016/j.cplett.2018.10.009 |
[5] |
Burgos MJC, Roa S, Cerda B, et al. (2021) Effects of PbS-NPs doping on the photovoltaic performance of natural dye-sensitized TiO2 photoelectrodes. Solid State Commun 340: 114523. https://doi.org/10.1016/j.ssc.2021.114523 doi: 10.1016/j.ssc.2021.114523
![]() |
[6] |
Erande KB, Hawaldar PY, Suryawanshi SR, et al. (2021) Extraction of natural dye (specifically anthocyanin) from pomegranate fruit source and their subsequent use in DSSC. Mater Today Proc 43: 2716-2720. https://doi.org/10.1016/j.matpr.2020.06.357 doi: 10.1016/j.matpr.2020.06.357
![]() |
[7] |
Siregar N, Panggabean JH, Sirait M, et al. (2021) Fabrication of dye-sensitized solar cells (DSSC) using Mg-doped ZnO as photoanode and extract of rose myrtle (Rhodomyrtus tomentosa) as natural dye. Int J Photoenergy 2021: 4033692. https://doi.org/10.1155/2021/4033692 doi: 10.1155/2021/4033692
![]() |
[8] |
Rani AJML, Shanmugasundaram K, Sundaramurthy D, et al. (2021) Correlation study on biopolymer-blended cobalt and iodine gel electrolyte to enhance the efficiency of natural dye-based DSSCs. Energy Fuels 35: 15033-15044. https://doi.org/10.1021/acs.energyfuels.1c02264 doi: 10.1021/acs.energyfuels.1c02264
![]() |
[9] |
González-Verjan VA, Trujillo-Navarrete B, Félix-Navarro RM, et al. (2020) Effect of TiO2 particle and pore size on DSSC efficiency. Mater Renew Sustain Energy 9: 13. https://doi.org/10.1007/s40243-020-00173-7 doi: 10.1007/s40243-020-00173-7
![]() |
[10] |
Subalakshmi K, Senthilselvan J (2018) Effect of fluorine-doped TiO2 photoanode on electron transport, recombination dynamics and improved DSSC efficiency. Sol Energy 171: 914-928. https://doi.org/10.1016/j.solener.2018.06.077 doi: 10.1016/j.solener.2018.06.077
![]() |
[11] |
Umale S, Sudhakar V, Sontakke SM, et al. (2019) Improved efficiency of DSSC using combustion synthesized TiO2. Mater Res Bull 109: 222-226. https://doi.org/10.1016/j.materresbull.2018.09.044 doi: 10.1016/j.materresbull.2018.09.044
![]() |
[12] |
Subalakshmi K, Kumar KA, Paul OP, et al. (2019) Platinum-free metal sulfide counter electrodes for DSSC applications: Structural, electrochemical and power conversion efficiency analyses. Sol Energy 193: 507-518. https://doi.org/10.1016/j.solener.2019.09.075 doi: 10.1016/j.solener.2019.09.075
![]() |
[13] |
Patil DS, Avhad KC, Sekar N (2018) Linear correlation between DSSC efficiency, intramolecular charge transfer characteristics, and NLO properties—DFT approach. Comput Theor Chem 1138: 75-83. https://doi.org/10.1016/j.comptc.2018.06.006 doi: 10.1016/j.comptc.2018.06.006
![]() |
[14] |
Yella A, Lee HW, Tsao HN, et al. (2011) Porphyrin-sensitized solar cells with cobalt (II/III)-based redox electrolyte exceed 12 percent efficiency. Science 334: 629-634. https://doi.org/10.1126/science.1209688 doi: 10.1126/science.1209688
![]() |
[15] |
Mathew S, Yella A, Gao P, et al. (2014) Dye-sensitized solar cells with 13% efficiency achieved through the molecular engineering of porphyrin sensitizers. Nat Chem 6: 242-247. https://doi.org/10.1038/nchem.1861 doi: 10.1038/nchem.1861
![]() |
[16] |
Kakiage K, Aoyama Y, Yano T, et al. (2015) Highly-efficient dye-sensitized solar cells with collaborative sensitization by silyl-anchor and carboxy-anchor dyes. Chem Commun 51: 15894-15897. https://doi.org/10.1039/C5CC06759F doi: 10.1039/C5CC06759F
![]() |
[17] |
Freitag M, Teuscher J, Saygili Y, et al. (2017) Dye-sensitized solar cells for efficient power generation under ambient lighting. Nat Photonics 11: 372-378. https://doi.org/10.1038/nphoton.2017.60 doi: 10.1038/nphoton.2017.60
![]() |
[18] |
Hara K, Wang ZS, Cui Y, et al. (2009) Long-term stability of organic-dye-sensitized solar cells based on an alkyl-functionalized carbazole dye. Energy Environ Sci 2: 1109-1114. https://doi.org/10.1039/b907486d doi: 10.1039/b907486d
![]() |
[19] | Kabir F, Sakib SN, Matin N (2019) Stability study of natural green dye based DSSC. Optik 181 458-464. https://doi.org/10.1016/j.ijleo.2018.12.077 |
[20] |
Kohn S, Wehlage D, Junger IJ, et al. (2019) Electrospinning a dye-sensitized solar cell. Catalysts 9: 975. https://doi.org/10.3390/catal9120975 doi: 10.3390/catal9120975
![]() |
[21] | Chen KF, Liu CH, Huang HK, et al. (2013) Polyvinyl butyral-based thin film polymeric electrolyte for dye-sensitized solar cell with long-term stability. Int J Electrochem Sci 8: 3524-3539. |
[22] |
Yoon JS, Kang Dk, Won JG, et al. (2012) Dye-sensitized solar cells using ion-gel electrolytes for long-term stability. J Power Sources 201: 395-401. https://doi.org/10.1016/j.jpowsour.2011.11.012 doi: 10.1016/j.jpowsour.2011.11.012
![]() |
[23] |
Iftikhar H, Sonai GG, Hashmi SG, et al. (2019) Progress on electrolytes development in dye-sensitized solar cells. Materials 12: 1998. https://doi.org/10.3390/ma12121998 doi: 10.3390/ma12121998
![]() |
[24] |
Dissanayake MAKL, Ekanayake EMBS, Bandara, LRAK, et al. (2016) Efficiency enhancement by mixed cation effect in polyethylene oxide (PEO)-based dye-sensitized solar cells. J Solid State Electrochem 20: 193-201. https://doi.org/10.1007/s10008-015-3018-1 doi: 10.1007/s10008-015-3018-1
![]() |
[25] |
Lei B, Li GR, Chen P, et al. (2019) A quasi-solid-state solar rechargeable battery with polyethylene oxide gel electrolyte. ACS Appl Energy Mater 2: 1000-1005. https://doi.org/10.1021/acsaem.8b02193 doi: 10.1021/acsaem.8b02193
![]() |
[26] |
Wang Y (2009) Recent research progress on polymer electrolytes for dye-sensitized solar cells. Sol Energy Mater Sol Cells 93: 1167-1175. https://doi.org/10.1016/j.solmat.2009.01.009 doi: 10.1016/j.solmat.2009.01.009
![]() |
[27] |
Shi Y, Zhan C, Wang L, et al. (2009) The electrically conductive function of high-molecular weight poly(ethylene oxide) in polymer gel electrolytes used for dye-sensitized solar cells. Phys Chem Chem Phys 11: 4230-4235. https://doi.org/10.1039/b901003c doi: 10.1039/b901003c
![]() |
[28] |
de Freitas JN, Nogueira AF, de Paoli MA (2009) New insights into dye-sensitized solar cells with polymer electrolytes. J Mater Chem 19: 5279. https://doi.org/10.1039/b900928k doi: 10.1039/b900928k
![]() |
[29] |
Tao L, Huo ZP, Ding Y, et al. (2015) High-efficiency and stable quasi-solid-state dye-sensitized solar cell based on low molecular mass organogelator electrolyte. J Mater Chem A 3: 2344-2352. https://doi.org/10.1039/C4TA06188H doi: 10.1039/C4TA06188H
![]() |
[30] |
Bella F, Popovic J, Lamberti A, et al. (2017) Interfacial effects in solid-liquid electrolytes for improved stability and performance of dye-sensitized solar cells. ACS Appl Mater Interfaces 9: 43. https://doi.org/10.1021/acsami.7b11899 doi: 10.1021/acsami.7b11899
![]() |
[31] |
Huo Z, Wang L, Tao L, et al. (2017) A supramolecular gel electrolyte formed from amide based co-gelator for quasi-solid-state dye-sensitized solar cell with boosted electron kinetic processes. J Power Sources 359: 80-87. https://doi.org/10.1016/j.jpowsour.2017.04.099 doi: 10.1016/j.jpowsour.2017.04.099
![]() |
[32] |
Bastianini M, Vivani R, Nocchetti M, et al. (2014) Effect of iodine intercalation in nanosized layered double hydroxides for the preparation of quasi-solid electrolyte in DSSC devices. Sol Energy 107: 692-699. https://doi.org/10.1016/j.solener.2014.06.014 doi: 10.1016/j.solener.2014.06.014
![]() |
[33] |
Bella F, Vlachopoulos N, Nonomura K, et al. (2015) Direct light-induced polymerization of cobalt-based redox shuttles: an ultrafast way towards stable dye-sensitized solar cells. Chem Commun 51: 16308-16311. https://doi.org/10.1039/C5CC05533D doi: 10.1039/C5CC05533D
![]() |
[34] |
Song D, Cho W, Lee JH, et al. (2014) Toward higher energy conversion efficiency for solid polymer electrolyte dye-sensitized solar cells: Ionic conductivity and TiO2 pore-filling. J Phys Chem Lett 5: 1249-1258. https://doi.org/10.1021/jz5002727 doi: 10.1021/jz5002727
![]() |
[35] |
Storck JL, Dotter M, Brockhagen B, et al. (2020) Evaluation of novel glycerol/PEO gel polymer electrolytes for non-toxic dye-sensitized solar cells with natural dyes regarding long-term stability and reproducibility. Crystals 10: 1158. https://doi.org/10.3390/cryst10121158 doi: 10.3390/cryst10121158
![]() |
[36] |
Storck JL, Dotter M, Adabra S, et al. (2020) Long-term stability improvement of non-toxic dye-sensitized solar cells via poly (ethylene oxide). Polymers 12: 3035. https://doi.org/10.3390/polym12123035 doi: 10.3390/polym12123035
![]() |
[37] |
Kobayashi K, Pagot G, Vezzù K, et al. (2021) Effect of plasticizer on the ion-conductive and dielectric behavior of poly(ethylene carbonate)-based Li electrolytes. Polym J 53: 149-155. https://doi.org/10.1038/s41428-020-00397-4 doi: 10.1038/s41428-020-00397-4
![]() |
[38] |
Mustafa MO, Ghareeb HB, Aziz S, et al. (2020) Electrochemical characteristics of glycerolized PEO-based polymer electrolytes. Membranes 10: 116. https://doi.org/10.3390/membranes10060116 doi: 10.3390/membranes10060116
![]() |
[39] |
Nagaraj P, Sasidharan A, David V, et al. (2017) Effect of cross-linking on the performances of starch-based biopolymer as gel electrolyte for dye-sensitized solar cell applications. Polymers 9: 667. https://doi.org/10.3390/polym9120667 doi: 10.3390/polym9120667
![]() |
[40] |
Li Q, Chen X, Tang Q, et al. (2013) Imbibition of polypyrrole into three-dimensional poly (hydroxyethyl methacrylate/glycerol) gel electrolyte for robust quasi-solid-state dye-sensitized solar cells. J Mater Chem A 1: 8055-8060. https://doi.org/10.1039/c3ta11166k doi: 10.1039/c3ta11166k
![]() |
[41] | Kumar M, Kumar A, Kumar S, et al. (2021) Fabrication and characterization of TiO2 based dye-sensitized solar cell. AIP Conf Proc 2352: 010043. |
[42] |
Gossen K, Ehrmann A (2020) Glycerin-based electrolyte for reduced drying of dye-sensitized solar cells. Optik 207: 163772. https://doi.org/10.1016/j.ijleo.2019.163772 doi: 10.1016/j.ijleo.2019.163772
![]() |
[43] |
Junger IJ, Homburg SV, Grethe T, et al. (2017) Examination of the sintering process-dependent properties of TiO2 on glass and textile substrates. J Photonics Energy 7: 015001. https://doi.org/10.1117/1.JPE.7.015001 doi: 10.1117/1.JPE.7.015001
![]() |
[44] |
Hö lscher F, Trümper PR, Junger IJ, et al. (2018) Raising reproducibility in dye-sensitized solar cells under laboratory conditions. J Renew Sustain Energy 10: 013506. https://doi.org/10.1063/1.5013181 doi: 10.1063/1.5013181
![]() |
[45] |
Udomrungkhajornchai S, Junger IJ, Ehrmann A (2020) Optimization of the TiO2 layer in DSSCs by a nonionic surfactant. Optik 203: 163945. https://doi.org/10.1016/j.ijleo.2019.163945 doi: 10.1016/j.ijleo.2019.163945
![]() |
[46] |
Hara K, Wang, ZS, Cui Y, et al. (2009) Long-term stability of organic-dye-sensitized solar cells based on an alkyl-functionalized carbazole dye. Energy Environ Sci 2: 1109-1114. https://doi.org/10.1039/b907486d doi: 10.1039/b907486d
![]() |
[47] |
Jung HS, Koo BK, Kim JY, et al. (2014) Enhanced photovoltaic properties and long-term stability in plasmonic dye-sensitized solar cells via noncorrosive redox mediator. ACS Appl Mater Interfaces 6: 21. https://doi.org/10.1021/am5051982 doi: 10.1021/am5051982
![]() |
[48] |
Sonai GG, Tiihonen A, Miettunen K, et al. (2017) Long-term stability of dye-sensitized solar cells assembled with cobalt polymer gel electrolyte. J Phys Chem C 121: 17577-17585. https://doi.org/10.1021/acs.jpcc.7b03865 doi: 10.1021/acs.jpcc.7b03865
![]() |
[49] |
Ma'alinia A, Moghaddam HA, Nouri E, et al. (2018) Long-term stability of dye-sensitized solar cells using a facile gel polymer electroly. New J Chem 42: 13256-13262. https://doi.org/10.1039/C8NJ02157K doi: 10.1039/C8NJ02157K
![]() |
[50] |
Griffini G, Bella F, Nisic F, et al. (2014) Multifunctional luminescent down-shifting fluoropolymer coatings: A straightforward strategy to improve the UV-light harvesting ability and long-term outdoor stability of organic dye-sensitized solar cells. Adv Energy Mater 5: 1401312. https://doi.org/10.1002/aenm.201401312 doi: 10.1002/aenm.201401312
![]() |
[51] |
Önen T, Karakus MÖ, Coskun R, et al. (2019) Reaching stability at DSSCs with new type gel electrolytes. J Photoch Photobio A 385: 112082. https://doi.org/10.1016/j.jphotochem.2019.112082 doi: 10.1016/j.jphotochem.2019.112082
![]() |
[52] |
Jang YJ, Thogiti S, Lee K, et al. (2019) Long-term stable solid-state dye-sensitized solar cells assembled with solid-state polymerized hole-transporting material. Crystals 9: 452. https://doi.org/10.3390/cryst9090452 doi: 10.3390/cryst9090452
![]() |
[53] |
Roy A, Ghosh A, Bhandari S, et al. (2019) Color comfort evaluation of dye-sensitized solar cell (DSSC) based building-integrated photovoltaic (BIPV) glazing after 2 years of ambient exposure. J Phys Chem C 123: 23834-23837. https://doi.org/10.1021/acs.jpcc.9b05591 doi: 10.1021/acs.jpcc.9b05591
![]() |
1. | Rami Ahmad El-Nabulsi, Waranont Anukool, Nonlocal thermal effects on biological tissues and tumors, 2022, 34, 24519049, 101424, 10.1016/j.tsep.2022.101424 | |
2. | Ivanna Andrusyak, Oksana Brodyak, Petro Pukach, Myroslava Vovk, Mathematical Modeling of Cell Growth via Inverse Problem and Computational Approach, 2024, 12, 2079-3197, 26, 10.3390/computation12020026 | |
3. | O. E. Adebayo, S. Urcun, G. Rolin, S. P. A. Bordas, D. Trucu, R. Eftimie, Mathematical investigation of normal and abnormal wound healing dynamics: local and non-local models, 2023, 20, 1551-0018, 17446, 10.3934/mbe.2023776 | |
4. | Johan Marguet, Raluca Eftimie, Alexei Lozinski, Numerical approaches for non-local transport-dominated PDE models with applications to biology, 2025, 44, 2238-3603, 10.1007/s40314-025-03146-6 |
Parameter | Value | Description | Reference |
D1 | 0.00675 | diffusion of primary tumour | [11] |
D2 | 0.00675 | diffusion of secondary tumour | [11] |
η1 | 2.85×10−2 | haptotaxis to ECM from c1 | [37] |
η2 | 2.85×10−2 | haptotaxis to ECM from c2 | [37] |
μc | 0.25 | proliferation of tumour cells c | [38] |
Kc | 1 | tissue carrying capacity | [16] |
ρ | 2 | ECM degradation coefficient | [39] |
μv | 0.40 | ECM remodelling coefficient | [16] |
Scc | (0.5000.3) | cell-cell adhesion function | [11] |
Scv | (0.1000.5) | cell matrix adhesion function | [11] |
t1,2 | 10 | time initiation for mutations | [38,40] |
ts | 3 | time-steepness coefficient | [38] |
δ0 | 0.3 | mutation from primary tumour | [38,40] |
Δx | 0.03125 | discretisation step size for GΩ | [16] |
Δt | 10−3 | time step size | [16] |
Δη | 0.0625 | mesh size used for G1M and G2M | Estimated |
Parameter | Value | Description | Reference |
D1 | 0.00675 | diffusion of primary tumour | [11] |
D2 | 0.00675 | diffusion of secondary tumour | [11] |
η1 | 2.85×10−2 | haptotaxis to ECM from c1 | [37] |
η2 | 2.85×10−2 | haptotaxis to ECM from c2 | [37] |
μc | 0.25 | proliferation of tumour cells c | [38] |
Kc | 1 | tissue carrying capacity | [16] |
ρ | 2 | ECM degradation coefficient | [39] |
μv | 0.40 | ECM remodelling coefficient | [16] |
Scc | (0.5000.3) | cell-cell adhesion function | [11] |
Scv | (0.1000.5) | cell matrix adhesion function | [11] |
t1,2 | 10 | time initiation for mutations | [38,40] |
ts | 3 | time-steepness coefficient | [38] |
δ0 | 0.3 | mutation from primary tumour | [38,40] |
Δx | 0.03125 | discretisation step size for GΩ | [16] |
Δt | 10−3 | time step size | [16] |
Δη | 0.0625 | mesh size used for G1M and G2M | Estimated |