Stability of gradient semigroups under perturbations

In this paper we prove that gradient-like semigroups (in the sense of Carvalho and Langa (2009 J. Diff. Eqns 246 2646-68)) are gradient semigroups (possess a Lyapunov function). This is primarily done to provide conditions under which gradient semigroups, in a general metric space, are stable under perturbation exploiting the known fact (see Carvalho and Langa (2009 J. Diff. Eqns 246 2646-68)) that gradient-like semigroups are stable under perturbation. The results presented here were motivated by the work carried out in Conley (1978 Isolated Invariant Sets and the Morse Index (CBMS Regional Conference Series in Mathematics vol 38) (RI: American Mathematical Society Providence)) for groups in compact metric spaces (see also Rybakowski (1987 The Homotopy Index and Partial Differential Equations (Universitext) (Berlin: Springer)) for the Morse decomposition of an invariant set for a semigroup on a compact metric space).

Strongly damped wave problems: Bootstrapping and regularity of solutions

The aim of the article is to present a unified approach to the existence, uniqueness and regularity of solutions to problems belonging to a class of second order in time semilinear partial differential equations in Banach spaces. Our results are applied next to a number of examples appearing in literature, which fall into the class of strongly damped semilinear wave equations. The present work essentially extends the results on the existence and regularity of solutions to such problems. Previously, these problems have been considered mostly within the Hilbert space setting and with the main part operators being selfadjoint. In this article we present a more general approach, involving sectorial operators in reflexive Banach spaces. (C) 2008 Elsevier Inc. All rights reserved.

Global analytic regularity for structures of co-rank one

We consider real analytic involutive structures V, of co-rank one, defined on a real analytic paracompact orientable manifold M. To each such structure we associate certain connected subsets of M which we call the level sets of V. We prove that analytic regularity propagates along them. With a further assumption on the level sets of V we characterize the global analytic hypoellipticity of a differential operator naturally associated to V. As an application we study a case of tube structures.

Existence of secondary bifurcations or isolas for PDEs

In this paper, we introduce a method to conclude about the existence of secondary bifurcations or isolas of steady state solutions for parameter dependent nonlinear partial differential equations. The technique combines the Global Bifurcation Theorem, knowledge about the non-existence of nontrivial steady state solutions at the zero parameter value and explicit information about the coexistence of multiple nontrivial steady states at a positive parameter value. We apply the method to the two-dimensional Swift-Hohenberg equation. (C) 2011 Elsevier Ltd. All rights reserved.

A discontinuous-Galerkin-based immersed boundary method

A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user-defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous-Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements. boundary locking is avoided and optimal-order convergence is achieved. This is shown through numerical experiments in reaction-diffusion problems. Copyright (c) 2008 John Wiley & Sons, Ltd.

Limit sets and the Poincare-Bendixson Theorem in impulsive semidynamical systems

We consider semidynamical systems with impulse effects at variable times and we discuss some properties of the limit sets of orbits of these systems such as invariancy, compactness and connectedness. As a consequence we obtain a version of the Poincare-Bendixson Theorem for impulsive semidynamical systems. (C) 2008 Elsevier Inc. All rights reserved.

Design of associative memories using cellular neural networks

Cellular neural networks (CNNs) have locally connected neurons. This characteristic makes CNNs adequate for hardware implementation and, consequently, for their employment on a variety of applications as real-time image processing and construction of efficient associative memories. Adjustments of CNN parameters is a complex problem involved in the configuration of CNN for associative memories. This paper reviews methods of associative memory design based on CNNs, and provides comparative performance analysis of these approaches.

Sound waves and solitons in hot and dense nuclear matter

Assuming that nuclear matter can be treated as a perfect fluid, we study the propagation of perturbations in the baryon density. The equation of state is derived from a relativistic mean field model, which is a variant of the non-linear Walecka model. The expansion of the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations leads to differential equations for the density perturbation. We solve them numerically for linear and spherical perturbations and follow the propagation of the initial pulses. For linear perturbations we find single soliton solutions and solutions with one or more solitons followed by ""radiation"". Depending on the equation of state a strong damping may occur. We consider also the evolution of perturbations in a medium without dispersive effects. In this case we observe the formation and breaking of shock waves. We study all these equations also for matter at finite temperature. Our results may be relevant for the analysis of RHIC data. They suggest that the shock waves formed in the quark gluon plasma phase may survive and propagate in the hadronic phase. (C) 2009 Elseiver. B.V. All rights reserved.

On the concept of electrical impedance for an electrolytic cell

The analysis of the electrical impedance of an electrolytic cell in the shape of a slab is performed. We have solved, numerically, the differential equations governing the phenomenon of the redistribution of the ions in the presence of an external electric field, and compared the results with the ones obtained by solving the linear approximation of these equations. The control parameters in our study are the amplitude and the frequency of the applied voltage, assumed a simple harmonic function of the time. We show that for the large amplitudes of the applied voltage, the actual current is no longer harmonic at low frequencies. From this result it follows that the concept of electrical impedance of a cell is a useful quantity only in the case where the linear approximation of the fundamental equations of problem work well.

Opinion dynamics of learning agents: does seeking consensus lead to disagreement?

