Quantum Hall effect in graphene - CiteSeerX

13 downloads 0 Views 703KB Size Report
Apr 27, 2007 - Z. Jianga,b, Y. Zhanga, Y.-W. Tana, H.L. Stormera,c, P. Kima,∗ a Department of Physics, Columbia University, New York, NY 10027, USA.
Solid State Communications 143 (2007) 14–19 www.elsevier.com/locate/ssc

Quantum Hall effect in graphene Z. Jiang a,b , Y. Zhang a , Y.-W. Tan a , H.L. Stormer a,c , P. Kim a,∗ a Department of Physics, Columbia University, New York, NY 10027, USA b National High Magnetic Field Laboratory, Tallahassee, FL 32310, USA c Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027, USA

Accepted 27 February 2007 by A. Geim Available online 27 April 2007

Abstract The quantum Hall (QH) effect in two-dimensional electron and hole gas is studied in high quality graphene samples. Graphene samples whose lateral size ∼10 µm were fabricated into mesoscopic devices for electrical transport measurement in magnetic fields. In an intermediate field range of up to 10 T, a distinctive half-integer QH effect is discovered with QH plateaus appearing at a filling factor sequence, ν = 4(n + 1/2), where n is the Landau level (LL) index. As the magnetic field increases to the extreme quantum limit, we observe additional QH plateaus at filling factors ν = 0, ±1, ±4. Further detailed investigations show that the presence of the ν = 0, ±1 QH states indicates the n = 0 LL at the charge neutral Dirac point splits into four sublevels. This lifts both the sublattice and the spin degeneracy, while the QH states at ν = ±4 can be attributed to lifting of the spin degeneracy of the LLs. Above 30 T of magnetic field, the large quasiparticle gaps between the n = 0 and n = ±1 LLs lead to the QH effect that can be observed even at room temperature. c 2007 Elsevier Ltd. All rights reserved.

PACS: 73.63.-b; 73.21.-b; 73.43.-f Keywords: A. Graphene; D. Quantum Hall effect; D. Landau level splitting

1. Introduction Over the last twenty years there has been a tremendous interest in the electronic properties of low dimensional condensed matter systems (e.g. two-dimensional (2D) electrons/holes confined in semiconductor devices; mesoscopic heterostructures incorporating thin metallic, superconducting, and/or ferromagnetic films). The interest in these systems mainly stems from the fact that, as device dimensions approach fundamental microscopic length scales, novel quantum effects manifest themselves. The quantum electronic properties of such systems are commonly described by the Schr¨odinger equation, in which quasiparticles behave non-relativistically with a finite effective mass. Recently, great interest has arisen in another class of condensed matter system, where the charge carrier dynamics appears effectively relativistic with zero effective carrier mass and

∗ Corresponding author.

E-mail address: [email protected] (P. Kim). c 2007 Elsevier Ltd. All rights reserved. 0038-1098/$ - see front matter doi:10.1016/j.ssc.2007.02.046

the transport properties are governed by the Dirac equation [1– 3]. This fascinating and unique system is graphene: a single atomic sheet of graphitic carbon atoms that are arranged into a honeycomb lattice. The hexagonal arrangement of carbon atoms in graphene can be decomposed into two interpenetrating sublattices of carbon atoms with inversion symmetry between them. This unique topology provides an unusual energy dispersion relation [5]. The energy dispersion near the charge neutrality points, termed “Dirac points”, is of particular interest, since there the 2D energy spectrum is linear, and thus the electrons always move at a constant speed, the band velocity v F ≈ 106 m/s. Many of the interesting electronic properties in graphene result from this dispersion relation which is analogous to that of relativistic, massless fermions. In fact, even before the experimental discovery of graphene, the distinct nature of charge carriers in a hexagonal carbon lattice has been speculated upon and was proposed to have important implications on the electronic transport in graphitic materials, including single-walled carbon nanotubes [6,7].

Z. Jiang et al. / Solid State Communications 143 (2007) 14–19

