|
|
||||||||
1Woolcock Institute of Medical Research, Camperdown, 2Cooperative Research Centre for Asthma, Camperdown, NSW 2050; and 3Department of Medicine, University of Sydney, NSW 2006, Australia
Submitted 4 December 2003 ; accepted in final form 8 July 2004
| ABSTRACT |
|---|
|
|
|---|
respiratory resistance; airway caliber; power law; velocity of shortening
In addition to the phasic variations in resistance related to the breath cycle, reductions in resistance caused by DI can persist for some time after the volume has returned to normal levels (11). To the extent that the measured resistance reflects airway caliber, this gradual increase in resistance may be said to represent delayed "renarrowing" of the airways after DI-induced stretch.1 Jensen et al. (19) suggest that the delayed renarrowing effect seen in some normal subjects, particularly after challenge, may be due to the temporal dynamics of airway smooth muscle (ASM). They propose that in asthmatic airways the ASM is in a "latch state" (15), which cannot be broken by the DI compared with normal airways where the stretch is able to break the actin-myosin cross bridges in the ASM, which then take some time to reform. In other experiments, Stephens and coworkers (10, 32) have shown that the velocity of shortening of ASM is comparable to the breathing rate and may be altered under certain conditions. Delayed renarrowing in some subjects after DI-induced stretch may therefore be a direct consequence of reduced ASM velocity of shortening restricting the ability of the airways to renarrow, and this could also be a determinant in the amount of shortening that occurs after tidal stretches (30). Indeed, in a recent study (29), our laboratory showed that the magnitude of the delayed renarrowing effect after DI was correlated with baseline tidal fluctuations in resistance, suggesting that similar mechanisms may underlie both effects.
The present study extends the work of both Salome et al. (29) and Jensen et al. (19) by modeling the time course of resistance change after DI and tidal inspirations by an explicit function of time. The underlying hypothesis is that the pattern of renarrowing is determined by a time-related constant that is the same for both small and large stretches. Provided that the slope of the renarrowing function is smaller than the equivalent rate of volume decrease during expiration, the delayed renarrowing effect will be the limiting factor in the corresponding resistance fluctuations. To test this hypothesis, the model was applied to recordings of resistance during tidal breathing and after DI obtained from normal and asthmatic subjects after methacholine challenge. The results show that the coefficient characterizing the post-DI recovery function is directly related to the equivalent coefficient representing resistance changes after each tidal inspiration. The mechanisms that determine the recovery from DI are therefore likely to be the same as those operating during quiet breathing and possibly directly related to dynamic properties of ASM such as its velocity of shortening.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Postmethacholine resistance was recorded during 1 min of quiet breathing, followed by a single, slow DI to total lung capacity followed by passive expiration to FRC, and a subsequent further minute of quiet breathing. Subjects were seated, supporting their cheeks with their hands, and wore nose clips.
Resistance calculation.
We measured respiratory system impedance at 6 Hz using an in-house-developed device that followed the general design of earlier forced oscillation systems where the subject can breathe through a resistive exhaust port while the oscillations are applied. Flow was determined by means of a 50-mm-diameter Fleisch pneumotach with the differential pressure drop measured by a ± 2.5 cmH2O silicon transducer (Sursense DCAL-4, Honeywell Sensing & Control). Mouth pressure was measured by a similar transducer, with a range of ± 12.5 cmH2O. The resistance of the apparatus itself was
1 cmH2O·l1·s1. The mouth pressure and flow signals were low-pass filtered at 25 Hz and sampled at 300 Hz with a 16-bit analog-to-digital converter (DT-9805, Data Translation). All subsequent filtering and processing was performed digitally by custom software written in the Matlab (The Mathworks) digital signal processing environment. The pressure and flow signals were filtered by a band-pass filter centered on the oscillation frequency (6 Hz) to produce oscillation-frequency pressure [Pos(t)] and flow [
os(t)], respectively. The bandwidth of this filter determines the temporal smoothing that is imposed on the flow and pressure signals and, hence, the ability to detect rapid changes in resistance. For the results presented here, a bandwidth of 2 Hz was used. The raw flow signal was separately low-pass filtered with a cutoff of 2 Hz to provide an estimate of the breathing flow signal [
b(t)]. The phase responses of the band-pass and low-pass filters were matched to preserve the temporal relationship between
os(t) and
b(t). The breathing volume signal [Vb(t)] was obtained by integration of
b(t). Drift correction was applied by subtracting from Vb(t) a straight-line fit through all the FRC points in the recording.
Following the approach previously described by Horowitz et al. (17), the time-varying respiratory impedance was computed from Pos(t) and
os(t) by division in the frequency domain of segments equal in length to the period of the oscillation signal (1/6 s). The impedance was calculated at a rate of 12 times/s by overlapping the segments. To provide an estimate of confidence in the impedance values, the coherence [
(t)] was also computed from three overlapping segments of Pos(t) and
os(t). Values of the time-varying respiratory impedance for which
2(t) was <0.95 were flagged as unreliable impedance estimates (12). For the analyses presented here, we retained only the resistance corresponding to the real part of the impedance.
The end-expiratory instants were automatically identified by locating the zero-crossings (tz) in the flow signal
b(t) for which the volume Vb(tz) was at local minima. To exclude possible false detection of breaths (i.e., zero flow points other than at breath end points), those zero crossings associated with volume excursions of less than one-eighth of the median tidal breath volume during the recording were excluded from the set of end-expiratory instants. The end-inspiratory instants were then defined at the points of local maxima in volume Vb(t) between every adjacent pair of end expirations.
Resistance model.
Following previous studies on volume and flow effects on airway resistance (27), we modeled the instantaneous intrabreath variations in resistance (see Fig. 1) by linear coefficients on breathing flow and volume
![]() | (1) |
b(tz) and Vb(t) are zero, Kfi and Kfe represent flow-related nonlinearities for inspiratory (
bi) and expiratory flow (
be), respectively, and Kv represents the volume-related decrease in resistance between end expiration and end inspiration.
|
![]() | (2) |
![]() | (3) |
|
Because the DI and subsequent renarrowing effects alter the direct relationship between resistance and volume, the coefficients for the instantaneous model were obtained from the breathing signal only up to and including the inspiratory portion of the DI.
Although we have defined the renarrowing function in terms of the slow increase in resistance often observed after a DI, an equivalent reasoning can be applied to the increase in resistance occurring after every breath (30) during tidal breathing. The implication of this reasoning is that the change in resistance from each end inspiration to the following end expiration, after the effects of flow nonlinearities are excluded, is primarily by a function of time rather than being determined by the change in volume. We therefore define a tidal renarrowing coefficient [KRN(QB)] based on the (non-flow related) resistance changes during expiration in tidal breathing. We calculated KRN(QB) by a least squares fit (in the log-log domain) of Eq. 3 to the resistance change during each individual expiration before the DI, after removing the expiratory flow effect Kfe described in Eq. 1. The values of KRN(QB) so obtained were compared with the equivalent values of KRN obtained from the post-DI data by calculating the intraclass correlation.
| RESULTS |
|---|
|
|
|---|
60 s) of post-DI data. Anthropometric data and lung function at baseline and after methacholine challenge are shown in Table 1. The asthmatic and nonasthmatic groups had similar anthropometric characteristics, except that there were more men in the asthmatic group (9:8 compared with 3:11). No subjects were current smokers, and only two asthmatic subjects had a significant smoking history (>10 pack-yr). Spirometric function was normal in the nonasthmatic subjects but significantly reduced (FEV1 = 82.9 ± 11.7% predicted) in the asthmatic subjects. None of the nonasthmatic subjects had a history of respiratory symptoms or regular respiratory medication use, whereas all asthmatic subjects had a history of symptoms and were receiving medication for their asthma.
|
Instantaneous model fit. Figure 1 shows an example of resistance recorded over a period of 1 min of quiet breathing, indicating the variability that typically occurs in conjunction with tidal breathing flow and volume. In this example, breath volume and resistance are shown as functions of time in Fig. 1, A and B, respectively, with Fig. 1, C and D, showing resistance as a function of breathing flow and volume, respectively. The instantaneous model represents these flow- and volume-related variations in resistance, indicated by the slopes of the straight lines drawn in Fig. 1, C and D. The coefficient values obtained at baseline and after methacholine challenge are detailed in Table 2. As indicated in Fig. 1C, the slope of the resistance nonlinearity differs between positive and negative flow (P < 0.01 at baseline and P < 0.001 post-methacholine challenge). The values of K0 and Kfi both increased significantly after methacholine challenge for both the nonasthmatic (P < 0.0001) and asthmatic (P < 0.005) groups. The change in Kfi after methacholine challenge was also significant (P < 0.01) for the nonasthmatic group after normalization by K0. The volume coefficient Kv (see Fig. 1D) increased slightly after methacholine challenge, but this was only significant in the asthmatic group (P < 0.001). The resistance predicted by this model is shown by the dotted line in Fig. 1B. On average, across the 14 nonasthmatic subjects, the instantaneous model accounted for 61 ± 21% of the variability in the resistance through each breath cycle during tidal breathing at baseline and 65 ± 17% of the variability after methacholine challenge. In the asthmatic subjects, the model accounted for 69 ± 14% of the variance at baseline and 68 ± 18% after methacholine challenge. There was a slightly worse fit (50 ± 8%) for the four nonasthmatic subjects with a Ros,mean of >140% predicted. Compared with the subjects with normal resistance, these subjects exhibited greater expiratory flow and volume effects, which contributed to the high mean resistance. For the asthmatic subjects, however, there was no difference in how well the model fit the subjects with low or high mean resistance.
|
Delayed renarrowing post-DI. Figure 2 shows an example of delayed renarrowing after DI-induced stretch in one of the subjects after inhalation of 100 µmol of methacholine. As indicated in Fig. 2A, the volume returned to the pre-DI levels within two breaths, but the resistance, shown in Fig. 2B, only gradually recovered over a period of about seven breaths and so no longer followed a simple volume and flow relationship, as indicated in Fig. 2, C and D. In addition, there appears to be a reduction of the intrabreath fluctuation in resistance during the recovery phase. This delayed renarrowing effect was modeled by including a time-dependent memory term in the resistance model so that the DI-induced resistance decrement corresponding to the volume stretch increases only slowly with respect to the subsequent breath pattern. The functional form of the renarrowing memory term then determines the overall pattern and slope of the post-DI resistance recovery. In Fig. 3, the resistance at each end-expiratory instant after the DI is shown on linear, log-linear, and log-log graphs with respect to the time after DI. As these three graphs indicate, the time course of resistance change after DI in this individual appeared to follow a power-law (log-time vs. log-amplitude) rather than an exponential (log-amplitude vs. linear-time) or linear relationship. Examples of the log-log fit for other subjects are shown in Figs. 4 and 5 for the nonasthmatic and asthmatic groups, respectively.
|
|
|
The mean value of KRN over the 10 subjects was 0.20 ± 0.10, and the intercept term t0 had an average value of 0.28 ± 0.42 s, but with a wide range from 0.004 to 1.3 s and a median of 0.09 s. For the four subjects in whom Ree returned to its pre-DI value within one or two breaths only, a renarrowing coefficient was computed by fixing the intercept t0 at 0.09 s, which was the geometric mean and median value of t0 for the other 10 subjects. The inclusion of these subjects changed the average value of KRN to 0.23 ± 0.11. There is a high correlation between KRN and log t0 (r2 = 0.74), but if t0 is fixed at 0.09 s for all 14 subjects the average value of KRN changed only slightly to 0.21 ± 0.09. The mean square error of fits (in effect including the DI value) was 0.06 ± 0.06, which was lower than the equivalent exponential model for all except one subject, and despite the wide range this difference is also significant (P < 0.05).
In the asthmatic subjects, the post-DI renarrowing coefficient KRN was 0.17 ± 0.011 in the 12 subjects for whom the post-DI renarrowing took more than two breaths (average 5.6 ± 2 breaths) until Ree reached its pre-DI level. However, due to apparent noise in the post-DI resistance estimates, some of these fits were poor, with an average r2 value of only 0.50 ± 0.29. By excluding outliers from the post-DI sequence of Ree, we were able to improve the fit of the model such that the average r2 value became 0.88 ± 0.11 over an average of 4.9 ± 2 breaths post-DI in 12 subjects. The average KRN was then 0.24 ± 0.08 with t0 = 0.59 ± 0.93 s (median 0.13 s, range 0.0013.2 s). These values were not, however, significantly different from corresponding values for nonasthmatic subjects.
Renarrowing during tidal breathing. Figure 6B shows an example of the resistance time course during tidal breathing, as indicated in Fig. 6A, with the flow effect (i.e., second 2 terms in Eq. 1) removed. In Fig. 6C, all end-inspiratory to end-expiratory increments are overlaid together with the renarrowing function obtained from the post-DI resistance, as defined in Eq. 3. In this example, it is clear that the time course of renarrowing after each tidal inspiration is quite similar to the time course after the DI, albeit with some degree of interbreath variability. The mean value of the tidal renarrowing quiet breathing coefficients [KRN(QB)] over all nonasthmatic subjects was 0.21 ± 0.14, with a range of 0.070.50. The corresponding quiet breathing values of t0 [t0(QB)] had a mean of 0.23 ± 0.10 s, with a range of 0.080.44 s. In the asthmatic subjects, the mean value of KRN(QB) was 0.17 ± 0.04, with corresponding values of t0(QB) being 0.22 ± 0.14 s and a range of 0.0080.47 s.
|
|
| DISCUSSION |
|---|
|
|
|---|
b(t) and Vb(t). Previous studies have shown that much of the periodic variation can be associated with breathing-related effects such as nonlinear flow resistance (1, 4, 7, 23), phasic changes in upper airway caliber between inspiration and expiration (8, 20), and the effects of volume stretch on the airway caliber (2, 6, 22). The model of flow and volume effects presented here explained a large proportion of the variation in oscillation resistance. However, such instantaneous effects cannot explain the slow increase over several breaths that can often be observed after DI, particularly in nonasthmatic subjects after a high-dose challenge (19). We modeled this delayed effect by a simple function of time and showed that the recovery followed a power-law (log time) rather than an exponential (linear time) type of decay. This result has important implications with respect to the possible mechanisms determining this slow renarrowing after stretch. Several explanations for the delayed effect of volume stretch have been proposed (35), including changes in the nature of the bridge dynamics in ASM after a larger than normal stretch (13) and the possibility that the velocity of shortening of the ASM may be slow enough in some cases that it takes several breaths for the muscle to shorten after a stretch (30, 32). In the perturbed equilibrium model (13), the force that the ASM is able to generate during cyclic stretching is greatly reduced from the maximum static force because a proportion of the actin-myosin cross bridges becomes detached by the applied stretches. Furthermore, after particularly large stretches, the generated force remains depressed for some time (13). Jensen et al. (19) suggest that, after a DI, the actin-myosin cross bridges that are detached by the stretch take some time to reattach, thus providing a reduced force on the airways and explaining the prolonged reduction in resistance after DI in normal subjects. In asthmatic subjects, however, the ASM is hypothesized to be in a "latch-bridge" state with stronger actin-myosin bonds, implying both that the ASM is stiffer and so not able to be stretched as much by the DI and also that the bonds are less easily broken, so that after the DI is complete, the airway returns immediately to its pre-DI state.
A power-law function implies that there is not a single time constant but rather that the response is produced by the combined effect of a multitude of interacting elements with a wide range of time constants. In the lung, such a mechanism could occur at several scales, ranging from interaction within the cytoskeleton of constituent cells (9), actin-myosin interactions within the ASM (14), hysteresis in the parenchymal forces acting on the airways (16), stress adaptation in lung tissue after a step-length change (25), or even heterogeneity in the airway tree arising from airway closure and reopening (33)
The dynamics of ASM may be well modeled by power-law functions because the active elements of the muscle effectively consist of loosely coupled molecular chains. The overall time course of narrowing is therefore determined by the combined effect of each of the individual molecular links, which can be modeled as probabilistic processes. It turns out that a combination of many such loosely coupled interacting elements gives rise to power-law behavior (9) in which there is no single time constant. Indeed, although our data are at a different level of structure in the organism, the value of our power-law coefficient KRN = 0.2 is quite similar to the equivalent coefficient (x 1 = 0.2) obtained by Fabry et al. (9) in their observations of the dynamic characteristics of cytoskeletal mechanics in human ASM cells. Peslin et al. (25) also obtained a similar power-law coefficient of 0.17 in their analysis of pressure change after a step increase in lung volume.
We did not find any significant differences in the renarrowing coefficients obtained from asthmatic and nonasthmatic subjects, even though the two groups differed significantly in their lung function. Indeed, both groups contained a wide (and overlapping) spread of values, with at least some asthmatic subjects (see Fig. 5, AD) showing slow renarrowing over five or more breaths and several nonasthmatic subjects renarrowing fully within only one or two breaths (see Fig. 4, E and I). This result suggests that there may not be a dichotomous difference between asthmatic and nonasthmatic ASM, contrary to theories such as those that contrast latch-bridge vs. cross-bridge states as explanations for the differences in post-DI behavior between asthmatic and nonasthmatic subjects (15, 19). Given the wide range of values for both groups, the dynamics of airway renarrowing are evidently not critical in determining whether an individual exhibits asthmatic symptoms. It is possible that the power-law pattern of renarrowing is a result not only of ASM behavior but also of other factors, such as stress adaptation in passive lung tissue (25). Indeed, Fabry et al. (9) found power-law renarrowing behavior in a wide range of different cell types.
Renarrowing during quiet breathing. Although the delayed effects of slow renarrowing are most easily seen after DI, our results show that the change in resistance (after accounting for nonlinear flow effects) after tidal inspirations is also well represented by the same coefficient that determines the post-DI shape. Several authors [e.g., Stephens et al. (32) and Solway and Fredberg (30)] have proposed that the finite velocity of shortening of activated ASM may limit the rate at which airways can narrow during tidal breathing. Indeed, measurements of the velocity of shortening in vitro suggest that it is roughly of the same time scale as breathing frequencies (32). According to this hypothesis, the same mechanism would operate after both small and large stretches, with the rate of subsequent renarrowing determined by the speed at which the ASM can shorten (if this is slower than the rate at which the applied volume stretch is removed). This suggests that tidal fluctuations in airway caliber may be limited by the magnitude of the velocity of shortening (30) in the same way that the recovery after a DI is limited at a relatively slow rate in some subjects. In other words, if the velocity of shortening is slower than the breathing rate, then the ASM does not have enough time to fully renarrow during the expiratory phase of the breathing cycle, and the tidal fluctuations are reduced from what they would be if there were no limit on the rate of renarrowing. Our results, at least in the nonasthmatic subjects, tend to support this hypothesis, because the rate at which the resistance increased after each tidal stretch could be well modeled by a time-dependent renarrowing function that was almost the same as the renarrowing that could be observed after a DI. Also, there was no correlation between the renarrowing coefficient KRN and the volume stretch coefficient Kv, implying that the renarrowing during expiration was unrelated to the volume-induced stretch during tidal inspiration.
In the asthmatic subjects, the renarrowing coefficient obtained during quiet breathing was not correlated to that after DI and was smaller (although not significantly so) than that obtained for the nonasthmatic subjects. This result does not seem to support the hypothesis that KRN is a measure of velocity of shortening, particularly if we expect that to be greater in asthmatic subjects (30). However, it is possible that our model of time-dependent narrowing during quiet breathing does not fit the asthmatic subjects so well. In particular, our calculation of KRN(QB) assumes that the immediate volume effect KvVb(t) in Eq. 1 does not affect the expiratory time course of resistance change. However, this is only true if volume decreases relatively fast after the inspiration (relative to the time over which we can measure the renarrowing). During a DI, this assumption is met, but during quiet breathing the time course of expiration is comparable to the time course over which we are trying to calculate KRN(QB). There is thus an upper limit on the value of KRN(QB), which we can measure, related to the rate of expiration. It may be necessary to measure the tidal renarrowing during artificially constrained breathing patterns, e.g., a short expiration followed by a pause.
Conclusions. In summary, we have presented a model that quantifies the time course of changes in airway resistance during periodic stretch applied to the airways during normal breathing as well as the change that occurs after a large-amplitude stretch due to DI. We have shown that the time course of slow renarrowing (over several breaths) after a DI can be well modeled by a power-law function of time and, furthermore, that the time course by which airways narrow after each tidal inspiration is similar (at least in nonasthmatic subjects) to the time course of narrowing after a DI. The analysis extends the earlier work presented by Salome et al. (29), who compared the relative changes in resistance between the DI and the subsequent breath to the tidal fluctuations during quiet breathing at baseline and found a high degree of correlation in both normal and asthmatic subjects.
The analysis and results presented here provide some evidence that the mechanism of delayed recovery after stretch is similar for both large and small inspirations, with a close similarity in the shape and time course of the renarrowing functions computed after DI and during quiet breathing. The results support the supposition of Solway and Fredberg (30) that temporal dynamics of airway mechanisms are important in maintaining airway caliber during breathing in normal subjects. They also suggest that it may be possible to make measurements in vivo of the temporal behavior of ASM. Such measurements could be made during quiet breathing and after DI, as in our data here, or by following a specific breathing pattern designed to better separate the instantaneous and temporal effects. Furthermore, it would be very worthwhile to utilize multiple-frequency forced oscillation technique to see whether there is any frequency dependence in the effect that could be associated with variations in airway tree heterogeneity (21).
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 In our description of resistance change after inspiration, it is convenient to assume that there is a direct relationship between the measured resistance and the underlying airway caliber through an inverse fourth power. Thus the decrease in resistance associated with inspiration may be termed "stretch," and the subsequent increase, whether during a single expiration or, in the case of gradual increase following DI, may be said to represent "renarrowing" of the airway caliber. Although we have chosen to present the data in terms of resistance, the description will, where necessary, make use of this terminology, which makes implicit assumptions about the relationship between measured resistance and the underlying airway caliber. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. G. Chapman, N. Berend, G. G. King, B. E. McParland, and C. M. Salome Deep inspirations protect against airway closure in nonasthmatic subjects J Appl Physiol, August 1, 2009; 107(2): 564 - 569. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Agrawal, S. K. Singh, V. P. Singh, E. Murphy, and I. Parikh Partitioning of nasal and pulmonary resistance changes during noninvasive plethysmography in mice J Appl Physiol, December 1, 2008; 105(6): 1975 - 1979. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. M. Slats, K. Janssen, A. van Schadewijk, D. T. van der Plas, R. Schot, J. G. van den Aardweg, J. C. de Jongste, P. S. Hiemstra, T. Mauad, K. F. Rabe, et al. Bronchial Inflammation and Airway Responses to Deep Inspiration in Asthma and Chronic Obstructive Pulmonary Disease Am. J. Respir. Crit. Care Med., July 15, 2007; 176(2): 121 - 128. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Diba, C. M. Salome, H. K. Reddel, C. W. Thorpe, B. Toelle, and G. G. King Short-term variability of airway caliber--a marker of asthma? J Appl Physiol, July 1, 2007; 103(1): 296 - 304. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. S. An, T. R. Bai, J. H. T. Bates, J. L. Black, R. H. Brown, V. Brusasco, P. Chitano, L. Deng, M. Dowell, D. H. Eidelman, et al. Airway smooth muscle dynamics: a common pathway of airway obstruction in asthma Eur. Respir. J., May 1, 2007; 29(5): 834 - 860. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Schweitzer, L. T. T. Vu, Y. T. Nguyen, C. Chone, B. Demoulin, and F. Marchal Estimation of the bronchodilatory effect of deep inhalation after a free run in children Eur. Respir. J., July 1, 2006; 28(1): 89 - 95. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. H. T. Bates, S. S. Wagers, R. J. Norton, L. M. Rinaldi, and C. G. Irvin Exaggerated airway narrowing in mice treated with intratracheal cationic protein J Appl Physiol, February 1, 2006; 100(2): 500 - 506. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |