7,036 research outputs found
Stability Estimates and Structural Spectral Properties of Saddle Point Problems
For a general class of saddle point problems sharp estimates for
Babu\v{s}ka's inf-sup stability constants are derived in terms of the constants
in Brezzi's theory. In the finite-dimensional Hermitian case more detailed
spectral properties of preconditioned saddle point matrices are presented,
which are helpful for the convergence analysis of common Krylov subspace
methods. The theoretical results are applied to two model problems from optimal
control with time-periodic state equations. Numerical experiments with the
preconditioned minimal residual method are reported
Shaping bursting by electrical coupling and noise
Gap-junctional coupling is an important way of communication between neurons
and other excitable cells. Strong electrical coupling synchronizes activity
across cell ensembles. Surprisingly, in the presence of noise synchronous
oscillations generated by an electrically coupled network may differ
qualitatively from the oscillations produced by uncoupled individual cells
forming the network. A prominent example of such behavior is the synchronized
bursting in islets of Langerhans formed by pancreatic \beta-cells, which in
isolation are known to exhibit irregular spiking. At the heart of this
intriguing phenomenon lies denoising, a remarkable ability of electrical
coupling to diminish the effects of noise acting on individual cells.
In this paper, we derive quantitative estimates characterizing denoising in
electrically coupled networks of conductance-based models of square wave
bursting cells. Our analysis reveals the interplay of the intrinsic properties
of the individual cells and network topology and their respective contributions
to this important effect. In particular, we show that networks on graphs with
large algebraic connectivity or small total effective resistance are better
equipped for implementing denoising. As a by-product of the analysis of
denoising, we analytically estimate the rate with which trajectories converge
to the synchronization subspace and the stability of the latter to random
perturbations. These estimates reveal the role of the network topology in
synchronization. The analysis is complemented by numerical simulations of
electrically coupled conductance-based networks. Taken together, these results
explain the mechanisms underlying synchronization and denoising in an important
class of biological models
The geometry of spontaneous spiking in neuronal networks
The mathematical theory of pattern formation in electrically coupled networks
of excitable neurons forced by small noise is presented in this work. Using the
Freidlin-Wentzell large deviation theory for randomly perturbed dynamical
systems and the elements of the algebraic graph theory, we identify and analyze
the main regimes in the network dynamics in terms of the key control
parameters: excitability, coupling strength, and network topology. The analysis
reveals the geometry of spontaneous dynamics in electrically coupled network.
Specifically, we show that the location of the minima of a certain continuous
function on the surface of the unit n-cube encodes the most likely activity
patterns generated by the network. By studying how the minima of this function
evolve under the variation of the coupling strength, we describe the principal
transformations in the network dynamics. The minimization problem is also used
for the quantitative description of the main dynamical regimes and transitions
between them. In particular, for the weak and strong coupling regimes, we
present asymptotic formulas for the network activity rate as a function of the
coupling strength and the degree of the network. The variational analysis is
complemented by the stability analysis of the synchronous state in the strong
coupling regime. The stability estimates reveal the contribution of the network
connectivity and the properties of the cycle subspace associated with the graph
of the network to its synchronization properties. This work is motivated by the
experimental and modeling studies of the ensemble of neurons in the Locus
Coeruleus, a nucleus in the brainstem involved in the regulation of cognitive
performance and behavior
A Bramble-Pasciak conjugate gradient method for discrete Stokes equations with random viscosity
We study the iterative solution of linear systems of equations arising from
stochastic Galerkin finite element discretizations of saddle point problems. We
focus on the Stokes model with random data parametrized by uniformly
distributed random variables and discuss well-posedness of the variational
formulations. We introduce a Bramble-Pasciak conjugate gradient method as a
linear solver. It builds on a non-standard inner product associated with a
block triangular preconditioner. The block triangular structure enables more
sophisticated preconditioners than the block diagonal structure usually applied
in MINRES methods. We show how the existence requirements of a conjugate
gradient method can be met in our setting. We analyze the performance of the
solvers depending on relevant physical and numerical parameters by means of
eigenvalue estimates. For this purpose, we derive bounds for the eigenvalues of
the relevant preconditioned sub-matrices. We illustrate our findings using the
flow in a driven cavity as a numerical test case, where the viscosity is given
by a truncated Karhunen-Lo\`eve expansion of a random field. In this example, a
Bramble-Pasciak conjugate gradient method with block triangular preconditioner
outperforms a MINRES method with block diagonal preconditioner in terms of
iteration numbers.Comment: 19 pages, 1 figure, submitted to SIAM JU
Gell-Mann - Low Function in QED for the arbitrary coupling constant
The Gell-Mann -- Low function \beta(g) in QED (g is the fine structure
constant) is reconstructed. At large g, it behaves as \beta_\infty g^\alpha
with \alpha\approx 1, \beta_\infty\approx 1.Comment: 5 pages, PD
- …