A sudden burst of experimental and theoretical work on graphene was initiated by the recent work that demonstrated experimentally accessible single and multi layered graphene samples by mechanical extraction [8–11,13] and chemical synthesis [12]. Ever since, numerous unique electrical, chemical, and mechanical properties of graphene have been discovered. In particular, a very unusual half-integer quantum Hall effect (QHE) and a non-zero Berry’s phase [14,15] were discovered in graphene, which provide the unambiguous evidence of the existence of Dirac fermions in graphene and distinguish it from other conventional 2D electronic systems with a finite carrier mass. In this paper, we review the QHE which appears in graphene based on the recent published work [14–16]. We first describe the sample fabrication and characterization method based on mechanical extraction. Then we discuss the experimental results on the QHE in high mobility graphene samples in the intermediate field regime (B < 10 T) followed by the QHE in the extreme quantum limit (B > 20 T). 2. Sample preparation The experimental discovery of an unusual half-integer QHE in graphene was possible because of the significant progress in device fabrication techniques for single layer graphene samples [8,13], and in particular the success in obtaining high mobility in such devices. The high mobility graphene samples used in our experiments were extracted from Kish graphite (Toshiba Ceramics Co.) or other natural graphite single crystals. Single layers of graphene were deposited on silicon wafers, following a method similar to the one described in Ref. [13]. This method involves a very simple, yet efficient method for mechanical extraction of mesoscopic graphite pieces. First, single crystals of graphite are fixed on an adhesive tape and cleaved to expose a fresh surface. Next, very thin flakes of graphite pieces are transferred from the crystals to the surface of silicon wafer by bringing the cleaved surface of the single crystals in contact with the wafer. Various thicknesses are deposited on the surface of SiO2 /Si substrate including single layers of graphite with a lateral size 1–100 µm. Optical microscope images show color contrast in samples of different thickness due to interference induced color shifts. These color shifts are sensitive enough to determine the number of graphene layers on top of the SiO2 /Si substrate. We also note that the thickness of SiO2 layer affects the contrast between different layers. 300 nm thick thermally grown SiO2 layer yields the best contrast for this purpose. After identifying suitable samples, we choose graphene pieces of lateral size 3–10 µm for device fabrication. Multiple electrodes, arranged in Hall bar (Fig. 1(a)) or van der Pauw geometry (Fig. 1(b)), are fabricated on the sample using electron beam lithography followed by Cr/Au (3/30 nm) evaporation and a lift-off process. For samples with non-ideal geometry, or for smaller single-layer pieces connected to larger pieces of thicker graphite, a second electron beam lithography step was performed to define an etch mask on the top of

15

Fig. 1. Optical microscope images of naturally Hall bar shaped (a), van der Pauw geometry (b), and an etched Hall bar shaped (c) graphene devices. (d) Schematic diagram for a graphene device and measurement schemes.

graphene. We then use an oxygen plasma etching process to remove unnecessary parts of the samples (Fig. 1(c)). For the longitudinal and Hall resistance measurement, we apply a small AC current (I . 10 nA) through the samples and measure the transverse voltage (Vx y ) and longitudinal voltage (Vx x ) using a lock-in amplifier. For the Hall bar geometry, the magnetoresistance Rx x and Hall resistance Rx y are then obtained from Rx x = Vx x /I and Rx y = Vx y /I , respectively (Fig. 1(d)). The degenerately doped silicon substrate serves as a gate electrode with silicon oxide acting as the gate dielectric. By applying a gate bias voltage, Vg , the charge density of the sample can be tuned. The induced charge density can be derived from a Hall measurement, where a magnetic field B is applied perpendicular to the sample surface, and a linear increase Rx y (B) is observed. In this measurement, we obtain the carrier density n e from the slope of these lines: n e = 1/e(dRx y /dB). The mobility of the sample is given by µ = (L/W )(en e Rx x )−1 , where W and L are the width and length of the Hall bar. The gate voltage corresponding to the charge neutrality point (i.e., n e = 0) is often found shifted to VDirac , indicating a natural charge doping from the environment to the samples. Although µ defined in this way diverges near n e ≈ 0 due to the failure of the semiclassical Drude model near the Dirac points of graphene, the limiting value of the mobility, µ L , defined in the large density limit (n e ≈ 5 × 1012 cm−2 ) serves as a good indicator of the sample quality. The majority of samples we measured so far exhibit µ L values in the 2000–20,000 cm2 /V s range. Empirically, we found that the samples with µ L > 5000 cm2 /V s show the QHE in the intermediate field regime (4 < B < 10 T). Only the highest quality samples, with quality (µ L > 15,000 cm2 /V s), exhibit structures beyond half-integer QHE in the extreme quantum limit (B > 20 T). 3. Half-integer quantum Hall effect in graphene Exceptionally high mobility graphene samples allow us to investigate transport phenomena in the magnetic quantum limit

16

Z. Jiang et al. / Solid State Communications 143 (2007) 14–19

Fig. 2. Quantized magnetoresistance and Hall resistance of a graphene device where n ≈ 1012 cm−2 and T = 1.6 K. The horizontal lines correspond to the inverse of the multiples e2 / h. The QHE in the electron gas is demonstrated by at least two quantized plateaus in R x y with vanishing R x x in the corresponding magnetic field regime.

where the QHE manifests itself. Fig. 2 shows Rx y and Rx x of a typical high mobility (µ > 10,000 cm2 /V s) graphene sample as a function of magnetic field B at a fixed gate voltage Vg > VDirac . The overall positive Rx y indicates that the contribution is mainly from electrons. At high magnetic field, Rx y (B) exhibits plateaus and Rx x is vanishing, which are the hallmark of the QHE. At least two well-defined plateaus with values (2e2 / h)−1 and (6e2 / h)−1 , followed by a developing (10e2 / h)−1 plateau, are observed before the QHE features transform into Shubnikov–de Haas (SdH) oscillations at lower magnetic field. We observed the equivalent QHE features for holes Vg < VDirac with negative Rx y (B) values. Alternatively we can access the QH plateaus by tuning the electron density by adjusting Vg at a fixed magnetic field. Fig. 3 shows Rx y of the sample of Fig. 2 as a function of gate voltage Vg at B = 9 T. A series of fully developed QH states, i.e., plateaus in h/(e2 ν) quantized to values with an integer filling factor ν, are observed, which are the hallmark of the QHE. Well-defined ν = ±2, ±6, ±10, ±14 QH states are clearly seen, with quantization according to   1 e2 Rx−1 = ±4 n + (1) y 2 h where n is a non-negative integer, and +/− stands for electrons and holes respectively. This quantization condition can be translated into the quantized filling factor ν = ±4(n + 1/2) in the usual QHE language. While the QHE has been observed in many 2D systems, the QHE observed in graphene is distinctively different from those ‘conventional’ QHE’s since the quantization condition Eq. (1) is shifted by a half integer. This so-called half-integer QHE is unique to graphene. It has been predicted by several theories which combine ‘relativistic’ Landau levels (LLs) with the particle–hole

Fig. 3. The Hall resistance as a function of gate voltage at fixed magnetic field B = 9 T, measured at 1.6 K. The horizontal lines correspond to the inverse of integer multiples of e2 / h values.

symmetry of graphene [1–3]. The experimental phenomena can be understood from the calculated LL spectrum in the Dirac spectrum. As in other 2D systems, application of a magnetic field B normal to the graphene plane quantizes the in-plane motion of charge carriers into LLs. The LL formation for electrons/holes in graphene has been studied theoretically using an analogy to 2 + 1 dimensional Quantum Electro Dynamics (QED) [4], in which the LL energy is given by q E n = sgn(n) 2eh¯ v 2F |n|B. (2) Here e and h¯ are electron charge and Planck’s constant divided by 2π , and the integer n represents an electron-like (n > 0) or a hole-like (n < 0) LL index. In particular, a single LL with n = 0 also occurs, where electrons and holes are degenerate. Note that in Eq. (2). we do not consider a spin degree of freedom, assuming the separation of E n is much larger than the Zeeman spin splitting. Therefore each LL has a degeneracy gs = 4, accounting for spin degeneracy and sublattice degeneracy. This assumption needs to be changed when the magnetic field becomes large as we will discuss in the next section. The observed QH sequence can be understood employing the symmetry argument for the Hall conductivity σx y = −Rx y /(Rx2y + (W/L)2 Rx2x ), where L and W are the length and width of the sample, respectively. With the given LL spectrum in Eq. (2), the corresponding Hall conductance σx y exhibits QH plateaus when an integer of LLs are fully occupied, and jumps by an amount of gs e2 / h when the Fermi energy, E F , crosses a LL. Time reversal invariance guarantees particle–hole symmetry and thus σx y is an odd function in energy across the Dirac point [4]. Here, in particular, the n = 0 LL is pinned at zero energy. Thus the first plateau for electrons (n = 1) and holes (n = −1) are situated exactly at

Z. Jiang et al. / Solid State Communications 143 (2007) 14–19

17

Fig. 5. Inset: R xmin x at ν = −4 (upper inset) and ν = +4 (lower inset) for 30 T (diamond), 37 T (down triangle) and 45 T (up triangle), applied normal to the graphene plane. The straight lines are linear fits to the data. Main figure: the activation energy, 1E, as a function of Btot for filling factor ν = −4 (filled circle) and ν = +4 (open circle). The straight lines are linear fits to the data (reproduced from Ref. [16]). Fig. 4. The Hall conductivity as a function of gate voltage at fixed magnetic field B = 25 T, measured at 1.6 K. The horizontal lines correspond to multiples of e2 / h values.

± 12 gs e2 / h. As E F crosses the next electron (hole) LL, Rx y increases (decreases) by an amount of gs e2 / h, which yields the quantization condition in Eq. (1). 4. Graphene quantum Hall effect in high magnetic fields As the magnetic field grows, the LL spectrum assumed in Eq. (2) is no longer valid due to increasing Zeeman spin splitting energy and many-body related correlation effect in the extreme quantum limit. Fig. 4 shows experimentally obtained σx y at B = 25 T as a function of Vg . A series of fully developed QH states with plateaus quantized to values gs (n + 1/2)e2 / h are observed. These well-defined QH states in σx y are in accordance with the previous low-magnetic field measurements (B < 20 T) shown in Fig. 3. In addition to these QH states, new QH states at ν = 0 and at ν = ±1, ±4 are clearly resolved in the high magnetic field data B > 20 T. These new QH states can be attributed to the magnetic field induced splitting of the n = 0 and n = ±1 LLs [16]. In lower magnetic fields, each LL at energy E n is assumed to be four-fold degenerate due to a two-fold spin degeneracy and a two-fold sublattice symmetry. The observed filling factor sequence ν = 0 and ±1 in higher fields (B > 20 T) implies that these degeneracies of n = 0 LL are fully lifted, such that σx y increases in steps of e2 / h as E F passes through the LLs. On the other hand, the fourfold degeneracy of the n = ±1 LLs is only partially resolved, leaving a two-fold degeneracy in each of these LLs and σx y increases in steps of 2e2 / h. The broken sublattice symmetry in graphene in high magnetic fields is quite surprising, because it has been known that the presence of a magnetic field alone does not break the inversion symmetry of the graphene lattice. Since the sublattice

symmetry is protected by the inversion symmetry of graphene, single-electron type of perturbations such as simple uniaxial strain does not lift the graphene sublattice degeneracy. This points to many-body electron correlations within LLs as an alternative origin for the lifting of the degeneracy at the Dirac point. There have been numerous theoretical discussion for the possible origin of this observed degeneracy lifting in the n = 0 LL [17–31]. The QH plateau located right at the Dirac point (Vg ≈ VDirac ) is particularly interesting. While σx y exhibits a clearly resolved ν = 0 plateau, the observed Rx x always shows a finite peak rather than zero for ν = 0 state. Thus this new type of QH state does not conform to the standard QH observation. Recent studies focused on this state reveals [32,33] that the ν = 0 QH state results from the splitting of the spin degeneracy, leaving two coexisting edge states circulating in opposite direction with opposite spin [32]. In order to further investigate the origin of the ν = ±4 QH states, we have carried out detailed temperature dependence measurements to determine the energy gaps between neighboring LLs, 1E [16]. In particular, we measured the temperature dependence of the associated Rx x minimum, Rxmin x , at ν = ±4. The inset of Fig. 5 shows the Arrhenius plots which reveal thermally activated behavior, Rxmin ∝ x exp[−1E/2kB T ], where kB is the Boltzmann constant. We find that 1E, derived from linear fits to these data, depends linearly on the magnitude of total magnetic field Btot (the main section of Fig. 5). Our experimental results suggest that the energy gap of the QH states can be written as [34]: 1E = g ∗ µ B Btot − 2Γ , where g ∗ is the effective g factor, µB is the Bohr magneton, and Γ is the half-width of the LL broadening at half-maximum. We assume that Γ is constant, then g ∗ and Γ can be extracted from the slope and y-intercept of the linear fit and find g ∗ = 2.0 and Γ = 18.2 K from Fig. 5 for the QH state at ν = −4 and g ∗ = 1.7 and Γ = 15.8 K for the ν = +4 QH state from a

