A Study of Early Afterdepolarizations in a Model for Human Ventricular Tissue

Sudden cardiac death is often caused by cardiac arrhythmias. Recently, special attention has been given to a certain arrhythmogenic condition, the long-QT syndrome, which occurs as a result of genetic mutations or drug toxicity. The underlying mechanisms of arrhythmias, caused by the long-QT syndrome, are not fully understood. However, arrhythmias are often connected to special excitations of cardiac cells, called early afterdepolarizations (EADs), which are depolarizations during the repolarizing phase of the action potential. So far, EADs have been studied mainly in isolated cardiac cells. However, the question on how EADs at the single-cell level can result in fibrillation at the tissue level, especially in human cell models, has not been widely studied yet. In this paper, we study wave patterns that result from single-cell EAD dynamics in a mathematical model for human ventricular cardiac tissue. We induce EADs by modeling experimental conditions which have been shown to evoke EADs at a single-cell level: by an increase of L-type Ca currents and a decrease of the delayed rectifier potassium currents. We show that, at the tissue level and depending on these parameters, three types of abnormal wave patterns emerge. We classify them into two types of spiral fibrillation and one type of oscillatory dynamics. Moreover, we find that the emergent wave patterns can be driven by calcium or sodium currents and we find phase waves in the oscillatory excitation regime. From our simulations we predict that arrhythmias caused by EADs can occur during normal wave propagation and do not require tissue heterogeneities. Experimental verification of our results is possible for experiments at the cell-culture level, where EADs can be induced by an increase of the L-type calcium conductance and by the application of I-Kr blockers, and the properties of the emergent patterns can be studied by optical mapping of the voltage and calcium.

Quantum Phase Diagram of One-Dimensional Spin and Hubbard Models with Transitions to Bond Order Wave Phases

Similar quantum phase diagrams and transitions are found for three classes of one-dimensional models with equally spaced sites, singlet ground states (GS), inversion symmetry at sites and a bond order wave (BOW) phase in some sectors. The models are frustrated spin-1/2 chains with variable range exchange, half-filled Hubbard models with spin-independent interactions and modified Hubbard models with site energies for describing organic charge transfer salts. In some range of parameters, the models have a first order quantum transition at which the GS expectation value of the sublattice spin < S-A(2)> of odd or even-numbered sites is discontinuous. There is an intermediate BOW phase for other model parameters that lead to two continuous quantum transitions with continuous < S-A(2)>. Exact diagonalization of finite systems and symmetry arguments provide a unified picture of familiar 1D models that have appeared separately in widely different contexts.

Quantum Phase Diagram of One-Dimensional Spin and Hubbard Models with Transitions to Bond Order Wave Phases

Similar quantum phase diagrams and transitions are found for three classes of one-dimensional models with equally spaced sites, singlet ground states (GS), inversion symmetry at sites and a bond order wave (BOW) phase in some sectors. The models are frustrated spin-1/2 chains with variable range exchange, half-filled Hubbard models with spin-independent interactions and modified Hubbard models with site energies for describing organic charge transfer salts. In some range of parameters, the models have a first order quantum transition at which the GS expectation value of the sublattice spin < S-A(2)> of odd or even-numbered sites is discontinuous. There is an intermediate BOW phase for other model parameters that lead to two continuous quantum transitions with continuous < S-A(2)>. Exact diagonalization of finite systems and symmetry arguments provide a unified picture of familiar 1D models that have appeared separately in widely different contexts.

Supersolid in a one-dimensional model of hard-core bosons

We study a system of hard-core boson on a one-dimensional lattice with frustrated next-nearest-neighbor hopping and nearest-neighbor interaction. At half filling, for equal magnitude of nearest- and next-nearest-neighbor hopping, the ground state of this system exhibits a first-order phase transition from a bond-ordered solid to a charge-density-wave solid as a function of the nearest- neighbor interaction. Moving away from half filling we investigate the system at incommensurate densities, where we find a supersolid phase which has concurrent off-diagonal long-range order and density-wave order which is unusual in a system of hard-core bosons in one dimension. Using the finite-size density-matrix renormalization group method, we obtain the complete phase diagram for this model.

Nonlinear gravity-wave interactions in stratified turbulence

