Scroll-Wave Dynamics in Human Cardiac Tissue: Lessons from a Mathematical Model with Inhomogeneities and Fiber Architecture

Cardiac arrhythmias, such as ventricular tachycardia (VT) and ventricular fibrillation (VF), are among the leading causes of death in the industrialized world. These are associated with the formation of spiral and scroll waves of electrical activation in cardiac tissue; single spiral and scroll waves are believed to be associated with VT whereas their turbulent analogs are associated with VF. Thus, the study of these waves is an important biophysical problem. We present a systematic study of the combined effects of muscle-fiber rotation and inhomogeneities on scroll-wave dynamics in the TNNP (ten Tusscher Noble Noble Panfilov) model for human cardiac tissue. In particular, we use the three-dimensional TNNP model with fiber rotation and consider both conduction and ionic inhomogeneities. We find that, in addition to displaying a sensitive dependence on the positions, sizes, and types of inhomogeneities, scroll-wave dynamics also depends delicately upon the degree of fiber rotation. We find that the tendency of scroll waves to anchor to cylindrical conduction inhomogeneities increases with the radius of the inhomogeneity. Furthermore, the filament of the scroll wave can exhibit drift or meandering, transmural bending, twisting, and break-up. If the scroll-wave filament exhibits weak meandering, then there is a fine balance between the anchoring of this wave at the inhomogeneity and a disruption of wave-pinning by fiber rotation. If this filament displays strong meandering, then again the anchoring is suppressed by fiber rotation; also, the scroll wave can be eliminated from most of the layers only to be regenerated by a seed wave. Ionic inhomogeneities can also lead to an anchoring of the scroll wave; scroll waves can now enter the region inside an ionic inhomogeneity and can display a coexistence of spatiotemporal chaos and quasi-periodic behavior in different parts of the simulation domain. We discuss the experimental implications of our study.

Two-Stage Extended Kalman Filters with Derivative-Free Local Linearizations

This paper proposes a derivative-free two-stage extended Kalman filter (2-EKF) especially suited for state and parameter identification of mechanical oscillators under Gaussian white noise. Two sources of modeling uncertainties are considered: (1) errors in linearization, and (2) an inadequate system model. The state vector is presently composed of the original dynamical/parameter states plus the so-called bias states accounting for the unmodeled dynamics. An extended Kalman estimation concept is applied within a framework predicated on explicit and derivative-free local linearizations (DLL) of nonlinear drift terms in the governing stochastic differential equations (SDEs). The original and bias states are estimated by two separate filters; the bias filter improves the estimates of the original states. Measurements are artificially generated by corrupting the numerical solutions of the SDEs with noise through an implicit form of a higher-order linearization. Numerical illustrations are provided for a few single- and multidegree-of-freedom nonlinear oscillators, demonstrating the remarkable promise that 2-EKF holds over its more conventional EKF-based counterparts. DOI: 10.1061/(ASCE)EM.1943-7889.0000255. (C) 2011 American Society of Civil Engineers.

Stream-Close: Fast mining of Closed Frequent Itemsets in high speed data streams

With the emergence of large-volume and high-speed streaming data, the recent techniques for stream mining of CFIpsilas (closed frequent itemsets) will become inefficient. When concept drift occurs at a slow rate in high speed data streams, the rate of change of information across different sliding windows will be negligible. So, the user wonpsilat be devoid of change in information if we slide window by multiple transactions at a time. Therefore, we propose a novel approach for mining CFIpsilas cumulatively by making sliding width(ges1) over high speed data streams. However, it is nontrivial to mine CFIpsilas cumulatively over stream, because such growth may lead to the generation of exponential number of candidates for closure checking. In this study, we develop an efficient algorithm, stream-close, for mining CFIpsilas over stream by exploring some interesting properties. Our performance study reveals that stream-close achieves good scalability and has promising results.

Role of the cesium antimonide layer in the Na2KSb/Cs3Sb photocathode