18

Z. Jiang et al. / Solid State Communications 143 (2007) 14–19

5. Summary Mesoscopic graphene devices are fabricated employing high mobility graphene samples mechanically extracted from graphite single crystals. Magentoresistance and Hall resistance measurements show the unusual half-integer quantum Hall effect in graphene. We find that the quantization condition in graphene is distinctively different from conventional 2D systems by a shift of a half-integer. In the extreme quantum limit, we observe that n = 0 LL splits into four sublevels, lifting both spin and sublattice degeneracy while n = ±1 LLs lift the spin degeneracy only. Detailed temperature dependence measurements show that the effective g factor at filling factors ν = ±4 is close to the bare electron g factor. Due to the large quasiparticle gaps between the n = 0 and n = ±1 LLs, a robust QHE can be observed at room temperature. Acknowledgements Fig. 6. The magnetoresistance and Hall resistance as a function of gate voltage at fixed magnetic field B = 45 T, measured at 300 K. The horizontal lines correspond to the inverse of the multiples of e2 / h values.

similar analysis. The fact that these measured effective g factors are very close to those of the bare electron (g ∗ = 2) suggests a spin-splitting origin of the ν = ±4 QH states. This conclusion has been further substantiated by the measurements in tilted magnetic fields [16], where a strong dependence of Rxmin x was found only for the total magnetic field applied. We also note that these values are in reasonable agreement with the estimation, Γ ∼ h¯ /τ , where τ ≈ 100 fs is the relaxation time deduced from the carrier mobility in the zero magnetic field limit. Finally, we discuss the robustness of QHE in graphene. As most quantum phenomena, the QHE has been usually associated with very low temperatures, typically below 4 K. Previous efforts to extend the QHE range by using semiconductors with small effective masses of charge carriers and ultra-high magnetic fields allow the observation of the QHE at T . 30 K [35]. Owing to the highly unusual nature of charge carriers that mimic massless relativistic particles, it has been recently demonstrated that the QHE can be observed even at elevated temperature T & 300 K. Fig. 6. shows Rx x and Rx y in graphene measured at room temperature at B = 45 T. One can see that Hall resistance Rx y exhibits clear plateaus at h/2e2 for both electrons and holes, while the longitudinal resistance Rx x shows a clear dip toward zero. This quantization is in agreement with the low-temperature QHE described in the previous section. The survival of the QHE to so high temperatures can be explained by a large cyclotron gap associated with massless quasiparticles in graphene within the experimentally attainable magnetic √ field range. Following Eq. (2), we have 1E ≈ 400 K B(T ) for the transition between the lowest LL (n = 0) and the next one (n = ±1). A quasiparticle energy gap as large as 2500 K is expected at B = 45 T, well exceeding the room temperature thermal energy kB T . Such a large gap is a result of the massless relativistic carrier dynamics unique to graphene, which makes the QHE robust enough to survive at elevated temperature.

This work is supported by the DOE (DE-AIO2-04ER46133 and DE-FG02-05ER46215), NSF (DMR-03-52738 and CHE0117752), NYSTAR, and the Keck Foundation. A portion of this work was performed at the National High Magnetic Field Laboratory, which is supported by NSF Cooperative Agreement No. DMR-0084173, by the State of Florida, and by the DOE. References [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12]

[13] [14] [15] [16]

[17] [18] [19] [20]

