SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Express

In-plane optical anisotropy of InAs/GaSb superlattices with alternate interfaces

Shujie Wu12, Yonghai Chen12*, Jinling Yu12, Hansong Gao12, Chongyun Jiang12, Jianliang Huang3, Yanhua Zhang3, Yang Wei3 and Wenquan Ma3

Author affiliations

1 Key Laboratory of Semiconductor Materials Science, Institute of Semiconductors, Chinese Academy of Sciences, P.O. Box 912, Beijing 100083, People’s Republic of China

2 Beijing Key Laboratory of Low Dimensional Semiconductor Materials and Devices, Institute of Semiconductors, Chinese Academy of Sciences, P.O. Box 912, Beijing 100083, People’s Republic of China

3 Laboratory of Nano-Optoelectronics, Institute of Semiconductors, Chinese Academy of Sciences, P.O. Box 912, Beijing 100083, People’s Republic of China

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2013, 8:298  doi:10.1186/1556-276X-8-298

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/8/1/298


Received:26 April 2013
Accepted:16 June 2013
Published:25 June 2013

© 2013 Wu et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The in-plane optical anisotropy (IPOA) in InAs/GaSb superlattices has been studied by reflectance difference spectroscopy (RDS) at different temperatures ranging from 80 to 300 K. We introduce alternate GaAs- and InSb-like interfaces (IFs), which cause the symmetry reduced from D2d to C2v. IPOA has been observed in the (001) plane along [110] and [1

1 ¯
0] axes. RDS measurement results show strong anisotropy resonance near critical point (CP) energies of InAs and GaSb. The energy positions show red shift and RDS intensity decreases with the increasing temperature. For the superlattice sample with the thicker InSb-like IFs, energy positions show red shift, and the spectra exhibit stronger IPOA. The excitonic effect is clearly observed by RDS at low temperatures. It demonstrates that biaxial strain results in the shift of the CP energies and IPOA is enhanced by the further localization of the carriers in InSb-like IFs.

Background

InAs/GaSb type-II superlattices (SLs) are a considerable interest in the application of middle and far infrared photodetection. These structures have broken-gap band alignment, which allows tuning optical and electronic properties by varying layer thickness [1,2]. As the InAs and GaSb share no common atoms (NCA) across the interface (IF), these IFs have to be controlled by both InAs-like, both GaSb-like or alternating InAs- and GaSb-like. Figure 1 illustrates a simplified ball-and-stick model of InAs/GaSb SL with lower GaAs-like and upper InSb-like IFs. This kind of CA/C’A’ zinc blende hetero-structures lost their ideal Td point-group symmetry along the [001] growth direction. C and A represent cation and anion, respectively. If SLs have only one type of IF such as C-A’ or C’-A, it exists a S4 rotation-reflection axis, the symmetry is described as D2d point-group symmetry. If SLs have both kinds of IFs alternately, the symmetry depends on the number of atomic monolayer (ML) of each components. SLs components with one or both odd numbers of atomic ML belong to the D2d point-group symmetry, with both even numbers of atomic ML are corresponding to C2v point-group symmetry. If the structure shares a common atom (CA) (A=A’ or C=C’), the IFs have a S4 rotation-reflection axis corresponding to the D2d point-group symmetry. It is supposed that C-A bonds lie in the (110) plane and A-C’ bonds are in the (1 1 ¯ 0) plane. When a beam of linear polarized light propagates along the [001] direction with its polarized direction parallel to the [110] or [1 1 ¯ 0] direction, it feels different chemical bonds. This kind of anisotropic-chemical-bond arrangement leads to in-plane optical anisotropy (IPOA) at the IFs, i.e., optical property of [110] and [1 1 ¯ 0] plane is different in the (001) plane. Exactly speaking, the IPOA of upper and lower IFs will cancel each other for the SLs with D2d symmetry. Although, it is hard to realize such perfect IFs by the growth process that has many uncontrollable factors, the weak IPOA is still well observed by reflectance difference spectroscopy (RDS) [3,4]. For the NCA SLs, it has been observed that the IPOA is very strong [5-8].

thumbnailFigure 1. Simple stick-and-ball model of InAs/GaSb SL with alternate GaAs and InSb IFs. The purple, blue, green, and brown balls denote In, As, Ga, and Sb atoms, respectively.

RDS is a very sensitive nondestructive optical detection technique for IPOA, which was invented by D.E Aspens [9]. This powerful tool is used to detect IPOA induced by strain, electric field, and atom segregation for bulk, surface, and IF. In this letter, we have measured the IPOA of (001) plane of InAs/GaSb SLs by RDS at different temperatures ranging from 80 to 300 K. In this experiment, two SL examples have different thickness of InSb-like IF. The spectra are ranging from 1.5 to 5.0 eV. In the spectra, the energies of main features are assigned to Γ (E0, E0+Δ0), Λ (E1, E1+Δ1), and other critical point (CP) interband transitions of InAs, GaSb, and the coupling of the components whereas the L, X, and Σ CP energies are complex and difficult to analyze. Table 1 shows a list of the CP energies of bulk InAs, GaSb, GaAs, and InSb [10]. Additional CP energies may be related to the IFs. The Λ CP energies are very sensitive to strain. The CP energies show red shift with the increasing temperature, which attributes to the enhancement of electron-phonon interaction and thermal expansion. The transitions show a clear exciton characteristic at low temperatures. Compared with sample A, the measured energies of Λ CPs show red shift for sample B and exhibit stronger IPOA. The red shift attributes to the increasing of average lattice constant. IPOA is enhanced by the further localization of carriers in InSb-like IFs.

Table 1. CP energies (in eV) of bulk InAs, GaSb, GaAs, and InSb measured by S.Adachi [[10]]

Methods

The InAs/GaSb SLs were grown on GaSb buffer layer, which is deposited on non-intentional doping GaSb (001) substrates by molecular beam epitaxy (MBE). The GaAs-like IFs were generated by employing As soaking after GaSb is deposited. The InSb-like IFs were formed by InSb deposition. Two samples have the same structure as 100 periods InAs (10 ML)/GaSb (8 ML) without capping layer. The difference of the two examples is only the thickness of InSb layer, 0.43 ML (sample A) and 1.29 ML (sample B), respectively. We used a Bede D1 high-resolution X-ray diffractometer to characterize structural quality of the samples. The lattice mismatch and one-period thickness can be predicted.

We measured the relative reflectance difference between [110] and [1 1 ¯ 0] in (001) plane, obtaining

Δr r = r 110 r 1 1 ¯ 0 r 110 + r 1 1 ¯ 0 , (1)

ranging from 80 to 300 K in a cryogenic Dewar bottle. In the RDS measurement, near-normal incidence reflectivity of two perpendicular directions was obtained in order to remove the influence of errors induced by optical components, averaging two spectra sample azimuth by 90°. The difference of dielectric functions ( Δε = ε 110 ε 1 1 ¯ 0 ) has a relation with Δr/r:

Δr r = ( α + ) Δε = ( αΔ ε r βΔ ε i ) + i ( αΔ ε i + βΔ ε r ) . (2)

Here, α and β are complicated functions of four refractive indices and the wavelength of light. Both the real and imaginary part of Δr/r are linear combinations of real and imaginary part of Δε[11]. The degree of polarization (DOP) is defined as P = M 110 M 1 1 ¯ 0 M 110 + M 1 1 ¯ 0 (M110 is the transition probability when light is polarized along [110] direction). Im(Δε) is proportional to ΔM, and Im(ε) is proportional to M. It can be deduced from the imaginary part of Δε and the imaginary part of ε: P = Im ( Δε ) 2 Im ( ε ) [12].

Results and discussion

Lattice constants of GaAs, InAs, GaSb, and InSb are 5.2430, 6.0173, 6.0959, and 6.8970 Å, respectively [13]. The lattice mismatch between InAs and GaSb is only 0.6%; however, that of GaAs/GaSb and InSb/GaSb are 8% and 6%, respectively. Inserting GaAs-like IFs equals to introduce compress strain for the SLs, while InSb-like IFs will result in tensile strain. Alternating GaAs- or InSb-like IF layers can compensate the lattice mismatch between InAs and GaSb by controlling the appropriate thickness of GaAs and InSb layers. If SLs are pseudomorphic-grown on GaSb substrate, the strains of GaAs, InAs, and InSb are determined by the substrate, which can be calculated by:

e z i = a i a i a i , e x i = e y i = 2 v i 1 v i e z i . (3)

e z i , e x i , and e y i are the strains of GaAs, InAs, and GaSb for directions parallel and perpendicular to the growth direction, respectively. asub, ai, and a i represent crystal constants of GaSb substrate, for each layer, and the layers of SLs after growth, respectively. vi is the Possion ratio. The band gap and energies of CPs will show blue or red shift for compress or tensile biaxial strain, respectively. The two SL samples have the same thickness of GaAs-like IFs and different thickness of InSb-like IFs. The average lattice constant of superlattice is increased as a result of red shift energies of the CPs.

X-ray diffraction (XRD) results indicate that the range of 0th peak of sample A and the substrate is 0.367° and 0.151° for sample B. The full width at half maximum (FWHM) of the first satellite peak is 34 arcsec for sample A and 43 aresec for sample B. Both of the samples show compression strain. The calculated strain is -0.0054 for sample A and -0.0023 for sample B. Increasing the thickness of InSb-like IF layers can reduce the average compression strain. We predicted one-period thickness from the spacing between the satellites. Each period thickness of sample A is 55.9 Å and 56.8 Å for sample B.

Figure 2a,b shows the real parts of the relative reflectance difference measured at 300 and 80 K, respectively. The resonances of two samples have the same lineshape. In the spectra, the sharp peak near 2.05 eV(CP1), which is related to E1energy of GaSb. The lineshape of real part is almost the derivative of the imaginary part. A small feature is observed at this region, which is coincidence that the InAs E1 and GaSb E1+Δ1energies are both near 2.50 eV(CP2). The InAs E1energy is a little larger than GaSb E1+Δ1 energy. Another feature is observed near 2.78 eV(CP3) corresponding to the critical point energy of InAs E1+Δ1. Two shoulder-like features were marked in Figure 2b on both sides of the sharp peak near 2.05 eV, which may be attributed to InSb-like IFs. The energy positions are near the E1 and E1+Δ1energies of bulk InSb, and it is more clearly shown in the 80-K measurement. However, the IPOA structures about GaAs are not observed. In comparison with sample A, it is observed that GaSb E1 and InAs E1+Δ1features show red shift for sample B, which attributes to the compensation of stress by increasing the thickness of InSb-like IF layer. It is anomalous that a blue shift peak is corresponding to InAs E1 and GaSb E1+Δ1. D. Behr et al. reported that it is complicated by inhomogeneity for E1 and transition of InAs and E1+Δ1 of GaSb [14].

thumbnailFigure 2. Real part of RD spectra of samples A and B measured at 300 and 80 K. (a) At 300 K. (b) At 80 K. The arrows indicate the CP energies.

For SL sample, reflectivity can be described by a three-phase model:

r = V 23 + V 12 exp ( i 2 φ 2 ) 1 + V 23 V 12 exp ( i 2 φ 2 ) ; (4)

with

V ij = n ~ j n ~ i n ~ j n ~ i , φ 2 = 2 π n ~ 2 d 2 λ , ()

where the indices i and j take the value 1, 2, and 3 for the substrate, SL layer, and air, respectively. n ~ i is complex refractive index of the ith layer, d2is the thickness of the SL layer, Λ is the wavelength of light in vacuum [15]. SL layer are treated as uniaxial medium, n ~ 2 is the weighted average refractive index of 100 periods of InAs (10 ML)/GaSb (8 ML) SL layer. We chose a simple three-phase model, with no capping layer:

Δr r = 4 π idΔε λ ( ε s 1 ) . (5)

εs is the dielectric function of GaSb substrate, d is the thickness of the superlattice, and Λ is the wavelength of light [16]. The εsdata of GaSb substrate is taken from Aspnes’ measurement [17]. Figure 3a,b shows the real and imaginary parts of anisotropy dielectric function Δε by Equation 5, respectively. The peaks and valleys in the imaginary anisotropic dielectric function spectra are corresponding to the CP energies. The imaginary part of ε is related to absorption, and the real part corresponds to transmittance properties.

thumbnailFigure 3. Calculated imaginary (a) and real (b) parts of Δ ε of samples A and B. The arrows indicate the CP energies.

