DEEP learning is a sub-discipline of artificial intelligence that uses machine learning algorithms based on artificial neural networks to create patterns and make predictions from large data sets [1] . We propose to solve PDEs by approximating the solution with a deep neural network (DNN) which is trained to satisfy the differential operator, initial condition, and boundary conditions [2] .

The increasing adoption of deep learning across healthcare domains together with the availability of highly characterized cancer datasets has accelerated research into the utility of deep learning in the analysis of the complex biology of cancer. In this work, the dynamical system of tumor cells driving in the human tissues are investigated. With the introduction of the deep learning framework, there have been numerous attempts to create more efficient approaches of numerical solutions of PDEs by using a graph kernel algorithm and a spectral decomposition of the Laplacian graph approach [2] . A cancel cell network of distribution into tissues can be represented as a dynamic graph. In this paper, we present a graph neural operator-based method to define mathematical problems to model the dynamical model of diffusion of tumor cells.

Cancer invasion and the ability of malignant tumor cells for directed migration and metastasis have remained a focus of research for many years. Numerous studies have confirmed the existence of two main patterns of cancer cell invasion: collective cell migration and individual cell migration, by which tumor cells overcome barriers of the extracellular matrix and spread into surrounding tissues. Each pattern of cell migration displays specific morphological features and the biochemical/molecular genetic mechanisms underlying cell migration. Cancer can also spread from where it first started to other parts of the body. This process is called metastasis [3] . Cancer cells can metastasize when they break away from the tumor and travel to a new location in the body through the blood or lymphatic system. Where cancer can spread and staging. Most cancers have a tendency to spread to certain areas of the body. This has helped doctors develop staging systems that are used to classify cancers based on information about where the cancer is in the body and if it has spread from where it started [4] .

Immunotherapies have revolutionized treatments for several types of cancers, but have been met with clinical trial failures in others. To improve immunotherapeutic approaches requires understanding how the immune system interacts around and within a tumor, allowing us to establish effective immunotherapeutic protocols. Mathematical modelling can help identify the mechanisms at the heart of immunotherapeutic efficacy and design successful therapeutic regimens.

The article provides a comprehensive investigation into the application of Artificial Neural Networks (ANNs) for solving boundary value problems (BVPs) associated with biomechanical systems, specifically tumor elasticity. The work bridges mathematical oncology and neural network methodologies, presenting an innovative approach to addressing challenges in cancer research.

Relevance and Innovation of this article is the focus on ANN-based solutions for PDEs in complex geometries demonstrates originality, offering potential advancements in personalized cancer treatment models. The methodological approach is based on the development of neural network architectures (Fourier neural operators) and its application in modeling tumor biomechanics.

The mechanics of cancer cells can be described using linear elasticity theory for biomechanical systems [5] . In contrast to the molecular dynamics simulations, finite element analysis requires less computation power and can be used to model larger portions or the entire length of cancer. Therefore, finite element analysis is a useful tool to approximate and estimate the mechanical properties diffuse interface models have gained a growing interest in cancer research for their ability to investigate the mechanic-biological features during tumor progression and to provide simulation tools for personalized anti-cancer strategies at an affordable computational cost. Here we propose a diffuse interface model for tumors evolution which accounts for an interfacial structure mimicking a finite elastic confinement at the tumor boundary, possibly due either to a localized elastic stress induced by host tissue displacements area. This model consists of the system of a partial differential equation of thermos-elasticity.

In contrast to the molecular dynamics simulations, finite element analysis requires less computation power and can be used to model larger portions or the entire length of microcells. Therefore, finite element analysis is a useful tool to approximate and estimate the mechanical properties of the microtubules. Fourier neural operators (FNOs) are used to create a neural operator pseudo oscillation system of PDEs. The relaxation dynamics of the cancer cells are described by the following system of partial differential equations:

$\forall x\in D^{+}L\left(\partial x,\tau \right)U\left(x,\tau \right)=F\left(x\right)$

$\forall z\in \partial D{\left[U\right]}^{+}=0$ (1)

where $U=\left(u,u_3\right)$, where $u=\left(u_1,u_2\right)$ is a displacement vector, $u_3$ is a temperature variation, $x=\left(x_j\right);y=\left(y_j\right);j=1,2$-points of $E^2$-two-dimensional Euclidean space, ( $D^{+}\subset E^2$, $D^{\left(0\right)}=E^2\backslash D^{+}\cup \partial D$) bounded by the close surface $S\in L_{\left(2\right)}\left(\alpha \right),\alpha >0$ of Holder class with outward positive normal vector.

