Phases of Na x CoO 2

Similar documents
A Hydrated Superconductor

Electrochemical synthesis and properties of CoO 2, the x = 0 phase of the A x CoO 2 systems (A = Li, Na)

Supporting Information

Organic Conductors and Superconductors: signatures of electronic correlations Martin Dressel 1. Physikalisches Institut der Universität Stuttgart

Role of the Octahedra Rotation on the Electronic Structures of 4d Transition Metal Oxides

Computational strongly correlated materials R. Torsten Clay Physics & Astronomy

Superconductivity, charge ordering and structural properties of

YBCO. CuO 2. the CuO 2. planes is controlled. from deviation from. neutron. , blue star for. Hg12011 (this work) for T c = 72

Research Highlights. Salient results from our group. Mixed phosphides in Sn-P and Sn-Mn-P systems

Study on Magnetic Properties of Vermiculite Intercalation compounds

Magnetism in correlated-electron materials

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION

SECOND PUBLIC EXAMINATION. Honour School of Physics Part C: 4 Year Course. Honour School of Physics and Philosophy Part C C3: CONDENSED MATTER PHYSICS

High temperature superconductivity

ARPES study of many-body effects and electronic reconstructions in misfit cobaltates

Coexistence of Ferromagnetic and Glassy States in Mechanically Milled GdAl 2

Magnetic Order versus superconductivity in the Iron-based

SUPPLEMENTARY INFORMATION

Harald Ibach Hans Lüth SOLID-STATE PHYSICS. An Introduction to Theory and Experiment

Arnab Pariari & Prabhat Mandal Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Calcutta , India

Origin of the anomalous low temperature upturn in resistivity in the electron-doped cuprates.

Linear temperature dependence of electron spin resonance linewidths in La 0.7 Ca 0.3 MnO 3 and YBaMn 2 O 6

Supplementary Figure S1: Number of Fermi surfaces. Electronic dispersion around Γ a = 0 and Γ b = π/a. In (a) the number of Fermi surfaces is even,

The Oxford Solid State Basics

Jim Freericks (Georgetown University) Veljko Zlatic (Institute of Physics, Zagreb)

Anisotropic Magnetic Structures in Iron-Based Superconductors

Material Science II. d Electron systems

Observation of topological surface state quantum Hall effect in an intrinsic three-dimensional topological insulator

1 G. Kotliar: Lecture 2

Transport and Magnetic Properties of La 0.85 Ca Mn 1-x Al x O 3 Manganites

Assessment of Variation in Zero Field Hall Constant of Colossal Magnetoresistive Manganites (Re1-x AxMnO3)

Self-compensating incorporation of Mn in Ga 1 x Mn x As

Heavy Fermion systems

Magnetic Properties and Scaling Behavior in Perovskite like La 0.7 (Ba 1-x Pb x ) 0.3 CoO 3 System

Key words: High Temperature Superconductors, ARPES, coherent quasiparticle, incoherent quasiparticle, isotope substitution.

Superconductivity at 41.0 K in the F-doped LaFeAsO 1-x F x

ARPES studies of cuprates. Inna Vishik Physics 250 (Special topics: spectroscopies of quantum materials) UC Davis, Fall 2016

Metal-insulator transitions

Superconductivity and Electron Correlations in Ruthenates

SUPPLEMENTARY INFORMATION

High-T c superconductors. Parent insulators Carrier doping Band structure and Fermi surface Pseudogap and superconducting gap Transport properties

SECOND PUBLIC EXAMINATION. Honour School of Physics Part C: 4 Year Course. Honour School of Physics and Philosophy Part C C3: CONDENSED MATTER PHYSICS

SUPPLEMENTARY INFORMATION

Name: (a) What core levels are responsible for the three photoelectron peaks in Fig. 1?

Supplementary Materials for

Fe Co Si. Fe Co Si. Ref. p. 59] d elements and C, Si, Ge, Sn or Pb Alloys and compounds with Ge

Intermediate valence in Yb Intermetallic compounds

Chemical Differences between K and Na in Alkali Cobaltates

Hole-concentration dependence of band structure in (Bi,Pb) 2 (Sr,La) 2 CuO 6+δ determined by the angle-resolved photoemission spectroscopy

Crystal Structure and Elementary Properties of Na x CoO 2 (x = 0.32, 0.5, 0.6, 0.75, and 0.92) in the Three-Layer NaCoO 2 Family

