[論文レビュー] Pandemic modeling and the renormalization group equations: Effect of contact matrices, fixed points and nonspecific vaccine waning

Michael McGuigan|arXiv (Cornell University)|Aug 3, 2020

COVID-19 epidemiological studies参考文献 21被引用数 5

ひとこと要約

本論文は、感染症モデルにおけるコンpartimental SIRDモデルと素性物理学におけるレノルミゼーション群（RG）方程式の類似性を明らかにし、接触行列によるオペレータ混合、フロー図、固定点といった共通する数学的構造を示している。これらの類似性を用いて、年齢層間の接触、国境を越えた伝播、および効果の持続期間が短い非特異的ワクチンの影響をモデル化した。10万人規模のモデルにおいて非特異的ワクチンを導入した場合、合計死亡者数は1,509人から1,369人に13%減少した。特に若年層ではワクチン効果の持続期間が長く、その影響が顕著であった。

ABSTRACT

In this paper we find common features between the equations that are used for pandemic or epidemic modeling and the renormalization group equations that are used in high energy physics. Some of these features include the relation of contact matrices in pandemic modeling and operator mixing in the renormalization group equations. Another common feature are the use of flow diagrams and the study of fixed points both in pandemic modeling and in evolution under renormalization group equations. We illustrate these relations through the study of some cases of interest to the current COVID-19 pandemic. These include pandemic modeling with mixing between different age groups and also contact matrices associated with contact between countries. For the final example we study the effect on mortality of waning from nonspecific vaccines which are designed to combat different pathogens but nevertheless may lessen the severity and mortality of COVID-19 infections.

研究の動機と目的

• 感染症のコンパートメントモデルと素性物理学におけるレノルミゼーション群（RG）方程式の間の共通する数学的構造を特定すること。
• フローダイアグラムや固定点解析といったRG技術を、感染症の伝播ダイナミクスの理解を深めるためにパンデミックモデルに適用すること。
• 免疫効果の持続期間が短い非特異的ワクチンの影響が、年齢層および国境を越えた接触の文脈で死亡率に与える影響をモデル化すること。
• 年齢層および国間の接触行列が、RGにインspiredされたフレームワークを用いてパンデミックの広がりに与える影響を検討すること。
• 伝播率βの変動に伴う死亡率予測のパラメータ感度と不確実性を評価すること。

提案手法

• 年齢別および国別接触行列を組み込んだ微分方程式を用いたSIRDコンパートメントモデルを用いて、感染、回復、死亡のダイナミクスをシミュレートする。
• スケールの進化、固定点、オペレータ混合といったレノルミゼーション群の概念を応用し、集団のコンパートメント間の遷移を量子場理論における結合定数の進化に類似したものとして解釈する。
• 感染率に影響を与えないが死亡率（μ）を低下させる非特異的ワクチン効果をモデル化し、その効果の時間的低下を、時間依存の有効性低下として表現する。
• Mathematica/MATLABを用いた数値シミュレーションにより、変動するβ（伝播率）、初期条件、ワクチンパラメータを用いたSIRD方程式を解く。
• 感染数および死亡者数の時間的推移をフローダイアグラムと時間的推移曲線で可視化し、非特異的ワクチンの有無によるシナリオを比較する。
• βを±10%変動させることで感度分析を行い、死亡率予測における不確実性を評価する。

実験結果

リサーチクエスチョン

• RQ1 パンデミックモデルにおける接触行列は、レノルミゼーション群方程式におけるオペレータ混合に対応するか？
• RQ2 RG方程式における固定点解析を用いて、感染症モデルにおける平衡状態を特定できるか？
• RQ3 持続期間が短い非特異的ワクチンの導入が、異なる年齢層における死亡率に与える影響は何か？
• RQ4 国境を越えた伝播ダイナミクスを接触行列でモデル化した場合、パンデミックの広がりに与える影響は何か？
• RQ5 伝播率βの不確実性に起因する死亡率予測の感度はどの程度か？

主な発見

• 非特異的ワクチンを導入したSIRDモデルでは、10万人規模の集団における合計死亡者数が1,509人から1,369人に13%減少した。
• 若年層（0–49歳）では死亡者数の割合的減少が顕著で、158人から86人に減少した。これはワクチン効果の持続期間が長く、若年層に効果的であることを示している。
• 高年齢層（50歳以上）では、ベースラインの死亡率が高く、ワクチン効果の持続期間が短いため、死亡者数の絶対的減少は小さく（1,351人から1,283人に）、効果が限定的であった。
• 非特異的ワクチン導入後も1日あたりの感染数に変化はなく、ワクチンが感染を防ぐのではなく、重症度と死亡率の低下に寄与していることが確認された。
• βを10%増加させた場合、ワクチン接種済み・未接種両シナリオで死亡者数が約5%増加した。これは伝播率βの不確実性に起因する死亡率予測の感度を示している。
• 本モデルは、フローダイアグラムや固定点解析といったRGにインspiredされた技術が、安定した均衡状態や伝播トレンドの特定に効果的に応用可能であることを示している。