Browse
Submit Data
Databases
API
Help

Dataset Information

15 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Minimum wave speeds in monostable reaction-diffusion equations: sharp bounds by polynomial optimization.

ABSTRACT: Many monostable reaction-diffusion equations admit one-dimensional travelling waves if and only if the wave speed is sufficiently high. The values of these minimum wave speeds are not known exactly, except in a few simple cases. We present methods for finding upper and lower bounds on minimum wave speed. They rely on constructing trapping boundaries for dynamical systems whose heteroclinic connections correspond to the travelling waves. Simple versions of this approach can be carried out analytically but often give overly conservative bounds on minimum wave speed. When the reaction-diffusion equations being studied have polynomial nonlinearities, our approach can be implemented computationally using polynomial optimization. For scalar reaction-diffusion equations, we present a general method and then apply it to examples from the literature where minimum wave speeds were unknown. The extension of our approach to multi-component reaction-diffusion systems is then illustrated using a cubic autocatalysis model from the literature. In all three examples and with many different parameter values, polynomial optimization computations give upper and lower bounds that are within 0.1% of each other and thus nearly sharp. Upper bounds are derived analytically as well for the scalar reaction-diffusion equations.

SUBMITTER: Bramburger JJ

PROVIDER: S-EPMC7544360 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Correction: sharp bounds and normalization of Wiener-type indices.

Project description:[This corrects the article DOI: 10.1371/journal.pone.0078448.].

| S-EPMC3988199 | biostudies-literature

Sharp bounds and normalization of Wiener-type indices.

Project description:Complex networks abound in physical, biological and social sciences. Quantifying a network's topological structure facilitates network exploration and analysis, and network comparison, clustering and classification. A number of Wiener type indices have recently been incorporated as distance-based descriptors of complex networks, such as the R package QuACN. Wiener type indices are known to depend both on the network's number of nodes and topology. To apply these indices to measure similarity of networks of different numbers of nodes, normalization of these indices is needed to correct the effect of the number of nodes in a network. This paper aims to fill this gap. Moreover, we introduce an f-Wiener index of network G, denoted by Wf(G). This notion generalizes the Wiener index to a very wide class of Wiener type indices including all known Wiener type indices. We identify the maximum and minimum of Wf(G) over a set of networks with n nodes. We then introduce our normalized-version of f-Wiener index. The normalized f-Wiener indices were demonstrated, in a number of experiments, to improve significantly the hierarchical clustering over the non-normalized counterparts.

| S-EPMC3832646 | biostudies-other

Poisson-Nernst-Planck equations for simulating biomolecular diffusion-reaction processes II: size effects on ionic distributions and diffusion-reaction rates.

Project description:The effects of finite particle size on electrostatics, density profiles, and diffusion have been a long existing topic in the study of ionic solution. The previous size-modified Poisson-Boltzmann and Poisson-Nernst-Planck models are revisited in this article. In contrast to many previous works that can only treat particle species with a single uniform size or two sizes, we generalize the Borukhov model to obtain a size-modified Poisson-Nernst-Planck (SMPNP) model that is able to treat nonuniform particle sizes. The numerical tractability of the model is demonstrated as well. The main contributions of this study are as follows. 1), We show that an (arbitrarily) size-modified PB model is indeed implied by the SMPNP equations under certain boundary/interface conditions, and can be reproduced through numerical solutions of the SMPNP. 2), The size effects in the SMPNP effectively reduce the densities of highly concentrated counterions around the biomolecule. 3), The SMPNP is applied to the diffusion-reaction process for the first time, to our knowledge. In the case of low substrate density near the enzyme reactive site, it is observed that the rate coefficients predicted by SMPNP model are considerably larger than those by the PNP model, suggesting both ions and substrates are subject to finite size effects. 4), An accurate finite element method and a convergent Gummel iteration are developed for the numerical solution of the completely coupled nonlinear system of SMPNP equations.

| S-EPMC3093557 | biostudies-literature

Sharp sensitivity bounds for mediation under unmeasured mediator-outcome confounding.

Project description:It is often of interest to decompose the total effect of an exposure into a component that acts on the outcome through some mediator and a component that acts independently through other pathways. Said another way, we are interested in the direct and indirect effects of the exposure on the outcome. Even if the exposure is randomly assigned, it is often infeasible to randomize the mediator, leaving the mediator-outcome confounding not fully controlled. We develop a sensitivity analysis technique that can bound the direct and indirect effects without parametric assumptions about the unmeasured mediator-outcome confounding.

| S-EPMC4890130 | biostudies-literature

Sharp bounds on sufficient-cause interactions under the assumption of no redundancy.

Project description:BackgroundSufficient-cause interaction is a type of interaction that has received much attention recently. The sufficient component cause model on which the sufficient-cause interaction is based is however a non-identifiable model. Estimating the interaction parameters from the model is mathematically impossible.MethodsIn this paper, I derive bounding formulae for sufficient-cause interactions under the assumption of no redundancy.ResultsTwo real data sets are used to demonstrate the method (R codes provided). The proposed bounds are sharp and sharper than previous bounds.ConclusionsSufficient-cause interactions can be quantified by setting bounds on them.

| S-EPMC5399844 | biostudies-literature

Exact Solutions of Coupled Multispecies Linear Reaction-Diffusion Equations on a Uniformly Growing Domain.

Project description:Embryonic development involves diffusion and proliferation of cells, as well as diffusion and reaction of molecules, within growing tissues. Mathematical models of these processes often involve reaction-diffusion equations on growing domains that have been primarily studied using approximate numerical solutions. Recently, we have shown how to obtain an exact solution to a single, uncoupled, linear reaction-diffusion equation on a growing domain, 0 < x < L(t), where L(t) is the domain length. The present work is an extension of our previous study, and we illustrate how to solve a system of coupled reaction-diffusion equations on a growing domain. This system of equations can be used to study the spatial and temporal distributions of different generations of cells within a population that diffuses and proliferates within a growing tissue. The exact solution is obtained by applying an uncoupling transformation, and the uncoupled equations are solved separately before applying the inverse uncoupling transformation to give the coupled solution. We present several example calculations to illustrate different types of behaviour. The first example calculation corresponds to a situation where the initially-confined population diffuses sufficiently slowly that it is unable to reach the moving boundary at x = L(t). In contrast, the second example calculation corresponds to a situation where the initially-confined population is able to overcome the domain growth and reach the moving boundary at x = L(t). In its basic format, the uncoupling transformation at first appears to be restricted to deal only with the case where each generation of cells has a distinct proliferation rate. However, we also demonstrate how the uncoupling transformation can be used when each generation has the same proliferation rate by evaluating the exact solutions as an appropriate limit.

| S-EPMC4583548 | biostudies-literature

Complete analytic solutions for convection-diffusion-reaction-source equations without using an inverse Laplace transform.

Project description:Transient mass-transfer phenomena occurring in natural and engineered systems consist of convection, diffusion, and reaction processes. The coupled phenomena can be described by using the unsteady convection-diffusion-reaction (CDR) equation, which is classified in mathematics as a linear, parabolic partial-differential equation. The availability of analytic solutions is limited to simple cases, e.g., unsteady diffusion and steady convective diffusion. The CDR equation has been considered analytically intractable, depending on the initial and boundary conditions. If spatial adsorption and desorption of matter are super-positioned in the CDR equation as sink and source functions, respectively, then the governing equation becomes an unsteady convection-diffusion-reaction-source (CDRS) equation, of which general solutions are unknown. In this study, a general 1D analytic solution of the CDRS equation is obtained by using a one-sided Laplace transform, by assuming constant diffusivity, velocity, and reactivity. This paper also provides a general formalism to derive 1D analytic solutions for Dirichlet/Dirichlet and Dirichlet/Neumann boundary conditions. Derivations of the analytic solutions are found to be straightforward if a combination of the source function and the initial concentration provide a non-zero singularity pole of inverse Laplace transform.

| S-EPMC7228951 | biostudies-literature

New mechanism of spiral wave initiation in a reaction-diffusion-mechanics system.

Project description:Spiral wave initiation in the heart muscle is a mechanism for the onset of dangerous cardiac arrhythmias. A standard protocol for spiral wave initiation is the application of a stimulus in the refractory tail of a propagating excitation wave, a region that we call the "classical vulnerable zone." Previous studies of vulnerability to spiral wave initiation did not take the influence of deformation into account, which has been shown to have a substantial effect on the excitation process of cardiomyocytes via the mechano-electrical feedback phenomenon. In this work we study the effect of deformation on the vulnerability of excitable media in a discrete reaction-diffusion-mechanics (dRDM) model. The dRDM model combines FitzHugh-Nagumo type equations for cardiac excitation with a discrete mechanical description of a finite-elastic isotropic material (Seth material) to model cardiac excitation-contraction coupling and stretch activated depolarizing current. We show that deformation alters the "classical," and forms a new vulnerable zone at longer coupling intervals. This mechanically caused vulnerable zone results in a new mechanism of spiral wave initiation, where unidirectional conduction block and rotation directions of the consequently initiated spiral waves are opposite compared to the mechanism of spiral wave initiation due to the "classical vulnerable zone." We show that this new mechanism of spiral wave initiation can naturally occur in situations that involve wave fronts with curvature, and discuss its relation to supernormal excitability of cardiac tissue. The concept of mechanically induced vulnerability may lead to a better understanding about the onset of dangerous heart arrhythmias via mechano-electrical feedback.

| S-EPMC3215707 | biostudies-literature

Wave-pinning and cell polarity from a bistable reaction-diffusion system.

Project description:Motile eukaryotic cells polarize in response to external signals. Numerous mechanisms have been suggested to account for this symmetry breaking and for the ensuing robust polarization. Implicated in this process are various proteins that are recruited to the plasma membrane and segregate at an emergent front or back of the polarizing cell. Among these are PI3K, PTEN, and members of the Rho family GTPases such as Cdc42, Rac, and Rho. Many such proteins, including the Rho GTPases, cycle between active membrane-bound forms and inactive cytosolic forms. In previous work, we have shown that this property, together with appropriate crosstalk, endows a biochemical circuit (Cdc42, Rac, and Rho) with the property of inherent polarizability. Here we show that this property is present in an even simpler system comprised of a single active/inactive protein pair with positive feedback to its own activation. The simplicity of this minimal system also allows us to explain the mechanism using insights from mathematical analysis. The basic idea resides in a well-known property of reaction-diffusion systems with bistable kinetics, namely, propagation of fronts. However, it crucially depends on exchange between active and inactive forms of the chemicals with unequal rates of diffusion, and overall conservation to pin the waves into a stable polar distribution. We refer to these dynamics as wave-pinning and we show that this phenomenon is distinct from Turing-instability-generated pattern formation that occurs in reaction-diffusion systems that appear to be very similar. We explain the mathematical basis of the phenomenon, relate it to spatial segregation of Rho GTPases, and show how it can account for spatial amplification and maintenance of polarity, as well as sensitivity to new stimuli typical in polarization of eukaryotic cells.

| S-EPMC2292363 | biostudies-other

Awake hippocampal sharp-wave ripples support spatial memory.

Project description:The hippocampus is critical for spatial learning and memory. Hippocampal neurons in awake animals exhibit place field activity that encodes current location, as well as sharp-wave ripple (SWR) activity during which representations based on past experiences are often replayed. The relationship between these patterns of activity and the memory functions of the hippocampus is poorly understood. We interrupted awake SWRs in animals learning a spatial alternation task. We observed a specific learning and performance deficit that persisted throughout training. This deficit was associated with awake SWR activity, as SWR interruption left place field activity and post-experience SWR reactivation intact. These results provide a link between awake SWRs and hippocampal memory processes, which suggests that awake replay of memory-related information during SWRs supports learning and memory-guided decision-making.

| S-EPMC4441285 | 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