Probing the Electronic Structure of Complex Systems by State-of-the-Art ARPES Andrea Damascelli

Structural and magnetic characterization of the new GdMn 1-x. O 3 perovskite material

arxiv:cond-mat/ v3 [cond-mat.supr-con] 23 May 2000

Chemical differences between K and Na in alkali cobaltates

Shigeki Yonezawa*, Yuji Muraoka and Zenji Hiroi Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba

Photoelectron Spectroscopy

SOLID STATE PHYSICS. Second Edition. John Wiley & Sons. J. R. Hook H. E. Hall. Department of Physics, University of Manchester

Electronic structure of correlated electron systems. Lecture 2

NEURON SCATTERING STUDIES OF THE MAGNETIC FLUCTUATIONS IN YBa 2 Cu 3 O 7¹d

Novel Magnetic Properties of Carbon Nanotubes. Abstract

New perspectives in superconductors. E. Bascones Instituto de Ciencia de Materiales de Madrid (ICMM-CSIC)

The Hubbard model in cold atoms and in the high-tc cuprates

First Experimental Evidence of Potassium Ordering in Layered K 4 Co 7 O 14

Emergent electronic matter : Fermi surfaces, quasiparticles and magnetism in manganites and pnictides de Jong, S.

Examples of Lifshitz topological transition in interacting fermionic systems

It is known that the chemical potential of lithium ion in the cathode material can be expressed as:

Superconducting Stripes

Effect of dimensionality in polymeric fullerenes and single-wall nanotubes

Supplementary Information for

Review of typical behaviours observed in strongly correlated systems. Charles Simon Laboratoire CRISMAT, CNRS and ENSICAEN, F14050 Caen.

Effect of randomness on anomalous Hall coefficient in antiferromagnet U 2 PdGa 3

A New Electronic Orbital Order Identified in Parent Compound of Fe-Based High-Temperature Superconductors

Iridium Containing Honeycomb Delafossites by Topotactic Cation Exchange. Supplemental Material

Dynamics of fluctuations in high temperature superconductors far from equilibrium. L. Perfetti, Laboratoire des Solides Irradiés, Ecole Polytechnique

Spin correlations in conducting and superconducting materials Collin Broholm Johns Hopkins University

High-Temperature Superconductivity in Lattice-Expanded C 60

High-T c superconductors

Supplementary Figures.

Ferromagnetic and spin-glass properties of single-crystalline U 2 NiSi 3

Photoemission and the electronic structure of magnetic oxides. Dan Dessau University of Colorado, Boulder Duane F625

Supplementary Figure 1 PtLuSb RHEED and sample structure before and after capping layer

2 B B D (E) Paramagnetic Susceptibility. m s probability. A) Bound Electrons in Atoms

Supplementary Figures

What's so unusual about high temperature superconductors? UBC 2005

ORBITAL SELECTIVITY AND HUND S PHYSICS IN IRON-BASED SC. Laura Fanfarillo

Dimensionality controlled Mott transition and correlation effects in

Magnetic Transition in the Kondo Lattice System CeRhSn 2. Z. Hossain 1, L.C. Gupta 2 and C. Geibel 1. Germany.

Discovery of spin-vortex-crystal magnetic order in Ni- and Co- doped CaKFe 4 As 4

Field-Temperature Evolution of Antiferromagnetic Phases in Ludvigites Ni 3-x Mn x BO 5

Chapter 7. Summary and Outlook

CDWs in ARPES. A momentum space picture of Fermi surface instabilities in crystalline solids. Physics 250, UC Davis Inna Vishik

LOCAL MOMENTS NEAR THE METAL-INSULATOR TRANSITION

shows the difference between observed (black) and calculated patterns (red). Vertical ticks indicate

Electromagnetism II. Instructor: Andrei Sirenko Spring 2013 Thursdays 1 pm 4 pm. Spring 2013, NJIT 1

A theoretical study of the single-molecule transistor

Tb 2 Hf 2 O 7 R 2 B 2 7 R B R 3+ T N

New Li-Ethylenediamine-Intercalated Superconductor Li x (C 2 H 8 N 2 ) y Fe 2-z Se 2 with T c = 45 K

Can superconductivity emerge out of a non Fermi liquid.

A BCS Bose-Einstein crossover theory and its application to the cuprates

Transcription:

Phases of Na x CoO 2 by Aakash Pushp (pushp@uiuc.edu) Abstract This paper deals with the various phases of Na x CoO 2 ranging from charge ordered insulator to Curie-Weiss metal to superconductor as the stoichiometry is changed and attempts to understand the related theory and experimental results. It starts with investigating the superconducting properties of this material, intercalated with H 2 O and goes on to describe the various phases that come into being as x is varied. It then concentrates on another remarkable phase at x = ½ which is an insulator. For all other doping, this material shows metallic behavior. Later, it tries to give a clear account of the various experiments performed including ARPES and discusses its anisotropic hole-type large hexagonal-like Fermi surface and its implications.

Introduction Since the discovery of high-t c superconductivity in layered copper oxides, many other metal oxides involving 3d-transition metals, such as cobalt and nickel, have been investigated for similar behavior. Such attempts have normally failed but in March 2003 Takada et al. [1] reported that Na x CoO 2.yH 2 O (x ~ 0.35, y ~ 1.3) is a superconductor with a T c of about 5 K. This compound consists of 2-dimensional CoO 2 separated by a thick insulating layer of Na + ions and H 2 O molecules. It shows a marked resemblance in the superconducting properties with high-t c copper oxides, suggesting that both the systems have similar underlying physics. X-Ray Diffraction data Takada et al. [1] obtained Na x CoO 2.yH 2 O sample through chemical oxidation of the parent compound Na 0.7 CoO 2 [4]. Before oxidation, Na 0.7 CoO 2 has a hexagonal layered structure (space group P6 3 /mmc) consisting of the 2D layers of CoO 2 and charge balancing Na + ions (Figure 1). Figure 2 shows an XRD pattern of the product obtained from oxidation. All the reflections were index matched on the basis of space Figure 1: Structural views of Na 0.7 CoO 2 (left) and Na x CoO 2 yh2o (right), where Na and H 2 O sites are partially occupied. [ref. 1] 2

group P6 3 /mmc. A marked increase in c was observed supporting the idea that the intercalation of certain guest molecules occurred in the oxidation process in addition to the de-intercalation of Na + ions. H 2 O is the most probable candidate for the guest molecule. The composition calculated by inductively coupled plasma atomic-emission spectroscopy (ICP-AES) based on this assumption was Na 0.35 CoO 2.1.3H 2 O. Further, infrared spectroscopy showed that the H 2 O content was 1.24 per formula unit. Figure 2: Rietveld refinement patterns for Na x CoO 2 yh 2 O. The observed diffraction intensities are represented by red plus signs, and the calculated pattern by the green solid line. Short light green vertical bars below the observed and calculated patterns indicate the positions of allowed Bragg reflections. The oxidation process slightly decreased a from 2.8292(3) Å for the parent compound to 2.8230(2) Å, where the parentheses indicate standard deviations, probably because Co 3+ ions were partially oxidized to smaller Co 4+ ions. On the other hand, c increased dramatically from 10.9628(8) Å to 19.6207(14) Å. [ref. 1] Rietveld analysis of the XRD data was carried out adopting a structural model shown in the right side of Figure 1. In this model, H 2 O molecules are located at 12k sites on the basis of electron-density distribution obtained by a maximum-entropy method. The Rietveld refinement patterns are shown in Figure 2. The H 2 O content resulting from these analyses came out to be 1.47 per formula unit, which is close to analytical values. 3

Magnetic Susceptibility and resistivity The magnetic susceptibility of Na x CoO 2.yH 2 O is plotted in Fig. 3 as a function of temperature. A steep decrease of susceptibility is observed in the measurements under the magnetic field H = 20 Oe at about 5 K both in zero-field cool (ZFC) and field cool (FC) processes. Figure 3: Magnetic susceptibility (χ) of Na x CoO 2 yh 2 O. Filled circles, zero-field cooling; open circles, field cooling. The susceptibility was measured in a magnetic field of 20 Oe. The inset shows χ measured under various magnetic fields by the zero-field cooling process. [ref. 1] From the susceptibility data, it can be concluded that the present phase undergoes superconducting transition at about 5 K, because only superconductivity can account for such a large diamagnetism. Downturn of the susceptibility was still observed near 5 K even at 70 koe, but the transition became quite broad and the susceptibility was not negative even at 2 K. 4

The electric resistivity (ρ) of Na x CoO 2 yh 2 O is shown in Fig. 4. Zero resistivity was not obtained down to 2 K in sample made by Takada et al. [1] because it was impossible to prepare a tightly sintered ceramic specimen. Figure 4: Resistivity (ρ) of Na x CoO 2 yh 2 O under zero magnetic field. The inset figure shows the resistivity measured under various magnetic fields. [ref. 1] Nevertheless, a sharp decrease of resistivity can be seen at around 4 K, supporting the superconducting transition. The onset temperature is lower in resistivity measurement than that of the susceptibility measurement because of variation of water content in the sample. The behavior of this material not only depends on the water content but also on the Na content. The phase change of this material ranges from that of a paramagnetic metal (x ~ 0.3) to a charge-ordered insulator (at x = 0.5) to a Curie-Weiss metal (x ~ 0.7), and finally to a weak-moment magnetically ordered state (x > 0.75). The next section discusses the various phases of Na x CoO 2. 5

Phases of Na x CoO 2 Research on oxide conductors has uncovered many interesting electronic states characterized by strong interaction, which include superconductivity, and charge or spin ordered states [5, 6]. Na x CoO 2 is one such oxide. At the doping x ~ 2/3, Na x CoO 2 exhibits an unusually large thermopower [7]. Although Na x CoO 2 shows metallic resistivity, its magnetic susceptibility shows a surprising Curie-Weiss profile [8], and its magnitude is comparable with that of antiferromagnetically coupled spin-1/2 local moments equal in number to the hole carriers [9]. Another interesting observation is that the thermopower at 2.5 K gets suppressed by an in-plane magnetic field [9], which implies that the enhanced thermopower is largely due to spin entropy carried by strongly correlated holes (Co 4+ sites) hopping on the triangular lattice. We have already seen in the previous section that when intercalated with water, Na x CoO 2.yH 2 O becomes superconducting for ¼ < x < 1/3 [10, 11, 12]. In order to find out if the Curie-Weiss state at 2/3 is continuous with the 1/3 state surrounding superconductivity, and if chargeordering effects are important, Cava et al. [2] have completed a phase diagram of nonhydrated Na x CoO 2. The phase diagram depicts a variety of electronic states as x is increased from 0.3 to 0.75. One of the most interesting among them is at x = ½ that involves charge ordering of the holes together with the Na + ions. Experimental details Cava et al. [2] accurately measured the variation of the c-axis lattice parameter variation as x is increased (Figure 5C). Figure 5A shows magnetic susceptibility χ vs. T as x is varied. For x ~ 2/3, χ vs. T follows the Curie-Weiss law χ = C/(T + θ) with θ ~ 70 K. Above 2/3, at x = 0.75, χ is slightly rounded below 20 K, consistent with the appearance of a weak magnetization M. For x > ¾, a spin-density-wave (SDW) has been suggested [13, 14]. When x is decreased below 2/3, the Curie-Weiss divergence in χ gets progressively reduced until it vanishes at x = ½. At x = ½, χ reveals two sharp cusps at T c1 = 88 K and T c2 = 53 K, which signal the onset of an insulating state. With further decrease of x below ½ to x = 0.31, the χ plot becomes featureless, but its magnitude remains large as compared with the Pauli susceptibility in conventional metals. The in-plane resistivity ρ shows a dramatic change in behavior across x = ½ (Figure 5B). Na x CoO 2 is metallic in its T dependence at all values away from x = ½. Insulating behavior abruptly appears at x = ½, in contrast to the metallic states at all other values of x. Initially, as the temperature is lowered towards T c1 and T c2, the resistivity increases gradually. Below T c2, however, ρ shoots up rapidly. At lower doping (x = 0.31), it behaves as a high-conductivity metallic state. Below 30 K, ρ follows a T 2 behavior and falls down to a value almost 5 times lower than at x = 0.71. 6

Figure 5: The susceptibility χ (A) and in-plane resistivity ρ (B) of single crystals of Na x CoO 2 with x determined by ICP (C). In Panel A, χ is measured in an in-plane field H = 5 T. In the crystal with x = 0.75, χ fits the Curie-Weiss form χ = C/(T + θ) (with θ ~ 150 K and C = 2.81 cm 3 K/mole). For x = 0.5, sharp transitions are observed at T c1 and T c2 (arrows). Panel B shows the T dependence of ρ at selected x. Insulating behavior is observed at x = 0.5 (data displayed at lower scale) in contrast to metallic behavior in the rest. At low T, ρ is T-linear for x = 0.71 but varies as T 2 for x = 0.3. In Panel C, the c-axis lattice parameter measured by XRD is plotted against the Na content x fixed by ICP in powder samples. [ref. 2] The susceptibility and resistivity measurements imply a phase diagram as shown in Figure 6. The interesting feature in the phase diagram is the narrow insulating state at x = ½ separating two distinct metallic states. It behaves as a paramagnetic material below ½ while shows a Curie-Weiss metallic state behavior above x = ½. 7

