Flexible single molecule simulation of reaction-diffusion processes (Q544570): Difference between revisions
From MaRDI portal
Created a new Item |
ReferenceBot (talk | contribs) Changed an Item |
||
(9 intermediate revisions by 6 users not shown) | |||
Property / author | |||
Property / author: Per Loetstedt / rank | |||
Property / author | |||
Property / author: Per Loetstedt / rank | |||
Normal rank | |||
Property / review text | |||
This article describes a computational methodology for the numerical simulation of the probability density of chemically-reactive Brownian particles at a mesoscopic scale. At the macroscopic scale, chemically-reactive systems are indeed well known to solve a system of reaction rate equations established by the chemical kinetics theory. But this does not take into account the specific spatial mobility of the various species. Various microscopic models exist that keep trace of the trajectory of each molecule. The molecules are typically Brownian ``hard spheres'' that move in space by diffusion in a solvent and then react together following the standard chemical kinetics laws as soon as in close contact. But these models are often computationally expensive. In contrast, the mesoscopic models take into account the specific spatial mobility of the various species and require only the probability density functional (PDF) of the whole system, which, of course, contains less information than the whole family of molecules trajectories. The PDF given by the microscopic model is a solution to the so-called chemical master equation (CME), a difficult high-dimensional problem to solve. Here, one simplifies the problem by considering only two-body reactions. This allows one to perform numerical simulations with an event-driven algorithm like the so-called Green's function reaction dynamics (GFRD). The time-evolution is computed at discrete times on-the-fly. At each discrete time, one forecasts the next time where reactions will occur; this gives the new time step. In between, one considers only the diffusion in space of the molecules using a splitting of the system into subsystems with two or one molecules. (The diffusion problem can be solved almost exactly with analytical solutions in those cases.) The algorithm tells you how to choose the new time step and the subsystems. The authors show a good performance of their method in a number of numerical benchmarks. They also give hints on the relative importance of the errors induced by the various numerical approximations and of the modelling assumptions (comparison with macroscopic and microscopic models in cases where they agree or not). | |||
Property / review text: This article describes a computational methodology for the numerical simulation of the probability density of chemically-reactive Brownian particles at a mesoscopic scale. At the macroscopic scale, chemically-reactive systems are indeed well known to solve a system of reaction rate equations established by the chemical kinetics theory. But this does not take into account the specific spatial mobility of the various species. Various microscopic models exist that keep trace of the trajectory of each molecule. The molecules are typically Brownian ``hard spheres'' that move in space by diffusion in a solvent and then react together following the standard chemical kinetics laws as soon as in close contact. But these models are often computationally expensive. In contrast, the mesoscopic models take into account the specific spatial mobility of the various species and require only the probability density functional (PDF) of the whole system, which, of course, contains less information than the whole family of molecules trajectories. The PDF given by the microscopic model is a solution to the so-called chemical master equation (CME), a difficult high-dimensional problem to solve. Here, one simplifies the problem by considering only two-body reactions. This allows one to perform numerical simulations with an event-driven algorithm like the so-called Green's function reaction dynamics (GFRD). The time-evolution is computed at discrete times on-the-fly. At each discrete time, one forecasts the next time where reactions will occur; this gives the new time step. In between, one considers only the diffusion in space of the molecules using a splitting of the system into subsystems with two or one molecules. (The diffusion problem can be solved almost exactly with analytical solutions in those cases.) The algorithm tells you how to choose the new time step and the subsystems. The authors show a good performance of their method in a number of numerical benchmarks. They also give hints on the relative importance of the errors induced by the various numerical approximations and of the modelling assumptions (comparison with macroscopic and microscopic models in cases where they agree or not). / rank | |||
Normal rank | |||
Property / Mathematics Subject Classification ID | |||
Property / Mathematics Subject Classification ID: 82C22 / rank | |||
Normal rank | |||
Property / Mathematics Subject Classification ID | |||
Property / Mathematics Subject Classification ID: 35K57 / rank | |||
Normal rank | |||
Property / zbMATH DE Number | |||
Property / zbMATH DE Number: 5908112 / rank | |||
Normal rank | |||
Property / zbMATH Keywords | |||
reaction-diffusion | |||
Property / zbMATH Keywords: reaction-diffusion / rank | |||
Normal rank | |||
Property / zbMATH Keywords | |||
operator splitting | |||
Property / zbMATH Keywords: operator splitting / rank | |||
Normal rank | |||
Property / zbMATH Keywords | |||
chemical-master equation | |||
Property / zbMATH Keywords: chemical-master equation / rank | |||
Normal rank | |||
Property / reviewed by | |||
Property / reviewed by: Sebastien Boyabal / rank | |||
Normal rank | |||
Property / describes a project that uses | |||
Property / describes a project that uses: MesoRD / rank | |||
Normal rank | |||
Property / describes a project that uses | |||
Property / describes a project that uses: Cell++ / rank | |||
Normal rank | |||
Property / describes a project that uses | |||
Property / describes a project that uses: Smoldyn / rank | |||
Normal rank | |||
Property / MaRDI profile type | |||
Property / MaRDI profile type: MaRDI publication profile / rank | |||
Normal rank | |||
Property / full work available at URL | |||
Property / full work available at URL: https://doi.org/10.1016/j.jcp.2011.02.020 / rank | |||
Normal rank | |||
Property / OpenAlex ID | |||
Property / OpenAlex ID: W2008057778 / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Q5785591 / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: A first-passage kinetic Monte Carlo algorithm for complex diffusion-reaction systems / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Simulation of Stochastic Reaction-Diffusion Processes on Unstructured Meshes / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Stochastic reaction-diffusion kinetics in the microscopic limit / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Handbook of stochastic methods for physics, chemistry and natural sciences. / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Q2758041 / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Finite Volume Methods for Hyperbolic Problems / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Q3126807 / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: What Makes Molecular Dynamics Work? / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: On the Construction and Comparison of Difference Schemes / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Stochastic simulation of coupled reaction-diffusion processes / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Q5342712 / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Computation of the Complex Error Function / rank | |||
Normal rank | |||
links / mardi / name | links / mardi / name | ||
Latest revision as of 04:40, 4 July 2024
scientific article
Language | Label | Description | Also known as |
---|---|---|---|
English | Flexible single molecule simulation of reaction-diffusion processes |
scientific article |
Statements
Flexible single molecule simulation of reaction-diffusion processes (English)
0 references
15 June 2011
0 references
This article describes a computational methodology for the numerical simulation of the probability density of chemically-reactive Brownian particles at a mesoscopic scale. At the macroscopic scale, chemically-reactive systems are indeed well known to solve a system of reaction rate equations established by the chemical kinetics theory. But this does not take into account the specific spatial mobility of the various species. Various microscopic models exist that keep trace of the trajectory of each molecule. The molecules are typically Brownian ``hard spheres'' that move in space by diffusion in a solvent and then react together following the standard chemical kinetics laws as soon as in close contact. But these models are often computationally expensive. In contrast, the mesoscopic models take into account the specific spatial mobility of the various species and require only the probability density functional (PDF) of the whole system, which, of course, contains less information than the whole family of molecules trajectories. The PDF given by the microscopic model is a solution to the so-called chemical master equation (CME), a difficult high-dimensional problem to solve. Here, one simplifies the problem by considering only two-body reactions. This allows one to perform numerical simulations with an event-driven algorithm like the so-called Green's function reaction dynamics (GFRD). The time-evolution is computed at discrete times on-the-fly. At each discrete time, one forecasts the next time where reactions will occur; this gives the new time step. In between, one considers only the diffusion in space of the molecules using a splitting of the system into subsystems with two or one molecules. (The diffusion problem can be solved almost exactly with analytical solutions in those cases.) The algorithm tells you how to choose the new time step and the subsystems. The authors show a good performance of their method in a number of numerical benchmarks. They also give hints on the relative importance of the errors induced by the various numerical approximations and of the modelling assumptions (comparison with macroscopic and microscopic models in cases where they agree or not).
0 references
reaction-diffusion
0 references
operator splitting
0 references
chemical-master equation
0 references