Kamran Arif , Guojiang Xi , Heng Wang , Weihua Deng
{"title":"具有厚尾lsamvy测度的高维偏微分方程的深度学习","authors":"Kamran Arif , Guojiang Xi , Heng Wang , Weihua Deng","doi":"10.1016/j.jcp.2025.114327","DOIUrl":null,"url":null,"abstract":"<div><div>The partial differential equations (PDEs) for jump process with Lévy measure have wide applications. When the measure has fat tails, it will bring big challenges for both computational cost and accuracy. In this work, we develop a deep learning method for high-dimensional PDEs related to fat-tailed Lévy measure, which can be naturally extended to the general case. Building on the theory of backward stochastic differential equations for Lévy processes, our deep learning method avoids the need for neural network differentiation and introduces a novel technique to address the singularity of fat-tailed Lévy measures. The developed method is used to solve four kinds of high-dimensional PDEs: the diffusion equation with fractional Laplacian; the advective diffusion equation with fractional Laplacian; the advective diffusion reaction equation with fractional Laplacian; and the nonlinear reaction diffusion equation with fractional Laplacian. The parameter <span><math><mi>β</mi></math></span> in fractional Laplacian is an indicator of the strength of the singularity of Lévy measure. Specifically, for <span><math><mrow><mi>β</mi><mo>∈</mo><mo>(</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>)</mo></mrow></math></span>, the model describes super-ballistic diffusion; while for <span><math><mrow><mi>β</mi><mo>∈</mo><mo>(</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>)</mo></mrow></math></span>, it characterizes super-diffusion. In addition, we experimentally verify that the developed algorithm can be easily extended to solve fractional PDEs with finite general Lévy measures. Our method achieves a relative error of <span><math><mrow><mi>O</mi><mo>(</mo><msup><mn>10</mn><mrow><mo>−</mo><mn>3</mn></mrow></msup><mo>)</mo></mrow></math></span> for low-dimensional problems and <span><math><mrow><mi>O</mi><mo>(</mo><msup><mn>10</mn><mrow><mo>−</mo><mn>2</mn></mrow></msup><mo>)</mo></mrow></math></span> for high-dimensional ones. We also investigate three factors that influence the algorithm’s performance: the number of hidden layers; the number of Monte Carlo samples; and the choice of activation functions. Furthermore, we test the efficiency of the algorithm in solving problems in 3D, 10D, 20D, 50D, and 100D. Our numerical results demonstrate that the algorithm achieves excellent performance with deeper hidden layers, a larger number of Monte Carlo samples, and the Softsign activation function.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"541 ","pages":"Article 114327"},"PeriodicalIF":3.8000,"publicationDate":"2025-08-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Deep learning for high-dimensional PDEs with fat-tailed Lévy measure\",\"authors\":\"Kamran Arif , Guojiang Xi , Heng Wang , Weihua Deng\",\"doi\":\"10.1016/j.jcp.2025.114327\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><div>The partial differential equations (PDEs) for jump process with Lévy measure have wide applications. When the measure has fat tails, it will bring big challenges for both computational cost and accuracy. In this work, we develop a deep learning method for high-dimensional PDEs related to fat-tailed Lévy measure, which can be naturally extended to the general case. Building on the theory of backward stochastic differential equations for Lévy processes, our deep learning method avoids the need for neural network differentiation and introduces a novel technique to address the singularity of fat-tailed Lévy measures. The developed method is used to solve four kinds of high-dimensional PDEs: the diffusion equation with fractional Laplacian; the advective diffusion equation with fractional Laplacian; the advective diffusion reaction equation with fractional Laplacian; and the nonlinear reaction diffusion equation with fractional Laplacian. The parameter <span><math><mi>β</mi></math></span> in fractional Laplacian is an indicator of the strength of the singularity of Lévy measure. Specifically, for <span><math><mrow><mi>β</mi><mo>∈</mo><mo>(</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>)</mo></mrow></math></span>, the model describes super-ballistic diffusion; while for <span><math><mrow><mi>β</mi><mo>∈</mo><mo>(</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>)</mo></mrow></math></span>, it characterizes super-diffusion. In addition, we experimentally verify that the developed algorithm can be easily extended to solve fractional PDEs with finite general Lévy measures. Our method achieves a relative error of <span><math><mrow><mi>O</mi><mo>(</mo><msup><mn>10</mn><mrow><mo>−</mo><mn>3</mn></mrow></msup><mo>)</mo></mrow></math></span> for low-dimensional problems and <span><math><mrow><mi>O</mi><mo>(</mo><msup><mn>10</mn><mrow><mo>−</mo><mn>2</mn></mrow></msup><mo>)</mo></mrow></math></span> for high-dimensional ones. We also investigate three factors that influence the algorithm’s performance: the number of hidden layers; the number of Monte Carlo samples; and the choice of activation functions. Furthermore, we test the efficiency of the algorithm in solving problems in 3D, 10D, 20D, 50D, and 100D. Our numerical results demonstrate that the algorithm achieves excellent performance with deeper hidden layers, a larger number of Monte Carlo samples, and the Softsign activation function.</div></div>\",\"PeriodicalId\":352,\"journal\":{\"name\":\"Journal of Computational Physics\",\"volume\":\"541 \",\"pages\":\"Article 114327\"},\"PeriodicalIF\":3.8000,\"publicationDate\":\"2025-08-29\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Journal of Computational Physics\",\"FirstCategoryId\":\"101\",\"ListUrlMain\":\"https://www.sciencedirect.com/science/article/pii/S0021999125006096\",\"RegionNum\":2,\"RegionCategory\":\"物理与天体物理\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q2\",\"JCRName\":\"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0021999125006096","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
Deep learning for high-dimensional PDEs with fat-tailed Lévy measure
The partial differential equations (PDEs) for jump process with Lévy measure have wide applications. When the measure has fat tails, it will bring big challenges for both computational cost and accuracy. In this work, we develop a deep learning method for high-dimensional PDEs related to fat-tailed Lévy measure, which can be naturally extended to the general case. Building on the theory of backward stochastic differential equations for Lévy processes, our deep learning method avoids the need for neural network differentiation and introduces a novel technique to address the singularity of fat-tailed Lévy measures. The developed method is used to solve four kinds of high-dimensional PDEs: the diffusion equation with fractional Laplacian; the advective diffusion equation with fractional Laplacian; the advective diffusion reaction equation with fractional Laplacian; and the nonlinear reaction diffusion equation with fractional Laplacian. The parameter in fractional Laplacian is an indicator of the strength of the singularity of Lévy measure. Specifically, for , the model describes super-ballistic diffusion; while for , it characterizes super-diffusion. In addition, we experimentally verify that the developed algorithm can be easily extended to solve fractional PDEs with finite general Lévy measures. Our method achieves a relative error of for low-dimensional problems and for high-dimensional ones. We also investigate three factors that influence the algorithm’s performance: the number of hidden layers; the number of Monte Carlo samples; and the choice of activation functions. Furthermore, we test the efficiency of the algorithm in solving problems in 3D, 10D, 20D, 50D, and 100D. Our numerical results demonstrate that the algorithm achieves excellent performance with deeper hidden layers, a larger number of Monte Carlo samples, and the Softsign activation function.
期刊介绍:
Journal of Computational Physics thoroughly treats the computational aspects of physical problems, presenting techniques for the numerical solution of mathematical equations arising in all areas of physics. The journal seeks to emphasize methods that cross disciplinary boundaries.
The Journal of Computational Physics also publishes short notes of 4 pages or less (including figures, tables, and references but excluding title pages). Letters to the Editor commenting on articles already published in this Journal will also be considered. Neither notes nor letters should have an abstract.