Figure 4a,b shows the measured IPOA of samples at different temperatures ranging from 80 to 300K. Figure 5a shows the temperature dependence of measured CP energy positions. Figure 5b shows the reflectance difference intensity of CP1 as a function of temperature. The energies of CPs show blue shift, and the amplitudes increase with the decreasing of measured temperature. There are no additional peaks observed. All the observed features are corresponding to CP energies. This kind of IPOA is stable and not caused by defects accumulated on the IF. The shoulder-like CP energy features about InSb clearly show character at low temperatures. Compared with sample A, all the spectra measured at different temperatures indicate that the CP energy are positioned on the red shift with a stronger RD intensity for sample B. J.S. Hang has reported that the GaSb critical point energies shift with temperature, as described by the Varshni expression [18], while J. Kim described the InAs CP energies and temperature dependence as Bose-Einstein statics [19]. We use the Varshni empirical formula to fit the temperature dependence:

E ( T ) = E 0 α T 2 T + β , (6)

thumbnailFigure 4. Real part ofΔ r/r of samples A and B measured ranging from 80 to 300 K.

thumbnailFigure 5. Measured CP energies of samples A and B as function of temperature and RD instensity of CP1. (a) Measured CP energies of samples A (squares) and B (circles) as a function of temperature. The lines are the Varshni empirical formula fitting. (b) Temperature-dependent RD intensity of CP1.

where β is a constant (K), Eois the width of semiconductor band gap, α is a fitting parameter (eVK−1), and T is the temperature. Table 2 lists the Varshni coefficients of samples A and B. It is found that excitonic transitions have important contributions to E1 and E1+Δ1 transitions. For this kind of transitions along eight equivalent Λ axes 〈111〉 direction of the Brillouin zone, the FWHM of the spectra decreases with the temperature decreasing. Since the spin orbit interaction in the valence band is large, the E1 transition split into E1and E1+Δ1 transitions. Δ1 is approximately 2/3 of Δ 0 at the Brillouin zone center [20]. The symmetry reduction remove the degeneracy of the four equivalent bands of two sets. As mentioned above, Δr/r is related to Δε; therefore, the line shape also depends on the symmetry of CP [21]. One electron approximation cannot explain the lifetime broadening; thus, it is suggested that Coulomb interaction should be taken into consideration [22]. The sharpening of spectra with reduction temperature indicates that excitons associate with the E1 transition [23].

Table 2. Varshni parameters for temperature-dependence fitting CPs of samples A and B

Samples A and B are both with GaAs-like and InSb-like alternate IFs and even number of InAs and GaSb MLs. The SLs possess C2vsymmetry in the ideal condition. At successive IFs, if In-Sb bonds lie in the (110) plane, while In-As bonds lie in the (1 1 ¯ 0) plane. Linearly polarized light propagates along the (001) direction. When the polarized direction is parallel to [110] and [1 1 ¯ 0] directions, it feels different chemical bonds at IFs. As a result, the optical properties along the [110] and [1 1 ¯ 0] directions are different. In the RDS spectra, InSb features were not observed clearly in room temperature, since the features of E0, E1, and E1+Δ1CPs are very broadening with few ML [24]. This effect is identified as the spread of carrier wave function of the ultra-thin IF to surrounding layers. Figure 6a shows the ΔEcand ΔEv of unstrained GaAs, InAs, InSb, and GaSb system at Γ point [25,26]. E1and E1+Δ1take place along the Λ directions of the Brillouin zone where the valence and conduction bands are nearly parallel. The energy gap of L and Λ are nearly equal. We have inferred the band alignment of L point in Figure 6b. The reflectance peaks of L transitions are not observed, since these transitions are too weak or hidden in the Λ transition structures [22]. In Figure 6b, the Λ 1conduction band offset between InAs and GaSb is 0.234 eV, and the Λ 3valence band offset is 0.544 eV. The staggered band alignment of bulk materials imply that in every InAs/GaSb SL, there is a InAs-like conduction band minimum and GaSb-like valence band maximum. The Λ 3valence band of InSb is much higher than GaSb, and the Λ 1conduction band is much higher than InAs. The Λ 3valence band splits into Λ (4,5)and Λ 6since the spin-orbital interaction. The red lines show the Λ 6energy positions. The Λ 6band of InSb is higher than Λ (4,5)band of InAs. As the thickness of InSb layers is increasing from 0.43 to 1.29 ML, compared to sample A, the effect of quantum well structures is enhanced. More holes are localized in InSb layers. However, there is no such effect for the GaAs layer. The IPOA intensities of CP1, CP2, and the shoulder-like CP about InSb are increased. While the IPOA intensities of CP3 are decreased and the transition energy position of CP2 are anomalous, blue shift may attribute to the coupling of these states.