We study opinion dynamics in a population of interacting adaptive agents voting on a set of issues represented by vectors. We consider agents who can classify issues into one of two categories and can arrive at their opinions using an adaptive algorithm. Adaptation comes from learning and the information for the learning process comes from interacting with other neighboring agents and trying to change the internal state in order to concur with their opinions. The change in the internal state is driven by the information contained in the issue and in the opinion of the other agent. We present results in a simple yet rich context where each agent uses a Boolean perceptron to state their opinion. If the update occurs with information asynchronously exchanged among pairs of agents, then the typical case, if the number of issues is kept small, is the evolution into a society torn by the emergence of factions with extreme opposite beliefs. This occurs even when seeking consensus with agents with opposite opinions. If the number of issues is large, the dynamics becomes trapped, the society does not evolve into factions and a distribution of moderate opinions is observed. The synchronous case is technically simpler and is studied by formulating the problem in terms of differential equations that describe the evolution of order parameters that measure the consensus between pairs of agents. We show that for a large number of issues and unidirectional information flow, global consensus is a fixed point; however, the approach to this consensus is glassy for large societies.

Canonical quantum potential scattering theory

A new formulation of potential scattering in quantum mechanics is developed using a close structural analogy between partial waves and the classical dynamics of many non-interacting fields. Using a canonical formalism we find nonlinear first-order differential equations for the low-energy scattering parameters such as scattering length and effective range. They significantly simplify typical calculations, as we illustrate for atom-atom and neutron-nucleus scattering systems. A generalization to charged particle scattering is also possible.

Self-dual hopfions

We construct static and time-dependent exact soliton solutions with nontrivial Hopf topological charge for a field theory in 3 + 1 dimensions with the target space being the two dimensional sphere S(2). The model considered is a reduction of the so-called extended Skyrme-Faddeev theory by the removal of the quadratic term in derivatives of the fields. The solutions are constructed using an ansatz based on the conformal and target space symmetries. The solutions are said self-dual because they solve first order differential equations which together with some conditions on the coupling constants, imply the second order equations of motion. The solutions belong to a sub-sector of the theory with an infinite number of local conserved currents. The equation for the profile function of the ansatz corresponds to the Bogomolny equation for the sine-Gordon model.

Static Hopfions in the extended Skyrme-Faddeev model

We construct static soliton solutions with non-zero Hopf topological charges to a theory which is an extension of the Skyrme-Faddeev model by the addition of a further quartic term in derivatives. We use an axially symmetric ansatz based on toroidal coordinates, and solve the resulting two coupled non-linear partial differential equations in two variables by a successive over-relaxation (SOR) method. We construct numerical solutions with Hopf charge up to four, and calculate their analytical behavior in some limiting cases. The solutions present an interesting behavior under the changes of a special combination of the coupling constants of the quartic terms. Their energies and sizes tend to zero as that combination approaches a particular special value. We calculate the equivalent of the Vakulenko and Kapitanskii energy bound for the theory and find that it vanishes at that same special value of the coupling constants. In addition, the model presents an integrable sector with an in finite number of local conserved currents which apparently are not related to symmetries of the action. In the intersection of those two special sectors the theory possesses exact vortex solutions (static and time dependent) which were constructed in a previous paper by one of the authors. It is believed that such model describes some aspects of the low energy limit of the pure SU(2) Yang-Mills theory, and our results may be important in identifying important structures in that strong coupling regime.

Axially symmetric soliton solutions in a Skyrme-Faddeev-type model with Gies`s extension

We construct static soliton solutions with non-zero Hopf topological charges to a theory which is the extended Skyrme-Faddeev model with a further quartic term in derivatives. We use an axially symmetric ansatz based on toroidal coordinates, and solve the resulting two coupled nonlinear partial differential equations in two variables by a successive over-relaxation method. We construct numerical solutions with the Hopf charge up to 4. The solutions present an interesting behavior under the changes of a special combination of the coupling constants of the quartic terms.

A multi-GPU algorithm for large-scale neuronal networks

Large-scale simulations of parts of the brain using detailed neuronal models to improve our understanding of brain functions are becoming a reality with the usage of supercomputers and large clusters. However, the high acquisition and maintenance cost of these computers, including the physical space, air conditioning, and electrical power, limits the number of simulations of this kind that scientists can perform. Modern commodity graphical cards, based on the CUDA platform, contain graphical processing units (GPUs) composed of hundreds of processors that can simultaneously execute thousands of threads and thus constitute a low-cost solution for many high-performance computing applications. In this work, we present a CUDA algorithm that enables the execution, on multiple GPUs, of simulations of large-scale networks composed of biologically realistic Hodgkin-Huxley neurons. The algorithm represents each neuron as a CUDA thread, which solves the set of coupled differential equations that model each neuron. Communication among neurons located in different GPUs is coordinated by the CPU. We obtained speedups of 40 for the simulation of 200k neurons that received random external input and speedups of 9 for a network with 200k neurons and 20M neuronal connections, in a single computer with two graphic boards with two GPUs each, when compared with a modern quad-core CPU. Copyright (C) 2010 John Wiley & Sons, Ltd.