Figure 6: The phase diagram of non-hydrated Na x CoO 2 (A), electron diffraction pattern of a crystal with x = ½ (B) and the Na superstructure lattice (C). In Panel A, the charge-ordered insulating state at ½ is sandwiched between the paramagnetic metal at 0.3 and the Curie-Weiss metallic state at 0.65-0.75. The superconducting state is obtained on intercalation with H 2 O [1, 10]. In Panel B, the diffraction pattern is taken along [001]. Bragg spots of the hexagonal lattice are encircled. The spots located at n 1 b 1 + n 2 b 2 correspond to the supercell with lattice vectors a 1 and a 2 drawn in Panel C. [ref. 2] Electron-diffraction (Figure 6B) reveals a very important feature of this interesting state at x = ½. Even at 300 K, the Na ions order as a superstructure with lattice vectors 3a x and 2a y (where a is the hexagonal parameter) as seen in the Bragg spots. At other doping levels, the superstructure Bragg spots are either much weaker or absent. The thermal conductivity κ measured parallel to the layers also clearly depicts the long-range nature of the Na superstructure. As phonons are strongly scattered by disorder in the Na sublattice, the phonon mean free path l ph is expected to be much longer for x = ½ due to the superstructure than at any other doping. This is confirmed by the data [2] of thermal conductivity κ vs. T for different doping (Figure 7). 8

Figure 7: The thermopower S and the Hall coefficient R H in Na x CoO 2 for x = 1/2 (A), and plots of the inplane thermal conductivity κ (B), and R H (C) for several values of x. In Panel A, both S and R H begin to increase steeply just below T c2 (88 K) and attain maxima at 40 and 30 K, respectively. Below the peak values, they fall towards zero as T 0, consistent with a particle-hole symmetric state. In Panel B, κ at x = ½ rises to a very large peak value, whereas it remains small at all other x. Panel C shows that the curve of R H vs. T at x = ½ is qualitatively different from all other x values. In the metallic samples, R H increases with T, consistent with strongly interacting carriers hopping on a triangular lattice. [ref. 2] In the metallic samples (x = 0.31 and 0.71), κ weakly depends on T and is also small in magnitude as compared to the κ for the x = ½ sample. κ also shows a very steep rise near the T c2 confirming the superstructure lattice idea. Hence both the diffraction and κ results provide strong evidence that the Na ions at x = ½, unlike other metallic state, order with very long range correlation even at room temperature. 9

Even the Hall coefficient R H at x = ½ shows a small peak before plunging down to large negative values. Above 200 K, R H is hole-like with a T-independent value. Near T c1, R H rises to a weak maximum and then plunges down to a large negative value. Below T c2, it increases dramatically in magnitude to a peak value 250 times greater than its magnitude above 200 K. This giant increase in R H indicates that the itinerant carrier density n decreases by ~2 orders of magnitude. Once again, the corresponding plots for the metallic states are featureless. However, R H in these metallic states is T-dependent above 200 K. The behavior of R H at low T has an unusual feature for the insulating state. R H falls steeply towards zero as T 0 as shown in Figure 7A. The usual behavior of semiconductors, for instance, is that it diverges as 1/n (as n 0). We know that in general, R H = σ H /Hσ 2, where σ H is the Hall conductivity and σ is the conductivity. The vanishing of Hall current at low temperature indicates that the charge excitations exhibit exact particle-hole symmetry, i.e., they are neither hole-like nor electron-like in their response to the applied magnetic field. From the Angle-Resolved Photoemission (ARPES) data [3], a large hole-type hexagonal-like Fermi surface (FS) is observed. Figure 8 reveals a large hole-pocket centered around the Γ-point. The FS assumes a hexagonal character with an average radius of 0.7 ± 0.05 Å -1. This is a bit larger than the size calculated for the related compound Na 0.5 CoO 2 using LDA, which may indicate that a fraction of the doped electrons are in the Fermi sea while the rest are presumably localized. The shape of the FS is rather hexagonal (anisotropic) in agreement with the LDA. Figure 8: Fermi surface: n(k) plot generated by integrating within 75 mev of Fermi level. A large holepocket is centered around the Γ-point. The Fermi surface, exhibiting some hexagonal anisotropy, is the inner edge of pocket as shown over the complete Brillouin zone. [ref. 3] 10

