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.

Fluidity, Permeability And Antioxidant Behavior Of Model Membranes Incorporated With Alpha-Tocopherol And Vitamin-E Acetate

Magnetic resonance studies reveal a marked difference between the binding of α-tocopherol and that of the corresponding acetate (vitamin E acetate) with dipalmitoylphosphatidylcholine (DPPC) vesicles. This is reflected in differences in the phase-transition curves of the DPPC vesicles incorporated with the two compounds, as well as in the 13C relaxation times and line widths. A model for the incorporation of these molecules in lipid bilayers has been suggested. α-Tocopherol binds strongly with the lipids, possibly through a hydrogen bond formation between the hydroxyl group of the former and one of the oxygen atoms of the latter. The possibility of such a hydrogen bond formation is excluded in vitamin E acetate, which binds loosely through the normal hydrophobic interaction. The model for lipid-vitamin interaction explains the in vitro decomposition of H2O2 by α-tocopherol. α-Tocopherol in conjuction with H2O2 can also act as a free-radical scavenger in the lipid phase. The incorporation of α-tocopherol and vitamin E acetate in DPPC vesicles enhances the permeability of lipid bilayers for small molecules such as sodium ascorbate.

Drosophila indirect flight muscle specific Act88F actin mutants as a model system for studying congenital myopathies of the human ACTA1 skeletal muscle actin gene

Most human ACTA1 skeletal actin gene mutations cause dominant, congenital myopathies often with severely reduced muscle function and neonatal mortality. High sequence conservation of actin means many mutated ACTA1 residues are identical to those in the Drosophila Act88F, an indirect flight muscle specific sarcomeric actin. Four known Act88F mutations occur at the same actin residues mutated in ten ACTA1 nemaline mutations, A138D/P, R256H/L, G268C/D/R/S and R372C/S. These Act88F mutants were examined for similar muscle phenotypes. Mutant homozygotes show phenotypes ranging from a lack of myofibrils to almost normal sarcomeres at eclosion. Aberrant Z-disc-like structures and serial Z-disc arrays, ‘zebra bodies’, are observed in homozygotes and heterozygotes of all four Act88F mutants. These electron-dense structures show homologies to human nemaline bodies/rods, but are much smaller than those typically found in the human myopathy. We conclude that the Drosophila indirect flight muscles provide a good model system for studying ACTA1 mutations.

Markov Chain Model for Cracking Behavior of Reinforced Concrete Beams

In this paper, nonhomogeneous Markov chains are proposed for modeling the cracking behavior of reinforced concrete beams subjected to monotonically increasing loads. The model facilitates prediction of the maximum crackwidth at a given load given the crackwidth at a lower load level, and thus leads to a better understanding of the cracking phenomenon. To illustrate the methodology developed, the results of three reinforced concrete beams tested in the laboratory are analyzed and presented.

Some experiments on human memory and a new neural model

An associative memory with parallel architecture is presented. The neurons are modelled by perceptrons having only binary, rather than continuous valued input. To store m elements each having n features, m neurons each with n connections are needed. The n features are coded as an n-bit binary vector. The weights of the n connections that store the n features of an element has only two values -1 and 1 corresponding to the absence or presence of a feature. This makes the learning very simple and straightforward. For an input corrupted by binary noise, the associative memory indicates the element that is closest (in terms of Hamming distance) to the noisy input. In the case where the noisy input is equidistant from two or more stored vectors, the associative memory indicates two or more elements simultaneously. From some simple experiments performed on the human memory and also on the associative memory, it can be concluded that the associative memory presented in this paper is in some respect more akin to a human memory than a Hopfield model.

Non-Fermi liquid behavior in a mean-field study of the t-J model: Role of zero-point spin fluctuations

We present analytic results to show that the Schwinger-boson hole-fermion mean-field state exhibits non-Fermi liquid behavior due to spin-charge separation. The physical electron Green's function consists of three additive components. (a) A Fermi-liquid component associated with the bose condensate. (b) A non-Fermi liquid component which has a logarithmic peak and a long tail that gives rise to a linear density of states that is symmetric about the Fermi level and a momentum distribution function with a logarithmic discontinuity at the Fermi surface. (c) A second non-Fermi liquid component associated with the thermal bosons which leads to a constant density of states. It is shown that zero-point fluctuations associated with the spin-degrees of freedom are responsible for the logarithmic instabilities and the restoration of particle-hole symmetry close to the Fermi surface.

Optimal management of beaver population using a reduced-order distributed parameter model and single network adaptive critics

Beavers are often found to be in conflict with human interests by creating nuisances like building dams on flowing water (leading to flooding), blocking irrigation canals, cutting down timbers, etc. At the same time they contribute to raising water tables, increased vegetation, etc. Consequently, maintaining an optimal beaver population is beneficial. Because of their diffusion externality (due to migratory nature), strategies based on lumped parameter models are often ineffective. Using a distributed parameter model for beaver population that accounts for their spatial and temporal behavior, an optimal control (trapping) strategy is presented in this paper that leads to a desired distribution of the animal density in a region in the long run. The optimal control solution presented, imbeds the solution for a large number of initial conditions (i.e., it has a feedback form), which is otherwise nontrivial to obtain. The solution obtained can be used in real-time by a nonexpert in control theory since it involves only using the neural networks trained offline. Proper orthogonal decomposition-based basis function design followed by their use in a Galerkin projection has been incorporated in the solution process as a model reduction technique. Optimal solutions are obtained through a "single network adaptive critic" (SNAC) neural-network architecture.

Chain folding and A:T pairing in human telomeric DNA: a model-building and molecular dynamics study

The various types of chain folding and possible intraloop as well as interloop base pairing in human telomeric DNA containing d(TTAG(3)) repeats have been investigated by model-building, molecular mechanics, and molecular dynamics techniques. Model-building and molecular mechanics studies indicate that it is possible to build a variety of energetically favorable folded-back structures with the two TTA loops on same side and the 5' end thymines in the two loops forming TATA tetrads involving a number of different intraloop as well as interloop A:T pairing schemes. In these folded-back structures, although both intraloop and interloop Watson-Crick pairing is feasible, no structure is possible with interloop Hoogsteen pairing. MD studies of representative structures indicate that the guanine-tetraplex stem is very rigid and, while the loop regions are relatively much more flexible, most of the hydrogen bonds remain intact throughout the 350-ps in vacuo simulation. The various possible TTA loop structures, although they are energetically similar, have characteristic inter proton distances, which could give rise to unique cross-peaks in two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiments. These folded-back structures with A:T pairings in the loop region help in rationalizing the data from chemical probing and other biochemical studies on human telomeric DNA.

WOMM: A Weak Operational Memory Model

Memory models of shared memory concurrent programs define the values a read of a shared memory location is allowed to see. Such memory models are typically weaker than the intuitive sequential consistency semantics to allow efficient execution. In this paper, we present WOMM (abbreviation for Weak Operational Memory Model) that formally unifies two sources of weak behavior in hardware memory models: reordering of instructions and weakly consistent memory. We show that a large number of optimizations are allowed by WOMM. We also show that WOMM is weaker than a number of hardware memory models. Consequently, if a program behaves correctly under WOMM, it will be correct with respect to those hardware memory models. Hence, WOMM can be used as a formally specified abstraction of the hardware memory models. Moreover; unlike most weak memory models, WOMM is described using operational semantics, making it easy to integrate into a model checker for concurrent programs. We further show that WOMM has an important property - it has sequential consistency semantics for datarace-free programs.

On the integrability and chaotic behavior of an ecological model

A three-species food chain model is studied analytically as well as numerically. Integrability of the model is studied using Painleve analysis while chaotic behavior is studied using numerical techniques, such as calculation of Lyapunov exponents, plotting the bifurcation diagram and phase plots. We correct and critically comment on the wrong results reported recently on this ecological model, in a paper by Rai [1995].

A model for tool-use traditions in primates: implications for the coevolution of culture and cognition

Inspired by the demonstration that tool-use variants among wild chimpanzees and orangutans qualify as traditions (or cultures), we developed a formal model to predict the incidence of these acquired specializations among wild primates and to examine the evolution of their underlying abilities. We assumed that the acquisition of the skill by an individual in a social unit is crucially controlled by three main factors, namely probability of innovation, probability of socially biased learning, and the prevailing social conditions (sociability, or number of potential experts at close proximity). The model reconfirms the restriction of customary tool use in wild primates to the most intelligent radiation, great apes; the greater incidence of tool use in more sociable populations of orangutans and chimpanzees; and tendencies toward tool manufacture among the most sociable monkeys. However, it also indicates that sociable gregariousness is far more likely to produce the maintenance of invented skills in a population than solitary life, where the mother is the only accessible expert. We therefore used the model to explore the evolution of the three key parameters. The most likely evolutionary scenario is that where complex skills contribute to fitness, sociability and/or the capacity for socially biased learning increase, whereas innovative abilities (i.e., intelligence) follow indirectly. We suggest that the evolution of high intelligence will often be a byproduct of selection on abilities for socially biased learning that are needed to acquire important skills, and hence that high intelligence should be most common in sociable rather than solitary organisms. Evidence for increased sociability during hominin evolution is consistent with this new hypothesis. (C) 2003 Elsevier Science Ltd. All rights reserved.

