AJP - Heart AJP: Lung Cellular and Molecular Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Heart Circ Physiol 277: H2065-H2082, 1999;
0363-6135/99 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Idriss, S. F.
Right arrow Articles by Ideker, R. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Idriss, S. F.
Right arrow Articles by Ideker, R. E.
Vol. 277, Issue 5, H2065-H2082, November 1999

Effect of pacing site on ventricular fibrillation initiation by shocks during the vulnerable period

Salim F. Idriss1,2, Patrick D. Wolf2, William M. Smith3, and Raymond E. Ideker3

1 Department of Pathology and School of Medicine, Duke University Medical Center, Durham 27710; 2 Department of Biomedical Engineering and the Engineering Research Center for Emerging Cardiovascular Technologies, School of Engineering, Duke University, Durham, North Carolina 27706; and 3 Departments of Medicine and Physiology, University of Alabama at Birmingham, Birmingham, Alabama 35233


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

The critical point hypothesis for the upper limit of vulnerability (ULV) states that the site of S1 pacing should not affect the ULV S2 shock strength for a single S2 shock electrode configuration but may affect the S1-S2 interval at which sub-ULV shocks induce ventricular fibrillation (VF). Furthermore, early post-S2 activations leading to VF should arise in areas with low potential gradients of similar magnitude, regardless of the S1 site. This hypothesis was tested in 10 pigs by determining ULVs for three S1 sites [left ventricular apex (LVA), LV base (LVB), and right ventricular outflow tract (RVOT)] with one S2 configuration (LVA patch to superior vena cava catheter). T-wave scanning was performed with biphasic S2 shocks incremented from 60 V in 40-V steps and stepped up or down in 20- and 10-V steps. Activations and S2 potential gradients were recorded at 528 epicardial sites. Although shocks just below the ULV induced VF significantly earlier in the T wave when the S1 site was the RVOT than when it was the LVA or LVB, ULVs were not significantly different for the three S1 pacing sites. Early post-S2 activations arose closer to the S2 electrode for weak S2s but moved to distant low potential gradient areas as the S2 strengthened. Just below the ULV, early post-S2 activations arose in the RVOT when the S1 site was the LVA or LVB but arose along the RV base when the S1 site was the RVOT. Early site potential gradients were not significantly different just below the ULV (LVA: 8.2 ± 4.1 V/cm; LVB: 8.6 ± 4.9 V/cm; RVOT: 8.7 ± 4.4 V/cm). At the ULV, early post-S2 activations arose from the same areas but did not induce VF. The results support the critical point hypothesis for the ULV. For this S2 configuration, no single point in the T wave could be used to determine the ULV for all S1 sites.

epicardial mapping; biphasic waveform; reentry; stimulation


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

A SINGLE SHOCK delivered during a portion of the T wave can induce ventricular fibrillation (VF) (17, 35). Thus ventricular vulnerability is bounded in terms of the timing of premature shocks that will induce VF. Upper and lower boundaries also exist for the strength of the shock. A shock must exceed the VF threshold to induce VF (34). If a shock is so strong that it exceeds the upper limit of ventricular vulnerability (ULV), then VF is not induced, regardless of the timing of the shock (9). This ULV exists in animals (9, 15, 28) and humans (2, 13).

Winfree (37) proposed a mechanism that explains both the induction of VF by a strong premature stimulus and the existence of the ULV. The mechanism, called the critical point mechanism by others (11), describes how a premature shock (S2) induces reentry in excitable tissue. According to this theory, the application of an S2 shock creates varying tissue effects that depend on the level of refractoriness and the strength of the S2 shock field at each tissue site. Reentry will form around a tissue site where a spatially varying S2 field encompasses a particular value, called the critical field strength, and where the tissue refractoriness also varies spatially and encompasses a particular level of recovery, called the critical recovery level. This tissue site, called the critical point, occurs where the critical level of the S2 potential gradient field intersects tissue that is critically refractory. Reentry should initiate wherever a critical point is formed as long as the surrounding area of myocardium is large enough to support the reentrant circuit.

The critical point hypothesis for the ULV states that S2 shocks stronger than the VF threshold but weaker than the ULV induce reentry, and subsequent VF, via the critical point mechanism (28, 36). Because the critical potential gradient areas are moved farther away from the S2 shocking electrodes as S2 strength is increased, the hypothesis predicts that the locations where VF is initiated will change with increasing shock strength. At or above a particular S2 shock strength, reentry will not occur because the critical point will be positioned very close to or beyond a heart border. This S2 shock strength is the ULV (36, 37).

Although the critical point mechanism for the induction of reentry and VF has been experimentally validated (11), the critical point hypothesis for the ULV has not. Shibata et al. (28) mapped the initiation of VF by shocks using 56 bipolar epicardial electrodes. Although they showed that increases in S2 shock strength move VF initiation sites away from the shocking electrodes until the ULV is reached, they did not map the S2 potential gradient field. A detailed study mapping the earliest post-S2 cardiac activations and the S2 shock potential gradient fields for shocks at and below the ULV has not previously been reported.

The purpose of this study was to examine the critical point hypothesis for the ULV by testing several predictions. The hypothesis predicts that, for a single S2 shock electrode configuration, the direction of activation and recovery preceding the shock (i.e., the S1 pacing site) should not affect the ULV shock strength but should affect the time points in the T wave at which VF is induced. The hypothesis predicts that the sites of earliest activation following S2 shocks should move away from the S2 electrodes as shock strength is increased. For shocks just below the ULV, the earliest postshock activation should appear where the potential gradient field is weakest in the ventricles. The potential gradients measured at the initiation points should be similar, regardless of location or S1 pacing site. In addition, the level of refractoriness at the time of the S2 should be similar at the initiation points, regardless of S1 pacing site. This study was accomplished by determining the ULV for three spatially distinct S1 pacing sites and a single S2 shock electrode configuration while mapping epicardial activations at 528 sites and shock potential gradients at 264 sites.


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Animal Preparation

All studies were conducted after protocol approval by the Institutional Animal Care and Use Committee. Ten pigs weighing 25-35 kg were premedicated with ketamine (22 mg/kg) and acepromazine (1.1 mg/kg). Anesthesia was induced with an intravenous bolus of pentobarbital sodium (10 mg/kg) and maintained with a continuous pentobarbital drip (0.05 mg · kg-1 · min-1). Morphine sulfate was used for analgesia (0.2 mg/kg im loading, followed by 0.1 mg/kg im every 2 h). The animals were intubated with a cuffed endotracheal tube and ventilated with a Harvard respirator (Harvard Apparatus, South Natick, MA). Central venous access was established by cut-down to the external jugular vein. A triple lumen catheter was inserted for infusion of anesthetic and fluids (0.9% saline at ~5-6 ml · kg-1 · h-1). A femoral or carotid arterial line was established for continuous hemodynamic monitoring and blood sampling. Arterial blood was sampled every 30-60 min, and blood gases and electrolyte levels were determined. Potassium chloride, sodium bicarbonate, and calcium chloride were given as indicated to maintain normal electrolyte values. Skeletal muscle paralysis was induced with intravenous succinylcholine (1 mg/kg) and maintained with a dosage of 0.25-0.50 mg/kg every hour as needed. Core body temperature was monitored and maintained at 37°C with a heating blanket and/or heat lamp.

The heart was exposed by median sternotomy, and a pericardial cradle was created. A stainless steel mesh defibrillation electrode (~3-cm2 area) was sewn to the epicardial surface of the left ventricular apex (LVA). The electrode was a cloverleaf shape with four leaflets to conform to the apical surface. A single-coil catheter defibrillation electrode (~3.5-cm length, ~3.5-mm diameter) was inserted into an internal or external jugular vein and positioned in the distal superior vena cava (SVC) with its tip at the entrance to the right atrium. Three stainless steel wires were sewn to the epicardium for S1 pacing: one near the LVA, one overlying the right ventricular outflow tract (RVOT) close to the pulmonic valve ring, and one posteriorly at the LV base (LVB). The return electrode for the S1 electrodes was connected to the chest retractor. The S1 and S2 electrode positions are diagrammed in Fig. 1. Two chlorided solid Ag beaded wires were inserted into the right and left thighs through small skin incisions to act as grounds for the mapping system.


View larger version (63K):
[in this window]
[in a new window]
 
Fig. 1.   Diagram showing anterior and posterior views of the heart with approximate positions of S1 and S2 electrodes and sock recording electrodes. S1 pacing was performed from each of the epicardial S1 positions referenced to the chest retractor. S2 shocks were delivered between left ventricular apex (LVA) patch electrode and superior vena cava (SVC) catheter electrode. , 1 of 528 unipolar recording sites. LAD, left anterior descending coronary artery, PDA, posterior descending coronary artery; RVOT, right ventricular outflow tract; LVB, left ventricular base.

S1 Pacing and S2 Shock Delivery

A computer-controlled programmable stimulator was used to deliver a train of S1 pacing stimuli followed by a pulse to trigger S2 shock delivery from a defibrillator. The S1 stimulus was a cathodal pulse set to twice the diastolic threshold for each individual S1 pacing site. The S1 pacing train consisted of 20 S1 stimuli with the S1-S1 interval set to 80% of the intrinsic R-R interval. The same S1-S1 interval was used for each pacing site. After the 20th S1 stimulus, the S2 shock was delivered after a specified delay (the S1-S2 coupling interval).