Along the Γ M energy crossing the quasiparticle band crosses the Fermi level from M toward Γ (Fig. 9(a)) as opposed to Γ to M, which gives the single-particle hopping (t) a negative sign. The total bandwidth of this system is estimated at 70(10) mev. For a tight binding hexagonal lattice total bandwidth W = 9t where t is the nearest neighbor single-particle hopping. Hence the effective nearest neighbor single particle hopping, t eff in this system can be estimated to be about 8(2) mev. It is interesting to note that this value is on the order of the exchange coupling J (~ 10 mev) [8], which suggests that the charge motion would be significantly affected by the spin dynamics of the system. Figure 9: Quasiparticle dispersion: (a) Γ M Fermi crossing. Color red reflects the highest intensity - yellow to green to blue is in the order of decreasing intensity. (b) EDCs corresponding to the image plot in (a). (c) A single EDC for k = 0.78 Å 1. To extract the peak position a background is subtracted. A constant step-like background is best seen in Fig 9(a). Based on the extracted peak positions E vs. k plots are made as shown in (d) for Γ M and in (e) for Γ K directions. [ref. 3] The quasiparticle spectral weight near the Fermi level as a function of temperature has been plotted in Figure 10. It is interesting to note that the quasiparticle weight decreases to almost zero on raising the temperature to around 120 K (Fig. 10(b)). This roughly coincides with the temperature region where the T-linear behavior in resisitivity gives way to stronger T-dependence. In other words, quasiparticles exist in this system only in the temperature regime where resisitivity is linear in T and transport behavior is non-fermiliquid like, in stark contrast with the conventional expectation that well defined quasiparticles are signatures of good Fermi liquid behavior. Further, a strong temperature dependence of the 0.7 ev feature in Na 0.7 CoO 2 is yet another signature of this system s strongly correlated behavior in addition to its highly flat band character which is typical of strong electron-electron correlation. 11

Figure 10: Temperature Dependence: (a) In-plane resistivity in Na 0.7 CoO 2 is linear up to 100 K [9] and then gradually crosses over to a stronger T-dependence. (b) Temperature dependence of quasiparticles near the Γ M Fermi crossing. The quasiparticle spectral weight ceases to exist above 120 K, close to the temperature where the T-linear behavior of the resistivity disappears. (c) Background subtracted integrated quasiparticle spectral weight is plotted as a function of temperature. [ref. 3] Conclusion The phase diagram in Figure 6 reveals the strong influence of Na ion on the electronic properties of Na x CoO 2. The insulating phase indicates a strong interaction between the ions and holes even though they occupy separate layers. For either charge species alone on the triangular lattice at x = ½, the liquid state seems to be more stable at finite temperature. If both populations are present, density fluctuations in one will strongly affect the other thereby giving interesting results in measurements of thermal conductivity, Hall coefficient, resistivity, magnetic susceptibility, etc. The superconducting state produced by intercalation with H 2 O evolves from the paramagnetic metal which is separated from the Curie-Weiss state by the charge-ordered insulator. In other words, this system is a very good example of strongly interacting particles and is an interesting state of matter that offers rich Physics. 12

Reference [1] K. Takada et al., Nature 422, 53 (2003) [2] M.L. Foo et al., cond-mat/0312174 [3] M.Z. Hasan et al., cond-mat/0308438 [4] C. Fouassier et al., J. Solid State Chem. 6, 532 (1973) [5] Y. Tokura, Physics Today 56, 50 (2003). [6] A. P. Mackenzie, and Y. Maeno, Rev. Mod Phys. 75, 657 (2003). [7] I. Terasaki et al., Phys. Rev. B 56, 12685 (1997). [8] R. Ray et al., Phys. Rev. B 59, 9454 (1999). [9] Y. Wang et al., Nature 423, 425 (2003). [10] R. E. Schaak et al., Nature 424, 527 (2003). [11] F. C. Chou et al., cond-mat/0306659 (2003). [12] R. Jin et al., Phys. Rev. Lett. 91, 217001 (2003). [13] T. Motohashi et al., Phys. Rev. B 67, 064406 (2003). [14] J. Sugiyama et al., cond-mat/0310516 (2003). 13