A Wavelet based Teager Energy Operator for Spike Detection in Microelectrode Array Recordings

Spike detection in neural recordings is the initial step in the creation of brain machine interfaces. The Teager energy operator (TEO) treats a spike as an increase in the `local' energy and detects this increase. The performance of TEO in detecting action potential spikes suffers due to its sensitivity to the frequency of spikes in the presence of noise which is present in microelectrode array (MEA) recordings. The multiresolution TEO (mTEO) method overcomes this shortcoming of the TEO by tuning the parameter k to an optimal value m so as to match to frequency of the spike. In this paper, we present an algorithm for the mTEO using the multiresolution structure of wavelets along with inbuilt lowpass filtering of the subband signals. The algorithm is efficient and can be implemented for real-time processing of neural signals for spike detection. The performance of the algorithm is tested on a simulated neural signal with 10 spike templates obtained from [14]. The background noise is modeled as a colored Gaussian random process. Using the noise standard deviation and autocorrelation functions obtained from recorded data, background noise was simulated by an autoregressive (AR(5)) filter. The simulations show a spike detection accuracy of 90%and above with less than 5% false positives at an SNR of 2.35 dB as compared to 80% accuracy and 10% false positives reported [6] on simulated neural signals.

Time Dependent Operator-split and Unsplit Schemes for One Dimensional Premixed Flames

This paper is concerned with a study of an operator split scheme and unsplit scheme for the computation of adiabatic freely propagating one-dimensional premixed flames. The study uses unsteady method for both split and unsplit schemes employing implicit chemistry and explicit diffusion, a combination which is stable and convergent. Solution scheme is not sensitive to the initial starting estimate and provides steady state even with straight line profiles (far from steady state) in small number of time steps. Two systems H2-Air and H2-NO (involving complex nitrogen chemistry) are considered in presentinvestigation. Careful comparison shows that the operator split approach is slightly superior than the unsplit when chemistry becomes complex. Comparison of computational times with those of existing steady and unsteady methods seems to suggest that the method employing implicit-explicit algorithm is very efficient and robust.

Negative differential resistance in Gd0.5Sr0.5MnO3: A consequence of Joule heating

Negative differential resistance (NDR) in current-voltage (I-V) characteristics and apparent colossal electroresistance were observed in Gd0.5Sr0.5MnO3 single crystals at low temperatures. The continuous dc I-V measurements showed a marked thermal drift. In addition, temperature of the sample surface was found to be significantly higher than that of the base at high applied currents. Two different strategies namely estimation and diminution of the Joule heating (pulsed I-V measurements) were employed to investigate its role in the electric transport properties. Our experiments reveal that the NDR in Gd0.5Sr0.5MnO3 is a consequence of Joule heating rather than the melting of charge order. (C) 2010 American Institute of Physics. doi:10.1063/1.3486221]

A maximal operator, Carleson's embedding, and tent spaces for vector-valued functions

This thesis is concerned with the area of vector-valued Harmonic Analysis, where the central theme is to determine how results from classical Harmonic Analysis generalize to functions with values in an infinite dimensional Banach space. The work consists of three articles and an introduction. The first article studies the Rademacher maximal function that was originally defined by T. Hytönen, A. McIntosh and P. Portal in 2008 in order to prove a vector-valued version of Carleson's embedding theorem. The boundedness of the corresponding maximal operator on Lebesgue-(Bochner) -spaces defines the RMF-property of the range space. It is shown that the RMF-property is equivalent to a weak type inequality, which does not depend for instance on the integrability exponent, hence providing more flexibility for the RMF-property. The second article, which is written in collaboration with T. Hytönen, studies a vector-valued Carleson's embedding theorem with respect to filtrations. An earlier proof of the dyadic version assumed that the range space satisfies a certain geometric type condition, which this article shows to be also necessary. The third article deals with a vector-valued generalizations of tent spaces, originally defined by R. R. Coifman, Y. Meyer and E. M. Stein in the 80's, and concerns especially the ones related to square functions. A natural assumption on the range space is then the UMD-property. The main result is an atomic decomposition for tent spaces with integrability exponent one. In order to suit the stochastic integrals appearing in the vector-valued formulation, the proof is based on a geometric lemma for cones and differs essentially from the classical proof. Vector-valued tent spaces have also found applications in functional calculi for bisectorial operators. In the introduction these three themes come together when studying paraproduct operators for vector-valued functions. The Rademacher maximal function and Carleson's embedding theorem were applied already by Hytönen, McIntosh and Portal in order to prove boundedness for the dyadic paraproduct operator on Lebesgue-Bochner -spaces assuming that the range space satisfies both UMD- and RMF-properties. Whether UMD implies RMF is thus an interesting question. Tent spaces, on the other hand, provide a method to study continuous time paraproduct operators, although the RMF-property is not yet understood in the framework of tent spaces.

Squeezed states, metaplectic group, and operator möbius transformations

We present an analysis, based on the metaplectic group Mp(2), of the recently introduced single-mode inverse creation and annihilation operators and of the associated eigenstates of different two-photon annihilation operators. We motivate and obtain a quantum operator form of the classical Mobius or fractional linear transformation. The subtle relation to the two unitary irreducible representations of Mp(2) is brought out. For problems involving inverse operators the usefulness of the Bargmann analytic function representation of quantum mechanics is demonstrated. Squeezing, bunching, and photon-number distributions of the four families of states that arise in this context are studied both analytically and numerically

Squeezed vacuum as an eigenstate of two-photon annihilation operator

We introduce the inverse annihilation and creation operators a-1 and a(dagger-1) by their actions on the number states. We show that the squeezed vacuum exp(1/2xia(dagger2)]\0] and squeezed first number state exp[1.2xia(dagger2)]\n = 1] are respectively the eigenstates of the operators (a(dagger-1)a) and (aa(dagger-1)) with the eigenvalue xi.

Polymorphism in Opto-Electronic Materials with a Benzothiazole-fluorene Core: A Consequence of High Conformational Flexibility of pi-Conjugated Backbone and Alkyl Side Chains

Fluorene and its derivatives are well-known organic semiconducting materials in the field of opto-electronic devices because of their charge transport properties. Three new organic semiconducting materials, namely, 2,2'-((9,9-butyl-9H-fluorene-2,7-diyl)bis(4,1 phenylene))bisbenzod]thiazole, C4; 2,2'-((octyl-9H-fluorene-2,7-diyl)bis(4,1 phenylene))bisbenzod]thiazole, C8; and 2,2'-((9,9-dodecayl-9H-fluorene-2,7-diyl)bis(4,1 phenylene))bisbenzod]thiazole, C12 with a benzothiazole-fluorene backbone, were synthesized and characterized for their photophysical properties. A phenomenon of concomitant polymorphism has been investigated in the first two derivatives (C4 and C8) and has been analyzed systematically in terms of the packing characteristics involving pi ... pi interactions. The conformational flexibility of the pi-conjugated 2,2'-(fluorene-2,7-diyl)bis(4,1 phenylene)bisbenzod]thiazole backbone coupled with orientational freedom of the terminal alkyl chains were found to be the key factors responsible for these polymorphic modifications. Attempts to grow suitable crystals for single crystal X-ray diffraction of compound C12 were unsuccessful.

The extended Enskog operator for simple fluids with continuous potentials: single particle and collective properties

We generalized the Enskog theory originally developed for the hard-sphere fluid to fluids with continuous potentials, such as the Lennard–Jones. We derived the expression for the k and ω dependent transport coefficient matrix which enables us to calculate the transport coefficients for arbitrary length and time scales. Our results reduce to the conventional Chapman–Enskog expression in the low density limit and to the conventional k dependent Enskog theory in the hard-sphere limit. As examples, the self-diffusion of a single atom, the vibrational energy relaxation, and the activated barrier crossing dynamics problem are discussed.

Wear of metals: a consequence of stable/unstable material response

Wear of metals in dry sliding is dictated by the material response to traction. This is demonstrated by considering the wear of aluminium and titanium alloys. In a regime of stable homogeneous deformation the material approaching the surface from the bulk passes through microprocessing zones of flow, fracture, comminution and compaction to generate a protective tribofilm that retains the interaction in the mild wear regime. If the response leads to microstructural instabilities such as adiabatic shear bands, the near-surface zone consists of stacks of 500 nm layers situated parallel to the sliding direction. Microcracks are generated below the surface to propagate normally away from the surface though microvoids situated in the layers, until it reaches a depth of 10-20 mum. A rectangular laminate debris consisting of a 20-40 layer stack is produced, The wear in this mode is severe.

An operator-splitting finite element method for the efficient parallel solution of multidimensional population balance systems

A finite element method for solving multidimensional population balance systems is proposed where the balance of fluid velocity, temperature and solute partial density is considered as a two-dimensional system and the balance of particle size distribution as a three-dimensional one. The method is based on a dimensional splitting into physical space and internal property variables. In addition, the operator splitting allows to decouple the equations for temperature, solute partial density and particle size distribution. Further, a nodal point based parallel finite element algorithm for multi-dimensional population balance systems is presented. The method is applied to study a crystallization process assuming, for simplicity, a size independent growth rate and neglecting agglomeration and breakage of particles. Simulations for different wall temperatures are performed to show the effect of cooling on the crystal growth. Although the method is described in detail only for the case of d=2 space and s=1 internal property variables it has the potential to be extendable to d+s variables, d=2, 3 and s >= 1. (C) 2011 Elsevier Ltd. All rights reserved.

Dilations of Gamma-contractions by solving operator equations

For a contraction P and a bounded commutant S of P. we seek a solution X of the operator equation S - S*P = (1 - P* P)(1/2) X (1 - P* P)(1/2) where X is a bounded operator on (Ran) over bar (1 - P* P)(1/2) with numerical radius of X being not greater than 1. A pair of bounded operators (S, P) which has the domain Gamma = {(z(1) + z(2), z(2)): vertical bar z(1)vertical bar < 1, vertical bar z(2)vertical bar <= 1} subset of C-2 as a spectral set, is called a P-contraction in the literature. We show the existence and uniqueness of solution to the operator equation above for a Gamma-contraction (S, P). This allows us to construct an explicit Gamma-isometric dilation of a Gamma-contraction (S, P). We prove the other way too, i.e., for a commuting pair (S, P) with parallel to P parallel to <= 1 and the spectral radius of S being not greater than 2, the existence of a solution to the above equation implies that (S, P) is a Gamma-contraction. We show that for a pure F-contraction (S, P), there is a bounded operator C with numerical radius not greater than 1, such that S = C + C* P. Any Gamma-isometry can be written in this form where P now is an isometry commuting with C and C. Any Gamma-unitary is of this form as well with P and C being commuting unitaries. Examples of Gamma-contractions on reproducing kernel Hilbert spaces and their Gamma-isometric dilations are discussed. (C) 2012 Elsevier Inc. All rights reserved.

An operator-splitting Galerkin/SUPG finite element method for population balance equations : stability and convergence

We present a heterogeneous finite element method for the solution of a high-dimensional population balance equation, which depends both the physical and the internal property coordinates. The proposed scheme tackles the two main difficulties in the finite element solution of population balance equation: (i) spatial discretization with the standard finite elements, when the dimension of the equation is more than three, (ii) spurious oscillations in the solution induced by standard Galerkin approximation due to pure advection in the internal property coordinates. The key idea is to split the high-dimensional population balance equation into two low-dimensional equations, and discretize the low-dimensional equations separately. In the proposed splitting scheme, the shape of the physical domain can be arbitrary, and different discretizations can be applied to the low-dimensional equations. In particular, we discretize the physical and internal spaces with the standard Galerkin and Streamline Upwind Petrov Galerkin (SUPG) finite elements, respectively. The stability and error estimates of the Galerkin/SUPG finite element discretization of the population balance equation are derived. It is shown that a slightly more regularity, i.e. the mixed partial derivatives of the solution has to be bounded, is necessary for the optimal order of convergence. Numerical results are presented to support the analysis.

Operator-splitting finite element algorithms for computations of high-dimensional parabolic problems

An operator-splitting finite element method for solving high-dimensional parabolic equations is presented. The stability and the error estimates are derived for the proposed numerical scheme. Furthermore, two variants of fully-practical operator-splitting finite element algorithms based on the quadrature points and the nodal points, respectively, are presented. Both the quadrature and the nodal point based operator-splitting algorithms are validated using a three-dimensional (3D) test problem. The numerical results obtained with the full 3D computations and the operator-split 2D + 1D computations are found to be in a good agreement with the analytical solution. Further, the optimal order of convergence is obtained in both variants of the operator-splitting algorithms. (C) 2012 Elsevier Inc. All rights reserved.

Quadrature approximation properties of the spiral-phase quadrature transform

The notion of the 1-D analytic signal is well understood and has found many applications. At the heart of the analytic signal concept is the Hilbert transform. The problem in extending the concept of analytic signal to higher dimensions is that there is no unique multidimensional definition of the Hilbert transform. Also, the notion of analyticity is not so well under stood in higher dimensions. Of the several 2-D extensions of the Hilbert transform, the spiral-phase quadrature transform or the Riesz transform seems to be the natural extension and has attracted a lot of attention mainly due to its isotropic properties. From the Riesz transform, Larkin et al. constructed a vortex operator, which approximates the quadratures based on asymptotic stationary-phase analysis. In this paper, we show an alternative proof for the quadrature approximation property by invoking the quasi-eigenfunction property of linear, shift-invariant systems. We show that the vortex operator comes up as a natural consequence of applying this property. We also characterize the quadrature approximation error in terms of its energy as well as the peak spatial-domain error. Such results are available for 1-D signals, but their counter part for 2-D signals have not been provided. We also provide simulation results to supplement the analytical calculations.

Riesz transforms and multipliers for the Grushin operator

We show that Riesz transforms associated to the Grushin operator G = -Delta - |x|(2 similar to) (t) (2) are bounded on L (p) (a''e (n+1)). We also establish an analogue of the Hormander-Mihlin Multiplier Theorem and study Bochner-Riesz means associated to the Grushin operator. The main tools used are Littlewood-Paley theory and an operator-valued Fourier multiplier theorem due to L. Weis.