When radiation of sufficiently high energy is incident on the surface of a semiconductor photocathode, electrons are excited from the valence band to the conduction band and these may contribute to the photocurrent. The photocurrent in a single-layer cathode is found to be small, because of collisions within the cathode material, the electron affinity condition, etc. It is observed that when a thin layer of n-type cesium antimonide (Cs3Sb) is deposited over a p-type layer of sodium potassium antimonide (Na2KSb), there occurs a sharp rise in the photocurrent. The causes for the dramatic increase in the photocurrent obtainable from a sodium potassium antimonide cathode, by depositing a thin layer of cesium antimonide are analyzed in this article. It has been shown that the interface between sodium potassium antimonide and cesium antimonide can result in lowering of the electron affinity to a level below the bottom of the conduction band of sodium potassium antimonide. The drift field that arises at the heterointerface enables the electrons to reach the surface, leading to the emission of almost all the photogenerated electrons within the cathode. The processes involved in photoemission from such a double-layer cathode are examined from a theoretical point of view. The spectral response of the two-layer cathode is also found to be better than that of a single-layer cathode.

Elerctron Transport Coefficients in N2 in Crossed Fields at Low E/p

The numerical solutions of Boltzmann transpott equation for the energy distribution of electrons moving in crossed fields in nitrogen have been obtained for 100 ÿ E/p ÿ 1000 V M-1 Torr-1 and for 0ÿ B/p ÿ 0.02 Tesla Torr-1 using the concept of energy dependent effective field intensity. From the derived distribution functions the electron mean energy, the tranaverse and perpendicular drift velocities and the averaged effective field intensity (Eavef) which signifies the average field intensity experienced by electron swarms in E àB field have been derived. The maximum difference between the electron mean energy for a given E ÃÂB field and that corresponding to Eavef/p (p is the gas pressure) is found to be within ñ3.5%.

First passage distributions for long memory processes

We study the distribution of first passage time for Levy type anomalous diffusion. A fractional Fokker-Planck equation framework is introduced.For the zero drift case, using fractional calculus an explicit analytic solution for the first passage time density function in terms of Fox or H-functions is given. The asymptotic behaviour of the density function is discussed. For the nonzero drift case, we obtain an expression for the Laplace transform of the first passage time density function, from which the mean first passage time and variance are derived.

A weighted average based external clock synchronization protocol for wireless sensor networks

Clock synchronization is an extremely important requirement of wireless sensor networks(WSNs). There are many application scenarios such as weather monitoring and forecasting etc. where external clock synchronization may be required because WSN itself may consists of components which are not connected to each other. A usual approach for external clock synchronization in WSNs is to synchronize the clock of a reference node with an external source such as UTC, and the remaining nodes synchronize with the reference node using an internal clock synchronization protocol. In order to provide highly accurate time, both the offset and the drift rate of each clock with respect to reference node are estimated from time to time, and these are used for getting correct time from local clock reading. A problem with this approach is that it is difficult to estimate the offset of a clock with respect to the reference node when drift rate of clocks varies over a period of time. In this paper, we first propose a novel internal clock synchronization protocol based on weighted averaging technique, which synchronizes all the clocks of a WSN to a reference node periodically. We call this protocol weighted average based internal clock synchronization(WICS) protocol. Based on this protocol, we then propose our weighted average based external clock synchronization(WECS) protocol. We have analyzed the proposed protocols for maximum synchronization error and shown that it is always upper bounded. Extensive simulation studies of the proposed protocols have been carried out using Castalia simulator. Simulation results validate our theoretical claim that the maximum synchronization error is always upper bounded and also show that the proposed protocols perform better in comparison to other protocols in terms of synchronization accuracy. A prototype implementation of the proposed internal clock synchronization protocol using a few TelosB motes also validates our claim.

Photocatalytic degradation with combustion synthesized WO3 and WO3-TiO2 mixed oxides under UV and visible light

The photocatalytic activity of commercial titanium dioxide under UV and visible radiation was improved by composites of tungsten trioxide (WO3) with TiO2. WO3 was prepared by solution combustion synthesis and the mixed oxides/composites of WO3-TiO2 were prepared in different weight ratios (0, 0.10, 0.15, 0.20, 0.25, 0.50, 0.75, and 1) by physical mixing. These catalysts were characterized by XRD, DRS, BET, SEM, TEM, pH drift method, TGA and photoluminescence. The photocatalytic activity varies with the WO3 loading in the composites. The optimum loading of WO3 in the composites was found to be 15 wt% for both UV and visible radiation. This loading showed faster dye degradation rate than commercial TiO2 (TiO2-C) and WO3 (WO3-C). The effect of initial concentrations of methylene blue (MB) and orange G (OG) and the effect of the functional group on dye degradation was studied with both anionic and cationic dyes with 15 wt% WO3-TiO2. (C) 2012 Elsevier B.V. All rights reserved.

