|
|
||||||||
EDITORIAL FOCUS
1California National Primate Research Center, 2Department of Psychology, and 3Department of Internal Medicine, Division of Pulmonary and Critical Care, University of California, Davis, Davis, California
Submitted 1 December 2006 ; accepted in final form 8 June 2007
| ABSTRACT |
|---|
|
|
|---|
stereology; parenchyma; alveolar ducts; postnatal development
36 wk of gestation to
1–2 yr of age and microvascular maturation by remodeling of interalveolar septa and restructuring of the capillary bed from birth to 2–3 yr of age (5, 40). The majority of alveoli are produced postnatally in humans to reach the adult number of
450 million alveoli (35). One morphometric study of lungs of children from 26 days to
5 yr of age identified two phases of postnatal lung development and growth: 1) from birth to
18 mo of age characterized by alterations in volumetric proportions of parenchymal compartments, and 2) 18 mo to adulthood with proportional growth of all lung compartments (49). In the first phase, there is a disproportionate increase in components that are concerned with gas exchange (air space and capillary volumes) with a proportional decrease in interstitial tissue mass. Since many human infants are exposed to a variety of inhaled infectious diseases and irritants such as air pollutants and environmental tobacco smoke that damage the lung, it is important that we understand the nature of these insults on the developmental process of the lung. This need necessitates the use of laboratory animals that represent good models of human lung development in which we can study normal lung development and its potential perturbation. Hence, we used precise, design-based stereological methods to sample whole rhesus monkey lungs to establish the normal parameters of alveolar growth. The time of postnatal developmental stages of lung parenchyma in the rhesus monkey is about a third that of the human. Thus we hypothesize that the most rapid phase of alveolar development in rhesus monkeys will be within the first year of life; however, we included animals from 4 to 2,675 days (7 yr) because somatic growth is complete in rhesus macaques by 6 yr of age (7, 42).
| MATERIALS AND METHODS |
|---|
|
|
|---|
All monkeys selected for these studies were California National Primate Research Center colony-born rhesus macaques (Macaca mulatta). All monkeys were given a comprehensive physical examination, including a chest radiograph and complete blood count, and were determined to be healthy monkeys. Care and housing of animals complied with the provisions of the Institute of Laboratory Animal Resources and conformed to practices established by the Association for Assessment and Accreditation of Laboratory Animal Care (AAALAC). Animal studies conformed to applicable provisions of the Animal Welfare Act and other federal statutes and regulations relating to animals (Guide for the Care and Use of Laboratory Animals; National Institutes of Health, revised 1985). Experimental protocols were reviewed and approved by the University of California, Davis Institutional Animal Care and Use Committee. Twenty-six rhesus monkeys (13 males ranging in age from 4 to 1,920 days, body wt from 0.408 to 12.66 kg, and lung volumes from 41.7 to 602 cm3; 13 females ranging in age from 22 to 2,675 days, body wt from 0.472 to 3.9 kg, and lung volumes from 43.5 to 380 cm3; Table 1) were killed with an overdose of pentobarbital after being sedated with Telazol (8 mg/kg im) and anesthetized with Diprivan (0.1–0.2 mg·kg–1·min–1 iv) with the dose adjusted as necessary by the attending veterinarian. The monkeys were necropsied following exsanguination through the posterior vena cava. The lungs were fixed via the airways through a tracheal cannula with 1% glutaraldehyde-1% paraformaldehyde in cacodylate buffer (330 mosM, pH 7.4) at 30-cm fluid pressure (8 h). After fixation, the trachea was tied off at the fixation pressure and held in the same fixative between 1–4 wk at room temperature. Subsequently, lungs were embedded in 4% agar, isotropically oriented using an orientator (31), sliced into 5-mm slabs, and sampled using a smooth fractionator (19) for 10–12 5 x 5 x 15-mm blocks for histology (26; Fig. 1 and Table 2). This sampling approach provides isotropic uniform random (IUR) sampling of lung tissue. Blocks were embedded in paraffin, cut in 5-µm serial sections, and stained with hematoxylin and eosin or Orcein stain for elastin.
|
|
|
Lung volume. Lung volume (VL) without the trachea and extrapulmonary bronchi was estimated by its buoyant weight in PBS (36). The right cranial lung lobes from three additional monkeys (7 mo, 2 yr 9 mo, and 5 yr 4 mo) were estimated by their buoyant weight in PBS and by the Cavalieri method, a volume estimate of the sectioned (5-mm slabs) fixed lobe using point counting to estimate slab areas that were multiplied by the slab thickness to estimate volume (34).
Estimation of the Volume of Parenchymal Components
We used the Intelligent Imaging Innovations software Slidebook 4.1 (Santa Monica, CA) at a magnification of x185 to evaluate 10–12 sampled blocks per lung. Using stratified random sampling with a spacing of 2,200 µm (
65 fields per section) and a double lattice test system of 25/100 points, we evaluated the volume densities (Vv) of parenchymal and nonparenchymal components using point counting. Points on alveolar air and interalveolar septa, alveolar duct core air, and other parenchymal compartments (terminal and respiratory bronchioles) were counted relative to VL. The absolute volumes of parenchyma (Vpar), alveoli (Valv), interalveolar septa (Vias), alveolar duct core air (Vad), and nonparenchyma (Vnp) were determined by multiplying their Vv by VL in units of cm3.
Estimation of Alveolar Number
The method for counting alveoli is based on the mathematical concept of the Euler characteristics of structures (26, 35). This provides an estimate of the number of holes in a two-dimensional (2-D) net (alveolar opening rings). To estimate the number of alveolar opening rings, the number of bridges (B) and islands (I) that appear on one section of a "disector" pair but not the other are counted (see Fig. 5; Ref. 26). Bridges (the more frequent event) are connections that appear between two separate profiles in one section but not the other. Islands (the rare event) are a new isolated profile. The Euler characteristic 
= (I-B/2). The total number of alveoli (Nalv) in a lung (alv,lung) was calculated using the fractionator principle
![]() |
![]() |
|
|
An unbiased estimate of the global volume change is the following factor for shrunken global volume (Fsgv), Fsgv = (
SAafter/
BAbefore)3/2, where the block area before processing, BAbefore, is estimated from the photographs, and section area after shrinkage, SAafter, is estimated by point counting; the summation is over all blocks for an animal (Table 2). All estimated absolute surfaces and volumes on histological sections were corrected for global shrinkage by dividing with the above factor under the explicit assumption that shrinkage is uniform across all tissue components (12).
Estimation of Alveolar Surface
The surface of interalveolar septa/volume of lung (SvIAS,Lung) was estimated as SvIAS,Lung = 2 x IIAS/LLung, where IIAS is the number of intersections with interalveolar septa by a linear probe, and LLung is the total probe length in sections of lung tissue (mm2/mm3) or mm–1 (38). This equation is valid for test lines that are isotropic, uniform, and random in 3-D space (see Fig. 1; Ref. 46). The total surface of interalveolar septa within the lung (SIAS) was estimated as SIAS = SvIAS,Lung x [1/(Fsgv)2/3] x VL, where (Fsgv)2/3 is the factor for shrunken global volume raised to the2/3 power, and the units are converted to m2.
Calculation of Number-Weighted Mean Alveolar Volume
Number-weighted mean alveolar volume (
n alv) was calculated from previously described values as
n alv = VL x Vv alv,lung/Nalv,lung, where the units are converted to µm3 (see Fig. 1; Ref. 20). In the number-weighted mean alveolar volume, each alveolus has an equal statistical "weight" regardless of its volume.
Calculation of Volume-Weighted Mean Alveolar Volume
Volume-weighted mean alveolar volume (
v alv) was calculated using the point-sampled intercept method that estimates the volume of structures provided that the tissue is sampled under IUR conditions (see Fig. 1; Ref. 20). A CAST-Grid system was used to make the calculation as
v alv = (
/3) x
IAS3 x (1/Fsgv), where
IAS3 is the mean of the cubed point-sampled intercepts of alveoli, Fsgv is the factor for shrunken global volume, and the units are in µm3. In the volume-weighted mean alveolar volume, each alveolus is "weighted" by its volume. Consequently, when there is any distribution of alveolar sizes in a lung, the volume-weighted mean alveolar volume will always be greater than the number-weighted mean alveolar volume. It becomes evident that the volume-weighted mean alveolar volume has size information embedded in it that can be estimated.
Calculation of the Coefficient of Variation of the Distribution of Number-Weighted Alveolar Volumes
The coefficient of variation of the distribution of number-weighted alveolar volumes (CVn alv) was calculated from previously described values (20) as
![]() |
Variance and Efficiency of Stereological Estimators
The observed biological variation among individuals is large for features of interest in biological tissues, and it is useful to know whether it is worth increasing the precision of the stereological sampling or including more animals in the study (22). It is possible to divide the observed variance (OCV) into its two components, the true biological variation (CV) and the average sampling variation of the stereological measurement (CE) in the following equation for the number of features, N: OCV2 (N) = CV2 (N) + CE2 (N).
For sensible fractionator sampling designs (19), the dominating component of sampling variation is the counting noise. Because of the very sparse sampling, the counting noise was calculated by the formula CE2 noise (N) = 1/(
B +
I), where the summation is over one animal (Table 3). Note that both bridges (B) and islands (I) contribute to the counting noise of alveolar number estimation. Contributions to stereological variation for ratio estimators like volume, number, surface, and length densities have been derived (9). However, we used some simple guidelines that usually suffice for stereological sample size within an animal (primary sampling unit): 100–200 probe interactions (e.g., point hits, intersections, or feature counts), 50 fields, and 10 blocks (25).
|
Statistical analyses were based on linear and nonlinear regression (SAS Institute, Cary, NC). A series of models were fitted to each outcome, considering age, VL, and weight as predictors. Five functions were considered for each outcome: linear, quadratic, two-piece linear spline, a two-parameter exponential function, and a three-parameter exponential function. For example, when considering an outcome as a function of age, the models fitted to each outcome were
Linear: y =
0 +
1Age + e
Quadratic: y =
0 +
1Age +
2Age2 + e
Two-piece linear spline: y =
0 +
1Age + e if Age <=
, and y =
0 +
1
+
2(Age –
) + e if Age >
Two-parameter exponential: y =
0exp(–
1Age) + e
Three-parameter exponential: y =
0 – (
0 –
1)exp(–
2Age) + e.
Assessment of model fit was based on the coefficient of determination (R2). The test for R2 is whether it is statistically different from 0. For the linear and two-parameter exponential functions, the minimum R2 value is 0.23. For the quadratic and three-parameter exponential functions, the minimum R2 value is 0.21. The best fitting function for each outcome is reported. Age, VL, and body weight were centered to the sample means of 743 days, 190.7 cm3, and 2.92 kg, respectively. Sex was treated as a moderator of model parameters to allow for sex differences in all regression coefficients. Statistical significance was defined as P
0.05.
| RESULTS |
|---|
|
|
|---|
n alv was best fit by quadratic functions with age (Adj R2 = 0.13) and body weight (Adj R2 = 0.39), whereas with VL (Adj R2 = 0.32), a three-parameter exponential function provided the best fit, indicating little consistency in change with age or VL (Table 4; Fig. 9). However, in females,
n alv and body weight showed a significantly increased nonlinear change rate and an increased alveolar volume at the mean value of 2.92 kg compared with males (Table 4). It should be noted that the mean values for
n alv, the number-weighted mean alveolar volume, were 12% and for
v alv, the volume-weighted alveolar volume, were 29% greater in females than males (Table 5). Furthermore, CVn alv had two-parameter exponential functions with age (Adj R2 = 0.48), VL (Adj R2 = 0.48), and body weight (Adj R2 = 0.42) (Table 4; Fig. 10), indicating a greater distribution of alveolar volumes with age, VL, and body weight. CVn alv was also 14% greater in females than males, an indication of a greater distribution of alveolar volumes in females than males (Table 5). Furthermore, for females, CVn alv and body weight showed a significantly increased linear change rate and an increased alveolar volume distribution at the mean value of 2.92 kg compared with males (Table 4). As expected, the global shrinkage because of dehydration and paraffin embedding was very pronounced: the lung tissue shrank to 51% of the fixed volume of the monkey (Table 2).
|
|
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
n alv showed little relationship with age, volume, or body weight, and in females,
n alv was larger, and alveoli showed a greater size distribution of volume than in males. Alveoli increase in number but not volume throughout all of the postnatal developmental/growth stages (infant, 1–12 mo; juvenile, 12–24 mo; adolescent, 2–4 yr; and young adult, 4–8 yr) in rhesus monkeys. Oxygen Diffusion and Lung Parenchyma
Oxygen exchange between air and blood occurs by diffusion across the air-blood barrier in the lung parenchyma. Oxygen flow by diffusion from air to blood is governed by Fick's law that has a permeability coefficient, surface, and thickness of the barrier between air and blood (16). A morphometric model of diffusion capacity that was proposed by Weibel (43) and subsequently refined by Weibel and colleagues (44) showed a direct correlation of diffusion capacity to alveolar surface and an indirect correlation to thickness of the blood air barrier. The morphometric estimate of diffusion capacity of oxygen in monkey lungs was very similar to dogs and scaled linearly with body weight (17). Additional studies using the morphometric estimate of diffusion capacity of oxygen in monkey lungs concur with these initial findings (1, 27). When Gehr and colleagues (17) compared the scaling of morphometric parameters (like alveolar surface area) with maximal oxygen consumption, they concluded that bigger animals required a larger pulmonary diffusion capacity to admit the flow of oxygen required by the organism. We would expect the increase in surface area in monkeys during the first 7 yr of life to show a proportional increase in pulmonary diffusion capacity even though there is no data on the maximal oxygen uptake in rhesus monkey lungs during postnatal development. It is noteworthy that membrane diffusion capacity and capillary blood volume showed an age-related increase consistent with alveolarization from birth to 8 wk of age in lambs (10).
Lung Growth Differences Between Sexes
In humans, female lungs tend to be smaller than male lungs throughout childhood and even in adolescence when girls attain greater height earlier than boys. As girls attain their maximum adult height in late puberty, their lung growth ceases, whereas that of boys continues longer, in some cases into early adulthood (2). In boys and girls, the growth of the lung parenchyma and its airways occurs independently, but this dysanapsis is more pronounced in boys (2). The configuration of the adult female lung is the result of proportional growth of its airways in relation to its parenchyma, but that of the adult male lung is the result of dysanaptic growth where growth of the airways has lagged behind that of the lung parenchyma. Our measurements of parenchymal volume in rhesus monkey lungs show very similar trends over the first 7 yr of life that correspond to growth into early adulthood in humans. In rodents, females have smaller alveoli and more alveoli and alveolar surface area per body weight than males (28, 29). Our observations in rhesus monkey lungs show similar results to those seen in rodents with more alveoli and alveolar surface area per body weight in females compared with males. However, rhesus monkeys have larger alveoli per body weight in females compared with males, a result opposite that observed in rodents. Comparison of body weights in rhesus monkeys between sexes is complicated by the greater musculoskeletal growth and mass in males compared with females that occurs between 4–9 yr of age (37). Furthermore, in our data, VL was the best independent variable for the prediction of alveolar number and surface area in both sexes (Table 4).
Estimation of Alveolar Number
Studies of Macaca fascicularis newborn (7 at 5 h), infant (2 at 14 days), and adult (6 mature females) lungs revealed that there was no increase in alveolar number during the first 14 days of life, and the alveolar number in adult lungs suggested that there was little alveolar multiplication after birth (24). Thus these results from the study of M. fascicularis are different than our study of M. mulatta as well as the majority of investigations of human postnatal lung development (40, 49). Studies of human lungs show a rapid alveolar growth phase from 36 wk of gestation to
1–2 yr of age with formation of secondary interalveolar septa (40, 49). A stage of late alveolarization in the human lung has been suggested (6). Although controversial in humans, early studies proposed that alveoli increase in number until 20 yr of age (14, 15), whereas more recent studies proposed an age of 8–11 yr (11, 13) or even 2 yr of age (40, 49). Therefore, the exact time at which alveolar multiplication ceases in humans is still obscure. Differences between studies may lie in the stereological approach used to estimate alveolar number. Traditionally, the estimation of alveolar number has been done by assumption of a specific geometric shape, a rotatory ellipsoid (45). Human acinar reconstructions showed that geometric assumption-based estimates of alveolar number provided underestimates of true numbers (23). Their acinar reconstructions identified six different geometric shapes for alveoli and thereby illustrate the difficulty in selecting one geometric shape as a mean of the representing all six shapes.
Recently, we (26) introduced an approach to alveolar counting that used the smooth fractionator (19), a rigorous design-based sampling approach to the lung and the Euler characteristic to estimate the number of alveoli in lung without bias. Estimation of the total number of any feature in an organ or any containing space with the fractionator is direct, and there is no need to know the reference volume. This method is unaffected by global and differential shrinkage, swelling, and distortion of the containing space during embedding and sectioning. Euler number estimation uses physical disectors that are true volume probes (39). Euler number estimation of alveoli makes no assumption regarding the size, shape, or orientation of the structures to be counted in contrast to 2-D analyses. In essence, the Euler number count is directed toward counting the rings of alveolar mouth openings. Another approach that has been applied to rodent lungs estimates the number of alveoli in the lung indirectly by selecting individual alveoli using serial sections and a disector and then making volume estimates by Cavalieri (33) or point-sample intercept (30) methods and dividing the mean alveolar volume into the volume of alveoli in the lung. This indirect approach is extraordinarily time consuming and is complicated by the need to define alveolar mouth openings in three dimensions as a (curved) wall of the alveolar air space in serial sections. The alveolar counting procedure used in this paper is one of counting discontinuities of a well-defined structure, the alveolar wall. Except for the necessity of a well-trained observer who can concentrate on the ends of the interalveolar septal walls, our experience is that the counting is straightforward and unproblematic (26).
Sampling
We previously optimized sampling for alveolar estimation in adult rhesus monkeys (26). One refinement in this study was that we used a thinner slab thickness of 5 mm to increase our sampling in infant lungs but maintained a 5-µm section thickness. This allowed us to achieve the appropriate sample total of 100–200 Euler counts per lung (Table 3) (21). Although the sampling and processing steps were monitored and corrected for shrinkage necessary for estimates of mean alveolar volumes, the fractionator design for the total number of alveoli is independent of shrinkage (26). Since only the counting noise is easily known in fractionator designs, it must be kept low by a sensible design (19). The estimated CE noise of 0.05 (Table 3) for both males and females easily meets this requirement. The total observed CV of 0.77 and 0.73 in males and females, respectively, indicate that the remaining biological variability of alveolar number during postnatal development is very high. Of course, this is not surprising since we are dealing with a developmental change of
11-fold in alveolar number in both males and females over the first 7 yr of life in rhesus monkeys. We can have confidence in our individual animal estimates because of a value of 0.05 for CE noise.
Number of Alveoli
The general description of rhesus monkey body weights and lung volumes in this study are similar to previously published values for these species (7). A precise morphometric study of regional differences in rat lungs fixed by intratracheal instillation at 20 cmH2O pressure showed significant decreases in the volume and surface densities of interalveolar septa in subpleural compared with central lung regions (48). These investigators recommended that for quantitative light microscopic analysis of lung tissue, the most appropriate sampling unit is at minimum the entire lobe (48). When monkey lung lobes are fixed by intratracheal instillation at 30 cmH2O pressure, more variation in alveolar number was observed between cranial and caudal lobes than between the entire left lung of adult monkeys (26). In this study, we analyzed all six lobes of the rhesus monkey lung that were fixed by intratracheal instillation at 30 cmH2O pressure, a fixation pressure that mimics total lung capacity. Our use of smooth fractionator sampling followed by stratified random sampling of disectors on sections guaranteed an alveolar number that was independent of inflation fixation variation and variations in alveolar size. Our finding of the greater increase in lung volumes and alveolar surface area in males compared with females over all ages is very similar to that reported for postnatal human lung growth (40). The greater numbers of alveoli in males compared with females over all ages is also very similar to that reported for postnatal human lung growth (40). At what age alveolar multiplication ends in humans is still open to question, primarily because of the stereological methods used and the variability between individuals. Our data in rhesus monkeys clearly documents alveolar addition in the first 7 yr that has not been previously reported in nonhuman primates. However, the best predictor of alveolar number in a rhesus monkey lung is VL.
Alveolar Volume
Mean alveolar volume showed a poor relationship with age and VL, implying that with increases in VL, alveoli are added and do not enlarge to any significant degree during postnatal growth in rhesus monkeys. It is noteworthy that females had mean alveolar volumes that were 8% greater than males. CVn alv, a measure of the size distribution of number-weighted alveolar volume, showed a steady increase with age for both males and females. Furthermore, CVn alv was 40% greater in females than males. It is possible from the greater size distribution that the greater mean volume in females is the result of a subpopulation of larger alveoli in females that are not found in males. Because little is known about the development of lung parenchyma, most review articles on sex differences in lung development focus on the airways (4). Since this study is the first to measure CVn alv, more studies are needed to estimate the effect of increasing age in both sexes into adulthood (7–18 yr) and with old age (
18 yr in rhesus monkeys). Does this greater CVn alv during postnatal development in females persist into adulthood and even advance with old age? If so, then perhaps there is an anatomical basis to the compelling evidence that women are more susceptible than men to the development of chronic airflow limitation (8). Furthermore, women who have chronic obstructive pulmonary disease (COPD) seem to have a greater risk of hospitalization than men. The greater severity of disease among women is consistent with evidence of a greater predisposition to develop COPD and the appearance of COPD at an earlier age than men (8).
The use of the fixed VL as estimated by the immersed buoyant weight in PBS may have overestimated our tissue and air space volumes. A comparison of right cranial lung lobe volumes from three monkeys spanning the range of this study showed an overestimation of 4.9–16% by the buoyant weight method compared with the Cavalieri method. There was no apparent relationship with age as the percentage of difference was least in the youngest (6.5%) and oldest (4.9%) and greatest in the middle age range monkey (16%). A comparison of lung volumes by these two methods in dog lungs showed that the volume of the intact fixed lung under positive pressure is systematically higher by 13–25% than that measured after sectioning and release of airway pressure (47). Only lung volumes estimated by the Cavalieri method were used in the morphometric calculations of dog lungs (47). The lower shrinkage from immersion to the slab volume estimate in monkey compared with dog lungs could be the result of a smaller size or fixative difference. Further studies are needed to investigate the shrinkage differences of these two fixatives (2.5% buffered glutaraldehyde and 1% glutaraldehyde-1% paraformaldehyde in cacodylate buffer), but we strongly recommend that investigators report both volume estimators for future morphometric studies.
We have shown that alveoli are added in rhesus monkey lungs proportional to age, VL, and body weight into young adulthood in both males and females. Furthermore, we recommend that previous reports of postnatal development in nonhuman primates and humans be reconsidered because of the potential bias of sample, tissue shrinkage, and geometric assumption of alveolar shape.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| 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.
| REFERENCES |
|---|
|
|
|---|
ata V, Van
atová MA, Chalyan VG, Meishvilli N. Longitudinal study of growth and body mass changes in ontogeny in captive rhesus macaques (Macaca mulatta) from the Institute of Medical Primatology, Sochi. Variability and Evolution 8: 51–81, 2000.This article has been cited by other articles:
![]() |
M. Avdalovic, L. Putney, N. Tyler, W. Finkbeiner, K. Pinkerton, and D. Hyde In Utero and Postnatal Exposure to Environmental Tobacco Smoke (ETS) Alters Alveolar and Respiratory Bronchiole (RB) Growth and Development in Infant Monkeys Toxicol Pathol, February 1, 2009; 37(2): 256 - 263. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Castro, M. I. Ramirez, J. E. Gern, G. Cutting, G. Redding, J. S. Hagood, J. Whitsett, S. Abman, J. U. Raj, R. Barst, et al. Strategic Plan for Pediatric Respiratory Diseases Research: An NHLBI Working Group Report Proceedings of the ATS, January 15, 2009; 6(1): 1 - 10. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. Schittny, S. I. Mund, and M. Stampanoni Evidence and structural mechanism for late lung alveolarization Am J Physiol Lung Cell Mol Physiol, February 1, 2008; 294(2): L246 - L254. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |