Browse
Submit Data
Databases
API
Help

Dataset Information

15 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Infinite product expansion of the Fokker-Planck equation with steady-state solution.

ABSTRACT: We present an analytical technique for solving Fokker-Planck equations that have a steady-state solution by representing the solution as an infinite product rather than, as usual, an infinite sum. This method has many advantages: automatically ensuring positivity of the resulting approximation, and by design exactly matching both the short- and long-term behaviour. The efficacy of the technique is demonstrated via comparisons with computations of typical examples.

SUBMITTER: Martin RJ

PROVIDER: S-EPMC4528656 | biostudies-literature | 2015 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Infinite product expansion of the Fokker-Planck equation with steady-state solution.

Martin R J RJ Craster R V RV Kearney M J MJ

Proceedings. Mathematical, physical, and engineering sciences 20150701 2179

We present an analytical technique for solving Fokker-Planck equations that have a steady-state solution by representing the solution as an infinite product rather than, as usual, an infinite sum. This method has many advantages: automatically ensuring positivity of the resulting approximation, and by design exactly matching both the short- and long-term behaviour. The efficacy of the technique is demonstrated via comparisons with computations of typical examples. ...[more]

PMID: 26346100

Similar Datasets

Applying the Fokker-Planck equation to grating-based x-ray phase and dark-field imaging.

Project description:X-ray imaging has conventionally relied upon attenuation to provide contrast. In recent years, two complementary modalities have been added; (a) phase contrast, which can capture low-density samples that are difficult to see using attenuation, and (b) dark-field x-ray imaging, which reveals the presence of sub-pixel sample structures. These three modalities can be accessed using a crystal analyser, a grating interferometer or by looking at a directly-resolved grid, grating or speckle pattern. Grating and grid-based methods extract a differential phase signal by measuring how far a feature in the illumination has been shifted transversely due to the presence of a sample. The dark-field signal is extracted by measuring how the visibility of the structured illumination is decreased, typically due to the presence of sub-pixel structures in a sample. The strength of the dark-field signal may depend on the grating period, the pixel size and the set-up distances, and additional dark-field signal contributions may be seen as a result of strong phase effects or other factors. In this paper we show that the finite-difference form of the Fokker-Planck equation can be applied to describe the drift (phase signal) and diffusion (dark-field signal) of the periodic or structured illumination used in phase contrast x-ray imaging with gratings, in order to better understand any cross-talk between attenuation, phase and dark-field x-ray signals. In future work, this mathematical description could be used as a basis for new approaches to the inverse problem of recovering both phase and dark-field information.

| S-EPMC6877582 | biostudies-literature

Beating the curse of dimension with accurate statistics for the Fokker-Planck equation in complex turbulent systems.

Project description:Solving the Fokker-Planck equation for high-dimensional complex dynamical systems is an important issue. Recently, the authors developed efficient statistically accurate algorithms for solving the Fokker-Planck equations associated with high-dimensional nonlinear turbulent dynamical systems with conditional Gaussian structures, which contain many strong non-Gaussian features such as intermittency and fat-tailed probability density functions (PDFs). The algorithms involve a hybrid strategy with a small number of samples [Formula: see text], where a conditional Gaussian mixture in a high-dimensional subspace via an extremely efficient parametric method is combined with a judicious Gaussian kernel density estimation in the remaining low-dimensional subspace. In this article, two effective strategies are developed and incorporated into these algorithms. The first strategy involves a judicious block decomposition of the conditional covariance matrix such that the evolutions of different blocks have no interactions, which allows an extremely efficient parallel computation due to the small size of each individual block. The second strategy exploits statistical symmetry for a further reduction of [Formula: see text] The resulting algorithms can efficiently solve the Fokker-Planck equation with strongly non-Gaussian PDFs in much higher dimensions even with orders in the millions and thus beat the curse of dimension. The algorithms are applied to a [Formula: see text]-dimensional stochastic coupled FitzHugh-Nagumo model for excitable media. An accurate recovery of both the transient and equilibrium non-Gaussian PDFs requires only [Formula: see text] samples! In addition, the block decomposition facilitates the algorithms to efficiently capture the distinct non-Gaussian features at different locations in a [Formula: see text]-dimensional two-layer inhomogeneous Lorenz 96 model, using only [Formula: see text] samples.