Towards formal verification of analog mixed signal designs using SPICE circuit simulation traces

An extension to a formal verification approach of hybrid systems is proposed to verify analog and mixed signal (AMS) designs. AMS designs can be formally modeled as hybrid systems and therefore lend themselves to the formal analysis and verification techniques applied to hybrid systems. The proposed approach employs simulation traces obtained from an actual design implementation of AMS circuit blocks (for example, in the form of SPICE netlists) to carry out formal analysis and verification. This enables the same platform used for formally validating an abstract model of an AMS design, to be also used for validating its different refinements and design implementation; thereby, providing a simple route to formal verification at different levels of implementation. The feasibility of the proposed approach is demonstrated with a case study based on a tunnel diode oscillator. Since the device characteristic of a tunnel diode is highly non-linear with a negative resistance region, dynamic behavior of circuits in which it is employed as an element is difficult to model, analyze and verify within a general hybrid system formal verification tool. In the case study presented the formal model and the proposed computational techniques have been incorporated into CheckMate, a formal verification tool based on MATLAB and Simulink-Stateflow Framework from MathWorks.

A rational approach to model the behavior of bonded soils

The paper presents a rational approach to model the behavior of bonded soils within the frame work of hardening plasticity. The approach is based on the premise that the resistance of bonded materials is a superposition of the two components of cement bond strength and soil frictional strength and that the deformation of the soil is associated with the frictional component of stresses just as in the case of a remoulded soil, the bonds offering additional resistance at any given strain level. This concept is similar to two stiffnesses acting in parallel for the same strain response. The proposed model considers the constitutive laws separately for the two components (bond and frictional) and adds the two to get the overall response. The unbonded soil component is described by the well known 'modified Cam clay' model. The response of the bond component is also described by a strain softening elasto-plastic model, considering the behavior to be elastic up to the yield surface and elasto-plastic beyond yield surface. To illustrate the capability of the proposed, model some laboratory test results of both compression and-extension shear tests are predicted. Despite the model being simple, several typical features of the behavior of bonded materials are well reproduced. The model parameters are well defined and easily determinable.

A Revisit to Capture the Entire Behavior of Ultrasonic Wave Dispersion Characteristics of Single Walled Carbon Nanotubes Based on Nonlocal Elasticity Theory and Flugge Shell Model

In this paper, an ultrasonic wave propagation analysis in single-walled carbon nanotube (SWCNT) is re-studied using nonlocal elasticity theory, to capture the whole behaviour. The SWCNT is modeled using Flugge's shell theory, with the wall having axial, circumferential and radial degrees of freedom and also including small scale effects. Nonlocal governing equations for this system are derived and wave propagation analysis is also carried out. The revisited nonlocal elasticity calculation shows that the wavenumber tends to infinite at certain frequencies and the corresponding wave velocity tends to zero at those frequencies indicating localization and stationary behavior. This frequency is termed as escape frequency. This behavior is observed only for axial and radial waves in SWCNT. It has been shown that the circumferential waves will propagate dispersively at higher frequencies in nonlocality. The magnitudes of wave velocities of circumferential waves are smaller in nonlocal elasticity as compared to local elasticity. We also show that the explicit expressions of cut-off frequency depend on the nonlocal scaling parameter and the axial wavenumber. The effect of axial wavenumber on the ultrasonic wave behavior in SWCNTs is also discussed. The present results are compared with the corresponding results (for first mode) obtained from ab initio and 3-D elastodynamic continuum models. The acoustic phonon dispersion relation predicted by the present model is in good agreement with that obtained from literature. The results are new and can provide useful guidance for the study and design of the next generation of nanodevices that make use of the wave propagation properties of single-walled carbon nanotubes.

Extracting Three Dimensional Surface Model of Human Kidney from the Visible Human Data Set using Free Software

Three dimensional digital model of a representative human kidney is needed for a surgical simulator that is capable of simulating a laparoscopic surgery involving kidney. Buying a three dimensional computer model of a representative human kidney, or reconstructing a human kidney from an image sequence using commercial software, both involve (sometimes significant amount of) money. In this paper, author has shown that one can obtain a three dimensional surface model of human kidney by making use of images from the Visible Human Data Set and a few free software packages (ImageJ, ITK-SNAP, and MeshLab in particular). Images from the Visible Human Data Set, and the software packages used here, both do not cost anything. Hence, the practice of extracting the geometry of a representative human kidney for free, as illustrated in the present work, could be a free alternative to the use of expensive commercial software or to the purchase of a digital model.