2025-11-22T04:13:15.983918

A Particle-Field Algorithm with Neural Interpolation for a Parabolic-Hyperbolic Chemotaxis System in 3D

Kim, Xin

Tumor angiogenesis involves a collection of tumor cells moving towards blood vessels for nutrients to grow. Angiogenesis, and in general chemo- taxis, systems have been modeled using partial differential equations (PDEs) and as such require numerical methods to approximate their solutions. Here we study a Parabolic-Hyperbolic Keller-Segel (PHKS) system in three space dimensions. The model arises in the angiogenesis literature. To compute solutions to the PHKS system, we develop a neural stochastic interacting particle-field (NSIPF) method where the density variable is represented as empirical measures of particles and the field variable (concentration of chemoattractant) approximated by a convolutional neural network (CNN). We discuss the performance of NSIPF in computing multi-bump solutions to the system.

academic

3次元放物線-双曲型化学走性系システムのための神経補間を用いた粒子場アルゴリズム

基本情報

論文ID: 2510.13199
タイトル: A Particle-Field Algorithm with Neural Interpolation for a Parabolic-Hyperbolic Chemotaxis System in 3D
著者: Jongwon David Kim, Jack Xin (カリフォルニア大学アーバイン校)
分類: math.NA cs.NA
発表日時: 2025年10月15日 (arXiv プレプリント)
論文リンク: https://arxiv.org/abs/2510.13199v1

要旨

腫瘍血管新生は、腫瘍細胞が栄養物質を獲得するために血管に向かって移動するプロセスを伴う。血管新生および一般的な化学走性系は偏微分方程式(PDE)によってモデル化されているため、数値解を近似するための数値方法が必要である。本論文は3次元空間における放物線-双曲型Keller-Segel (PHKS)系を研究し、このモデルは血管新生文献に由来する。PHKS系を解くために、著者らは神経確率相互作用粒子場(NSIPF)法を開発した。この方法では、密度変数は粒子の経験測度として表現され、場変数(化学走性因子濃度)は畳み込みニューラルネットワーク(CNN)によって近似される。本論文はNSIPFが系の多峰解を計算する際の性能について論じている。

研究背景と動機

解決すべき問題:
- 腫瘍血管新生における細胞の血管への移動の数学的モデル化と数値解法
- 3次元空間における放物線-双曲型Keller-Segel系の効率的な数値計算
問題の重要性:
- 血管新生は腫瘍転移の重要な生物学的プロセスであり、その機構の理解は癌治療に重要な意義を持つ
- 化学走性系は細菌採食、組織形成など複数の生物学的プロセスで重要な役割を果たす
既存方法の限界:
- 従来の有限差分法(FDM)は3次元高分解能計算において計算時間が長い
- 古典的なスプライン補間法は3次元空間での実装が複雑で計算コストが高い
- 放物線-双曲型系に特化した効率的な粒子法が不足している
研究動機:
- ニューラルネットワークの高次元近似能力を利用して補間効率を改善する
- 粒子法のメッシュフリーの利点とCNNの補間性能を組み合わせる
- 3次元化学走性系に対してより効率的な数値解法を提供する

核心的貢献

神経確率相互作用粒子場(NSIPF)法を提案し、CNN補間器を確率粒子場アルゴリズムに統合した
PHKS系の混沌伝播理論基礎を確立し、適度な相互作用を通じて神経補間に理論的支援を提供した
効率的なCNNアーキテクチャを開発し、放射状対称解の訓練データを使用して、2次元および3次元系の補間に成功した
著しい計算効率の向上を実現し、従来の方法と比較して解の品質を保ちながら計算コストを大幅に削減した

方法の詳細説明

タスク定義

3次元空間における放物線-双曲型Keller-Segel系を解く:

ρₜ = ∇ · (γ∇ρ - χρ∇c)  (1)
cₜ = -cρ                  (2)

ここで、ρは細菌密度、cは化学走性因子濃度、γは拡散係数、χは化学走性感受性である。

モデルアーキテクチャ

1. 粒子表現

密度関数は粒子の経験測度によって近似される:

ρₜ ≈ (M₀/P) Σⱼ₌₁ᴾ δ(x - Xⱼₜ)

ここで、Pは粒子数、M₀は保存質量である。

2. SDE表現

粒子の運動は確率微分方程式に従う:

dXᵢₜ = χ∇c(Xᵢₜ,t)dt + √(2γ)dBᵢₜ

3. CNN補間器アーキテクチャ

入力層: 低分解能濃度場データ
畳み込み層: 5層の3次元畳み込み、チャネル数1→16→32→32→32→16→1
活性化関数: ReLU
スキップ接続: 勾配伝播を強化
出力: 高分解能補間場

4. NSIPFアルゴリズムフロー

アルゴリズム1: 神経SIPF
1. ρ₀に基づいて粒子位置を初期化
2. 各時間ステップについて:
   a. 粒子をビニングしてρₙを得る
   b. 濃度を更新: cⁿ⁺¹ = cⁿ - δt·cⁿ·ρⁿ
   c. CNNで補間してc(x,tₙ)を得る
   d. ∇cを計算して粒子位置を更新

技術的革新点

古典的補間を神経補間に置き換え:
- スプライン補間の代わりにCNNを使用して濃度場勾配を計算
- 訓練データは1次元放射状系の廉価な解から得られる
- 非放射状初期条件への一般化をサポート
混沌伝播理論による支援:
- mollifier正則化を通じて理論基礎を確立
- 適度な相互作用が粒子系の収束性を保証
- CNN補間の合理性に数学的根拠を提供
効率的な訓練戦略:
- 放射状解のデータ拡張(ダウンサンプリング、平行移動、ぼかし)を使用
- パッチベースの訓練により3次元計算負担を削減
- 100エポックの訓練で良好な補間効果を達成