Invariants of binary differential equations

In this paper, we study binary differential equations a(x, y)dy (2) + 2b(x, y) dx dy + c(x, y)dx (2) = 0, where a, b, and c are real analytic functions. Following the geometric approach of Bruce and Tari in their work on multiplicity of implicit differential equations, we introduce a definition of the index for this class of equations that coincides with the classical Hopf`s definition for positive binary differential equations. Our results also apply to implicit differential equations F(x, y, p) = 0, where F is an analytic function, p = dy/dx, F (p) = 0, and F (pp) not equal aEuro parts per thousand 0 at the singular point. For these equations, we relate the index of the equation at the singular point with the index of the gradient of F and index of the 1-form omega = dy -aEuro parts per thousand pdx defined on the singular surface F = 0.

Recording from Two Neurons: Second-Order Stimulus Reconstruction from Spike Trains and Population Coding

We study the reconstruction of visual stimuli from spike trains, representing the reconstructed stimulus by a Volterra series up to second order. We illustrate this procedure in a prominent example of spiking neurons, recording simultaneously from the two H1 neurons located in the lobula plate of the fly Chrysomya megacephala. The fly views two types of stimuli, corresponding to rotational and translational displacements. Second-order reconstructions require the manipulation of potentially very large matrices, which obstructs the use of this approach when there are many neurons. We avoid the computation and inversion of these matrices using a convenient set of basis functions to expand our variables in. This requires approximating the spike train four-point functions by combinations of two-point functions similar to relations, which would be true for gaussian stochastic processes. In our test case, this approximation does not reduce the quality of the reconstruction. The overall contribution to stimulus reconstruction of the second-order kernels, measured by the mean squared error, is only about 5% of the first-order contribution. Yet at specific stimulus-dependent instants, the addition of second-order kernels represents up to 100% improvement, but only for rotational stimuli. We present a perturbative scheme to facilitate the application of our method to weakly correlated neurons.

Second-order negative-curvature methods for box-constrained and general constrained optimization

A Nonlinear Programming algorithm that converges to second-order stationary points is introduced in this paper. The main tool is a second-order negative-curvature method for box-constrained minimization of a certain class of functions that do not possess continuous second derivatives. This method is used to define an Augmented Lagrangian algorithm of PHR (Powell-Hestenes-Rockafellar) type. Convergence proofs under weak constraint qualifications are given. Numerical examples showing that the new method converges to second-order stationary points in situations in which first-order methods fail are exhibited.

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.

Métodos numéricos para resolução de equações diferenciais ordinárias lineares baseados em interpolação por spline

In this work we have elaborated a spline-based method of solution of inicial value problems involving ordinary differential equations, with emphasis on linear equations. The method can be seen as an alternative for the traditional solvers such as Runge-Kutta, and avoids root calculations in the linear time invariant case. The method is then applied on a central problem of control theory, namely, the step response problem for linear EDOs with possibly varying coefficients, where root calculations do not apply. We have implemented an efficient algorithm which uses exclusively matrix-vector operations. The working interval (till the settling time) was determined through a calculation of the least stable mode using a modified power method. Several variants of the method have been compared by simulation. For general linear problems with fine grid, the proposed method compares favorably with the Euler method. In the time invariant case, where the alternative is root calculation, we have indications that the proposed method is competitive for equations of sifficiently high order.

Implicit differential equations with impasse singularities and singular perturbation problems

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

Periodic solutions of abstract neutral functional differential equations

We characterize the existence of periodic solutions of some abstract neutral functional differential equations with finite and infinite delay when the underlying space is a UMD space. (C) 2011 Elsevier B.V. All rights reserved.

Bose-Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods

We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross-Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose-Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x, y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge-Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank-Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.

Second order gauge theory

A gauge theory of second order in the derivatives of the auxiliary field is constructed following Utiyama's program. A novel field strength G = partial derivative F + fAF arises besides the one of the first order treatment, F = partial derivative A - partial derivative A + fAA. The associated conserved current is obtained. It has a new feature: topological terms are determined from local invariance requirements. Podolsky Generalized Eletrodynamics is derived as a particular case in which the Lagrangian of the gauge field is L-P alpha G(2). In this application the photon mass is estimated. The SU(N) infrared regime is analysed by means of Alekseev-Arbuzov-Baikov's Lagrangian. (c) 2006 Elsevier B.V. All rights reserved.

Cosmic acceleration from second order gauge gravity

We construct a phenomenological theory of gravitation based on a second order gauge formulation for the Lorentz group. The model presents a long-range modification for the gravitational field leading to a cosmological model provided with an accelerated expansion at recent times. We estimate the model parameters using observational data and verify that our estimative for the age of the Universe is of the same magnitude than the one predicted by the standard model. The transition from the decelerated expansion regime to the accelerated one occurs recently (at similar to 9.3 Gyr).

Superspace approach to the renormalization of the O'Raifeartaigh model up to the second order in the linear delta expansion parameter

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Stability of differential equations with piecewise constant argument via discrete equations

Lyapunov stability for a class of differential equation with piecewise constant argument (EPCA) is considered by means of the stability of a discrete equation. Applications to some nonlinear autonomous equations are given improving some linear known cases.

Stability results for impulsive functional differential equations with infinite delay

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

Second-order statistics for equal gain and maximal ratio diversity-combining reception

Exact and closed-form expressions for the level crossing rate and average fade duration are presented for equal gain combining and maximal ratio combining schemes, assuming an arbitrary number of independent branches in a Rayleigh environment. The analytical results are thoroughly validated by simulation.

Second-order statistics for diversity-combining, techniques in Nakagami-fading channels

Exact and closed-form expressions for the level crossing rate and average fade duration are presented for the M branch pure selection combining (PSC), equal gain combining (EGC), and maximal ratio combining (MRC) techniques, assuming independent branches in a Nakagami environment. The analytical results are thoroughly validated by reducing the general case to some special cases, for which the solutions are known, and by means of simulation for the more general case. The model developed here is general and can be easily applied to other fading statistics (e.g., Rice).