$L\left(\partial x,\tau \right)={\left|L_{jk}\left(\partial x,\tau \right)\right|}_{3\times 3}$ is neural operator for PDEs of $L\left(\partial x,\partial t\right)$-elliptic operator of PDEs of biomechanical system described as follows:

$L\left(\partial x,\tau \right)=\Vert \begin{array}{cc}\left(\lambda +2\mu \right)\Delta -\varsigma {\tau }^2& -{\gamma }_{\tau }\Delta \\ -\eta \tau & \Delta -\frac{\tau }{{\wp }_{\tau }}\end{array}\Vert$ (2)

where:

$\varsigma >0,\mu >0,3\lambda +2\mu >0,\wp >0,\frac{\gamma }{\eta }>0$ are the constants of elasticity, temperature and diffusion parameters in D⁺ domain ${\gamma }_{\tau }={\gamma }_{\tau }\left(1+{\tau }_1\tau \right)$, $\frac{1}{{\Re }_{\tau }}=\frac{1}{{\Re }_{\tau }}\left(1+{\tau }_0\tau \right)$ ${\tau }_1\ge {\tau }_0>0$ constants of relaxation, Δ is a two-dimensional Laplacian operator; $\tau =\sigma +iq,\sigma >0$ (corresponds to the general dynamical problems) [6] [7] , $F=\left(F_1,F_2,F_3\right)\subset C^{0,\alpha }\left(D^{+}\right),\alpha >0$ is a given vectors.

Constructing the analytic Green’s function for a given domain is in general a difficult matter, but accordingly theory of elasticity the (1) problem has regular solutions and Green’s function can be represented in quadrature [9] .

In order (2) neural Green’s function is derived as follows:

$G\left(x,y;\tau ,D^{+}\right)={\displaystyle {\sum }_{j=1}^{\infty }{\omega }_j{\phi }_j\left(x\right){\phi }_j\left(y\right)}$ (4)

where $\left({\omega }_j,{\phi }_j\left(x\right)\right)$ is accountable sequence of Eigen pairs.

As a result of (2) and (3) when utilizing Green’s function to obtain the solution, the exact solution and its numerical approximation can be represented through

$U\left(x,\tau \right)=-\frac{1}{2}{\displaystyle {\sum }_{j=1}^{\infty }{\omega }_j}\left({\displaystyle {\int }_{D^{+}}{\phi }_j\left(y\right)F\left(y\right)\text{d}y}\right){\phi }_j\left(x\right)$ (5)

$U^N\left(x,\tau \right)=-\frac{1}{2}{\displaystyle {\sum }_{j=1}^N{\omega }_j\left({\displaystyle {\int }_{D^{+}}{\phi }_j\left(y\right)F\left(y\right)\text{d}y}\right){\phi }_j\left(x\right)}$ (6)

In order to assumptions of Green’s tensor [9] and applying the Cauchy-Schwarz inequality [4] , we have

$U\left(x,\tau \right)=\underset{N\to \infty }{\lim }{U}^N\left(x,\tau \right)$ (7)

Let us consider particular case of boundary conditions for (1) in finite space $\left(0\le \left|x\right|\le 1\right)$, with respect to ${\phi }_j\left(x\right)$ we have: ${\phi }_j\left(0\right)={\phi }_j\left(1\right)=0$.

Let us introduce the following notation:

$G\left(x,y\right)=\hat{G}\left(x,y,r\left(x,y\right)\right),\forall x\in D^{+},\forall y\in D^{+}\cup \partial D,r=\left|x-y\right|$

In order to investigate the spectral bias in the training dynamics related to (2) expand the approximation error during the training process in terms of the orthonormal eigenfunctions in the following form [10] :

$\left(x,y;\tau ,D^{+}\right)-G\left(x,y;\tau ,D^{+}\right)={\displaystyle {\sum }_{j=1}^{\infty }{\theta }_j{\phi }_j\left(x\right){\phi }_j\left(y\right)}$ (8)

where

${\theta }_j={\displaystyle {\int }_0^1{\displaystyle {\int }_0^1\left(\hat{G}\left(x,y\right)-G\left(x,y\right)\right){\phi }_j\left(x\right){\phi }_j\left(y\right)\text{d}y\text{d}x}}\text{in}\left(0,1\right)$

Here we only focus on the interior approximation error, as the Dirichlet boundary condition can be enforced by defining the $L^2$ loss over high frequencies $j\ge \eta$ for our neural approximation