To investigate the dynamics of gravity waves in stratified Boussinesq flows, a model is derived that consists of all three-gravity-wave-mode interactions (the GGG model), excluding interactions involving the vortical mode. The GGG model is a natural extension of weak turbulence theory that accounts for exact three-gravity-wave resonances. The model is examined numerically by means of random, large-scale, high-frequency forcing. An immediate observation is a robust growth of the so-called vertically sheared horizontal flow (VSHF). In addition, there is a forward transfer of energy and equilibration of the nonzero-frequency (sometimes called ``fast'') gravity-wave modes. These results show that gravity-wave-mode interactions by themselves are capable of systematic interscale energy transfer in a stratified fluid. Comparing numerical simulations of the GGG model and the full Boussinesq system, for the range of Froude numbers (Fr) considered (0.05 a parts per thousand currency sign Fr a parts per thousand currency sign 1), in both systems the VSHF is hardest to resolve. When adequately resolved, VSHF growth is more vigorous in the GGG model. Furthermore, a VSHF is observed to form in milder stratification scenarios in the GGG model than the full Boussinesq system. Finally, fully three-dimensional nonzero-frequency gravity-wave modes equilibrate in both systems and their scaling with vertical wavenumber follows similar power-laws. The slopes of the power-laws obtained depend on Fr and approach -2 (from above) at Fr = 0.05, which is the strongest stratification that can be properly resolved with our computational resources.

Nonequilibrium ``Melting'' of a Charge Density Wave Insulator via an Ultrafast Laser Pulse

We employ an exact solution of the simplest model for pump-probe time-resolved photoemission spectroscopy in charge-density-wave systems to show how, in nonequilibrium, the gap in the density of states disappears while the charge density remains modulated, and then the gap reforms after the pulse has passed. This nonequilibrium scenario qualitatively describes the common short-time experimental features in TaS2 and TbTe3, indicating a quasiuniversality for nonequilibrium ``melting'' with qualitative features that can be easily understood within a simple picture.

Ultrasonic guided wave sensing properties of PVDF thin film with inter digital electrodes

Ultrasonic strain sensing performance of the large area PVDF with Inter Digital Electrodes (IDE) is studied in this work. Procedure to obtain IDE on a beta-phase PVDF is explained. PVDF film with IDE is bonded on a plate structure and is characterized for its directional sensitivity at different frequencies. Guided waves are induced on the IDE-PVDF sensor from different directions by placing a piezoelectric wafer actuator at different angles. Strain induced on the IDE-PVDF sensor by the guided waves in estimated by using a Laser Doppler Vibrometer (LDV) and a wave propagation model. Using measured voltage response from IDE-PVDF sensor and the strain measurements from LDV the piezoelectric coefficient is estimated in various directions. The variation of 11 e at different angles shows directional sensitivity of the IDE-PVDF sensor to the incident guided waves. The present study provides an effective technique to characterize thin film piezoelectric sensors for ultrasonic strain sensing at very high frequencies of 200 kHz. Often frequency of the guided wave is changed to alter the wavelength to interrogate damages of different sizes in Structural Health Monitoring (SHM) applications. The unique property of directional sensitivity combined with frequency tunability makes the IDE-PVDF sensor most suitable for SHM of structures.

Spectral finite element based on an efficient layerwise theory for wave propagation analysis of composite and sandwich beams

In this paper, we present a spectral finite element model (SFEM) using an efficient and accurate layerwise (zigzag) theory, which is applicable for wave propagation analysis of highly inhomogeneous laminated composite and sandwich beams. The theory assumes a layerwise linear variation superimposed with a global third-order variation across the thickness for the axial displacement. The conditions of zero transverse shear stress at the top and bottom and its continuity at the layer interfaces are subsequently enforced to make the number of primary unknowns independent of the number of layers, thereby making the theory as efficient as the first-order shear deformation theory (FSDT). The spectral element developed is validated by comparing the present results with those available in the literature. A comparison of the natural frequencies of simply supported composite and sandwich beams obtained by the present spectral element with the exact two-dimensional elasticity and FSDT solutions reveals that the FSDT yields highly inaccurate results for the inhomogeneous sandwich beams and thick composite beams, whereas the present element based on the zigzag theory agrees very well with the exact elasticity solution for both thick and thin, composite and sandwich beams. A significant deviation in the dispersion relations obtained using the accurate zigzag theory and the FSDT is also observed for composite beams at high frequencies. It is shown that the pure shear rotation mode remains always evanescent, contrary to what has been reported earlier. The SFEM is subsequently used to study wavenumber dispersion, free vibration and wave propagation time history in soft-core sandwich beams with composite faces for the first time in the literature. (C) 2014 Elsevier Ltd. All rights reserved.

Wave propagation analysis in laminated composite plates with transverse cracks using the wavelet spectral finite element method

