Source author record

Mauricio J. del Razo

Mauricio J. del Razo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.stat-mech math-ph math.MP math.PR physics.chem-ph Biological Physics math.NA physics.comp-ph physics.med-ph Subcellular Processes

Catalog footprint

What is connected

5works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A probabilistic framework for particle-based reaction-diffusion dynamics using classical Fock space representations

The modeling and simulation of stochastic reaction-diffusion processes is a topic of steady interest that is approached with a wide range of methods. \rev{At the level of particle-resolved descriptions, where chemical reactions are coupled to the spatial diffusion of individual particles, there exist comprehensive numerical simulation schemes, while the corresponding mathematical formalization is relatively underdeveloped. The aim of this paper is to provide a framework to systematically formulate the probabilistic evolution equation, termed chemical diffusion master equation (CDME), that governs particle-based stochastic reaction-diffusion processes. To account for the non-conserved and unbounded particle number of this type of open systems, we employ a classical analogue of the quantum mechanical Fock space that contains the symmetrized probability densities of the many-particle configurations in space. Following field-theoretical ideas of second quantization, we introduce creation and annihilation operators that act on single-particle densities and provide natural representations of symmetrized probability densities as well as of reaction and diffusion operators. These operators allow us to consistently and systematically formulate the CDME for arbitrary reaction schemes. The resulting form of the CDME further serves as the foundation to derive more coarse-grained descriptions of reaction-diffusion dynamics. In this regard, we show that a discretization of the evolution equation by projection onto a Fock subspace generated by a finite set of single-particle densities leads to a generalized form of the well-known reaction-diffusion master equation, which supports non-local reactions between grid cells and which converges properly in the continuum limit.

preprint2022arXiv

Chemical diffusion master equation: formulations of reaction--diffusion processes on the molecular level

The chemical diffusion master equation (CDME) describes the probabilistic dynamics of reaction--diffusion systems at the molecular level [del Razo et al., Lett. Math. Phys. 112:49, 2022]; it can be considered the master equation for reaction--diffusion processes. The CDME consists of an infinite ordered family of Fokker--Planck equations, where each level of the ordered family corresponds to a certain number of particles and each particle represents a molecule. The equations at each level describe the spatial diffusion of the corresponding set of particles, and they are coupled to each other via reaction operators --linear operators representing chemical reactions. These operators change the number of particles in the system, and thus transport probability between different levels in the family. In this work, we present three approaches to formulate the CDME and show the relations between them. We further deduce the non-trivial combinatorial factors contained in the reaction operators, and we elucidate the relation to the original formulation of the CDME, which is based on creation and annihilation operators acting on many-particle probability density functions. Finally we discuss applications to multiscale simulations of biochemical systems among other future prospects.

preprint2016arXiv

A Discrete Stochastic Formulation for Reversible Bimolecular Reactions via Diffusion Encounter

The classical models for irreversible diffusion-influenced reactions can be derived by introducing absorbing boundary conditions to over-damped continuous Brownian motion (BM) theory. As there is a clear corresponding stochastic process, the mathematical description takes both Kolmogorov forward equation for the evolution of the probability distribution function and the stochastic sample trajectories. This dual description is a fundamental characteristic of stochastic processes and allows simple particle based simulations to accurately match the expected statistical behavior. However, in the traditional theory using the back-reaction boundary condition to model reversible reactions with geminate recombinations, several subtleties arise: it is unclear what the underlying stochastic process is, which causes complications in producing accurate simulations; and it is non-trivial how to perform an appropriate discretization for numerical computations. In this work, we derive a discrete stochastic model that recovers the classical models and their boundary conditions in the continuous limit. In the case of reversible reactions, we recover the back-reaction boundary condition, unifying the back-reaction approach with those of current simulation packages. Furthermore, all the complications encountered in the continuous models become trivial in the discrete model. Our formulation brings to attention the question: With computations in mind, can we develop a discrete reaction kinetics model that is more fundamental than its continuous counterpart?

preprint2015arXiv

Computational and in vitro studies of blast-induced blood-brain barrier disruption

There is growing concern that blast-exposed individuals are at risk of developing neurological disorders later in life. Therefore, it is important to understand the dynamic properties of blast forces on brain cells, including the endothelial cells that maintain the blood-brain barrier (BBB), which regulates the passage of nutrients into the brain and protects it from toxins in the blood. To better understand the effect of shock waves on the BBB we have investigated an {\em in vitro} model in which BBB endothelial cells are grown in transwell vessels and exposed in a shock tube, confirming that BBB integrity is directly related to shock wave intensity. It is difficult to directly measure the forces acting on these cells in the transwell container during the experiments, and so a computational tool has been developed and presented in this paper. Two-dimensional axisymmetric Euler equations with the Tammann equation of state were used to model the transwell materials, and a high-resolution finite volume method based on Riemann solvers and the Clawpack software was used to solve these equations in a mixed Eulerian/Lagrangian frame. Results indicated that the geometry of the transwell plays a significant role in the observed pressure time series in these experiments. We also found that pressures can fall below vapor pressure due to the interaction of reflecting and diffracting shock waves, suggesting that cavitation bubbles could be a damage mechanism. Computations that include a simulated hydrophone inserted in the transwell suggest that the instrument itself could significantly alter blast wave properties. These findings illustrate the need for further computational modeling studies aimed at understanding possible blast-induced BBB damage.

preprint2014arXiv

Fluorescence Correlation Spectroscopy and Nonlinear Stochastic Reaction-Diffusion

The currently existing theory of fluorescence correlation spectroscopy(FCS) is based on the linear fluctuation theory originally developed by Einstein, Onsager, Lax, and others as a phenomenological approach to equilibrium fluctuations in bulk solutions. For mesoscopic reaction-diffusion systems with nonlinear chemical reactions among a small number of molecules, a situation often encountered in single-cell biochemistry, it is expected that FCS time correlation functions of a reaction-diffusion system can deviate from the classic results of Elson and Magde. We first discuss this nonlinear effect for reaction systems without diffusion. For nonlinear stochastic reaction-diffusion systems here are no closed solutions; therefore, stochastic Monte-Carlo simulations are carried out. We show that the deviation is small for a simple bimolecular reaction; the most significant deviations occur when the number of molecules is small and of the same order. Our results show that current linear FCS theory could be adequate for measurements on biological systems that contain many other sources of uncertainties. At the same time it provides a framework for future measurements of nonlinear, fluctuating chemical reactions with high-precision FCS. Extending Delbrück-Gillespie's theory for stochastic nonlinear reactions with rapidly stirring to reaction-diffusion systems provides a mesoscopic model for chemical and biochemical reactions at nanometric and mesoscopic level such as a single biological cell.