thumbnailFigure 6. Band alignments of InAs, GaAs, GaSb and InSb binary system. (a) At Γ point of Brillouin zone. (b) At L point of Brillouin zone. The red lines are the spin-orbital splitting energies at L point.

Conclusions

The IPOA of InAs/GaSb SLs with InAs-like and GaSb-like alternate IFs were observed by RDS. The main mechanism can attribute to the symmetry reduction to C2v. The increasing of InSb IFs’ thickness release the mismatch between the SL layer and substrate. The red shift of CP energies was observed. Meanwhile, the holes are further localized in the InSb IFs, leading to the intensities of IPOA further increased.

Abbreviations

IPOA: In-plane optical anisotropy; SLs: Superlattices; RDS: Reflectance difference spectroscopy; CP: Critical point; IFs: Interfaces; NCA: No common atom; ML: Monolayer; MBE: Molecular beam epitaxy; DOP: Degree of polarization.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SW carried out the analysis, did the measurements, and drafted the manuscript. YC conceived of the study and participated in its design and coordination. JY and HG participated in the design of the study. JY and CJ participated in the revision of the manuscript and discussed the analysis. JH, YZ, and YW prepared the samples and measured the quality by XRD. WM designed the structure and supervised the preparation of samples. All authors read and approved the final manuscript.

Acknowledgements

This work was supported by the 973 Program (2013CB632805, 2012CB921304 and 2010CB327602) and the National Natural Science Foundation of China (No. 60990313, No. 61176014, and No. 61290303).