| S-EPMC5724285 | biostudies-literature

Solving the inverse problem of time independent Fokker-Planck equation with a self supervised neural network method.

Project description:The Fokker-Planck equation (FPE) has been used in many important applications to study stochastic processes with the evolution of the probability density function (pdf). Previous studies on FPE mainly focus on solving the forward problem which is to predict the time-evolution of the pdf from the underlying FPE terms. However, in many applications the FPE terms are usually unknown and roughly estimated, and solving the forward problem becomes more challenging. In this work, we take a different approach of starting with the observed pdfs to recover the FPE terms using a self-supervised machine learning method. This approach, known as the inverse problem, has the advantage of requiring minimal assumptions on the FPE terms and allows data-driven scientific discovery of unknown FPE mechanisms. Specifically, we propose an FPE-based neural network (FPE-NN) which directly incorporates the FPE terms as neural network weights. By training the network on observed pdfs, we recover the FPE terms. Additionally, to account for noise in real-world observations, FPE-NN is able to denoise the observed pdfs by training the pdfs alongside the network weights. Our experimental results on various forms of FPE show that FPE-NN can accurately recover FPE terms and denoising the pdf plays an essential role.

| S-EPMC8324819 | biostudies-literature

Pedestrian motion modelled by Fokker-Planck Nash games.

Project description:A new approach to modelling pedestrians' avoidance dynamics based on a Fokker-Planck (FP) Nash game framework is presented. In this framework, two interacting pedestrians are considered, whose motion variability is modelled through the corresponding probability density functions (PDFs) governed by FP equations. Based on these equations, a Nash differential game is formulated where the game strategies represent controls aiming at avoidance by minimizing appropriate collision cost functionals. The existence of Nash equilibria solutions is proved and characterized as a solution to an optimal control problem that is solved numerically. Results of numerical experiments are presented that successfully compare the computed Nash equilibria to the output of real experiments (conducted with humans) for four test cases.

| S-EPMC5627107 | biostudies-literature

Fokker-Planck analysis of optical near-field traps.

Project description:The motion of a nanoparticle in the vicinity of a near-field optical trap is modeled using the Fokker-Planck equation. A plasmonic C-shaped engraving on a gold film is considered as the optical trap. The time evolution of the position probability density of the nanoparticle is calculated to analyze the trapping dynamics. A spatially varying diffusion tensor is used in the formulation to take into account the hydrodynamic interactions. The steady-state position distribution obtained from the Fokker-Planck equation is compared with experimental results and found to be in good agreement. Computational cost of the proposed method is compared with the conventionally used Langevin equation based approach. The proposed method is found to be computationally efficient (requiring 35 times less computation time) and scalable to more complex lab-on-a-chip systems.

| S-EPMC6606609 | biostudies-literature

Finite element solution of the steady-state Smoluchowski equation for rate constant calculations.

Project description:This article describes the development and implementation of algorithms to study diffusion in biomolecular systems using continuum mechanics equations. Specifically, finite element methods have been developed to solve the steady-state Smoluchowski equation to calculate ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to mouse acetylcholinesterase. Rates for inhibitor binding to mAChE were calculated at various ionic strengths with several different reaction criteria. The calculated rates were compared with experimental data and show very good agreement when the correct reaction criterion is used. Additionally, these finite element methods require significantly less computational resources than existing particle-based Brownian dynamics methods.

| S-EPMC1304055 | biostudies-literature

Fokker-Planck dynamics of the El Nino-Southern Oscillation.