Y. Zheng, T. Ando, Phys. Rev. B 65 (2002) 245420. V.P. Gusynin, S.G. Sharapov, Phys. Rev. Lett. 95 (2005) 146801. N.M.R. Peres, F. Guinea, A.H.C. Neto, Phys. Rev. B 73 (2006) 125411. F.D.M. Haldane, Phys. Rev. Lett. 61 (1988) 2015. P.R. Wallace, Phys. Rev. 71 (1947) 622. T. Ando, T. Nakaishi, R. Saito, J. Phys. Soc. Japan 67 (1998) 2857. P.L. McEuen, M. Bockrath, D.H. Cobden, Y. Yoon, S.G. Louie, Phys. Rev. Lett. 83 (1999) 5098. K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, Y. Zhang, S.V. Dubonos, I.V. Grigorieva, A.A. Firsov, Science 306 (2004) 666–669. Y. Zhang, J.P. Small, W.V. Pontius, P. Kim, Appl. Phys. Lett. 86 (2005) 073104. Y. Zhang, J.P. Small, E.S. Amori, P. Kim, Phys. Rev. Lett. 94 (2005) 176803. J.S. Bunch, Y. Yaish, M. Brink, K. Bolotin, P.L. McEuen, Nano Lett. 5 (2005) 287–290. C. Berger, Z.M. Song, T.B. Li, X.B. Li, A.Y. Ogbazghi, R. Feng, Z.T. Dai, A.N. Marchenkov, E.H. Conrad, P.N. First, W.A. de Heer, J. Phys. Chem. B 108 (2004) 19912. K.S. Novoselov, D. Jiang, F. Schedin, T.J. Booth, V.V. Khotkevich, S.V. Morozov, A.K. Geim, Proc. Natl. Acad. Sci. USA 102 (2005) 10451. K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, M.I. Katsnelson, I.V. Grigorieva, S.V. Dubonos, A.A. Firsov, Nature 438 (2005) 197. Y. Zhang, Y.-W. Tan, H.L. Stormer, P. Kim, Nature 438 (2005) 201. Y. Zhang, Z. Jiang, J.P. Small, M.S. Purewal, Y.W. Tan, M. Fazlollahi, J.D. Chudow, J.A. Jaszczak, H.L. Stormer, P. Kim, Phys. Rev. Lett. 96 (2006) 136806. D.V. Khveshchenko, Phys. Rev. Lett. 87 (2001) 206401; D.V. Khveshchenko, Phys. Rev. Lett. 87 (2001) 246802. E.V. Gorbar, V.P. Gusynin, V.A. Miransky, I.A. Shovkovy, Phys. Rev. B 66 (2002) 045108. C.L. Kane, E.J. Mele, Phys. Rev. Lett. 95 (2005) 146802. A.H. Castro Neto, F. Guinea, N.M.R. Peres, Phys. Rev. B 73 (2006) 205408.

Z. Jiang et al. / Solid State Communications 143 (2007) 14–19 [21] [22] [23] [24] [25] [26] [27] [28] [29]

D.A. Abanin, P.A. Lee, L.S. Levitov, Phys. Rev. Lett. 96 (2006) 176803. K Nomura, A.H. MacDonald, Phys. Rev. Lett. 96 (2006) 256602. K. Yang, S. Das Sarma, A.H. MacDonald, Phys. Rev. B 74 (2006) 075423. H.A. Fertig, L. Brey, Phys. Rev. Lett. 97 (2006) 116805. J. Alicea, M.P.A. Fisher, Phys. Rev. B 74 (2006) 075422. M.O. Goerbig, R. Moessner, B. Doucot, Phys. Rev. B 74 (2006) 161407. V.M. Apalkov, T. Chakraborty, Phys. Rev. Lett. 97 (2006) 126801. I.F. Herbut, Phys. Rev. Lett. 97 (2006) 146401. V.P. Gusynin, V.A. Miransky, S.G. Sharapov, I.A. Shovkovy, Phys. Rev.

19

B 74 (2006) 195429. [30] J.-N. Fuchs, P. Lederer, Phys. Rev. Lett. 98 (2007) 016803. [31] Vinu Lukose, R. Shankar, G. Baskaran, cond-mat/0603594. [32] D.A. Abanin, K.S. Novoselov, U. Zeitler, P.A. Lee, A.K. Geim, L.S. Levitov, Phys. Rev. Lett. 98 (2007) 196806. [33] Z. Jiang, Y. Zhang, H.L. Stormer, P. Kim, cond-mat/07051102. Preprint. [34] T. Ando, A.B. Fowler, F. Stern, Rev. Modern Phys. 54 (1982) 437. [35] S.Q. Murphy, et al., Physica E 6 (2000) 293; G. Landwehr, et al., Physica E 6 (2000) 713.