Abstract:The Fokker-Planck (FP) equation governing the evolution of the probability density function (PDF) is applicable to many disciplines but it requires specification of the coefficients for each case, which can be functions of space-time and not just constants, hence requiring the development of a data-driven modeling approach. When the data available is directly on the PDF, then there exist methods for inverse problems that can be employed to infer the coefficients and thus determine the FP equation and subsequently obtain its solution. Herein, we address a more realistic scenario, where only sparse data are given on the particles' positions at a few time instants, which are not sufficient to accurately construct directly the PDF even at those times from existing methods, e.g., kernel estimation algorithms. To this end, we develop a general framework based on physics-informed neural networks (PINNs) that introduces a new loss function using the Kullback-Leibler divergence to connect the stochastic samples with the FP equation, to simultaneously learn the equation and infer the multi-dimensional PDF at all times. In particular, we consider two types of inverse problems, type I where the FP equation is known but the initial PDF is unknown, and type II in which, in addition to unknown initial PDF, the drift and diffusion terms are also unknown. In both cases, we investigate problems with either Brownian or Levy noise or a combination of both. We demonstrate the new PINN framework in detail in the one-dimensional case (1D) but we also provide results for up to 5D demonstrating that we can infer both the FP equation and} dynamics simultaneously at all times with high accuracy using only very few discrete observations of the particles.

Time-dependent parameter identification in a Fokker-Planck equation based magnetization model of large ensembles of nanoparticles

Approaches for modeling magnetic nanoparticle dynamics

High-dimensional Bayesian Parameter Estimation: Case Study for a Model of JAK2/STAT5 Signaling.

Equilibrium Model with Anisotropy for Model-Based Reconstruction in Magnetic Particle Imaging

Simulations of magnetic nanoparticle Brownian motion

Modeling of 2D Diffusion Processes Based on Microscopy Data: Parameter Estimation and Practical Identifiability Analysis

Parameter Identification for Partial Differential Equations with Spatiotemporal Varying Coefficients

Simultaneous estimation of magnetic moment and Brownian relaxation time distributions of magnetic nanoparticles based on magnetic particle spectroscopy for biosensing application

The Fokker–Planck equation for coupled Brown–Néel-rotation

Stochastic parameterization identification using ensemble Kalman filtering combined with expectation-maximization and Newton-Raphson maximum likelihood methods

A time-varying approach to single particle tracking with a nonlinear observation model

Efficient numerical approximation of a non-regular Fokker-Planck equation associated with first-passage time distributions

Trajectory dependent particle response for anisotropic mono domain particles in magnetic particle imaging

Nonequilibrium response of magnetic nanoparticles to time-varying magnetic fields: Contributions from Brownian and Néel processes

Dynamic image reconstruction with motion priors in application to 3D magnetic particle imaging

Numerical Solution of the Fokker-Planck Equation by Spectral Collocation and Finite-Element Methods for Stochastic Magnetization Dynamics

Analysis and systematic discretization of a Fokker-Planck equation with Lorentz force

Solving Inverse Stochastic Problems from Discrete Particle Observations Using the Fokker-Planck Equation and Physics-informed Neural Networks

Joint Super-Resolution Image Reconstruction and Parameter Identification in Imaging Operator: Analysis of Bilinear Operator Equations, Numerical Solution, and Application to Magnetic Particle Imaging

Dynamic Programming for Finite Ensembles of Nanomagnetic Particles

Comparative Modeling of Frequency Mixing Measurements of Magnetic Nanoparticles Using Micromagnetic Simulations and Langevin Theory