Modeling the effects of anesthesia on the electroencephalogram

Changes to the electroencephalogram (EEG) observed during general anesthesia are modeled with a physiological mean field theory of electrocortical activity. To this end a parametrization of the postsynaptic impulse response is introduced which takes into account pharmacological effects of anesthetic agents on neuronal ligand-gated ionic channels. Parameter sets for this improved theory are then identified which respect known anatomical constraints and predict mean firing rates and power spectra typically encountered in human subjects. Through parallelized simulations of the eight nonlinear, two-dimensional partial differential equations on a grid representing an entire human cortex, it is demonstrated that linear approximations are sufficient for the prediction of a range of quantitative EEG variables. More than 70 000 plausible parameter sets are finally selected and subjected to a simulated induction with the stereotypical inhaled general anesthetic isoflurane. Thereby 86 parameter sets are identified that exhibit a strong “biphasic” rise in total power, a feature often observed in experiments. A sensitivity study suggests that this “biphasic” behavior is distinguishable even at low agent concentrations. Finally, our results are briefly compared with previous work by other groups and an outlook on future fits to experimental data is provided.

Electrorhythmogenesis and anaesthesia in a physiological mean field theory

We solve eight partial-differential, two-dimensional, nonlinear mean field equations, which describe the dynamics of large populations of cortical neurons. Linearized versions of these equations have been used to generate the strong resonances observed in the human EEG, in particular the α-rhythm (8–), with physiologically plausible parameters. We extend these results here by numerically solving the full equations on a cortex of realistic size, which receives appropriately “colored” noise as extra-cortical input. A brief summary of the numerical methods is provided. As an outlook to future applications, we explain how the effects of GABA-enhancing general anaesthetics can be simulated and present first results.

On wave action and phase in the non-canonical Hamiltonian formulation

The long time–evolution of disturbances to slowly–varying solutions of partial differential equations is subject to the adiabatic invariance of the wave action. Generally, this approximate conservation law is obtained under the assumption that the partial differential equations are derived from a variational principle or have a canonical Hamiltonian structure. Here, the wave action conservation is examined for equations that possess a non–canonical (Poisson) Hamiltonian structure. The linear evolution of disturbances in the form of slowly varying wavetrains is studied using a WKB expansion. The properties of the original Hamiltonian system strongly constrain the linear equations that are derived, and this is shown to lead to the adiabatic invariance of a wave action. The connection between this (approximate) invariance and the (exact) conservation laws of pseudo–energy and pseudomomentum that exist when the basic solution is exactly time and space independent is discussed. An evolution equation for the slowly varying phase of the wavetrain is also derived and related to Berry's phase.

Maximum principles for vectorial approximate minimisers on non-convex functionals

We establish Maximum Principles which apply to vectorial approximate minimizers of the general integral functional of Calculus of Variations. Our main result is a version of the Convex Hull Property. The primary advance compared to results already existing in the literature is that we have dropped the quasiconvexity assumption of the integrand in the gradient term. The lack of weak Lower semicontinuity is compensated by introducing a nonlinear convergence technique, based on the approximation of the projection onto a convex set by reflections and on the invariance of the integrand in the gradient term under the Orthogonal Group. Maximum Principles are implied for the relaxed solution in the case of non-existence of minimizers and for minimizing solutions of the Euler–Lagrange system of PDE.

A dynamical systems framework for intermittent data assimilation

We consider the problem of discrete time filtering (intermittent data assimilation) for differential equation models and discuss methods for its numerical approximation. The focus is on methods based on ensemble/particle techniques and on the ensemble Kalman filter technique in particular. We summarize as well as extend recent work on continuous ensemble Kalman filter formulations, which provide a concise dynamical systems formulation of the combined dynamics-assimilation problem. Possible extensions to fully nonlinear ensemble/particle based filters are also outlined using the framework of optimal transportation theory.

Dynamics of a class of ODEs more general than almost periodic

A temporally global solution, if it exists, of a nonautonomous ordinary differential equation need not be periodic, almost periodic or almost automorphic when the forcing term is periodic, almost periodic or almost automorphic, respectively. An alternative class of functions extending periodic and almost periodic functions which has the property that a bounded temporally global solution solution of a nonautonomous ordinary differential equation belongs to this class when the forcing term does is introduced here. Specifically, the class of functions consists of uniformly continuous functions, defined on the real line and taking values in a Banach space, which have pre-compact ranges. Besides periodic and almost periodic functions, this class also includes many nonrecurrent functions. Assuming a hyperbolic structure for the unperturbed linear equation and certain properties for the linear and nonlinear parts, the existence of a special bounded entire solution, as well the existence of stable and unstable manifolds of this solution are established. Moreover, it is shown that this solution and these manifolds inherit the temporal behaviour of the vector field equation. In the stable case it is shown that this special solution is the pullback attractor of the system. A class of infinite dimensional examples involving a linear operator consisting of a time independent part which generates a C(0)-semigroup plus a small time dependent part is presented and applied to systems of coupled heat and beam equations. (C) 2010 Elsevier Ltd. All rights reserved.