High precision Mg isotope measurements of meteoritic samples by secondary ion mass spectrometry

The possibility of establishing an accurate relative chronology of the early solar system events based on the decay of short-lived Al-26 to Mg-26 (half-life of 0.72 Myr) depends on the level of homogeneity (or heterogeneity) of Al-26 and Mg isotopes. However, this level is difficult. to constrain precisely because of the very high precision needed for the determination of isotopic ratios, typically of +/- 5 ppm. In this study, we report for the first time a detailed analytical protocol developed for high precision in situ Mg isotopic measurements ((25)mg/(24)mg and (26)mg/Mg-24 ratios, as well as Mg-26 excess) by MC-SIMS. As the data reduction process is critical for both accuracy and precision of the final isotopic results, factors such as the Faraday cup (FC) background drift and matrix effects on instrumental fractionation have been investigated. Indeed these instrumental effects impacting the measured Mg-isotope ratios can be as large or larger than the variations we are looking for to constrain the initial distribution of Al-26 and Mg isotopes in the early solar system. Our results show that they definitely are limiting factors regarding the precision of Mg isotopic compositions, and that an under- or over-correction of both FC background instabilities and instrumental isotopic fractionation leads to important bias on delta Mg-25, delta(26)mg and Delta Mg-26 values (for example, olivines not corrected for FC background drifts display Delta Mg-26 values that can differ by as much as 10 ppm from the truly corrected value). The new data reduction process described here can then be applied to meteoritic samples (components of chondritic meteorites for instance) to accurately establish their relative chronology of formation.