$L_k={\displaystyle {\sum }_{j=k}^{\infty }{\theta }_j^2}$ (9)

Let us consider numerical experiments to show effectiveness of our approach in solving the one-dimensional indefinite boundary value problems, throughout the training process, we employ a fully-connected neural network with a tanh activation function [10] [11] .

As a result, the approximate solution $\hat{U}$ of problem (1) can be immediately obtained using our neural Green’s function ( Figure 1 and Figure 2 ). The numerical solution using our neural Green’s function has an approximation error ( Figure 3 ).

Let us consider the particular case of boundary conditions with respect to (1).

Problem. It is required to find in D the regular solution $U={\left(u,u_3\right)}^{\text{T}}={\Vert u_k\Vert }_{3\times 1}$, of system (1), which satisfies the boundary condition; the displacement vector and temperature variations are given on the $x_2=0$ line:

$u\left(x_1,0\right)={\psi }^{\left(1\right)}\left(x_1\right),u_3\left(x_1,0\right)={\psi }_3\left(x_1\right)$.

Therefore the solutions are described as follows:

$u_1=\frac{1}{\pi }{\displaystyle \underset{-\infty }{\overset{+\infty }{\int }}\frac{\partial }{\partial n\left(y\right)}H\left(i\tau \sqrt{\frac{\varsigma +\eta }{\lambda +2\mu }}\left|x-y\right|\right){\psi }_1\left(y\right)\text{d}y}-\frac{1}{\pi }{\displaystyle \underset{x_2>0}{\iint }{G}_1\left(x,z\right)F_1\left(z\right)\text{d}z}$

$u_2=\frac{1}{\pi }{\displaystyle \underset{-\infty }{\overset{+\infty }{\int }}\frac{\partial }{\partial n\left(y\right)}H\left(i\tau \sqrt{\frac{\varsigma +\eta }{\lambda +2\mu }}\left|x-y\right|\right){\psi }_2\left(y\right)\text{d}y}-\frac{1}{\pi }{\displaystyle \underset{x_2>0}{\iint }{G}_2\left(x,z\right)F_2\left(z\right)\text{d}z}$

$u_3=\frac{1}{\pi }{\displaystyle \underset{-\infty }{\overset{+\infty }{\int }}\frac{\partial }{\partial n\left(y\right)}H\left(i\tau \sqrt{\frac{\tau }{{\gamma }_{\tau }{\wp }_{\tau }}}\left|x-y\right|\right){\psi }_3\left(y\right)\text{d}y}-\frac{1}{\pi }{\displaystyle \underset{x_2>0}{\iint }{G}_3\left(x,z\right)F_3\left(z\right)\text{d}z}$

where

$G_1\left(x,z\right)=H_0^{\left(1\right)}\left(i\tau \sqrt{\frac{\varsigma +\eta }{\lambda +2\mu }}\left|x-z\right|\right)+H_0^{\left(1\right)}\left(i\tau \sqrt{\frac{\varsigma +\eta }{\lambda +2\mu }}\left|x^{*}-z\right|\right)$

$G_2\left(x,z\right)=H_0^{\left(1\right)}\left(i\tau \sqrt{\frac{\varsigma +\eta }{\lambda +2\mu }}\left|x-z\right|\right)-H_0^{\left(1\right)}\left(i\tau \sqrt{\frac{\varsigma +\eta }{\lambda +2\mu }}\left|x^{*}-z\right|\right)$

$G_3\left(x,z\right)=H_0^{\left(1\right)}\left(i\tau \sqrt{\frac{\tau }{{\gamma }_{\tau }{\wp }_{\tau }}}\left|x-y\right|\right)+H_0^{\left(1\right)}\left(i\tau \sqrt{\frac{\tau }{{\gamma }_{\tau }{\wp }_{\tau }}}\left|x^{*}-z\right|\right)$

where $x=\left(x_1,x_2\right),x^{*}=\left(x_1,-x_2\right)$, $H_0^{\left(1\right)}\left(\right)$-is Hankel Function of the First Kind (the zero order) [11] .

In this work, a novel framework for solving regular PDEs using Deep Neural Networks is developed, extended, and analyzed. A novel data-enrichment algorithm into neural Green’s function framework is presented, Benefits of graph embedding in treatment being able to represent data using graph embedding offers great benefits, including:

The method to modify the boundary problem through a deep learning algorithm to perform the long-time simulation for the rogue wave by the superior numerical errors is developed. The proposed method gives rise to better numerical results with MI in comparison with the ones obtained by traditional methods.