A collocation method for solving singular integro-differential equations

This paper describes a collocation method for numerically solving Cauchy-type linear singular integro-differential equations. The numerical method is based on the transformation of the integro-differential equation into an integral equation, and then applying a collocation method to solve the latter. The collocation points are chosen as the Chebyshev nodes. Uniform convergence of the resulting method is then discussed. Numerical examples are presented and solved by the numerical techniques.

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.

BREAKING WAVES IN AN IDEAL QUARK GLUON PLASMA

We study the propagation of perturbations in the energy density in a quark gluon plasma. Expanding the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations we obtain a nonlinear differential equation called the breaking wave equation. We solve it numerically and follow the time-evolution of initially localized pulses. We find that, quite unexpectedly, these pulses live for a very long time (compared to the reaction time-scales) before breaking. In practice, they mimick the Korteweg-de Vries solitons. Their existence may have some observable consequences.

A Direct Numerov Sixth-order Numerical Scheme to Accurately Solve the Unidimensional Poisson Equation with Dirichlet Boundary Conditions

In this article, we present an analytical direct method, based on a Numerov three-point scheme, which is sixth order accurate and has a linear execution time on the grid dimension, to solve the discrete one-dimensional Poisson equation with Dirichlet boundary conditions. Our results should improve numerical codes used mainly in self-consistent calculations in solid state physics.

Controle inteligente de sistemas eletroidráulicos utilizando redes neurais artificiais

This work describes the development of a nonlinear control strategy for an electro-hydraulic actuated system. The system to be controlled is represented by a third order ordinary differential equation subject to a dead-zone input. The control strategy is based on a nonlinear control scheme, combined with an artificial intelligence algorithm, namely, the method of feedback linearization and an artificial neural network. It is shown that, when such a hard nonlinearity and modeling inaccuracies are considered, the nonlinear technique alone is not enough to ensure a good performance of the controller. Therefore, a compensation strategy based on artificial neural networks, which have been notoriously used in systems that require the simulation of the process of human inference, is used. The multilayer perceptron network and the radial basis functions network as well are adopted and mathematically implemented within the control law. On this basis, the compensation ability considering both networks is compared. Furthermore, the application of new intelligent control strategies for nonlinear and uncertain mechanical systems are proposed, showing that the combination of a nonlinear control methodology and artificial neural networks improves the overall control system performance. Numerical results are presented to demonstrate the efficacy of the proposed control system

Stability of neutral functional differential equations in terms of measures

In this paper we investigate the relationships between different concepts of stability in measure for the solutions of an autonomous or periodic neutral functional differential equation.

On retarded differential equations with impulses

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Force and displacement transmissibility of a nonlinear isolator with high-static-low-dynamic-stiffness

Engineers often face the challenge of reducing the level of vibrations experienced by a given payload or those transmitted to the support structure to which a vibrating source is attached. In order to increase the range over which vibrations are isolated, soft mounts are often used in practice. The drawback of this approach is the static displacement may be too large for reasons of available space for example. Ideally, a vibration isolator should have a high-static stiffness, to withstand static loads without too large a displacement, and at the same time, a low dynamic stiffness so that the natural frequency of the system is as low as possible which will result in an increased isolation region. These two effects are mutually exclusive in linear isolators but can be overcome if properly configured nonlinear isolators are used. This paper is concerned with the characterisation of such a nonlinear isolator comprising three springs, two of which are configured to reduce the dynamic stiffness of the isolator. The dynamic behaviour of the isolator supporting a lumped mass is investigated using force and displacement transmissibility, which are derived by modelling the dynamic system as a single-degree-of-freedom system. This results in the system dynamics being approximately described by the Duffing equation. For a linear isolator, the dynamics of the system are the same regardless if the source of the excitation is a harmonic force acting on the payload (force transmissibility) or a harmonic motion of the base (displacement transmissibility) on which the payload is mounted. In this paper these two expressions are compared for the nonlinear isolator and it is shown that they differ. A particular feature of the displacement transmissibility is that the response is unbounded at the nonlinear resonance frequency unless the damping in the isolator is greater than some threshold value, which is not the case for force transmissibility. An explanation for this is offered in the paper. (C) 2011 Elsevier Ltd. All rights reserved.

Particular solution for anomalous diffusion equation with source term

In literature the phenomenon of diffusion has been widely studied, however for nonextensive systems which are governed by a nonlinear stochastic dynamic, there are a few soluble models. The purpose of this study is to present the solution of the nonlinear Fokker-Planck equation for a model of potential with barrier considering a term of absorption. Systems of this nature can be observed in various chemical or biological processes and their solution enriches the studies of existing nonextensive systems.