We describe in this chapter a promising system and method to search for novel superconductors by investigating appropriately chosen antiferromagnets and creating carriers by doping using a high throughput spatial composition spread (combinatorial) approach. The method has been applied to the cuprate superconductors and has the potential to enhance our understanding of these materials and push the boundaries of the field by quickly exploring novel ones.
Finding novel superconducting materials, which superconduct at much higher temperatures, now seems to be a realistic goal because of three recent developments: the 2008 discovery of superconductivity in iron arsenide based materials; the observation that a number of superconductors are doped antiferromagnets; and the tremendous progress researchers have made over the past 20 years in understanding the physical properties of existing superconductors. These developments suggest a path to novel superconductors - explore the electrical transport properties of doped antiferromagnets. While a number of applications of existing superconductors have been realized, their widespread use depends on raising the transition temperatures substantially above the current world record Tc of 138 K.
Expanding the number of known systems which exhibit superconductivity also allows researchers to identify its essential elements. These observations help reduce the number of models which purport to explain the mechanism of pair formation, and allow researchers to ignore irrelevant peculiarities of a particular system. In general, the preparation of novel phases of matter increases the likelihood of the discovery of novel material properties. The combinatorial approach to materials discovery (Xiang et al, 1995) allows one to realize these goals at an unprecedented rate.
Discoveries by groups in Japan and China over the past few years have added to the class of systems for which antiferromagnetic order exists in close proximity to superconductivity, and in some cases may even coexist. Iron arsenide-based, cuprate, fulleride and heavy fermion superconductors [detailed references below] populate this class, and the diversity of hosts highlight how useful it is to search for new examples of superconductors in aid of an empirical identification of the important parameters on which to build a correct theory. As the number of examples has grown, the importance of spin fluctuations has emerged. This chapter describes a route to discover novel superconductors in doped antiferromagnets to enhance our understanding of superconductivity. To accomplish these goals we use a combinatorial approach to materials discovery, which we have recently demonstrated allows us to map the superconducting properties of the La214 cuprate superconductor and search for new superconductors by quickly and efficiently exploring phase space in a chosen system. The single layer K2CuF4 and double layer K3Cu2F7 perovskites will be shown to be likely candidates to exhibit superconductivity on the border of antiferromagnetic or ferromagnetic phases. Rapid characterization using high throughput resistivity apparatus, such as the one (Hewitt et al, 2005) developed in our lab, allows one to identify superconducting phases at an accelerated pace.
The spatial composition spread approach is described in section 2, its application to cuprate superconductors in section 3 and section 4 proposes likely superconductors on the boundary of antiferromagnetic and ferromagnetic phases.
2. Spatial composition spread approach
Combinatorial materials science (CMS) methods represent a powerful technique to produce a large number of compositions on a spatially addressable substrate. In particular, one may produce a continuous variation in composition across a substrate using physical vapor deposition techniques. For example, one can locate three evaporation targets at the vertices of a triangle to prepare a ternary phase diagram (Kennedy et al, 1965); or four targets at the corners of a square to prepare a quaternary phase diagram (van Dover et al, 1998) by sputter deposition. The drawback of this approach is that a non-linear variation in composition is produced as a function of position on the substrate. In order to obtain a linear relationship which would allow for a much easier interpretation of the data, we altered the flux produced by sputtering targets through design of masks which intercept the flux to create a linear or constant deposition as a function of position. Although CMS methods have been used to show that particular superconducting phases can be prepared (Xiang et al, 1995), to the author’s knowledge the composition spread approach has not been used for this purpose until recently by our group (Sanderson and Hewitt, 2005, 2007) and subsequently by a Brookhaven group (Logvenov et al, 2007) using molecular beam epitaxy.
To implement a linear composition spread approach, the target deposition profile needs to be altered to provide the appropriate linear (or constant) variation across the substrate. One way to produce this variation is to interrupt the flux with a physical mask placed over the target. The process of determining the target mask profiles, for the Corona Vacuum Coaters V-3T sputtering system used in our lab is described in detail here, as outlined elsewhere (Dahn et al, 2002).
The sputtering flux from a circular magnetron can be determined by sputtering a target in front of a stationary substrate and measuring the mass as a function of position. The result for an aluminum target sputtered in an argon atmosphere at 5.6 mTorr is shown in Figure 2. A Gaussian profile of the form fits the data well, with the full width at half-max, = 10 cm.
To intimately mix the elements the table must be rotated at a high rate, so one must calculate the deposition a rotating substrate. To obtain a specific deposition profile (linear out, linear in, constant) one must be able to calculate the amount of material deposited upon a rotating substrate, through a target mask. Figure 3 illustrates the geometry of the problem in which vectors rT (table centre to target centre location on the table), r1 (from table centre to start of the deposited film) and r2 (distance from the table centre to the film end) are drawn. If r1 and r2 are equidistant from rT, then the deposition would be the same at r1 and r2 when the substrate table is stationary. However, when the table is rotating, a point at r1 has a smaller tangential velocity than a point at r2, so more flux needs to reach r2 to produce the same deposition. For the Corona Vacuum Coaters V-3T sputtering system rT = 13.33 cm, r1 = 9.5 cm and r2 = 17.1 cm.
To calculate the mask profile needed consider the geometry shown in Fig. 3B. R is a vector from the centre of the substrate table to a point below the centre of the target on the table, while q is a vector from the centre of the target to the point on the substrate where the flux is being calculated. By integrating the flux along the heavy arc shown in Figure 3B the deposition D(s) can be found.
It can be shown easily that,
and the deposition as a function of radial position (s) is,
The expression of D(s) can be used to determine the desired deposition profile, whether linearly increasing or decreasing with s, or simply a constant independent of s. Finally, the mask shape can be numerically determined by solving Equations 1 and 2 for max and s. The mask shape is found by converting the max and s values to x and y coordinates. Masks have been designed for a constant deposition, as well as depositions that vary linearly, increasing inwards and outwards. With these three masks a wide variation of film compositions can be made. Figure 4 shows an image of each mask as well as masses of the film deposited onto a rotating substrate through the corresponding mask.
The results presented in the next section usually employ the set-up depicted in Figure 5. Using targets of composition A (placed behind the linear-in mask, blue in Fig. 5) and B (placed behind the linear out mask, red in Fig. 5), results in compositions varying linearly as A1-xBx with radial position (s).
3. Densely mapping the phase diagram of the cuprate superconductors
We first applied the linear composition spread approach to the cuprate superconductor Bi2212 (Sanderson and Hewitt, 2005, 2007), and then La214 (Saadat and Hewitt, 2010). The results for La214 are described here.
The cuprates are doped antiferromagnetic Mott-Hubbard insulators, becoming Fermi liquid metals at large hole concentrations. (They can also be electron doped, although the Tc’s are not as high, reflecting some type of electron-hole asymmetry.) The limiting regimes are well described by the two-dimensional Hubbard model (Hubbard, 1963). When the onsite coulomb repulsion (U) dominates the kinetic energy of hopping (t) (U>>t), the solutions to the Hubbard Hamiltonian produce highly localized electron wavefunctions which gives rise to an insulator, and super-exchange interaction produces the Néel state. When the hopping term (t) dominates (t>>U), the solution is a modified Bloch wave function which gives rise to the Fermi liquid state. At intermediate hole concentrations where t and U are comparable, the transition from localized to itinerant electrons produces interesting physics and remains enigmatic. It is within this region of the phase diagram where superconductivity occurs. Thus it is important to understand the normal state properties of metal-insulator transitions (Mott, 1968), in order to decipher the nature of the superconducting state. Also, in this region of phase space there exists a partial suppression of low energy excitations, a Pseudogap (PG) (Timusk and Statt, 1999), which appears below a temperature T* > Tc. Its exact description is a subject of intense study because it may hold the key to understanding the transition to the superconducting state. Theories of the PG can be divided into two broad categories: ones which identify it as a precursor to the onset of superconductivity, or others which classify it as a competing phase. It is thought that a dense map of the doping dependence of T* is the key to deciphering the nature of the PG (Norman et al, 2005). Whether T* merges with Tc in the overdoped regime or ends at optimal doping determines whether the PG is a “friend” or “foe” of superconductivity.
We have recently obtained data showing that one can densely map the temperature-hole concentration phase diagram of La2-xSrxCuO4+δ (0 ≤ x ≤ 0.18) using the spatial composition spread approach, as described in section 2 of this chapter, and presented in our recent article (Saadat et al, 2010). First applied to the cuprate superconductor Bi2212 (Sanderson and Hewitt, 2005, 2007), it was used to synthesize a La2-xSrxCuO4+δ (0 ≤ x ≤ 0.18) library. In this approach, targets of La2CuO4 and La1.82Sr0.18CuO4 were co-sputtered with specially designed target masks, which ultimately produce a linear composition gradient varying from x = 0 (at 0 mm) to x = 0.18 (at 75 mm) on a set of eight single crystal substrates. The libraries’ structures are characterized by X-ray diffraction (Fig. 6) and cation composition by Energy and Wavelength Dispersive Spectroscopy (EDS/WDS) (Fig. 7).
While we have shown it is possible to prepare single phase films in this manner, they are polycrystalline. Epitaxial films are required to separate contributions to T* from c-axis and ab-plane transport. Polycrystalline films are sufficient, however, to measure Tc. Therefore, a high-throughput resistivity apparatus (Hewitt et al, 2005) was used to measure the DC resistivity of the 52 member library.
Tc and T* were determined and plotted versus Sr content as shown in Fig. 9. We found that Tc is suppressed near 1/8 (x = 0.125) doping, consistent with the formation of a stripe phase (Tranquada et al, 1995). The lowest Sr content (x) at which superconductivity appears is 0.03, not at the expected value of 0.05 (Ando et al, 2004). Independent measurements of the hole
concentration can be obtained by measuring the intensity of a feature located 2 eV below the O Kedge (530 eV) by X-ray absorption spectroscopy (XAS) (Kuiper et al, 1988). These measurements would not only allow us to evaluate whether this unexpected result is a truly novel feature and not simply a consequence of oxygen non-stoichiometry and/or film strain, it would also test a theory that predicts the existence of charge 2e bosons in the PG state (Leigh et al, 2007; Choy et al, 2008).
4. Prospective 2D antiferromagnetic and ferromagnetic superconductors
In this chapter we also propose promising hosts for superconductivity and describe the use of the combinatorial approach to rapidly and efficiently scan the horribly large phase space of possible dopings. In the field of superconductivity there are very few “right” substitutions which produce superconductivity, and then only over a rather limited range of dopings.
Until recently, the cuprates and heavy fermion systems have presented a rather unique example of superconductivity in doped antiferromagnets. The phase diagram of the recently discovered FeAs-based (e.g. RE(O,F)FeAs, RE = La, Sm, or Ce (Luetkens et al, 2009; Drew et al, 2009; Zhao et al, 2008) and (Ba,K)Fe2As2 (Chen et al, 2009) superconductors share many features with the cuprates: a) they are doped antiferromagnets, b) superconductivity occurs
in 2D planes, though corrugated in FeAs-based materials, c) the planes are doped by adjacent charge reservoir layers, and d) the maximum superconducting transition temperature is a similar fraction (~1/3) of the maximum Néel temperature. Recent re-examination of the fulleride (A3C60; A = K, Cs, and Rb) superconductors has revealed that they are also doped antiferromagnetic insulators (Arvanitidis, 2007; Takabayashi, 2009). The phase diagram of cobaltate (NaxCoO2.yH2O) (Foo et al, 2004; Takada et al, 2003) and heavy fermion (e.g. CeCoIn5) (Petrovic et al, 2001) materials have this property as well.
A consistent picture is emerging (Uemura, 2009) that unconventional superconductivity is intimately related to antiferromagnetism (Fig. 10). Imai and coworkers have shown recently (Imai, 2009) that the strength of antiferromagnetic spin fluctuations in FeSe are correlated with the pressure-induced increase in Tc, suggesting a link between spin fluctuations and the mechanism of superconductivity, and have prompted some to suggest that it is a common thread linking organic, heavy-fermion, actinide, cuprate and Fe superconductors (Scalapino, 2009; Uemura, 2009). Spin fluctuation mediated pairing has always been a strong candidate for the mechanism of superconductivity in the cuprates, among the more than twenty candidates (Cho, 2006).
Monthoux and Lonzarich have proposed (Monthoux & Lonzarich, 1999; Monthoux & Lonzarich, 2001; Monthoux et al, 2007) a spin fluctuation based mechanism for superconductivity in systems close to a ferromagnetic or antiferromagnetic instability, making the convincing argument that on the border of long-range magnetic order the dominant interaction channel must be of magnetic origin and depend on the relative spin orientations of the interacting quasiparticles. Superconductivity is predicted to be more robust in doped antiferromagnets vs ferromagnets, and the amplitude of the oscillations in the interaction is enhanced by low dimensionality. For example, the range of temperature and pressure over which superconductivity is observed was increased by about one order of
magnitude in going from cubic CeIn3 to its tetragonal analogues CeMIn5 (M = Rh, Ir or Co) as anticipated by the magnetic interaction model (Monthoux & Lonzarich, 2002). Thus in the search for higher temperature superconductors one should explore the border of antiferromagnetism in a quasi two-dimensional tetragonal system with high characteristic spin fluctuation frequencies.
The conditions favourable for magnetic pairing include: (i) strong quasi two-dimensional antiferromagnetic correlations (large J) for spin singlet pairing and for large amplitude oscillations of the spin-spin interaction (gives small correlation length ξ which is inversely proportional to Tc), (ii) a single band of relatively high characteristic energy scale, and (iii) a crystal structure that enables the repulsive regions of the pairing potential to be optimally neutralized. Favourable Tc‘s can be achieved in layered d-electron systems of moderate electron densities (n) and bandwidths (t) and can be controlled by chemical doping or hydrostatic pressure (Monthoux et al, 2007).
One system which satisfies most of these requirements is the perovskite-type single layer compounds of composition A2MX4 and double layer compounds of composition A3M2X7, where A1+ = K, Rb, Cs, M2+ = Mg, Mn, Fe, Co, Cu, Cd and X = F, Cl, or Br (see Geick, 2001 for a review). In these perovskite-type layer structures the dominant magnetic interaction is the nearest-neighbour Heisenberg exchange within the layers which causes their 2D character. These compounds have a metal ion (M) surrounded by 6 halides (X) in an octahedral arrangement. The magnetic properties depend on the intralayer superexchange interaction (J) mediated by the halide (X) between two M ions. Theory predicts an exponential dependence of J on the nearest neighbour distance (aMXM) and experiments find a power law dependence J(ann) = J(ann,0)(ann/ann,0)-12 for small ann.
A classic perovskite layer compound is La2CuO4 (X = O2-, M = Cu2+, and A = La2+ ) which when appropriately doped (e.g. hole doped by replacing Sr2+ for La3+ or electron doped by replacing La3+ with Nd3+and doping with Ce3+,4+) forms a high temperature superconductor. It must be noted that
In looking for promising hosts, the single layer K2CuF4 and double layer K3Cu2F7 compounds seem to have the right structure, and the d-band of Cu is the highest partially filled band. However, the intralayer interaction is small (J/kB = 11K), produces ferromagnetic order (Feldkemper et al, 1995) and the Cu2+ ions exhibit alternating occupation of z2 - x2 and z2 - y2 hole states unlike the x2 - y2 ordering in La2CuO4 (Fig. 11). However, by inducing distortive changes at pressures larger than 9.5 GPa in the basal plane of the CuF6 octahedra (Ishizuka et al, 1996; Ishizuka et al, 1998) was able to obtain (Fig. 12) antiferromagnetic order in K2CuF4 with x2 - y2 hole orbital overlap, exactly as found in the prototype cuprate superconductor La2CuO4.
SCAN PHASE SPACE: Since the high pressure phase of K2CuF4 is so similar to La2CuO4 in its orbital ordering, structure and magnetic properties, it satisfies the conditions set out by Monthoux and Lonzarich and should become a superconductor when appropriately doped. To obtain the high pressure phase of K2CuF4 one may attempt (a) pseudomorphic growth of
films onto substrates which produce compressive strain. The M-X-M distance in K2CuF4 is ~ 4.124 Å. Therefore, SrLaAlO4 (3.756 Å, -10%), SrTiO3 (3.905 Å, -6%), LaAlO3 (3.821 Å, -8%) and SrLaGaO4 (3.843 Å, -8%) substrates should all produce compressive strain, while MgO (a = 4.212 Å, +1%) should produce tensile strain in epitaxial films. It must be noted that film stresses of more than 10 GPa have been achieved in pseudomorphic Fe layers (Sander, 1999). Epitaxial films of the cuprate superconductors sometimes show enhanced Tc perhaps because of the increase in J. One may also (b) dope smaller cations (Na+, Li+) to create pressure by cation substitution. Substrate or cation-induced decreases in the nearest neighbour M-X-M distances should also exponentially enhance the intra-layer interaction. In addition, we plan to replace the fluorine anion with other halides (Cl, Br, I ) to enhance our understanding of the effect of the exchange interaction on the appearance of superconductivity in the films.
As pointed out earlier, it is difficult to predict
New materials form the basis of new products which drive economic development. Superconducting materials have held great promise for some time because they pass a current without resistance and expel magnetic fields. These properties make them the most sensitive magnetic sensors, best source of large magnetic fields (e.g. for use in medical imaging - MRI), most efficient transmisson lines; and are a leading candidate for high speed quantum computers (B. G. Levi, 2009). However, they have not found widespread application mainly because the materials require cooling to at least -136 degrees C. Finding materials that superconduct at much higher temperatures is now thought to be a realistic goal with the recent discovery of superconductivityin iron arsenide based materials, the observation that a number of superconductors are doped antiferromagnets, and the tremendous progress researchers have made in understanding the physical properties of existing superconductors. These developments have re-ignited the field by offering a path to novel superconductors - explore the transport properties of doped antiferromagnets.
To explore the properties of a large number of samples, a spatial composition spread approach has been developed at Dalhousie to more quickly and efficiently prepare new materials. In a single experiment hundreds of compositions can be studied, whereas a serial preparation approach would take several years. We have described in this chapter the feasibility of the approach to densely map the physical properties of an existing superconductor, La2-xSrxCuO4. To identify novel superconductors, we have proposed that layered fluoride perovskite compounds be screened using the high-throughput resistivity apparatus developed in our labs. To enhance our understanding of existing superconductors we have also shown that the phase diagram can be mapped at very high density to deduce the doping dependence of a feature, the pseudogap onset temperature, which helps determine the class of theories that apply to the cuprate superconductors.
We acknowledge the financial support of the Natural Science and Enginnering Research council of Canada, and use of the facilities of the Institute for Research in Materials. We also acknowledge recent fruitful discussion with Paul Canfield during a recent visit to Dalhousie.