Isolation of a High Affinity Neutralizing Monoclonal Antibody against 2009 Pandemic H1N1 Virus That Binds at the `Sa' Antigenic Site

Influenza virus evades host immunity through antigenic drift and shift, and continues to circulate in the human population causing periodic outbreaks including the recent 2009 pandemic. A large segment of the population was potentially susceptible to this novel strain of virus. Historically, monoclonal antibodies (MAbs) have been fundamental tools for diagnosis and epitope mapping of influenza viruses and their importance as an alternate treatment option is also being realized. The current study describes isolation of a high affinity (K-D = 2.1 +/- 0.4 pM) murine MAb, MA2077 that binds specifically to the hemagglutinin (HA) surface glycoprotein of the pandemic virus. The antibody neutralized the 2009 pandemic H1N1 virus in an in vitro microneutralization assay (IC50 = 0.08 mu g/ml). MA2077 also showed hemagglutination inhibition activity (HI titre of 0.50 mu g/ml) against the pandemic virus. In a competition ELISA, MA2077 competed with the binding site of the human MAb, 2D1 (isolated from a survivor of the 1918 Spanish flu pandemic) on pandemic H1N1 HA. Epitope mapping studies using yeast cell-surface display of a stable HA1 fragment, wherein `Sa' and `Sb' sites were independently mutated, localized the binding site of MA2077 within the `Sa' antigenic site. These studies will facilitate our understanding of antigen antibody interaction in the context of neutralization of the pandemic influenza virus.

A novel filtering framework through Girsanov correction for the identification of nonlinear dynamical systems

Using a Girsanov change of measures, we propose novel variations within a particle-filtering algorithm, as applied to the inverse problem of state and parameter estimations of nonlinear dynamical systems of engineering interest, toward weakly correcting for the linearization or integration errors that almost invariably occur whilst numerically propagating the process dynamics, typically governed by nonlinear stochastic differential equations (SDEs). Specifically, the correction for linearization, provided by the likelihood or the Radon-Nikodym derivative, is incorporated within the evolving flow in two steps. Once the likelihood, an exponential martingale, is split into a product of two factors, correction owing to the first factor is implemented via rejection sampling in the first step. The second factor, which is directly computable, is accounted for via two different schemes, one employing resampling and the other using a gain-weighted innovation term added to the drift field of the process dynamics thereby overcoming the problem of sample dispersion posed by resampling. The proposed strategies, employed as add-ons to existing particle filters, the bootstrap and auxiliary SIR filters in this work, are found to non-trivially improve the convergence and accuracy of the estimates and also yield reduced mean square errors of such estimates vis-a-vis those obtained through the parent-filtering schemes.

Multi-resistance states in the electrical switching behavior of amorphous Si15Te75Ge10 thin films: Possibility of multi-bit storage

Electrical switching studies on amorphous Si15Te75Ge10 thin film devices reveal the existence of two distinct, stable low-resistance, SET states, achieved by varying the electrical input to the device. The multiple resistance levels can be attributed to multi-stage crystallization, as observed from temperature dependant resistance studies. The devices are tested for their ability to be RESET with minimal resistance degradation; further, they exhibit a minimal drift in the SET resistance value even after several months of switching. (c) 2013 Elsevier B.V. All rights reserved.

Temperature dependence of resistance and crystallization in amorphous Si15Te85-xGex thin films

Amorphous thin chalcogenide Si15Te85-xGex films (x: 5, 9, 10, 11, 12) are prepared by flash evaporation and the temperature dependence of resistance of these films has been studied in the temperature range 25-250 degrees C. All the compositions show a linear variation of resistance in this temperature range. Apart from the linear variation, a sharp reduction in resistance at one or at two distinct temperatures (T-TR1/T-TR2) is seen. Thin films annealed at these temperatures, when subjected to X-ray diffraction studies suggest that the dominant crystalline phase at T-TR1 and at T-TR2 is the same and the two dips are associated with varying levels of crystallization. This is also reflected in the atomic force microscopic (AFM) study. Further, the resistance of these two phases shows no drift when the films are annealed for varying lengths of time (10 min to 120 min) suggesting the stability of the phases.

Study of electron focusing in thick GEM based photon detectors using semitransparent photocathodes

The detection efficiency of a GEM based UV sensitive gaseous photomultiplier (GPM) depends on the focusing of electrons from the drift gap to the GEM aperture. We have studied the effect of drift parameters on the efficiency of electron focusing into Thick GEM (THGEM) holes in a GPM with semitransparent UV photoconverter. This study comprises simulation of electron focusing into THGEM holes using GARFIELD for different Ar and Ne based gas mixtures and experimental investigations of the same with P10 gas mixture. (c) 2013 Elsevier B.V. All rights reserved.

A Nearly Exact Reformulation of the Girsanov Linearization for Stochastically Driven Nonlinear Oscillators

The Girsanov linearization method (GLM), proposed earlier in Saha, N., and Roy, D., 2007, ``The Girsanov Linearisation Method for Stochastically Driven Nonlinear Oscillators,'' J. Appl. Mech., 74, pp. 885-897, is reformulated to arrive at a nearly exact, semianalytical, weak and explicit scheme for nonlinear mechanical oscillators under additive stochastic excitations. At the heart of the reformulated linearization is a temporally localized rejection sampling strategy that, combined with a resampling scheme, enables selecting from and appropriately modifying an ensemble of locally linearized trajectories while weakly applying the Girsanov correction (the Radon-Nikodym derivative) for the linearization errors. The semianalyticity is due to an explicit linearization of the nonlinear drift terms and it plays a crucial role in keeping the Radon-Nikodym derivative ``nearly bounded'' above by the inverse of the linearization time step (which means that only a subset of linearized trajectories with low, yet finite, probability exceeds this bound). Drift linearization is conveniently accomplished via the first few (lower order) terms in the associated stochastic (Ito) Taylor expansion to exclude (multiple) stochastic integrals from the numerical treatment. Similarly, the Radon-Nikodym derivative, which is a strictly positive, exponential (super-) martingale, is converted to a canonical form and evaluated over each time step without directly computing the stochastic integrals appearing in its argument. Through their numeric implementations for a few low-dimensional nonlinear oscillators, the proposed variants of the scheme, presently referred to as the Girsanov corrected linearization method (GCLM), are shown to exhibit remarkably higher numerical accuracy over a much larger range of the time step size than is possible with the local drift-linearization schemes on their own.