Project description:The asymmetric nature of the El Niño-Southern Oscillation (ENSO) is explored by using a probabilistic model (PROM) for ENSO. Based on a Fokker-Planck Equation (FPE), PROM describes the dynamics of a nonlinear stochastic ENSO recharge oscillator model for eastern equatorial Pacific temperature anomalies and equatorial Pacific basin-averaged thermocline depth changes. Eigen analyses of PROM provide new insights into the stationary and oscillatory solutions of the stochastic dynamical system. The first probabilistic eigenmode represents a stationary mode, which exhibits the asymmetric features of ENSO, in case deterministic nonlinearities or multiplicative noises are included. The second mode is linked to the oscillatory nature of ENSO and represents a cyclic asymmetric probability distribution, which emerges from the key dynamical processes. Other eigenmodes are associated with the temporal evolution of higher order statistical moments of the ENSO system. The model solutions demonstrate that the deterministic nonlinearity plays a stronger role in establishing the observed asymmetry of ENSO as compared to the multiplicative stochastic part.

| S-EPMC7529818 | biostudies-literature

A Fokker-Planck reaction model for the epitaxial growth and shape transition of quantum dots.

Project description:We construct a Fokker-Planck reaction (FPR) model to investigate the dynamics of the coupled epitaxial growth and shape transition process of an array of quantum dots (QDs). The FPR model is based on a coupled system of Fokker-Planck equations wherein the distribution of each island type is governed by its own Fokker-Planck equation for growth, with reaction terms describing the shape transitions between islands of different types including asymmetric shapes. The reaction terms for the shape transitions depend on the island size and are determined from explicit calculations of the lowest barrier pathway for each shape transition. This mean-field model enables us to consider the kinetics of asymmetric shape transitions and study the evolution of island shape distributions during the coupled growth and transition process. Asymmetric metastable shapes play a crucial role in the dynamics, with asymmetric QDs comprising up to 10% of the population, and with up to 100% of the shape transitions passing through asymmetric shapes. Moreover, we find that the characteristic multimodal distribution of pyramid/dome QD coarsening can be eliminated at sufficiently high temperature and deposition rate.

| S-EPMC5666231 | biostudies-literature

Exact and efficient solution of the radiative transport equation for the semi-infinite medium.

Project description:An accurate and efficient solution of the radiative transport equation is proposed for modeling the propagation of photons in the three-dimensional anisotropically scattering half-space medium. The exact refractive index mismatched boundary condition is considered and arbitrary rotationally invariant scattering functions can be applied. The obtained equations are verified with Monte Carlo simulations in the steady-state, temporal frequency, and time domains resulting in an excellent agreement.

| S-EPMC3684809 | biostudies-other

Differential equation methods for simulation of GFP kinetics in non-steady state experiments.

Project description:Genetically encoded fluorescent proteins, combined with fluorescence microscopy, are widely used in cell biology to collect kinetic data on intracellular trafficking. Methods for extraction of quantitative information from these data are based on the mathematics of diffusion and tracer kinetics. Current methods, although useful and powerful, depend on the assumption that the cellular system being studied is in a steady state, that is, the assumption that all the molecular concentrations and fluxes are constant for the duration of the experiment. Here, we derive new tracer kinetic analytical methods for non-steady state biological systems by constructing mechanistic nonlinear differential equation models of the underlying cell biological processes and linking them to a separate set of differential equations governing the kinetics of the fluorescent tracer. Linking the two sets of equations is based on a new application of the fundamental tracer principle of indistinguishability and, unlike current methods, supports correct dependence of tracer kinetics on cellular dynamics. This approach thus provides a general mathematical framework for applications of GFP fluorescence microscopy (including photobleaching [FRAP, FLIP] and photoactivation to frequently encountered experimental protocols involving physiological or pharmacological perturbations (e.g., growth factors, neurotransmitters, acute knockouts, inhibitors, hormones, cytokines, and metabolites) that initiate mechanistically informative intracellular transients. When a new steady state is achieved, these methods automatically reduce to classical steady state tracer kinetic analysis.

| S-EPMC6003217 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data