References

  1. Sai-Halasz GA, Tsu R, Esaki L: A new semiconductor superlattice.

    Appl Phys Lett 1997, 30:651-653. OpenURL

  2. Smith DL, Mailhiot C: Proposal for strained type II superlattice infrared detectors.

    J Appl Phys 1987, 62:2545-2548. Publisher Full Text OpenURL

  3. Koopmans B, Richards B, Santos P, Eberl K, Cardona M: In-plane optical anisotropy of GaAs/AlAs multiple quantum wells probed by microscopic reflectance difference spectroscopy.

    Appl Phys Lett 1996, 69:782-784. Publisher Full Text OpenURL

  4. Chen YH, Yang Z, Wang ZG, Bo Xu, Liang JB: Quantum-well anisotropic forbidden transitions induced by a common-atom interface potential.

    Phys Rev B 1999, 60:1783-1786. Publisher Full Text OpenURL

  5. Krebs O, Voisin P: Giant optical anisotropy of semiconductor heterostructures with no common atom and the quantum-confined pockels effect.

    Phys Rev Lett 1996, 77:1829-1832. PubMed Abstract | Publisher Full Text OpenURL

  6. Krebs O, Rondi D, Gentner JL, Goldstein L, Voisin P: Inversion asymmetry in heterostructures of zinc-blende semiconductors: interface and external potential versus bulk effects.

    Phys Rev Lett 1998, 80:5770-5773. Publisher Full Text OpenURL

  7. Ivchenko EL, Toropov AA, Voisin P: Interface optical anisotropy in a heterostructure with different cations and anions.

    Phys Solid State 1998, 40:1748-1753. Publisher Full Text OpenURL

  8. Krebs O, Voisin P: Light-heavy hole mixing and in-plane optical anisotropy of InP−AlxIn1−xAs type-II multiquantum wells.

    Phys Rev B 2000, 61:7265-7268. Publisher Full Text OpenURL

  9. Aspnes DE, Harbison JP, Studna AA, Florez LT: Application of reflectance difference spectroscopy to molecular-beam epitaxy growth of GaAs and AlAs.

    J Vac Sci Technol A-Vac Surf Films 1988, 6:1327-1332. Publisher Full Text OpenURL

  10. Adachi S: Optical dispersion relations for GaP, GaAs, GaSb, InP, InAs, InSb, Alx, Ga1−x As, and In1−x Gax Asy P1−y.

    J Appl Phys 1989, 66:6030-6040. Publisher Full Text OpenURL

  11. Ye X-L, Chen YH, Wang JZ, Wang ZG, Yang Z: Determination of the values of hole-mixing coefficients due to interface and electric field in GaAs/Alx, Ga1−x As superlattices.

    Phys Rev B 2001, 63:115317. OpenURL

  12. Chen YH, Ye XL, Xu B, Wang ZG: Strong in-plane optical anisotropy of asymmetric (001) quantum wells.

    J Appl Phys 2006, 99:096102. Publisher Full Text OpenURL

  13. Vurgaftman I, Meyer JR, Ram-Mohan LR: Band parameters for III–V compound semiconductors and their alloys.

    J Appl Phys 2001, 89:5815-5875. Publisher Full Text OpenURL

  14. Behr D, Wagner J, Schmitz J, Herres N, Ralston JD, Koidl P, Ramsteiner M, Schrottke L, Jungk G: Resonant Raman scattering and spectral ellipsometry on InAs/GaSb superlattices with different interfaces.

    Appl Phys Lett 1994, 65:2972-2974. Publisher Full Text OpenURL

  15. McIntyre JDE, Aspnes DE: Differential reflection spectroscopy of very thin surface films.

    Surf Sci 1971, 24:417-434. Publisher Full Text OpenURL

  16. Kamiya I, Aspnes DE, Florez LT, Harbison JP: Reflectance-difference spectroscopy of (001) GaAs surfaces in ultrahigh vacuum.

    Phys Rev B 1992, 46:15894-15904. Publisher Full Text OpenURL

  17. Aspnes DE, Studna AA: Dielectric functions and optical parameters of Si, Ge, GaP, GaAs, GaSb, InP, InAs, and InSb from 1.5 to 6.0 eV.

    Phys Rev B 1983, 27:985-1009. Publisher Full Text OpenURL

  18. Hwang JS, Tyan SL, Lin MJ, Su YK: Studies of interband transitions and thermal annealing effects on ion-implented (100) GaSb by photoreflectance and Raman spectra.

    Solid State Commun 1991, 80:891-896. Publisher Full Text OpenURL

  19. Kim TJ, Hwang SY, Byun JS, Barange NS, Kim JY, Kim YD: Temperature dependence of the dielectric function and critical-point energies of InAs.

    J Korean Phys Soc 2012, 61:97-101. Publisher Full Text OpenURL

  20. Cardona M, Christensen NE, Fasol G: Relativistic band structure and spin-orbit splitting of zinc-blende-type semiconductors.

    Phys Rev B 1988, 38:1806-1827. Publisher Full Text OpenURL

  21. Welkowsky M, Braunstein R: Interband transitions and exciton effects in semiconductors.

    Phys Rev B 1972, 5:497-509. Publisher Full Text OpenURL

  22. Zucca RRL, Shen YR: Wavelength-modulation spectra of some semiconductors.

    Phys Rev B 1970, 1:2668-2676. Publisher Full Text OpenURL

  23. Lautenschlager P, Garriga M, Vina L, Cardona M: Temperature dependence of the dielectric function and interband critical points in silicon.

    Phys Rev B 1987, 36:4821-4830. Publisher Full Text OpenURL

  24. Weimar U, Wagner J, Gaymann A, Köhler K: Broadening of interband resonances in thin AlAs barriers embedded in GaAs.

    Appl Phys Lett 1996, 68:3293-3295. Publisher Full Text OpenURL

  25. Wei S-H, Zunger A: Calculated natural band offsets of all II–VI and III–V semiconductors: chemical trends and the role of cation d orbitals.

    Appl Phys Lett 1998, 72:2011-2013. Publisher Full Text OpenURL

  26. Magri R, Zunger A: Effects of interfacial atomic segregation on optical properties of InAs/GaSb superlattices.

    Phys Rev B 2001, 64:081305. OpenURL