This paper presents a newly developed wavelet spectral finite element (WFSE) model to analyze wave propagation in anisotropic composite laminate with a transverse surface crack penetrating part-through the thickness. The WSFE formulation of the composite laminate, which is based on the first-order shear deformation theory, produces accurate and computationally efficient results for high frequency wave motion. Transverse crack is modeled in wavenumber-frequency domain by introducing bending flexibility of the plate along crack edge. Results for tone burst and impulse excitations show excellent agreement with conventional finite element analysis in Abaqus (R). Problems with multiple cracks are modeled by assembling a number of spectral elements with cracks in frequency-wavenumber domain. Results show partial reflection of the excited wave due to crack at time instances consistent with crack locations. (C) 2014 Elsevier B.V. All rights reserved.

A spectral multiscale method for wave propagation analysis: Atomistic-continuum coupled simulation

In this paper, we present a new multiscale method which is capable of coupling atomistic and continuum domains for high frequency wave propagation analysis. The problem of non-physical wave reflection, which occurs due to the change in system description across the interface between two scales, can be satisfactorily overcome by the proposed method. We propose an efficient spectral domain decomposition of the total fine scale displacement along with a potent macroscale equation in the Laplace domain to eliminate the spurious interfacial reflection. We use Laplace transform based spectral finite element method to model the macroscale, which provides the optimum approximations for required dynamic responses of the outer atoms of the simulated microscale region very accurately. This new method shows excellent agreement between the proposed multiscale model and the full molecular dynamics (MD) results. Numerical experiments of wave propagation in a 1D harmonic lattice, a 1D lattice with Lennard-Jones potential, a 2D square Bravais lattice, and a 2D triangular lattice with microcrack demonstrate the accuracy and the robustness of the method. In addition, under certain conditions, this method can simulate complex dynamics of crystalline solids involving different spatial and/or temporal scales with sufficient accuracy and efficiency. (C) 2014 Elsevier B.V. All rights reserved.

Dynamic Site Characterization and Correlation of Shear Wave Velocity with Standard Penetration Test `N' Values for the City of Agartala, Tripura State, India

Seismic site characterization is the basic requirement for seismic microzonation and site response studies of an area. Site characterization helps to gauge the average dynamic properties of soil deposits and thus helps to evaluate the surface level response. This paper presents a seismic site characterization of Agartala city, the capital of Tripura state, in the northeast of India. Seismically, Agartala city is situated in the Bengal Basin zone which is classified as a highly active seismic zone, assigned by Indian seismic code BIS-1893, Indian Standard Criteria for Earthquake Resistant Design of Structures, Part-1 General Provisions and Buildings. According to the Bureau of Indian Standards, New Delhi (2002), it is the highest seismic level (zone-V) in the country. The city is very close to the Sylhet fault (Bangladesh) where two major earthquakes (M (w) > 7) have occurred in the past and affected severely this city and the whole of northeast India. In order to perform site response evaluation, a series of geophysical tests at 27 locations were conducted using the multichannel analysis of surface waves (MASW) technique, which is an advanced method for obtaining shear wave velocity (V (s)) profiles from in situ measurements. Similarly, standard penetration test (SPT-N) bore log data sets have been obtained from the Urban Development Department, Govt. of Tripura. In the collected data sets, out of 50 bore logs, 27 were selected which are close to the MASW test locations and used for further study. Both the data sets (V (s) profiles with depth and SPT-N bore log profiles) have been used to calculate the average shear wave velocity (V (s)30) and average SPT-N values for the upper 30 m depth of the subsurface soil profiles. These were used for site classification of the study area recommended by the National Earthquake Hazard Reduction Program (NEHRP) manual. The average V (s)30 and SPT-N classified the study area as seismic site class D and E categories, indicating that the city is susceptible to site effects and liquefaction. Further, the different data set combinations between V (s) and SPT-N (corrected and uncorrected) values have been used to develop site-specific correlation equations by statistical regression, as `V (s)' is a function of SPT-N value (corrected and uncorrected), considered with or without depth. However, after considering the data set pairs, a probabilistic approach has also been presented to develop a correlation using a quantile-quantile (Q-Q) plot. A comparison has also been made with the well known published correlations (for all soils) available in the literature. The present correlations closely agree with the other equations, but, comparatively, the correlation of shear wave velocity with the variation of depth and uncorrected SPT-N values provides a more suitable predicting model. Also the Q-Q plot agrees with all the other equations. In the absence of in situ measurements, the present correlations could be used to measure V (s) profiles of the study area for site response studies.

Comparative study of different nonconserving time integrators for wave propagation in hyperelastic waveguides