All S2 and defibrillation shocks were delivered from an external defibrillator (HVS-02, Ventritex, Sunnyvale, CA). A biphasic truncated exponential waveform was used in which both phases were 6 ms in duration. The LVA patch electrode was the first-phase anode. To mimic a single capacitor discharge, the leading-edge voltage of the second phase was set equal to the trailing edge of the first phase. For each shock, the voltage and current waveforms were acquired at a 20-kHz sampling rate with a waveform analyzer (model 6100, Data Precision, Danvers, MA). Shock parameters were then measured (voltage and current) or calculated (energy, impedance, and time constant) and recorded.

Epicardial Recording Electrodes

Forty-four multicontact electrode plaques were manually constructed (14). Each plaque had a solder layer, a circuit-board backing layer, and a recording layer (Fig. 2A). The recording layer was formed by embedding 24 solid Ag wires (254-µm diameter, California Fine Wire, Grover City, CA) in casting resin (Castin'craft, ETI, Fields Landing, CA). The wires were positioned as a 4 × 6 rectangular grid with 2.54-mm center-center spacing. Before each study, the recording layer surface was polished. Each Ag wire tip was then electrolytically chlorided (12, 14) to form an Ag-AgCl recording electrode.


View larger version (47K):
[in this window]
[in a new window]
 
Fig. 2.   A: top, side, and bottom views of a multicontact electrode plaque. Side view shows the 3 layers of the plaque: solder layer, circuit-board layer, and recording layer with embedded Ag wire recording electrodes. B: recording face of a plaque with 24 Ag-AgCl recording sites. Exploded view shows group of 4 sites oriented as a square with 2.54-mm side lengths. Within each square, cardiac activations were measured at 2 sites in unipolar mode from + electrodes. Bipolar recordings were made during shocks between + and - electrodes across each diagonal for potential gradient calculations.

An elastic recording sock (Xspan, Baxter Healthcare, Deerfield, IL) was constructed to which the 44 multicontact electrode plaques were secured with small snaps. The recording sock was stretched to cover both ventricles with the plaques against the epicardial surface. Multiple sutures were placed along the atrioventricular (AV) groove to hold the sock in place. At the end of each study, the plaques were removed from the sock and cleaned. A new sock was used for each animal study. Care was taken to maintain consistent plaque positioning and orientation across animals.

Recording Methods

Each multicontact plaque contained six squares of four recording electrodes (Fig. 2B). Unipolar cardiac activations were recorded using two electrodes per square referenced to the left leg for a total of 528 sites. During S2 shocks, bipolar recordings were made with electrodes on opposing corners of the square. Therefore, the recorded shock waveforms represented the potential differences across each diagonal. The potential gradient across each diagonal was calculated by dividing the peak potential difference, measured at the leading edge of the S2 shock (16), by the measured distance between recording electrodes. The two-dimensional epicardial potential gradient magnitude at each of the 264 squares was calculated as the square root of the sum of the squared diagonal potential gradients (23). The component of the potential gradient magnitude directed normal to the epicardium was not measured.

All electrode recordings were made at a sampling rate of 2 kHz using a 528-channel mapping system (38). Unipolar cardiac activations were recorded until 10 ms before shock delivery. The mapping system was then switched to bipolar recording until 5 ms after the end of the shock, when unipolar recording was resumed. Unipolar cardiac signals were band-pass filtered from 0.05 to 500 Hz except in the 50-ms window after a shock. During that postshock period, the high-pass filter was set to 5 Hz to rapidly reduce any postshock offset and to return the signal to baseline. The bipolar signals recorded during the shock were band-pass filtered from direct current to 4 kHz.

The gain settings for recording shocks were determined at the beginning of each study by delivering a 180- and a 380-V S2 shock after S1 pacing trains with the S1-S2 coupling interval set to the S1-S1 interval to prevent VF. The S2 peak voltages were determined for each channel for each shock (16). Optimal gains for any S2 strength were estimated from a linear equation relating peak voltage and S2 strength for each channel.

ULV Testing

ULV testing was performed by systematically varying the S1-S2 coupling interval to "scan" the T wave with S2 shocks. The S1-S2 coupling intervals were determined at the beginning of the study individually for the three S1 pacing sites. To determine the scanning interval, the last S1 stimulus and its ventricular response were recorded on the waveform analyzer using electrocardiogram (ECG) signals recorded through the catheter-patch S2 electrodes. This was done because the surface ECG signals were of poor quality due to the median sternotomy and the presence of the electrode sock. The time from the beginning of the S1 stimulus to the beginning (S1-Tb), peak (S1-Tpk), and end (S1-Te) of the T wave was then measured. This sequence was repeated four more times. The shortest S1-S2 interval used during scanning was the average of the five S1-Tb times minus 10 ms. Similarly, the longest S1-S2 coupling interval used was the average of the five S1-Te times plus 10 ms. These intervals were rechecked after every six VF inductions during ULV testing.

Initially, the S2 shock was set to 60 V of leading-edge voltage (LEV) and the S1-S2 interval was set to the S1-Tpk interval (Fig. 3). If the S2 test shock did not induce VF, the S1-S2 interval was incremented by 10 ms (S1-Tpk + 10 ms). If the S2 shock again did not induce VF, the S1-S2 interval was decremented by 20 ms (S1-Tpk - 10 ms). Testing continued by alternating between scanning later and earlier than Tpk in 10-ms steps until VF was induced or until the minimum and maximum scanning intervals were reached. When either the maximum or minimum scanning interval was obtained and more than one step was needed to obtain the other limit, successive incrementing or decrementing of the S1-S2 interval was used to complete the T-wave scan. At least 15 s were allowed to elapse between S2 delivery and restarting the S1 train when VF was not induced (33). If VF was induced, the animal was defibrillated with one or more rescue shocks (equal to or 50-400 V higher than the S2) delivered through the S2 electrodes. After a recovery period of at least 4 min following defibrillation, testing was resumed.


View larger version (19K):
[in this window]
[in a new window]
 
Fig. 3.   Scanning protocol for upper limit of variability (ULV) testing. A portion of T wave is shown at bottom. In this case, testing is started with S2 = 60 V. The first S2 shock (shock 1) is timed to occur at peak of T wave (Tpk). Ventricular fibrillation (VF) was not induced, and the S1-S2 interval was incremented by 10 ms for the next shock (shock 2). VF was not induced and S1-S2 was decremented by 20 ms (3 ). After shock 3 failed to induce VF, S1-S2 was incremented by 30 ms. VF was induced when shock 4 was delivered. After defibrillation, S2 was set to 100 V and scanning was restarted (top). The first shock at the new S2 strength (shock 5) was timed to intersect the T wave where VF was previously induced. Back-and-forth scanning was continued from this point until VF was again induced (shock 7).

If VF was not induced at any S1-S2 coupling interval, the S2 shock strength was decremented by 40 V, and the scanning sequence was restarted. If VF was induced, testing at that S2 strength was halted, the S2 strength was incremented by 40 V, and the scanning sequence was restarted. The first S1-S2 coupling interval to be tested at the new S2 strength was the interval that first induced VF at the previous S2 shock strength (Fig. 3). The ULV was determined when a single increment or decrement of the S2 shock strength changed the result of the T-wave scan compared with the preceding scan. For example, ULV testing was halted when VF was induced at one S2 shock strength but not at the next higher S2 shock strength. The ULV was considered the lowest S2 shock strength that did not induce VF at any point in the T-wave scanning interval.

The ULV was determined for all three S1 pacing sites. The ULV was measured three times for each S1 site, first using a 40-V step size as described above, next with a 20-V step size, and last with a 10-V step size for the S2 shock. The starting voltage for the second and third ULV determinations was either the ULV or one step below the preceding ULV. The first ULV is called ULV40 (40-V step size), the second ULV20 (20-V step size), and the third ULV10 (10-V step size). A randomized Latin-square design was used to determine the testing order of the three S1 pacing sites for each animal. Within an animal, the S1 pacing site tested was changed after each T-wave scanning sequence. For the second and third ULV determinations, the S2 shock strength was not incremented after VF was induced, but scanning was continued to determine all S1-S2 intervals that induced VF at that S2 strength.

At the end of the study, the animal was killed with electrically induced VF. The catheter electrode position was verified, and the heart was excised, rinsed, weighed, tagged, and placed in 10% Formalin.

Data Analysis

Statistical analysis was performed using SAS software (25). A P value <0.05 was considered statistically significant.

Global shock variables. The ULVs for each S1 site were compared on the basis of the first-phase LEV and leading-edge current and on the total energy of the two phases combined. Data were analyzed by repeated-measures ANOVA to control for interanimal variability.

Timing intervals. Timing intervals measured from the S2 electrodes were compared for shocks that induced VF from the ULV40 group for the three S1 sites.

Early site locations. To determine the earliest sites of activation following the S2 shocks, cardiac activations were visualized and animated on a three-dimensional representation of the electrode sock obtained from a magnetic resonance image of the sock on a fixed pig heart (Fig. 1). All VF induction episodes that occurred during determination of ULV40 were analyzed. For ULV20 and ULV10, VF episodes for S2 shocks whose strength was one step below the ULV were analyzed. Earliest sites of postshock activation were also determined for each shock of ULV strength that was given at the same S1-S2 interval as a lower-strength shock that induced VF. Although these ULV strength shocks did not induce VF, they usually induced at least one ectopic beat.

The 528 unipolar recordings were passed through a 60-Hz filter for noise reduction (1). The first derivatives of the filtered data were calculated using a five-point algorithm (21). The derivatives were animated on the three-dimensional display to track the activation sequences. Each electrode was represented by a small square that was colored with one of eight different hues, each representing a different range of the value of the first derivative. Activation at a single electrode was represented by the first color change, which was set to occur when the derivative became more negative than approximately -0.2 to -0.5 mV/ms (3).

A single site of early activation following the S2 was considered to be present if this site was the origination point for a propagating wave front that activated the remainder of the ventricles. However, if other propagating wave fronts developed before the first wave front activated the rest of the epicardium, then multiple early sites were considered to exist. Each early site was determined by visually identifying the region where the activation wave front started. The electrograms from the electrodes in this region were then displayed one at a time, and the electrode with the earliest activation was identified.

Each early site location was specified in terms of its polar distance from the tip of the LVA and its rotational angle from the anterior midline. Each electrode position was represented by rho  (distance from tip of LVA), theta  (angle from midline in transverse plane), and phi (angle from central vertical axis). A single view of the entire sock was obtained by plotting rho  versus theta  at which the early sites were visualized and quantified (Fig. 4). Early site location data were analyzed by ANOVA.


View larger version (30K):
[in this window]
[in a new window]
 
Fig. 4.   Sock electrode positions plotted as theta  on horizontal axis and rho  along vertical axis. theta  is the rotational angle around the heart. The "line" at which theta  = 0 divides the sock into right and left halves along the anterior midline. At top, the middle of the right ventricular outflow tract (RVOT) is positioned along the theta  = 0 line. rho  is the distance from the tip of the LVA. Each square represents 1 of 528 unipolar recording sites. The approximate position of the SVC catheter is shown at top left. The LVA patch is positioned along the lower border of the sock. The approximate positions of the three S1 pacing sites are shown as filled circles. RCA, right coronary artery.

Early site activations. The local electrogram recorded at each early site was displayed from the next-to-last (19th) S1 paced beat to several beats following the S2 shock. The times of local activation following the last S1 stimulus (VR) and after the S2 shock (V1) were determined for each early site. The activation time was considered to be the time when the derivative was maximally negative. In addition, the time of the maximum positive derivative of the local recorded T wave following the last S1 stimulus was found at each early site. This time (RS1) was used as a reference point for recovery of the underlying tissue (22). On the basis of these measured time points, several intervals were calculated and compared. Early site activation data were analyzed by repeated-measures ANOVA and analysis of covariance.

Early site potential gradients. The two-dimensional potential gradient magnitudes were calculated for every shock. For shocks that were one step weaker than the ULV, the potential gradients measured at the early activation sites were compared for the three S1 pacing sites. Early site potential gradient data were analyzed by ANOVA and analysis of covariance.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

A total of 301 S2 shock episodes were analyzed. In five animals, all three ULVs were determined. In four animals, only the ULV40 was obtained. In one animal, only ULV40 and ULV20 were determined. In one animal, the mapping system gain bank used immediately after S2 shock delivery was set improperly so that the sites of early activations following shock delivery could not be determined. This animal was excluded from early site analysis. Thus, for early site determinations, 294 S2 shock episodes were analyzed.

Global Shock Variables

The global shock strength data for shocks at ULV are shown in Table 1 for each S1 site. For the ULV40 data, the individual results are tabulated. Data for ULV40, ULV20, and ULV10 are summarized at the bottom of Table 1. For most animals, the ULV40 voltage either did not differ or differed by one step in the test shock strength. There was a two-step difference in the ULV40 between S1 sites in two animals. For pigs 4 and 10, there was a three-step difference in the ULV40 voltage. However, in both animals, the ULVs converged closer with further testing (pig 4 ULV20: 380, 340, and 400 V; pig 10: 380, 420, and 400 V for LVA, LVB, and RVOT S1 sites, respectively). For ULV40, changing the S1 pacing did not significantly affect the measured ULV for voltage (P = 0.48), current (P = 0.83), or energy (P = 0.92) when repeated-measures ANOVA was used to account for interanimal variation. The same was true for ULV20 and ULV10. Also, for each S1 pacing site, ULV40, ULV20, and ULV10 were not significantly different in terms of voltage, current, or energy.

                              
View this table:
[in this window]
[in a new window]
 
Table 1.   ULV values for the three S1 sites

Timing Intervals

S1-S1 pacing intervals ranged from 280 to 500 ms (mean ± SD: 366 ± 67 ms). The S1-S1 pacing intervals remained constant for the duration of the study in all but one animal. In that case, the pacing interval was decreased after the first two sets of shocks to maintain overdrive pacing. The pacing interval then remained constant for the duration of the study. Data obtained before the change in pacing interval were not used for analysis.

The timing of VF inductions relative to the peak of the T wave (Tpk-S2 interval) as a function of the S2 shock LEV during determination of ULV40 is shown in Fig. 5. A positive Tpk-S2 interval indicates that the S2 shock was delivered after the peak. At the lowest shock level (60 V), VF was frequently induced by the first (Tpk) or second (Tpk + 10) S2 shock given. As the S2 shock strength increased, there was a clear separation of the timing of VF inductions as a function of the S1 pacing site. At one step below the ULV for the LVA and LVB pacing sites, S2 shocks induced VF at later points in the T wave than the peak (Tpk-S2 LVA: 32.4 ± 17 ms; LVB: 20.6 ± 21 ms). However, for the RVOT S1 site, S2 shocks induced VF at points significantly earlier in the T wave than the peak (Tpk-S2 -19.1 ± 19 ms).


View larger version (21K):
[in this window]
[in a new window]
 
Fig. 5.   Plot of mean Tpk-S2 interval at which VF was induced vs. S2 leading-edge voltage (LEV) for all animals combined. Tpk-S2 was measured using the electrogram recorded between the LVA patch and the SVC catheter. The Tpk-S2 interval was the time between Tpk and an S2 shock that induced VF at a specific S2 strength. The Tpk-S2 = 0 line represents the peak of the T wave. Initial S2 shocks were 60 V and were delivered at Tpk. Lines show change in Tpk-S2 interval with increasing S2 shock strength for LVA (squares), LVB (triangles), and RVOT (circles) S1 pacing sites. Filled symbols represent mean Tpk-S2 intervals; open symbols represent 1 SD above or below mean. Values without SD points represent single-animal data.

During determination of ULV20 and ULV10, full T-wave scanning was performed. VF was induced at only a single time point in 27 of the 33 full T-wave scans that were analyzed one step (20 or 10 V) below the ULV. In the remaining six T-wave scans, VF was induced at only two times in the T wave. For five of the six cases, these two VF induction points were either 10 or 20 ms apart. In the remaining case, the VF induction points were separated by 50 ms.

Early Site Locations

The mean number of early sites following a VF-inducing S2 shock was 1.9 ± 0.94 (range: 1-6). For all VF inductions, one, two, or three early sites were identified in 38%, 40%, and 17% of the episodes, respectively. The number of early sites per VF episode did not vary significantly for different S1 pacing sites (LVA: 1.9 ± 1.0; LVB: 1.8 ± 0.8; RVOT: 2.0 ± 1.1; P = 0.50, ANOVA). For S2 shocks that were weaker than 120 V below the ULV, two or three early sites were identified in most cases, regardless of S1 site. As S2 shock strength increased, the proportion of episodes with a single early site increased. At one step below the ULV (-10, -20, and -40 V), a single early site was identified for 60% of the shocks that induced VF. Two early sites accounted for 27% of the cases at these shock levels.

Examples of activation sequences for S1 paced beats from the LVA, LVB, and RVOT are shown in Fig. 6. For the LVA pacing site, activation proceeded along the anterior and posterior surfaces simultaneously. The entire epicardium was activated in 74 ms in this example, with the last portion activated in the high RVOT. For LVB S1 pacing, activation began in the high posterior base and propagated symmetrically around the right and left lateral epicardium. The last area to be activated was in the high RVOT, similar to that for LVA. However, for LVB pacing, it took 97 ms to activate the entire epicardium. For RVOT S1 pacing, the activation sequence was almost exactly opposite to that for LVB pacing and took nearly the same amount of time to activate the epicardium (98 ms).


View larger version (69K):
[in this window]
[in a new window]
 
Fig. 6.   Examples of activation sequences for a single paced beat for all three S1 pacing sites. Isochronal maps show representative activation sequences following an S1 stimulus from LVA (A), LVB (B), and RVOT sites (C). Anterior (left) and posterior (right) views are shown; shaded areas represent the atria. The LAD (anterior view) and circumflex with posterior descending coronary arteries (posterior view) are indicated by thick shaded lines. The initial site of recorded activation is shown as time 0 (arrow) in each panel. Isochrone lines are shown at intervals of 10 ms and indicate times at which underlying tissue became activated relative to the earliest recorded activation time. For S1 pacing from the LVA site (A), the last recorded activation occurred in the high RVOT 74 ms after the earliest activation. For S1 pacing from the LVB site (B), the last recorded activation occurred in the high RVOT 97 ms after the earliest activation. For S1 pacing from the RVOT site (C), the last recorded activation occurred in the posterior LVB 98 ms after the earliest activation.

Figure 7 shows an example of the activation sequences for the first two postshock cycles leading to VF following a 140-V S2 shock after LVB S1 pacing in the same animal as in Fig. 6 (views shifted 90° with respect to Fig. 6B). After the shock, activations were initiated from two early sites. Activation first appeared on the epicardium in the right lateral wall along the base 22 ms after the start of the shock. A second wave front appeared on the opposite side of the heart in the lateral LVB 30 ms after the start of the S2. The potential gradients measured at first and second early sites were 2.3 and 3.7 V/cm, respectively. The time between the last preshock activation following the 20th S1 at the early site and the time of S2 delivery (VR-S2, see below) was 157 ms for the first site and 141 ms for the second site. The activation wave front in the RV proceeded clockwise around a line of block while the LV wave front proceeded counterclockwise around a similar line of block. In both the RV and the LV, as the activation fronts completed propagation around the ends of the lines of block, wave fronts appeared in the nearby previously activated tissue at the end of the block, indicating reentry into these areas. The second cycle wave fronts appeared in the RV and LV at 107 and 121 ms, respectively, after the S2 while the first cycle wave fronts completed activation of nearby previously unactivated myocardium. The maps for the second cycle show that the lines of block extended further, but the activation fronts still propagated around the block and reentered the previously activated tissue. In this case, after the second cycle, wave fronts propagated from the original sites, but these wave fronts broke into numerous smaller segments that activated the remainder of the epicardium. Subsequently, other small wave fronts appeared on the epicardial surface in other areas that could not be directly linked on the epicardium to the original propagating activation front. Multiple activation wavelets were present on the epicardial surface, and VF ensued. Two bilateral reentrant rotors were seen in this example, with the plane of rotation appearing parallel with the epicardial surface. In other cases with identifiable reentry near the early site, the circuit was more distorted, with the plane of rotation apparently angled into the myocardium. Frequently, activation was followed by wave front propagation away from the early site; however, one or more wave fronts could not be identified as reentering. When this occurred, activations originated at the early site repeatedly until multiple small wave fronts were present and VF ensued. It could not be distinguished whether the earliest wave fronts represented activations propagating from an automatic focus or from a rotor deep within the myocardium.


View larger version (91K):
[in this window]
[in a new window]
 
Fig. 7.   Anterolateral views of RV (left) and LV (right) showing beats 1 (A) and 2 (B) after a 140-V shock that induced VF in the same animal as in Fig. 6. Shaded areas represent the atria. The LAD can be seen at bottom right of anterolateral RV view and at top left of anterolateral LV view (thick shaded lines). S1 pacing was from the LVB site, located at top right of anterolateral LV view. Isochrone lines are spaced in 10-ms increments; numbers indicate activation times of underlying tissue after the beginning of S2 shock. Arrows indicate initial sites of activation for beats 1 and 2. At 22 ms after the S2, an activation front initiated in the lateral RV base (beat 1). A 2nd front initiated at 30 ms after the S2 in the lateral LVB. After initial activation, the 2 bilateral wave fronts propagated around lines of block (thick lines) and reentered tissue that was previously activated (curving arrows). At 107 and 121 ms after the S2 (beat 2), the 2nd postshock cycle began on the lateral RV and LV walls, respectively. Wave fronts for 2nd cycle proceeded in the same direction around extended lines of block in a reentrant fashion. After beat 2, wave fronts broke down into smaller fragments and VF occurred.

Figure 8 shows all early site positions following the S2 for LVA, LVB, and RVOT S1 pacing sites, respectively, in a single animal. For LVA pacing at lower- strength S2 shocks (400 V below ULV), three early sites of activation were observed about midway between the apex and the base, one in the RV and two in the LV. As the S2 shock strength increased, the sites moved to more basal positions and the number of early sites decreased. With pacing from the posterior LVB, bilateral early sites were slightly more posterior. However, as the shock strength increased, those sites moved to anterior and basal positions similar to that with LVA pacing. With the S1 pacing site at the RVOT, however, the sites originated anteriorly. The early sites moved to more lateral and basal positions with increasing S2 shock strength. At one step below the ULV, two early sites were located at the edge of the AV groove.


View larger version (29K):
[in this window]
[in a new window]
 
Fig. 8.   Plots of all early site locations for a single animal in terms of (rho , theta ) coordinates (see Fig. 4). Early site locations are shown for LVA (A), LVB (B), and RVOT pacing sites (C). Numbers represent S2 shock voltage relative to ULV. Arrows indicate movement of early site locations as S2 strength is increased closer to ULV level. As S2 shock strength increased, the number of early sites decreased. At 40 V below ULV, there was a single early site in the RVOT area for LVA (A) and LVB pacing (B). For RVOT pacing with S2 40 V below ULV, there were 2 early sites near each other at the right lateral RV base (C). Movements and final locations of early sites at strengths just below ULV were consistent across animals.

The progression of early sites was similar across all animals for each S1 pacing site. The distance from the tip of the LVA to the epicardial early site is plotted in Fig. 9 as a function of the S2 voltage relative to the ULV for all animals and each S1 site. There was a consistent increase in the distance from the S2 shocking electrode to the epicardial early site location as S2 shock strength increased, regardless of S1 pacing site.


View larger version (21K):
[in this window]
[in a new window]
 
Fig. 9.   Plot of polar distance from LVA to epicardial early site (y-axis) vs. S2 shock strength relative to ULV (x-axis). Polar distance (rho ) was measured as the distance in spherical coordinates from the tip of the LVA to the epicardial electrode. The y-axis range from 30 to 100 mm represents the range of epicardial electrode distances. Filled symbols represent mean values across animals for each S1 pacing site; open symbols represent 1 SD above or below mean. For each S1 pacing site as S2 strength increased, the early sites moved progressively away from the LVA S2 shock electrode.

The early site positions in all animals for shocks just below the ULV are shown in Fig. 10. In most cases, the early sites of activation for LVA and LVB S1 pacing arose in the RVOT. On the other hand, the early sites were located in slightly different locations across animals for pacing from the RVOT. Whereas some early sites were present over the RVOT, the majority were located along the AV groove at various points around the heart. Note that in the two animals in which the VF induction time corresponded to the earliest scanning time for RVOT S1 pacing, the early sites at one step below the ULV were located in the RVOT closer to the pacing site. In the remaining seven animals, the early sites were not present in the RVOT.


View larger version (48K):
[in this window]
[in a new window]
 
Fig. 10.   Positions on heart of all early sites for S2 shocks 1 step below ULV strength (ULV at -10, -20, and -40 V) in all animals. Left, anterior view; right, posterior view (see Fig. 1). , S1 pacing locations. Early sites for RVOT tend to group along lateral RV and posterolateral RV base. Sites for LVA and LVB group in RVOT.

Early Site Potential Gradients

Figure 11 shows early site locations and potential gradient values for shock strengths of 60, 220, and 380 V in the same animal as in Fig. 8. The S2 shock electrode configuration produces a field that is strong at the apex and decreases toward the base. The field is weakest along the lateral LV and RV along the AV groove and in the RVOT. At low S2 shock levels, the total epicardial surface area that exceeds 10 V/cm is small. As the shock strength is raised, however, the epicardial potential gradient field intensifies. As this occurs, the earliest sites of postshock activation move to more basal positions, remaining in low potential gradient areas.


View larger version (103K):
[in this window]
[in a new window]
 
Fig. 11.   Four views of electrode sock showing sites of earliest postshock activation leading to VF and the potential gradient fields for 380-V (A), 220-V (B), and 60-V (C) S2 shocks delivered between LVA patch and SVC catheter electrodes for same animal as in Fig. 8. Pacing was from LVA S1 site. Each potential gradient measurement site is represented by a shaded patch. Scale (bottom) shows corresponding potential gradient level (in V/cm) for each level of shading: low potential gradient regions are represented by lighter shades; darkest shade represents values greater than ~10.5 V/cm. White "holes" represent sites with missing values. For S2 = 60 V (C), the potential gradient field was <10 V/cm at most measurement sites. In this case there were 3 early sites (star ) all in the midtransverse region. For S2 = 220 V (B), the potential gradient field exceeded ~10.5 V/cm in almost the entire lower one-half of the heart. Two early sites were present at more basal positions compared with the early site locations for S2 = 60 V. For S2 = 380 V (A), the potential gradient field was strong in all but a few areas of the heart. Weak potential gradient areas were present in RVOT and along RV and LV base. At this shock level, there were 2 early sites, 1 in high RVOT and 1 in posterolateral RV base in top row of electrodes. In this case, the ULV was achieved when S2 was increased by another 80 V (S2 = 420 V).

The average potential gradient magnitudes at the early sites are plotted in Fig. 12 as a function of S2 LEV for all three S1 pacing sites. The potential gradients at the early sites increase with increasing S2 shock strength for each S1 site (P < 0.001 for each slope vs. zero, analysis of covariance, r = 0.4). The gradients for each S1 site increase at the same rate (P > 0.10 for slope comparisons, analysis of covariance).


View larger version (31K):
[in this window]
[in a new window]
 
Fig. 12.   Plot of early site potential gradient magnitudes vs. S2 shock strength relative to ULV level. Each filled symbol represents mean value across animals for a specific S2 shock level and S1 pacing site; open symbols represent 1 SD for each S1 pacing site.

Figure 13 shows an example of the location of the early sites relative to the potential gradient field after shocks one step below the ULV for all three S1 pacing sites in the same animal as in Figs. 8 and 10. The early sites at this shock level arose in regions of the heart remote from the S2 electrodes, where the potential gradient field was weak. The early site potential gradient data measured at one step below the ULV are summarized in Table 2. Among S1 pacing sites, the early site potential gradient magnitudes were not significantly different for ULV40, ULV20, or ULV10. The same was true when the data for ULV40, ULV20, and ULV10 were combined (LVA: 7.0 ± 3.9 V/cm; LVB: 6.9 ± 4.3 V/cm; RVOT: 7.9 ± 4.1 V/cm; P = 0.53 between S1 sites, ANOVA). For the ULV40 group, the maximum measured potential gradients on the heart were 264 ± 128, 266 ± 142, and 243 ± 142 V/cm for the LVA, LVB, and RVOT pacing sites, respectively. The respective minimum measured potential gradients were 1.0 ± 0.8, 1.7 ± 1.0, and 1.4 ± 1.0. Thus the early sites were located in low potential gradient regions but were not the minimum values on the epicardium.


View larger version (91K):
[in this window]
[in a new window]
 
Fig. 13.   Comparison of early site locations and potential gradient fields for each S1 site when S2 strength was 40 V less than ULV in same animal as in Figs. 8 and 11. Each potential gradient measurement site is represented by a shaded patch. Scale (right) shows corresponding potential gradient level (in V/cm) for each level of shading: low potential gradient regions are represented by lighter shades; darkest shade represents values greater than ~10.5 V/cm. White holes represent sites with missing values. Early site locations are indicated by an arrow and asterisk. Regardless of pacing site, distribution of potential gradient values was similar in each case. For LVA (A) and LVB (B), single early sites were located in RVOT. However, for RVOT (C), a single early site was present along lateral RV base.


                              
View this table:
[in this window]
[in a new window]
 
Table 2.   Potential gradients at earliest postshock activation sites for shocks one step lower than the ULV

Early Site Activations

To determine whether there was a relationship between early site activation and recovery following the last S1 stimulus and the timing of the S2 shock that induced VF, timing intervals were calculated and compared for S2 shocks one step below the ULV in the ULV40 group. The VR-S2 interval was calculated as the time between preshock activation after the last S1 at each early site and the time of S2 delivery. There were no significant differences in the early site VR-S2 intervals between S1 pacing sites (LVA: 192 ± 52 ms; LVB: 212 ± 38 ms; RVOT: 209 ± 41 ms; P = 0.48, repeated-measures ANOVA). This indicates that, regardless of pacing site, the earliest postshock activations arose in areas that had been activated at approximately the same time preceding S2 delivery.

To further assess the timing of S2 delivery relative to the timing of local recovery, the RS1-S2 interval was calculated as the time between the maximum positive slope of the T wave following the last S1 stimulus at each early site and the time of S2 delivery. A negative RS1-S2 interval indicates that the S2 shock was delivered earlier than the maximum positive slope of the local T wave. Across animals, the RS1-S2 interval was not significantly (P > 0.14) different from zero for any S1 pacing site (LVA: -8.0 ± 19 ms; LVB: -7.6 ± 17 ms; RVOT: -2.0 ± 12 ms). Additionally, there were no significant differences (P = 0.41, repeated-measures ANOVA) among pacing sites for this interval. Because the RS1 point for unipolar electrograms has been shown to correlate with the refractory period of the underlying tissue, (22) these data suggest that early sites were located where the tissue was in a particular state of refractoriness at the time of the S2.

The S2-V1 interval is the time delay between the beginning of the S2 shock and the time of the first recorded post-S2 activation at the early site. The S2-V1 interval is plotted in Fig. 14 as a function of shock strength relative to the ULV for all S1 pacing sites. As S2 shock strength increased, the S2-V1 interval increased for all three S1 pacing sites. However, the S2-V1 interval increased at a faster rate for the RVOT S1 pacing site compared with the LVA and LVB pacing sites (P < 0.01, analysis of covariance). At one S2 step below the ULV (-40, -20, and -10 V combined), the S2-V1 interval was significantly (P < 0.05) longer for RVOT (82.2 ± 27 ms) than for LVA (64.1 ± 21 ms) and LVB (61.8 ± 21 ms). S2-V1 intervals were not significantly different for LVA versus LVB. For all S1 pacing sites, the early site potential gradient and early site S2-V1 interval showed slight positive correlation (r = 0.44).


View larger version (30K):
[in this window]
[in a new window]
 
Fig. 14.   Plot of time delay between beginning of S2 shock and time of first recorded post-S2 activation at early site (S2-V1 interval) vs. S2 shock strength relative to ULV. Each data point (filled symbol) is average value across animals at a specific relative S2 shock strength; open symbols are 1 SD above or below mean. As S2 shock strength increased to ULV, average S2-V1 interval increased for all three S1 pacing sites. Interval for RVOT increased at a significantly faster rate than that for LVA and LVB (P < 0.01, analysis of covariance).

The potential gradient measured at each early site is plotted in Fig. 15 against the RS1-S2 interval at the same site. Early sites where the S2-V1 interval was <25 ms are shown with a different marker. As seen in Fig. 15, the earliest sites of activation with the shortest S2-V1 intervals had a smaller divergence of potential gradient (mean 2.8 ± 1.6 V/cm) and RS1-S2 (mean 12.6 ± 11.1 ms) values. These points of earliest activation arose where, simultaneously, there was a similar degree of relative refractoriness and similar potential gradient at the time of the S2 shock.


View larger version (19K):
[in this window]
[in a new window]
 
Fig. 15.   Plot of epicardial potential gradient magnitude at each early site vs. time of maximum positive derivative of local recorded T wave following last S1 stimulus at each early site (RS1-S2 interval) at same site for all pacing sites in all animals. Small shaded points represent early sites at which S2-V1 interval exceeded 25 ms. black-triangle, data for early sites at which S2-V1 interval was 25 ms or earlier. Clustering of data points for sites with short latency (S2-V1 <=  25 ms) indicates that across animals and pacing configurations, the earliest postshock activations arose where there was a similar potential gradient and similar state of relative refractoriness at time of S2 shock.

ULV Shocks

Early sites of activation were analyzed for 67 ULV shock sequences. In 73% of the cases, a single beat occurred after the S2. This was followed by a long pause (>400 ms) during which there was no detectable electrical activity. Sinus rhythm or an "escape" rhythm ensued. In 14% of the cases, two to four early repetitive activations were observed before a pause. A period of initial post-S2 asystole followed by a delayed escape or "sinus" beat was the first activation in only 13% of the cases. The early sites following ULV-level shocks that did not induce VF were located in areas that were similar to those for shocks just below the ULV that did induce VF (Fig. 10). In most cases, activation initiated as if the shock were going to induce VF. An activation wave front commenced and rotated as if it were going to reenter. In some cases, a single rotation occurred but the wave front ceased. In other cases of repetitive activity, a rotating wave front was present that extinguished after two cycles. In one case, activation from the earliest site rotated around a point; however, a wave front from a nearby second early site appeared 58 ms after the first and appeared to collide with the first wave front. In other examples, a wave front appeared from a central location and propagated outward, circumferentially. It could not be determined whether the wave front initiated from a reentrant site deep beneath the epicardium or from a focal source. Activation ceased after a single event in most cases.

The potential gradients measured at these early sites for ULV level shocks are summarized in Table 3 for each of the S1 pacing sites. There were no significant differences in the early site potential gradients at the ULV between S1 pacing sites (P = 0.67). However, the early site potential gradients measured at the ULV were significantly (P < 0.01) higher than those measured for shocks that were one step below the ULV (ULV: 8.9 ± 4.8 V/cm vs. ULV-step: 7.3 ± 4.1 V/cm).

                              
View this table:
[in this window]
[in a new window]
 
Table 3.   Potential gradients at earliest postshock activation sites for shocks at the ULV


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Global Shock Variables and Shock Timing

The relationship between S2 shock strength and the myocardial potential gradient field is determined by the S2 electrode configuration (6, 27), cardiac geometry and blood volume (15, 16), and myocardial conductivity (27). If these factors are held constant, the spatial distribution of potential gradient values should not vary for repeated shocks at a given strength. Under these conditions, the critical point hypothesis predicts that the ULV shock strength should be a constant, and the ULV should not be affected by the direction of activation and recovery preceding the shock. In other words, the shock strength at the ULV should be the same, regardless of the S1 pacing site.

In the current study, the ULV was determined using a single S2 shock electrode configuration and three S1 pacing sites that were placed in distinct areas of the heart. Regardless of this disparity in the origin of ventricular activation, the measured ULVs at resolutions of 10, 20, and 40 V were not significantly different with the use of repeated-measures analysis to adjust for interanimal variation. These findings are consistent with the predictions of the critical point hypothesis for the ULV.

In a few animals, there was a variation in the measured ULVs during the step-up testing for the ULV40 (Table 1). No particular S1 pacing site appeared to result in a different ULV. These results may be reflective of the dose-response relationship between shock strength and the ULV (29). Step-up testing will determine a point in a lower and potentially flatter portion of the curve. Thus there may be more variation in the measured ULV shock strength. Up/down testing with sequentially smaller step sizes will determine a more accurate measure of the 50% point on the dose-response curve. At that point, the slope is steeper and there is less potential for random variation of the ULV (32). When we performed further testing for the ULV in the ULV20 and ULV10 sets, the ULVs converged closer between S1 pacing sites in those animals.

Whereas the global shock strength at the ULV did not appear to be affected by the choice of S1 pacing site, the timing of VF inductions for shocks below the ULV was altered by this choice. For all three S1 sites, S2 shocks were delivered starting at the same point in the T wave. However, as shock strength increased, the VF induction points in the T wave diverged on the basis of the origin of ventricular activation (Fig. 5). S2 shocks had to be delivered earlier than the peak of the T wave for RVOT pacing but later than the peak for LVA and LVB pacing.

On the basis of the critical point theory, the time points in the electrical cycle at which an S2 induces VF should depend on the S1 pacing site. The critical point hypothesis predicts that properly timed S2 shocks slightly weaker than the ULV initiate VF because the critical potential gradient level is present in one or more areas of the heart. The critical potential gradient should be located in these same areas of the heart as long as the S2 shock strength, S2 configuration, and intrinsic factors remain unaltered. As a result, the critical point hypothesis predicts that after shocks slightly weaker than the ULV, VF will be initiated in one of these areas regardless of the S1 pacing site. However, the S1 pacing site will affect the timing of S2 shocks that induce VF because the time at which one of the regions becomes activated and critically recovered depends on the sequence of ventricular activation. The critical potential gradient area will activate soon after the S1 and pass through the critical degree of recovery early in the T wave if the S1 site is nearby. However, pacing from a distant S1 site will cause the area to activate later after the S1 and pass through the critical degree of recovery later in the T wave. Therefore, the critical point hypothesis predicts that VF will be induced for S2 shocks delivered early in the T wave for S1 sites close to a critical area and later in the T wave for sites distant from a critical area.

These findings have important implications concerning testing methods used to determine the ULV. In this study, the entire T wave was scanned with S2 shocks at a particular shock strength. This method requires the delivery of a large number of shocks, thus decreasing its clinical practicality. In clinical practice, others have employed techniques using a single point in the T wave or a limited scanning interval to determine the ULV (5, 13). The use of such techniques in the current study would probably have altered the results. If only a small segment on the upstroke of the T wave were scanned (see Fig. 5), the ULV for S1 pacing from the RVOT site would probably have been similar to that found with full T-wave scanning. However, for the LVA or LVB S1 sites, VF induction times at later points in the T wave would have been untested. Therefore, the measured ULV would have been lower than the absolute ULV. Thus, whereas a point in the upstroke of the T wave may be optimal for testing with RVOT S1 pacing, the same point would be much less desirable for LVA or LVB S1 pacing.

Early Sites

On the basis of the critical point hypothesis, low- strength S2 shocks should initiate activations leading to VF at one or more locations in low potential gradient regions that contain the critical potential gradient. The specific locations will depend on where the critical degree of refractoriness intersects the critical potential gradient. However, as shock strength increases, the critical potential gradient will be exceeded in increasingly larger areas, encompassing tissue farther from the shock electrodes. Consequently, activations will be initiated at an increasing distance from the shock electrodes. For shocks just below the ULV strength, the earliest sites of activation should arise from the few areas of the heart that contain the critical potential gradient. Although the earliest sites of activation do not represent the critical point locations themselves, because activation will not be occurring at the singularity point, the measured potential gradients at those sites should be similar. In addition, the level of tissue refractoriness at the early sites at the time of the shock should be similar. These measured values at the early sites should not be affected by the choice of S1 pacing site.

Early site locations. In this study, as predicted by the critical point hypothesis, shocks well below the ULV induced activations leading to VF at multiple apical sites on the heart, closer to the LVA shocking electrode. As shock strength increased, the earliest sites of activation moved away from the S2 shocking electrode and decreased in number. The movement of early sites for increasing shock strengths is similar to that shown by Shibata et al. (28) using a similar S2 configuration but a different waveform in dogs. As in their study, we observed oppositely rotating bilateral reentrant wave fronts for some lower-strength shocks. We have shown further that the early sites for initiating these wave fronts originated where the potential gradients were low and of similar value. The early sites also originated where the tissue was in a similar state of recovery following prior S1 pacing (similar VR-S2 interval).

As predicted by the critical point hypothesis, as shock strength continued to increase, the early sites converged to specific regions of the heart. At shock levels just below the ULV, activation fronts first appeared on the epicardium primarily in one region for the LVA and LVB pacing sites but in a separate region for the RVOT pacing site. The early sites were grouped in the high RVOT when pacing was performed from the LVA or LVB. When S1 pacing was performed from the RVOT, the early activations arose primarily in the posterolateral RV base near the AV groove. These locations at one shock step below the ULV were consistent across animals, even though the shocks that were delivered ranged from 250 to 500 V (LEV). Therefore, the global shock strength was not a good predictor of the early site locations for different animals.

There are several possible reasons for the finding that the early sites for S1 pacing from the RVOT differ from the early sites for S1 pacing from the LVA and the LVB. Winfree (37) describes the formation of reentry via a critical point as the result of the intersection of an isorefractory contour shell with an isogradient contour shell. For reentry to form, there must be some finite angle between these two contours. It is possible that with S1 pacing from the RVOT, the angles of intersection between these two contours in the RVOT region were too shallow to initiate reentry via the critical point mechanism. Additionally, prior pacing close to a "critical" region and/or the direction of activation and recovery through the region may have influenced the formation of early activation sites. The direction of activation and recovery relative to myocardial fiber orientation in the region also may have played a role in this finding (4, 24). It is also possible that pacing from the RVOT S1 site resulted in activations that were initiated deeper within the myocardium. This may be reflected in the significantly longer S2-V1 interval for RVOT compared with that for LVA and LVB for shocks just below the ULV (Fig. 14).

Early sites at the ULV. Repetitive activity following ULV-strength S2 shocks has been noted previously (28). However, the origin of these postshock activations after ULV-level shocks has not been documented. In this study, we found that most activations following ULV- level S2 shocks arose in the same areas of the heart where the early sites occurred for S2 shocks just weaker than the ULV. These activations following ULV- level shocks either self-extinguished or were halted by collisions with other early site wave fronts. The results indicate that whereas a single mechanism may be responsible for initiating the postshock activations, there may be several mechanisms that prevent their continuation. As we have shown, increases in the S2 shock strength moved the early sites away from the S2 shocking electrodes. When a single rotational wave front was observed but failed to reenter, the increased S2 strength may have moved the early site just close enough to a tissue border that an activation front was induced but reentry could not proceed (36). In cases with multiple beats from a single early site, the reentrant circuit may have drifted (8) into a border and extinguished, although active "drifting" was not observed, perhaps due to the mapping resolution in this study. Recent modeling combines the critical point mechanism with bidomain theory. In that model, the critical potential gradient remains present, regardless of shock strength. However, due to anisotropy and its spatial effects on stimulation, higher-strength shocks may create rotors that are too close together to survive (24). With any of these mechanisms, small local variations in potential gradient or refractoriness may result in continuing reentrant activity rather than cessation after one or two beats. Thus a probability factor exists that may explain the dose-response relationship for the ULV (29). High-resolution three-dimensional mapping close to the tissue borders needs to be performed to determine the exact mechanism for the ULV.

Early site activations. Two intervals from the early site electrograms were used to assess the timing of the S2 shock relative to the local recovery time for S2 shocks just below the ULV. The interval from the time of local activation to the time of the S2 (VR-S2 interval) was not significantly affected by the S1 pacing site. If it is assumed that the refractory periods are the same at all early sites, then the results would indicate that the earliest sites of activation always arose in areas where the tissue was at a similar degree of relative refractoriness.

The time of the maximum positive slope of the T wave in unipolar electrograms has been shown to correlate well with the downstroke of the action potential in underlying tissue (39). The time of this maximum positive derivative corresponds to times that are slightly earlier than the effective refractory period (22). Thus the timing of the S2 relative to this point in the early site T wave (RS1-S2) would give an indication of the timing of the S2 relative to the effective refractory period. This interval could then be used to compare the early sites for the different S1 pacing sites. There was a substantial amount of variation in the measured values, which could reflect measurements for some earliest sites of activation that originated deeper in the myocardium and propagated to the surface. Overall, the intervals were not significantly different from zero for early sites from each S1 pacing site. In addition, there were no significant differences across pacing sites. These results indicate that the level of tissue recovery at the time of the S2 was approximately the same for the early sites from each pacing origin. The data indicate that, on average, the earliest sites of activation arose where the tissue was close to its effective refractory period at the time of the S2 shock.

Interestingly, the interval (S2-V1) from the beginning of the S2 shock to the earliest recorded activation on the epicardium increased with increasing shock strength (Fig. 14). This interval, termed "latency" (4), has been shown to increase with increasing shock prematurity (28). Although this may explain the longer latency for the RVOT site versus the other S1 sites at shocks just below the ULV level, because the RVOT shocks that induced VF were delivered at shorter coupling intervals, it does not explain the overall increase in latency with shock strength. Furthermore, the degree of prematurity relative to the refractoriness of the tissue should be the same for the RVOT shocks as for the other shocks. One explanation for the latency may be that the early sites of activation originated in tissue below the epicardium. The latency would therefore represent a conduction delay from the deeper tissue to the epicardial surface.

Early site potential gradients. The potential gradients measured at the earliest sites of activation for shocks just below the ULV averaged 5-9 V/cm. As predicted by the critical point hypothesis for the ULV, the S1 pacing site did not have an effect on the early site potential gradients, even though the early sites arose in a different region of the heart for the RVOT pacing site. Frazier et al. (11) reported potential gradients of ~5 V/cm at the critical point of reentry induced in a small region of the RVOT. A different shock waveform, S1 site, and S2 electrode configuration were used in that study. However, the potential gradients measured at the early sites in this study would be expected to be slightly different. With the critical point mechanism, activations would initiate in an area adjacent to the singularity point, where the potential gradient field is slightly weaker than at the critical point itself.

For shocks just below the ULV, the early sites were located in low potential gradient regions on the heart (Fig. 13) but were not located at the minimum potential gradient areas. The formation of a critical point requires the presence of a spatial dispersion of potential gradient that spans the critical value. In addition, there must be enough distance between a boundary and the critical point for the reentrant circuit to circulate (37). Therefore, potential gradients weaker than the critical potential gradient must be present.

Although the measured potential gradient magnitudes at the early sites for shocks just below the ULV averaged 5-9 V/cm, the early site gradients for lower- strength shocks were weaker (~3 V/cm for shocks at 440 V below the ULV). As shock strengths increased and the early site locations moved to more basal positions at higher shock strengths, there was a consistent increase in the early site potential gradients, regardless of S1 pacing site (Fig. 12).

If early sites at the base are well below the epicardium, as supported by the increased latency for these sites, the shock potential gradient may be different at these intramural early sites from those at the mapped early site on the epicardium. If so, the increase in epicardial potential gradient with an increase in shock strength at the earliest site of activation on the epicardium (Fig. 12) may not be true intramurally at the actual early site. The clustering of potential gradients for short-latency early sites (Fig. 15) would therefore suggest that the potential gradients close to a critical point were ~3 V/cm. It is also possible that the results indicate a difference in the critical potential gradient based on heart location. It is known that tissue anisotropy influences local wave front propagation and potential gradient distribution (10). Anisotropy can also influence the field requirements for stimulation (10, 20) and for refractory period prolongation (18). As suggested by the bidomain formulation, other factors such as fiber curvature may also be influential (30). Therefore, differences in tissue structure at different levels of the heart may have altered the potential gradient field requirements for critical point formation. Besides tissue structure, bidomain theory predicts that the changes in transmembrane potential are proportional to the derivative of the potential gradient field rather than to the potential gradient itself (26, 31). Because the derivative of the potential gradient decreases with distance from the S2 electrode, for this derivative to remain constant the potential gradient would have to be smaller for early sites close to the S2 electrode than for more distant early sites.

The results may also simply indicate that the two-dimensional potential gradient magnitude on the epicardial surface was not the optimal quantity to measure. It is possible that the results would have been different if measurement of the third vector component of the potential gradient were included. The vector component in the third direction would have measured the potential gradient within the tissue and normal to the epicardium. It is possible that this component was more important in areas that were closer to the apex shocking electrode. Thus inclusion of this component may have increased the calculated potential gradient magnitude for early sites located at lower positions on the heart. Furthermore, the leading-edge potential gradient magnitude may not be the optimal measurement point. It is possible that trailing-edge potential gradient or mean potential gradient may have yielded a more constant value.

Limitations

A major limitation of the current study was the use of mapping from the epicardial surface alone. This restricted the measurement of activations to those close to or just beneath the surface of the recording electrodes. It is possible that activations were initiated from deep within the myocardium by processes other than the critical point and then simply propagated to the surface. Unrecorded deep early sites could have been responsible for the initiation of VF. However, given this limitation, the results in this study have fulfilled multiple predictions based on the critical point hypothesis for the ULV. In addition, new information has been obtained with regard to the location of early nonsustained repetitive activity following ULV strength shocks. Three-dimensional mapping with electrodes inserted into the myocardium would have been optimal. However, a large number of intramural electrodes could potentially damage the myocardium as well as cause distortions of the shock potential gradient field.

We used multielectrode plaques to record from the epicardial surface. This was performed to reduce potential gradient measurement error caused by errors in interelectrode distance measurements. Although this optimized potential gradient measurements, it restricted our ability to map evenly across the epicardial surface. Our ability to track activations back to the earliest sites would have been improved by reducing the interplaque distances and increasing the number of recording electrodes; however, this was not possible due to technical limitations.

In this study, we used the leading-edge S2 potential gradient magnitude as the parameter to describe the local shock "force" affecting the tissue. In theory, the formation of critical point reentry is most dependent on the cellular state of the tissue at the instant the S2 stimulus is removed. As long as the regions of varying cellular states exist in the proper spatial orientation with an intervening singularity point, reentry will ensue spontaneously. The use of this single parameter was based on previous studies involving the critical point and potential gradient measurements during S2 shock delivery (6, 11). This single parameter measured at the beginning of the S2 may not be the best predictor of cellular response. Differing cell responses based on S2 intensity, duration, polarity, and orientation have been documented (19). Other parameters such as the mean potential gradient, trailing-edge potential gradient, or potential gradient combined with shock duration and/or polarity may be better predictors of cellular response and reentry formation.

Interestingly, although we used a biphasic waveform in this study, it is unknown whether the results would have been different had a monophasic S2 shock been utilized. The similarity of results between our study and that of Shibata et al. (28) in terms of early site location and movement with S2 shock strength suggests that waveform polarity may not play a significant role, because they used a monophasic S2 shock. However, it is unknown whether the early site potential gradients would have been different had a monophasic waveform been used.

In summary, in this study the ULV was determined by using a single S2 shock electrode configuration and three separate S1 pacing sites. The measured ULVs were not affected by the S1 pacing site in terms of global shock strength variables (LEV, current, and total shock energy). However, the choice of S1 pacing site changed the time points in the T wave at which VF was induced for shocks below the ULV. These results were predicted by the critical point hypothesis for the ULV.

Epicardial mapping of cardiac activations revealed that shocks at strengths far below the ULV induced activations in areas closer to the S2 shocking electrode. As shock strength was increased, the early sites were moved away from the shocking electrode to more remote areas where the potential gradient field remained weak. When slightly stronger shocks were delivered, one or several early activations originated from the same remote areas but halted and did not result in VF.

For this S2 configuration, shocks at strengths just weaker than the ULV initiated activations in the RVOT when pacing was from the LVA or LVB. However, the early sites were located at the lateral aspect of the RV when pacing was performed from the RVOT. Regardless of the pacing origin, shocks just below the ULV level induced activations for which the potential gradient magnitude was ~5-9 V/cm, and the level of tissue recovery was similar at the time of the shock. However, for shocks that were much weaker than the ULV, the early site potential gradient magnitudes were lower.

On the basis of the critical point hypothesis, predictions were made regarding the effect of S1 pacing site on the ULV shock strength, the timing of VF inductions in the T wave, the locations of earliest post-S2 epicardial activations, and the potential gradients and level of tissue recovery at the early sites. The results support most of these predictions. The few predictions not verified in this study, such as the prediction of a constant potential gradient at the early sites, indicate either that the critical point hypothesis is correct but the lack of intramural electrodes in this study did not allow us to identify the true early site or that other factors in addition to the critical point, such as bidomain effects and the direction of activation and recovery, are influential. The results indicate that the critical value of the potential gradient may not be a constant but may change with location in the heart. The ULV appears to result from movement of the critical potential gradient to remote areas where other factors, such as the proximity of heart borders, restrict the formation of reentry around the critical point.


    ACKNOWLEDGEMENTS

This study was supported in part by National Heart, Lung, and Blood Institute Grants HL-42760 and HL-44066, American Heart Association-North Carolina Affiliate Grant NC-93-SA-05, and National Science Foundation Engineering Research Center Grant CDR-8622201.


    FOOTNOTES

Address for reprint requests and other correspondence: R. E. Ideker, Univ. of Alabama at Birmingham, 1670 Univ. Blvd., Volker Hall B140, Birmingham, AL 35294 (E-mail: rei{at}crml.uab.edu).

Received 10 October 1997; accepted in final form 30 June 1999.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

1.   Ahlstrom, M. L., and W. J. Tompkins. Digital filters for real-time ECG signal processing using microprocessors. IEEE Trans. Biomed. Eng. 32: 708-713, 1985[Medline].

2.   Bacon, M. E., R. N. Vitullo, J. H. Kasell, T. Blitchington, J. C. Pressley, P. K. Smith, J. E. Lowe, R. E. Ideker, and J. M. Wheaton. Upper limit of vulnerability and its relationship to defibrillation threshold in humans (Abstract). Pacing Clin. Electrophysiol. 15: 530, 1992.

3.   Blanchard, S. M., R. J. Damiano, Jr., T. Asano, W. M. Smith, R. E. Ideker, and J. E. Lowe. The effects of distant cardiac electrical events on local activation in unipolar epicardial electrograms. IEEE Trans. Biomed. Eng. 34: 539-546, 1987[Medline].

4.   Chen, P. S., Y. M. Cha, B. B. Peters, and L. S. Chen. Effects of myocardial fiber orientation on the electrical induction of ventricular fibrillation. Am. J. Physiol. 264 (Heart Circ. Physiol. 33): H1760-H1773, 1993[Abstract/Free Full Text].

5.   Chen, P. S., G. K. Feld, J. M. Kriett, M. M. Mower, R. Y. Tarazi, R. P. Fleck, C. D. Swerdlow, E. S. Gang, and R. M. Kass. Relation between upper limit of vulnerability and defibrillation threshold in humans. Circulation 88: 186-192, 1993[Abstract/Free Full Text].

6.   Chen, P. S., P. D. Wolf, F. Claydon, E. G. Dixon, H. J. Vidaillet, Jr., N. D. Danieley, T. C. Pilkington, and R. E. Ideker. The potential gradient field created by epicardial defibrillation electrodes in dogs. Circulation 74: 626-636, 1986[Abstract/Free Full Text].

8.   Davidenko, J. M., A. V. Pertsov, R. Salomonsz, W. Baxter, and J. Jose. Stationary and drifting spiral waves of excitation in isolated cardiac muscle. Nature. 355: 349-351, 1992[Medline].

9.   Fabiato, A., P. Coumel, R. Gourgon, and R. Saumont. Le seuil de reponse synchrone des fibres myocardiques. Application a la comparaison experimentale de l'efficacite des differences formes de chocs electriques de defibrillation. Arch. Mal. Coeur Vaiss. 60: 527-544, 1967[Medline].

10.   Frazier, D. W., W. Krassowska, P.-S. Chen, P. D. Wolf, E. G. Dixon, W. M. Smith, and R. E. Ideker. Extracellular field required for excitation in three-dimensional anisotropic canine myocardium. Circ. Res. 63: 147-164, 1988[Abstract/Free Full Text].

11.   Frazier, D. W., P. D. Wolf, J. M. Wharton, A. S. L. Tang, W. M. Smith, and R. E. Ideker. Stimulus-induced critical point. Mechanism for electrical initiation of reentry in normal canine myocardium. J. Clin. Invest. 83: 1039-1052, 1989.

12.   Geddes, L. A., L. E. Baker, and A. G. Moore. Optimum electrolytic chloriding of silver electrodes. Med. Biol. Eng. 7: 49-56, 1969[Medline].

13.   Hwang, C., C. D. Swerdlow, R. M. Kass, E. S. Gang, W. J. Mandel, C. T. Peter, and P.-S. Chen. Upper limit of vulnerability reliably predicts the defibrillation threshold in humans. Circulation 90: 2308-2314, 1994[Abstract/Free Full Text].

14.   Idriss, S. F. Characterization of the Upper Limit of Ventricular Vulnerability (PhD thesis). Durham, NC: Duke University, 1995.

15.   Idriss, S. F., M. P. Anstadt, G. L. Anstadt, and R. E. Ideker. The effect of cardiac compression on defibrillation efficacy and the upper limit of vulnerability. J. Cardiovasc. Electrophysiol. 6: 368-378, 1995[Medline].

16.   Idriss, S. F., S. B. Melnick, P. D. Wolf, W. M. Smith, and R. E. Ideker. Predicting the potential gradient field in ventricular fibrillation from shocks delivered in paced rhythm. Am. J. Physiol. 268 (Heart Circ. Physiol. 37): H2336-H2344, 1995[Abstract/Free Full Text].

17.   King, B. G. The Effect of Electric Shock on Heart Action With Special Reference to Varying Susceptibility in Different Parts of the Cardiac Cycle (PhD thesis). New York: Columbia University, 1934.

18.   Knisley, S. B., W. M. Smith, and R. E. Ideker. Effect of field stimulation on cellular repolarization in rabbit myocardium: implications for reentry induction. Circ. Res. 70: 707-715, 1992[Abstract/Free Full Text].

19.   Knisley, S. B., W. M. Smith, and R. E. Ideker. Effect of intrastimulus polarity reversal on electric field stimulation thresholds in frog and rabbit myocardium. J. Cardiovasc. Electrophysiol. 3: 239-254, 1992.

20.   Krassowska, W., D. W. Frazier, T. C. Pilkington, and R. E. Ideker. Potential distribution in three-dimensional periodic myocardium---Part II: Application to extracellular stimulation. IEEE Trans. Biomed. Eng. 37: 267-284, 1990[Medline].

21.   Marble, A. E., M. McIntyre, R. Hastings-James, and C. W. Hor. A comparison of digital algorithms used in computing the derivative of left ventricular pressure. IEEE Trans. Biomed. Eng. 28: 524-529, 1981[Medline].

22.   Millar, C. K., F. A. Kralios, and R. L. Lux. Correlation between refractory periods and activation-recovery intervals from electrograms: effects of rate and adrenergic interventions. Circulation 72: 1372-1379, 1985[Abstract/Free Full Text].

23.   Plonsey, R., and R. C. Barr. Bioelectricity: A Quantitative Approach (1st ed.). New York: Plenum, 1988.

24.   Roth, B. J. The pinwheel experiment revisited. J. Theor. Biol. 190: 389-393, 1998[Medline].

25.   SAS Institute. SAS/STAT User's Guide, Release 6.03. Cary, NC: SAS Institute, 1988.

26.   Sepulveda, N. G., B. J. Roth, and J. P. Wikswo, Jr. Current injection into a two-dimensional anisotropic bidomain. Biophys. J. 55: 987-999, 1989[Medline].

27.   Sepulveda, N. G., J. P. Wikswo, Jr., and D. S. Echt. Finite element analysis of cardiac defibrillation current distributions. IEEE Trans. Biomed. Eng. 37: 354-365, 1990[Medline].

28.   Shibata, N., P.-S. Chen, E. G. Dixon, P. D. Wolf, N. D. Danieley, W. M. Smith, and R. E. Ideker. Influence of shock strength and timing on induction of ventricular arrhythmias in dogs. Am. J. Physiol. 255 (Heart Circ. Physiol. 24): H891-H901, 1988[Abstract/Free Full Text].

29.   Souza, J. J., R. A. Malkin, and R. E. Ideker. Comparison of upper limit of vulnerability and defibrillation probability of success curves using a nonthoracotomy lead system. Circulation 91: 1247-1252, 1995[Abstract/Free Full Text].

30.   Trayanova, N. A., B. J. Roth, and L. J. Malden. The response of a spherical heart to a uniform electric field: a bidomain analysis of cardiac stimulation. IEEE Trans. Biomed. Eng. 40: 899-908, 1993[Medline].

31.   Tung, L. A Bi-Domain Model for Describing Ischemic Myocardial DC Potentials (PhD thesis). Cambridge, MA: MIT, 1978.

32.   Wetherill, G. B. Sequential estimation of quantal response curves. J. Royal Stat. Soc. B25: 1-43, 1963.

33.   Wharton, J. M., V. J. Richard, C. E. Murry, Jr., E. G. Dixon, K. A. Reimer, J. Meador, W. M. Smith, and R. E. Ideker. Electrophysiological effects of monophasic and biphasic stimuli in normal and infarcted dogs. Pacing Clin. Electrophysiol. 13: 1158-1172, 1990[Medline].

34.   Wiggers, C. J. The mechanism and nature of ventricular defibrillation. Am. Heart J. 20: 399-412, 1940.

35.   Wiggers, C. J., and R. Wegria. Ventricular fibrillation due to single, localized induction and condenser shocks applied during the vulnerable phase of ventricular systole. Am. J. Physiol. 128: 500-505, 1940.

36.   Winfree, A. T. The electrical thresholds of ventricular myocardium. J. Cardiovasc. Electrophysiol. 1: 393-410, 1990.

37.  Winfree, A. T. Sudden cardiac death: a problem in topology. Sci. Am. 248: 144-149, 152-157, 160-161, 1983.

38.   Wolf, P. D., D. L. Rollins, T. F. Blitchington, R. E. Ideker, and W. M. Smith. Design for a 512 channel cardiac mapping system. In: Biomedical Engineering: Opening New Doors. Proceedings of the Fall 1990 Annual Meeting of the Biomedical Engineering Society. New York: New York Univ. Press, 1990, p. 5-13.

39.  Wyatt, R. F. Comparison of activation and recovery times from bipolar and unipolar electrograms to in vivo transmembrane action potential durations.Proc. IEEE Eng. Med. Biol. Soc. 2nd Ann. Conf. 1980, vol. 1, p. 22-25.


Am J Physiol Heart Circ Physiol 277(5):H2065-H2082
0002-9513/99 $5.00 Copyright © 1999 the American Physiological Society



This article has been cited by other articles:


Home page
CirculationHome page
N. Chattipakorn, P. C. Fotuhi, X. Zheng, and R. E. Ideker
Left Ventricular Apex Ablation Decreases the Upper Limit of Vulnerability
Circulation, May 30, 2000; 101(21): 2458 - 2460.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Idriss, S. F.
Right arrow Articles by Ideker, R. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Idriss, S. F.
Right arrow Articles by Ideker, R. E.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online