This paper addresses the formulation and numerical efficiency of various numerical models of different nonconserving time integrators for studying wave propagation in nonlinear hyperelastic waveguides. The study includes different nonlinear finite element formulations based on standard Galerkin finite element model, time domain spectral finite element model, Taylor-Galerkin finite element model, generalized Galerkin finite element model and frequency domain spectral finite element model. A comparative study on the computational efficiency of these different models is made using a hyperelastic rod model, and the optimal computational scheme is identified. The identified scheme is then used to study the propagation of transverse and longitudinal waves in a Timoshenko beam with Murnaghan material nonlinearity.

Turbulent electrical activity at sharp-edged inexcitable obstacles in a model for human cardiac tissue

Wave propagation around various geometric expansions, structures, and obstacles in cardiac tissue may result in the formation of unidirectional block of wave propagation and the onset of reentrant arrhythmias in the heart. Therefore, we investigated the conditions under which reentrant spiral waves can be generated by high-frequency stimulation at sharp-edged obstacles in the ten Tusscher-Noble-Noble-Panfilov (TNNP) ionic model for human cardiac tissue. We show that, in a large range of parameters that account for the conductance of major inward and outward ionic currents of the model fast inward Na+ current (INa), L-type slow inward Ca2+ current (I-CaL), slow delayed-rectifier current (I-Ks), rapid delayed-rectifier current (I-Kr), inward rectifier K+ current (I-K1)], the critical period necessary for spiral formation is close to the period of a spiral wave rotating in the same tissue. We also show that there is a minimal size of the obstacle for which formation of spirals is possible; this size is similar to 2.5 cm and decreases with a decrease in the excitability of cardiac tissue. We show that other factors, such as the obstacle thickness and direction of wave propagation in relation to the obstacle, are of secondary importance and affect the conditions for spiral wave initiation only slightly. We also perform studies for obstacle shapes derived from experimental measurements of infarction scars and show that the formation of spiral waves there is facilitated by tissue remodeling around it. Overall, we demonstrate that the formation of reentrant sources around inexcitable obstacles is a potential mechanism for the onset of cardiac arrhythmias in the presence of a fast heart rate.

Spiral-wave dynamics in ionically realistic mathematical models for human ventricular tissue: the effects of periodic deformation

We carry out an extensive numerical study of the dynamics of spiral waves of electrical activation, in the presence of periodic deformation (PD) in two-dimensional simulation domains, in the biophysically realistic mathematical models of human ventricular tissue due to (a) ten-Tusscher and Panfilov (the TP06 model) and (b) ten-Tusscher, Noble, Noble, and Panfilov (the TNNPO4 model). We first consider simulations in cable-type domains, in which we calculate the conduction velocity theta and the wavelength lambda of a plane wave; we show that PD leads to a periodic, spatial modulation of theta and a temporally periodic modulation of lambda; both these modulations depend on the amplitude and frequency of the PD. We then examine three types of initial conditions for both TP06 and TNNPO4 models and show that the imposition of PD leads to a rich variety of spatiotemporal patterns in the transmembrane potential including states with a single rotating spiral (RS) wave, a spiral-turbulence (ST) state with a single meandering spiral, an ST state with multiple broken spirals, and a state SA in which all spirals are absorbed at the boundaries of our simulation domain. We find, for both TP06 and TNNPO4 models, that spiral-wave dynamics depends sensitively on the amplitude and frequency of PD and the initial condition. We examine how these different types of spiral-wave states can be eliminated in the presence of PD by the application of low-amplitude pulses by square- and rectangular-mesh suppression techniques. We suggest specific experiments that can test the results of our simulations.

Wave propagation analysis in adhesively bonded composite joints using the wavelet spectral finite element method

A wavelet spectral finite element (WSFE) model is developed for studying transient dynamics and wave propagation in adhesively bonded composite joints. The adherands are formulated as shear deformable beams using the first order shear deformation theory (FSDT) to obtain accurate results for high frequency wave propagation. Equations of motion governing wave motion in the bonded beams are derived using Hamilton's principle. The adhesive layer is modeled as a line of continuously distributed tension/compression and shear springs. Daubechies compactly supported wavelet scaling functions are used to transform the governing partial differential equations from time domain to frequency domain. The dynamic stiffness matrix is derived under the spectral finite element framework relating the nodal forces and displacements in the transformed frequency domain. Time domain results for wave propagation in a lap joint are validated with conventional finite element simulations using Abaqus. Frequency domain spectrum and dispersion relation results are presented and discussed. The developed WSFE model yields efficient and accurate analysis of wave propagation in adhesively-bonded composite joints. (C) 2014 Elsevier Ltd. All rights reserved.