Papers in Physics, vol. 5, art. 050007 (2013) Received: 25 January 2013, Accepted: 10 October 2013 Edited by: S. A. Grigera Reviewed by: E. M. Forgan, School of Physics & Astronomy, University of Birmingham, U.K. Licence: Creative Commons Attribution 3.0 DOI: http://dx.doi.org/10.4279/PIP.050007 www.papersinphysics.org ISSN 1852-4249 Invited review: Graphite and its hidden superconductivity P. Esquinazi1∗ We review experimental results, from transport to magnetization measurements, on dif- ferent graphite samples, from bulk oriented graphite, thin graphite films to transmission electron microscope lamellae, that indicate the existence of granular superconductivity at temperatures above 100 K. The accumulated evidence speaks for a localization of the superconducting phase(s) at certain interfaces embedded in semiconducting crystalline regions with Bernal stacking order. I. Introduction Over the past decade, our interpretation of the magnetic and transport properties of ordered graphite bulk samples has experienced a change respect to the partially accepted general descrip- tion of their intrinsic properties. The description of graphite one finds in the not-so-old literature tells us that it is a kind of (semi)metal with a fi- nite Fermi energy and carrier (electron plus hole) densities per graphene layer at low temperatures n0 ∼ 1010 . . . 1012 cm−2, see e.g., [1–3]. However, real samples are not necessarily ideal, we mean defect-free, and therefore those carrier densities are not necessarily intrinsic of ideal graphite. The exhaustive experience accumulated in gapless and narrow band semiconductors [4] already indicates us how important defects and impurities (not nec- essarily magnetic ones but, for example, hydrogen) are in determining some of the measured properties. Therefore, taking experimental data of real samples as intrinsic, without knowing their microstructure ∗E-mail: esquin@physik.uni-leipzig.de 1 Division of Superconductivity and Magnetism, Institute for Experimental Physics II, Fakultät für Physik und Ge- owissenschaften, Universität Leipzig, Linnéstrasse 5, D- 04103 Leipzig, Germany. and/or defect concentration, was indeed a mislead- ing assumption in the past. This assumption has drastically influenced the description of the band structure of graphite we found nowadays in several books and publications. For example, if graphite has a finite Fermi energy EF (whatever the ma- jority carriers are), as assumed everywhere, up to seven free parameters have to be introduced [2, 5, 6] to describe the apparently ideal band structure of Bernal graphite with the well-known ABAB stack- ing order of the graphene layers. The impact of well defined two-dimensional in- terfaces inside graphite samples [7, 8] had not been realized until recent studies of the transport proper- ties as a function of thickness of the graphite sample provided a link to the microstructure of the sam- ples obtained by transmission electron microscope (TEM) studies. We also have to add the sensi- tivity of the graphite transport properties to very small amount of defects [9]. Those results [7, 9] do not only indicate us that at least a relevant part of the carrier densities measured in graphite is not intrinsic but also that the metallic-like be- havior of the electrical resistance does not reflect ideal, defect-free graphite [10]. An anomalous van- ishing of the amplitude of the Shubnikov-de Haas (SdH) oscillations decreasing the thickness of the graphite samples was published, more than 10 years 050007-1 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi ago, [11] without attracting the necessary atten- tion, although those results already suggested that the SdH oscillations are probably not intrinsic of the graphite structure. These results are supported by the absence of SdH oscillations, i.e., no evi- dence for the existence of a Fermi surface, found recently in bulk oriented samples of high grade and high purity but without internal interfaces [12]. All these results indicate that the internal microstruc- ture of the graphite samples play an important role, a microstructure that was neither characterized nor considered in the discussion of the measured prop- erties of different graphite samples, from highly ori- ented pyrolytic graphite (HOPG) to Kish or natu- ral graphite, even in nowadays literature [6, 13, 14]. What does this have to do with superconduc- tivity? If we start searching for superconductiv- ity in graphite by measuring the behavior of the electrical resistance (R) with temperature (T) and magnetic field (H), for example, it should be clear that the knowledge of the intrinsic, normal state dependence is needed. Otherwise, we may mislead- ingly interpret an anomalous behavior due to, for example, the influence of non-percolative, granu- lar superconducting regions embedded in a (nor- mal state) graphite matrix, as intrinsic of the ma- terial, clearly missing an interesting aspect of the sample. A reader with expertise in superconduc- tivity might not be convinced that such a mistake could be ever made. However, the ballistic trans- port characteristics of the graphene layers in ideal graphite with their huge mobility and mean free path [15–18] provide a high conductivity path in parallel; such that it is not at all straightforward by simple experiments to realize and prove the ex- istence of superconductivity at certain regions in some, not all, graphite samples. One needs indeed to do systematic experiments decreasing the size of the graphite samples (but not too much) to ob- tain clear evidence for the embedded or “hidden” superconductivity. A note on samples: The internal ordering or mosaicity of the graphite crystalline regions in- side commercial HOPG samples is given usually by the grade. For example, the highest ordered pyrolytic graphite samples have is a grade “A”, which means a rocking curve width ∆ ∼ 0.4◦±0.2◦ (“B”, ∆ ∼ 0.8◦, etc.). Interestingly, and due to the contribution of two dimensional highly conduct- ing internal interfaces between crystalline regions [7, 10], the highest grade, i.e., smaller rocking curve width, does not always mean that the used sam- ple provides the intrinsic transport properties of ideal graphite. The characterization of the internal structure of usual HOPG samples, as well as the thickness dependence of R(T) to understand the transport and the magnetic properties of graphite, indicate that these two dimensional interfaces are of importance. The existence of rhombohedral in- clusions [19, 20] (stacking order ABCABC instead of ABAB of the usual Bernal graphite structure) in HOPG as well as in Kish graphite samples can also have a relationship with the hidden superconduc- tivity in graphite, following the theoretical work in Ref. [21]. According to literature (see e.g., Fig. 2- 2 in Ref. [8]), the density of interfaces parallel to the graphene layers in Kish graphite, in regions of several microns length, is notable. Therefore, quantifying the perfection of any graphite sample through the resistivity ratio between 300 K and 4.2 K [8] is not necessarily the best criterion to be used if we are interested on the intrinsic properties of the graphene layers in graphite, because of the high conductivity of the interfaces in parallel to the graphene layers of the sample [10]. Two examples of the interfaces we are referring to can be seen in Fig. 1. On the other hand, commercial HOPG bulk samples are of high purity with average total impu- rity concentrations below 20 ppm. Especially the existence of magnetic impurities are of importance if the Defect-Induced Magnetism (DIM) is the main research issue. Their concentration remains below a few ppm for high grade HOPG samples [22]. The graphite flakes discussed in this work were obtained by exfoliation of HOPG samples of differ- ent batches, by careful mechanical press and rub- bing the initial material on a previously cleaned substrate. As substrate, we used p-doped Si with a 150 nm SiN layer on top. We selected the flakes using microscopic and micro-Raman techniques to check their quality. More details on the preparation can be taken from Ref. [7] and other publications cited below. This review is organized as follows. In the next section we discuss the experimental data for R(T,H) from different graphite samples published in the last 12 years and argue that the first hints on unusual superconducting contribution can be al- ready found in those measurements. In section III, we discuss the anomalous hysteresis in the magne- 050007-2 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi toresistance, a first clear indication for embedded granular superconductivity. Section IV, deals with the Josephson behavior measured in TEM lamellae whereas section V deals with the granular super- conducting behavior found in the magnetization of water-treated graphite powder as well as in bulk HOPG samples for fields normal to the interfaces found inside those samples. In the last section, sec- tion VI, before the conclusion, we discuss possible origins for the superconducting signals on the basis of earlier and recent experimental and theoretical work. II. The behavior of the resistance vs. temperature at different applied magnetic fields In this section, we discuss the behavior of the resis- tance R(T,H) of different HOPG samples includ- ing Kish graphite. The data we present here were taken from [7, 23, 24] and a quick search in litera- ture demonstrates that these data are reproducible and can be found in different publications, see, e.g., [2, 25–27]. Figure 2(a) shows the R(T) for different bulk graphite samples of different grades (rocking curve width) and for one sample (HOPG-1) at zero and under a magnetic field applied normal to the main area, i.e., normal to the graphene planes of the sample. This figure reveals a general behav- ior, namely that the lower the resistivity ρ of the HOPG sample, the more metallic-like its temper- ature dependence. It is appealing to assume that these characteristics, low ρ, low ∆ and the metal- lic behavior are clear signs for more ideal graphite. Therefore, from the measured R(T), we may con- clude that sample HOPG-3 is more ideal than sam- ple HOPG-1 and the latter being more ideal than sample HOPG-2, see Fig. 2(a). This is indeed the usual interpretation found in several reviews in the literature, see e.g., [8, 28]. From a quick look at all the curves in Fig. 2, however, one recognizes a striking similarity be- tween them, although we are comparing different samples with different thickness and some of the curves were measured under a magnetic field ap- plied normal to the graphene layers of the samples; also normal to the interfaces commonly found in some ordered samples [7, 8]. Let us start discussing the metallic-like behavior Figure 1: Transmission electron microscope pic- tures of two different kinds of interfaces and their distribution in HOPG samples. The TEM pictures were taken from two different lamellae, each about 300 nm thick and with the electron beam nearly parallel to the graphene planes of the samples. (a) The interfaces are recognized at the borders of crys- talline regions of different gray colors. Taken from [7]. (b) Interfaces found in a HOPG sample used for magnetization measurements (see section V)that reveals hysteretic behavior in field and tempera- ture. Taken from [29]. of R(T) of sample HOPG-3 in Fig. 2(a). This sam- ple behaves as the HOPG-UC sample shown in (c) at zero magnetic field, having also a maximum at T ∼ 150 K. A “better” metallic character shows the HOPG sample in (b) or the Kish graphite sample in (d), without any maximum in the shown tem- perature range. Is this metallic-like behavior really intrinsic of ideal graphite? 050007-3 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi The following experimental evidence does not support such interpretation: First, for samples from the same batch, the metal- lic character of R(T) vanishes in the whole tem- perature range when the sample thickness is below ∼ 50 nm [7, 11, 27], see Fig. 2(b). Second, the metallic-like behavior vanishes in the whole temperature range after applying a magnetic field of the order of 1 to 2 kOe, an interesting behav- ior known as the Metal-Insulator Transition (MIT) [23, 25, 26], see Figs. 2(a), (c) and (d). Note that such magnetic field strength influences mainly the metallic-like region, see e.g., the change of sample HOPG-1 in Fig. 2(a) at zero and at 1 kOe field, an interesting behavior noted first in [30] and in- terpreted as due to superconducting instabilities. At those field strengths, i.e., H ∼ 1 kOe, the obtained R(T) curves, for samples showing at zero field a metallic-like behavior, resemble the semiconducting-like curves obtained for sam- ple HOPG-2 (Fig. 2(a)) or for samples with small thickness (Fig. 2(b)). At fields higher than a few kOe, the rather large magnetoresistance of graphite starts to play the main role and the R(T) curve in- creases in the whole temperature range. Finally, all these results added to the existence of well defined interfaces in the metallic-like HOPG samples as well as in Kish graphite, with distances in the c−axis direction usually larger than ∼ 30 nm, indicate that the metallic-like behavior is due to the contribution of these interfaces and it is not intrin- sic of the graphene layers of ideal graphite [7, 10]. Therefore, explanations of the MIT based on ideal graphite band models with a large number of free parameters [25, 26] are certainly not the appropri- ate ones. All the different R(T) curves for different sam- ples shown in Fig. 2 and at zero field can be very well understood assuming the parallel contribution of semiconducting graphite paths in parallel to the one from the highly conducting interfaces [10]. The saturation of the resistance at T → 0 K is inter- preted as due to the finite resistance of the sample surfaces (the free one and the one on the substrate) short circuiting the intrinsic behavior of the bulk graphene layers at low enough temperatures. The question is now whether parts of these inter- faces hide superconducting regions. It is certainly appealing to suggest that the huge MIT at fields normal to the interfaces (and below ∼ 2 kOe) is related to Josephson-coupled superconducting re- gions embedded in some of the interfaces. Note that the huge anisotropy of the MIT (parallel fields to the interfaces main plain do not affect the elec- trical transport) already implies that the regions responsible for the MIT must be laying parallel to the graphene layers [31]. Without the knowledge on the existence of these interfaces, an interpretation of the low field MIT based on the influence of su- perconductivity has been discussed in detail in the reviews [24, 32]. In those reviews, one can recog- nize the remarkable similarity between the scaling approaches used to characterize the magnetic-field- induced superconductor-insulator quantum phase transition [33] or of the field-driven MIT in 2D elec- tron (hole) systems [34] and the one obtained for the MIT observed in graphite. As we will see in the next sections, the experimental evidence obtained in the last recent years indicates that granular su- perconductivity exists within some of those inter- faces, indeed. If superconducting patches exist embedded in parts of the interfaces or in other two dimensional regions of the bulk ordered samples, one expects to measure some signs of granular superconductivity, as for example nonlinear I − V curves or hystere- sis in the magnetoresistance. However, this is not really observed in bulk large samples. There are at least two reasons for the apparent absence of these expected phenomena. One is the distribution of the input current between the ballistic channels given by the graphene layers [17, 18], the metallic, normal conducting parts of the interfaces and the regions where the superconducting patches exist. In other words, the usual maximum currents used in transport experiments reported in bulk samples may have been small enough so that the current through the superconducting regions remained be- low the critical Josephson one. The other reason is the experimental voltage sensitivity to measure the possible irreversibility in the magnetoresistance due to the existence of pinned vortices or fluxons. We will see in the next section that part of these prob- lems can be overcome decreasing the sample size; in this way, one obtains the voltage signals from the regions of interest with enough sensitivity. Apart from the large magnetic field sensitiv- ity of the metallic-like resistance measure in bulk graphite samples with interfaces, is there any fur- ther hint for the existence of granular supercon- 050007-4 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi 0 50 100 150 200 250 300 0.4 0.6 0.8 1.0 1.2 1.4 HOPG-3, H = 0 HOPG-1, H = 0 HOPG-1, H = 1 kOe HOPG-2, H = 0 R (T )/ R (3 0 0 K ) TEMPERATURE, T (K) 1 10 100 10 -3 10 -2 10 -1 1 10 100 10 -3 10 -2 10 -1 HOPG-UC TEMPERATURE (K) 1 T 5 kOe 2 kOe 1 kOe 500 Oe 300 Oe 100 Oe H = 0 B A S A L -P L A N E R E S IS T A N C E ( ) K-1 1 T 5 kOe 2.5 kOe 1 kOe 500 Oe 300 Oe H = 0 (c) (a) (b) (d) Figure 2: (a) Normalized resistance vs. temperature for three different HOPG bulk samples. The bottom, metallic-like curve corresponds to the sample HOPG-3, the curves above correspond to HOPG- 1 (H = 0), HOPG-1 (H = 1 kOe), and HOPG-2. The grade and resistivity values are HOPG-1 (∆ = 1.4◦, resistivity at 300 K ρ(300K) = 45 µΩcm), HOPG-2 (∆ = 1.2◦, ρ(300K) = 135 µΩcm) and HOPG-3 (∆ = 0.5◦, ρ(300K) = 5 µΩcm). Taken from [23]. (b) Similar to (a) but for HOPG samples from the same batch but with different size, namely (thickness × length × width) L5: 12 ± 3 nm, 27 µm, 14 µm, L2A: 20 ± 5 nm, 5 µm, 10 µm, L8A: 13 ± 2 nm, 14 µm, 10 µm, L8B: 45 ± 5 nm, 3 µm, 3 µm, L7: 75 ± 5 nm, 17 µm, 17 µm, HOPG: 17 ± 2 µm, 4.4 mm, 1.1 mm, taken from [7]. (c) and (d) Resistance of bulk graphite samples vs. temperature at different applied fields normal to the graphene layers. The sample in (c) is a HOPG bulk sample from Union Carbide of grade A and the sample in (d) is Kish graphite. Taken from [24]. ductivity in those R(T) curves? Yes, this hint is related to the thermally activated function (∝ exp(−Ea/kBT) with Ea, a sample dependent effec- tive thermal barrier ∼ 30 K) one needs in order to fit the metallic-like contribution below T ∼ 200 K [10]. This function is relevant in spite of only a fac- tor five increase of the resistance between low and high temperatures, see Fig. 2. Skeptical readers can convince themselves about its relevance tak- ing a similar example, as the exponential function 050007-5 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi used to fit the increase, by a similar factor, of the ultrasonic attenuation with temperature below Tc in conventional superconductors. We note that this exponential function has already been used to de- scribe the increasing resistance of bulk graphite samples with temperature and it was speculated to be related to some superconducting-like behav- ior in graphite [32]. It is clear that this func- tion is not the usual one, one expected for met- als or semimetals and that cannot be understood within the usual electron-phonon interaction mech- anisms, nor in two dimensions. A similar depen- dence has been observed in granular AlGe [35], which shows for a particular Al concentration a superconductor-semiconductor transition similar to that reported in Ref. [10] or, after an appropri- ate scaling in temperature, to some of the curves shown in Fig. 2. The observed thermally acti- vated behavior might be understood on the ba- sis of the Langer-Ambegaokar-McCumber-Halperin (LAMH) model [36, 37] that applies to narrow su- perconducting channels in which thermal fluctua- tions can cause phase slips. This interpretation gets further support from the evidence we discuss in the following sections. III. Hysteresis in the magnetoresis- tance In order to reveal by transport measurements the existence of granular superconductivity in some re- gions of the graphite samples, we need to increase the sensitivity of the measured voltage to those re- gions. To achieve this, we decrease the size of the sample enhancing in this way the probability to get some measurable influence of this phenomenon in the voltage. The work in Ref. [38] reported the first observations of an anomalous irreversible behavior in the magnetoresistance (MR) in a few tens of nm thick and several micrometer large multigraphene samples. Hysteresis in the magnetoresistance is a key evidence on the existence of either magnetic order (domains with their walls, for example) or vortices/fluxons and therefore on the existence of superconductivity. Because defects as well as hy- drogen can trigger magnetic order in graphite, a first attempt would be to relate the measured hys- teresis in the MR with the existence of magnetic order and magnetic domains, for example. How- ever, the data exhibited anomalous hysteresis loops in the MR [38], similar to those observed in granu- lar superconductors with Josephson-coupled grains [39–41]. The anomalous hysteresis was observed only for magnetic fields perpendicular to the planes, whereas in the parallel to the planes direction, the MR remains negligible. This fact already points out to a remarkable large anisotropic response of the su- perconducting phase(s) in agreement with the hy- pothesis that these superconducting regions might be embedded in some of the interfaces found inside some bulk graphite samples [7, 8]. The amplitude of the hysteresis in the MR reported in Ref. [38] vanishes at temperatures T ∼ 10 K, clearly below the temperature at which the resistance shows a maximum, as it is the case for samples HOPG-1 in Fig. 2(a) or sample L2A in Fig. 2(b). It is clear that thermal fluctuations can prevent the establishment of a coherent superconducting state in parts of the sample and therefore zero re- sistance state is not so simple to be achieved if the superconducting distribution is a mixture of super- conducting patches at the interfaces and these are embedded in a multigraphene semiconducting ma- trix. Moreover, we should take also into account that the voltage electrodes are usually connected at the top surface of the graphite samples picking the voltage difference coming from a non-negligible normal conducting path. One possibility to increase the sensitivity of the measured voltage to the field hysteresis these re- gions produce is to make a constriction in the middle of the two voltage electrodes, see inset in Fig. 3(a). In this case, we expect a locally nar- rower distribution of superconducting and normal regions at the constriction such that averaging ef- fects should be less important. Simultaneously, through the constriction the main part of the volt- age drop depends mostly on the region at the con- striction, see Fig. 2(c) in Ref. [16]. Then, a higher sensitivity to the superconducting paths can be achieved in case they remain at or near the con- striction. This idea has been successfully realized in [42] and its main results will be reviewed in this section. Let us take two slightly different samples, 1 and 2, with R(T)−curves as shown in Fig. 3(a). The aim of the experiment is to study the hysteresis in the MR those samples might show below the temperature at which a maximum in the resistance 050007-6 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi -200 -100 0 100 200 0 1 2 3 R (B )- R (0 )/ R (0 ) (1 0 -3 ) Sample 2 W = 3 µm T= 10 K Applied field H(Oe) (d) 0 10 20 30 40 50 0.00 0.05 0.10 0.15 0.20 0.25 0.30 without constriction without constriction 6 µ m 4 µ m ∆ M R ( 1 0 -3 ) Temperature T(K) (c) -200 -100 0 100 200 0 4 8 12 16 -400 -200 0 200 400 -3 -2 -1 0 1 2 3 ∆ M R ( 1 0 -3 ) Sample 1 W= 4 µm T= 2 K M R = R (B )- R (0 )/ R (0 ) (1 0 -3 ) Applied field H(Oe) (b) 0 50 100 150 200 250 300 44 46 48 50 52 54 56 58 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 R e s is ta n c e R (Ω ) Temperature T(K) Sample 2 Sample 1 (a) Figure 3: (a) Resistance vs. temperature, without constrictions and at zero applied field, for two graphite flakes of size (distance between voltage electrodes × width × thickness) for sample 1 (2): 13×16×0.015 (2.6×6×0.040) µm3. The observed temperature dependence remains for all constrictions widths. The inset shows a scanning electron microscope picture of sample 1 with a constriction width of 4.3 µm between the two voltage electrodes. The scale bar is 5 µm. (b) Magnetoresistance (MR) vs. applied magnetic field for sample 1 with a 4 µm constriction width and at 2 K. The input current was 1 µA. Note the clear hysteresis in the MR when the field is swept from |Hmax| = 1000 Oe. The inset shows the difference ∆MR between the curve obtained starting from Hmax = +1000 Oe and the return curve measured from Hmin = −1000 Oe. (c) The absolute difference between the two MR curves of the hysteresis loop obtained at a fixed magnetic field of 16.6 Oe for sample 1 without constriction (?) and for two different constriction widths. The figure also shows the corresponding data for another graphite sample without constrictions (�) from [38]. (d) Magnetoresistance measured from a starting maximum field of 1.4 kOe at 10 K for sample 2 with a constriction width of 3 µm. Some of the figures and the data were taken from [42]. is measured, in case that maximum is related to the Josephson coupling between superconducting regions. Figure 3(b) shows one example of the anomalous hysteresis loop in the MR. The going down curve (from high, positive to low, negative fields, red arrow), for example, runs below the go- ing up curve (green arrow) in the same quadrant as the field sweep was started, showing a minimum at positive fields of the order of 20 Oe, see also sim- ilar curves in Ref. [38]. To present the anomalous behavior clearly, the inset in Fig. 3(b) shows the difference between the two curves. This difference 050007-7 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi is in clear contrast to the usual hysteresis in super- conductors as well as ferromagnets [38, 39], where the minimum (or maximum) in the MR is observed in the opposite field quadrant, and the increasing field resistance curve is usually below the decreas- ing field one. Figure 3(c) shows the temperature dependence of the difference in the MR between the decreasing and increasing field curves at a fixed magnetic field for sample 1, without and with two constrictions. The results show that the smaller the constriction width, the higher the temperature at which the anomalous hysteresis is observed, decreasing below the sensitivity limit at T > 50 K for a constriction width of 4 µm, whereas the maximum in the R(T) curve is at ∼ 70 K, see Fig. 3(a). The absence of any hysteresis in the MR for sample 2 with a constriction width of 3 µm and at T = 10 K indi- cates that the hysteresis does not come from some artifact due to the used focused ion beam method [42, 43] or due to an artifact in the measurement of the real field applied to the sample. As expected from the R(T) curve, see Fig. 3(a), sample 2 shows the anomalous hysteresis in the MR at lower tem- peratures than sample 1 [42]. Summarizing this section, the observation of the anomalous hysteresis in the MR – together with the MIT and the relatively large MR at tempera- tures below the maximum in R(T) – provides al- ready striking hints that granular superconductiv- ity is at work in some regions of these samples. The increase in the temperature region where the hysteresis is observed, decreasing the constriction width, demonstrates the problem of current averag- ing and voltage sensitivity limits usual experiments with large samples have. IV. Direct evidence for Josephson behavior in the transport prop- erties of graphite: Measure- ments in TEM lamellae If the embedded interfaces (or some other quasi two dimensional regions) inside the measured graphite samples have superconducting properties, the best way to check them would be contacting electrodes as near as possible to those interfaces or interface regions and study the behavior as a function of any useful parameter one can take to influence their re- sponse. It should be clear that one cannot simply open the graphite sample at the interface and put voltage electrodes at the open surfaces of the inter- face, simply because it will not remain anymore. A tentative approach to put the contacts as near as possible to an interface has been done in Ref. [46]. Indeed, the observed behavior at low temperatures and as a function of magnetic field appeared to be superconducting-like. A better and appealing evidence for the super- conducting behavior embedded in some graphite or- dered samples can be obtained by trying to locate the voltage electrodes directly at the inner edges of the interfaces. The work in Ref. [45] prepared TEM lamellae from bulk HOPG samples and using lithography and focused ion beam techniques, cur- rent and voltages electrodes at different positions of the samples were prepared. In this way, one tries to contact several of those interface edges simulta- neously, as shown in Fig. 4(a). We note, however, that a thin surface layer of disordered graphite ex- ists due to the Ga+ ion irradiation used to cut the lamella from the bulk HOPG sample. This layer has a much larger resistance than the one of the graphene layers or of the interfaces [43] and there- fore the input current goes through the lowest re- sistance path as well as the voltage electrodes pick up the response of the graphite sample with its in- terfaces. One can see this comparing first the R(T) curves obtained at large enough currents in the lamellae (Fig. 4(c)) with those of graphite samples with top electrodes (Fig. 2). The fact that a zero re- sistance state (minimum voltage noise ± 5 nV upon sample) is obtained at low currents with I−V char- acteristic curves that resemble the one expects for Josephson coupled grains leaves little doubt about the origin of the obtained signals. In the TEM pic- ture of Fig. 4(b), one can see the graphite single crystalline regions (different gray colors) oriented differently between them about the common c−axis and having well defined two dimensional interfaces, as high resolution TEM studies revealed [47]. Figure 4(c) shows the voltage vs. temperature measured in a TEM lamella of oriented graphite [45] at different input DC currents, from 100 nA to 10 µA. The clear sharp transition, observed at ∼ 150 K at the lowest current, shifts to lower tem- peratures increasing the input DC current. For the largest input currents, the temperature dependence of the resistance of the contacted lamella shows a 050007-8 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi o (c) (d) Figure 4: (a) Scanning Electron Microscopy (SEM) image of a lamella of 300 nm thickness on a Si/SiN substrate where the yellowish colored areas are the electrodes. A four-point configuration has been prepared with the outer electrodes used to apply current and the inner ones to measure the voltage drop. The c−axis runs parallel to the substrate surface and normal to the current direction. (b) Transmission Electron Microscopy (TEM) image of a HOPG lamella. The different brightness corresponds to a different orientation within the a − b plane of the crystalline regions with thickness > 30 nm. (c) Voltage vs. temperature at different input currents for a lamella of ∼ 800 nm thickness and with Van der Pauw contact configuration. (d) Current-Voltage characteristics at different temperatures for a lamella of ∼ 300 nm thickness in reduced coordinates, where R is the normal state resistance, I the input current, and Ic the critical Josephson current. The continuous curves are fitted to the model proposed in Ref. [44] with Ic(T) as the only free parameter. Figures taken from [45]. maximum or follows the intrinsic semiconducting behavior of the graphene layers. This behavior al- ready suggests the existence of high temperature granular superconductivity at some parts of the sample. The study reported in Ref. [45] shows that the transition temperature depends on the pre- pared sample. This indicates a sample dependent distribution of the superconducting regions and/or some influence of the preparation process or sam- ple size on the superconductivity [47]. We also note that the observed sharp decrease in the mea- sured voltage does not necessarily indicate the criti- 050007-9 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi cal temperature of the superconducting regions but the temperature below which a percolative granu- lar system shows negligible resistance due to the Josephson coupling at the used input current. Current-voltage characteristic curves at differ- ent temperatures and in different lamellae obtained from different HOPG samples have been studied in Ref. [45]. An example of this I − V curves at three temperatures is shown in Fig. 4(d) obtained for a different lamella. The curves follow the ex- pected dependence for a Josephson junction [44] with a temperature dependent critical current, the only free parameter in the fit. Further evidence that speaks for a superconduct- ing origin of the I − V curves is given by the ex- pected detrimental effect of a magnetic field on the superconducting state. This effect can be due to an orbital depairing effect or due to the alignment of the electron spins at much higher fields, in case of singlet coupling. The effect of a magnetic field ap- plied normal and parallel to the interfaces has been studied in detail for thick and thin lamellae in Ref. [45]. Upon sample size (thickness, i.e., width of the graphene planes inside the lamella) the observed effects are from the usual vanishing of the zero re- sistance state or no effect at all for thin lamellae. A magnetic field of a few kOe applied normal to the interfaces is enough to destroy the Josephson coupling at low temperatures, an effect compatible with the MIT observed in several graphite samples, see section II. Whereas a field applied parallel to them does not influence the I −V curves at all, a fact that speaks for the two dimensionality of the superconducting regions. Nevertheless, the influence of a magnetic field in HOPG samples with interfaces is not as “simple” as in conventional superconductors. For high fields applied normal to the interfaces, the I −V curves show a recovery to the zero resistance state. The observed reentrance appears to be related to the magnetic-field driven reentrance observed at low temperatures in the longitudinal resistance at high enough magnetic fields [48]. This interesting behav- ior as well as the insensitivity of the I −V curves to magnetic fields in very thin lamellae [45] deserve further studies. We would like to note here that the possible ef- fects of a magnetic field on the superconducting state of quasi two-dimensional superconductors, or in case the coupling does not correspond to a sin- glet state, are not that clear as in conventional superconductors. For example, results in two dif- ferent two-dimensional superconductors, including one produced at the interfaces between non super- conducting regions [49], show that superconductiv- ity can even be enhanced by a parallel magnetic field. In case the pairing is p−type [50], the influ- ence of a magnetic field is expected to be qualita- tively different from the conventional, singlet cou- pling behavior [51,52] with even an enhancement of the superconducting state at intermediate fields. In case the London penetration depth is much larger than the size of the superconducting regions at the interfaces of our lamellae or if the superconducting coherence length is of the order or larger than the thickness of the lamella, the influence of a magnetic field should be less detrimental. Through these studies, and taking into account that in samples without these interfaces no signa- ture of superconducting or metallic-like behavior has been observed (see also section V) it is ap- pealing to suggest that superconductivity is some- where hidden at some of those interfaces or inter- face regions. It should be also clear that not all those interfaces have superconducting regions with similar critical parameters. Those interfaces are formed during the preparation of the HOPG sam- ples based on treatments at very high temperatures (T > 3400◦C) and high pressures (P ∼10 kg/cm3) and in a non-systematic way [8]. Actually, they are not at all an aim of the production but actually the opposite, they should be avoided in order to enhance the crystal perfection of the bulk HOPG material. It is even possible that, upon the proce- dure used to control the structure and texture of the graphite sample, the near surface region, for example, can have a different degree of graphiti- zation as inside the bulk HOPG sample [8]. This means that one may obtain different results from different parts of the same HOPG sample. There- fore, disconcerting situations and an apparent lack of reproducibility are preprogrammed in case the research studies are done without taking care of the internal microstructure of the studied samples. V. Magnetization measurements In this section, we present and discuss magneti- zation measurements done in bulk HOPG samples 050007-10 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi with and without embedded interfaces and in water treated graphite powders. One of the main prob- lems in interpreting magnetization data for fields applied parallel to the c−axis of the graphite struc- ture, i.e., normal to the graphene layers and in- terfaces, is the need of subtraction of a large dia- magnetic background. Due to the small amplitude of the superconducting-like signals in the studied samples, the subtraction of this linear in field back- ground is not so simple, because it is not known with enough certainty to obtain the true field hys- teresis after its subtraction. That means that we always have a certain arbitrariness in the shape of the obtained field hysteresis, a situation that will improve with the increase of the amount of material responsible for those superconducting-like signals. The small SQUID signals of interest im- ply that one should take additional efforts to min- imize or rule out possible SQUID artifacts [53–55]. Therefore, systematic studies of samples of differ- ent or equal geometry and magnetic background, with and without interfaces, are necessary. Taking into account: the overall shape of the hys- teresis, the slope of the virgin curve at low fields where the subtraction does not affect too much, and the overall experience with ferromagnetic graphite [22, 56], one can rely to a certain extent on the ob- tained hysteresis shape. Certainly, not only the field hysteresis but also other evidence one gets from magnetization measurements as, e.g., the re- manence at zero field as a function of the maxi- mal field applied (see for example measurements for YBa2Cu3O7 in Ref. [57]) and the hysteresis in temperature dependent measurements helps to convince oneself about the existence of some kind of granular superconductivity. The hysteresis be- tween the field cooled (FC) and zero-field cooled (ZFC) curves can help to discern between a su- perconducting or ferromagnetic-like behavior. The most obvious evidence that speaks against a simple ferromagnetic order of the hysteresis observed as a function of temperature and field is the two dimen- sionality of the obtained hysteretic signals [29], i.e., the superconducting-like signals are mainly mea- sured for fields normal to the interfaces. This fact is not compatible with any kind of magnetic order including shape or magneto crystalline anisotropy, whatever large they might be. We note that the ferromagnetic response of graphite due to DIM is mostly measured for fields parallel to the graphene -60 -40 -20 0 20 40 60 -3 -2 -1 0 1 2 3 Applied field µ 0 H (mT) M ag ne tiz at io n M ( 10 - 4 e m u/ g) HOPG-1 HOPG-2 T = 300 K (a) -5 0 5 WTGP M (10 - 5 em u/g) 0 100 200 300 400 -2 0 2 4 6 8 10 12 (b) 0.5 T (b.a.) 0.5 T (a.a.) m F C -m Z F C ( 10 -6 e m u) Temperature T(K) 0 20 40 60 80 4 T (b.a.) 4 T (a.a.) Figure 5: (a) Magnetization of two HOPG bulk samples (HOPG-1 and HOPG-2) after subtraction of a diamagnetic background and of water treated graphite powder (WTGP, right y−axis) at 300 K. The HOPG-2 sample shows no hysteresis in con- trast to the other two samples. (b) Temperature dependence of the difference between FC and ZFC magnetic moments of the HOPG-1 sample before (b.a.) and after (a.a.) warming the sample up to ' 600 K, at two constant applied fields, 0.5 T (left y−axis) and 4 T (right y-axis). The field was ap- plied always normal to the interfaces or graphene planes of the samples. Data taken from [29]. layers, parallel to the main area of the samples [22]. i. Bulk graphite samples Figure 5(a) shows the field hysteresis, after sub- traction of the corresponding diamagnetic linear background, at 300 K of two bulk HOPG sam- ples, HOPG-1 and HOPG-2 and a water treated graphite powder (WTGP) (right y−axis). A TEM characterization of the internal microstructure of 050007-11 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi the HOPG-1 sample shows clear evidence for well defined interfaces running parallel to the graphene layers, in contrast to the HOPG-2 sample [29], see Fig. 1(b). These results clearly indicate that the origin of the hysteresis is related to the exis- tence of the interfaces in the HOPG-1 sample. The absence of the hysteresis in the HOPG-2 sample, which has a similar diamagnetic background and overall geometry as the HOPG-1 sample, also indi- cates that the hysteresis is not due to an obvious SQUID artifact or an artifact in the background subtraction. The field hysteresis is similar to that of WTGP. The narrowing of the hysteresis observed at high fields is expected for granular supercon- ductors [58, 60–62]. From the hysteresis, as well as measuring the remanent magnetic moment as a function of the applied field [29], one obtains the characteristic Josephson critical fields hJc1(T) and hJc2(T) with values similar to the WTGP [58] and a similar ratio hJc2(T)/h J c1(T) ∼ 3 [29]. Figure 5(b) shows the magnetic moment hystere- sis in temperature (FC minus the ZFC curve) for the HOPG-1 sample as received (b.a.) and after sweeping the temperature up to 500 K (a.a.) [29], at two applied fields. We would like to stress the fol- lowing features: The hysteresis for the as-received sample starts from the turning point (390 K) and it is positive. The hysteresis in temperature at both applied fields are qualitatively similar, showing a crossing to negative values at low temperatures. Larger ZFC values (smaller in absolute value) than FC ones in the magnetic moment are usually not observed, neither in superconductors nor in ferro- magnets and it does appear to be a SQUID arti- fact [29]. This negative hysteresis in temperature would suggest that the superconducting properties can be enhanced to some extent under a magnetic field, an effect that might be related to the reen- trance we have shortly mentioned in section IV. A slight annealing of the HOPG-1 sample of less than one hour at ∼ 500 K changes drastically the ob- served hysteresis for both fields (open symbols in Fig. 5(b)). The hysteresis appears to be shifted to lower temperatures but with negative values at high temperatures and high fields. We note that anneal- ing at similar temperatures for several hours pro- duced a decrease in the overall hysteresis observed in WTGP (see supporting information of Ref. [58]). At the state of this research, it is unclear whether pinning properties of vortices and/or of fluxons or 0 100 200 300 -2 0 2 4 6 8 10 0 10 20 30 m F C - m Z F C ( 10 -5 e m u) Temperature T(K) 10 mT 20 mT 30 mT 50 mT 0.1 T 0.15 T 0.2 T 0.3 T 40 (b) 0.5 T 1 T -40 -20 0 20 40 -4 -2 0 2 4 M ag ne tiz at io n M ( 10 -5 e m u/ g) Applied Field µ 0 H (mT) S1 S2 S3 5 K (a) Figure 6: (a) Field hysteresis at 5 K for a maxi- mum applied field of 40 mT for the water treated graphite powder (S1), the same powder but after pressing it in a pellet with a pressure of 18±5 MPa (S2) and after pressing it again with a pressure of 60±20 MPa (S3). The corresponding diamagnetic linear backgrounds were subtracted from the mea- sured data. (b) Difference between the FC and ZFC curve at different applied fields for a water treated graphite powder. Data taken from Ref. [58]. the existence of different superconducting phases play a main role in the hysteresis that is observed for fields applied normal to the interfaces. ii. Water treated graphite powder The work of Ref. [58] reports on the magnetic re- sponse of WTGPs. The main message of that work is that the WTGP shows a hysteretic behavior in field and temperature compatible with granu- lar superconductivity. As an example, we show in Fig. 6(a) the field hysteresis at 5 K of a WTGP (S1, lose powder without applying significant pres- 050007-12 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi sure) and the same WTGP but after pressing it into a pellet with two different pressures (S2,S3). After the diamagnetic background subtraction, the field hysteresis is similar to that obtained for bulk HOPG sample with interfaces, see Fig. 5(a) for sim- ilar data but at 300 K. The fact that the hystere- sis vanishes after applying pressure to the powder rules out simple SQUID artifacts (the diamagnet background does not diminish after making a pel- let from the graphite powder, but the contrary) and also it rules out that the hysteresis is due to a fer- romagnetic response due to impurities. Figure 6(b) shows the difference in the magnetic moment between the ZFC and FC curves, as in Fig. 5(b). The behavior of this difference as a func- tion of the applied field appears to be compatible with the one expected for granular superconductors [58]. Note the following features: The hysteresis in- creases at all T for fields µ0H . 50 mT, showing a maximum near the turning point of 300 K, similar to the HOPG-1 sample in the as-received state, see Fig. 5(b). At fields 0.1 T . µ0H . 0.2 T the dif- ference decreases at all T and remains rather field independent. At higher fields, however, it increases showing a shift of the crossing point (from nega- tive to positive values) to higher T . This behavior is at odds to the one expected for ferromagnets, even for ferromagnetic nanoparticles [63] as well as for superconductors with a pinning force that de- creases with applied field in the shown field range. From the results in [58], and using basic concepts of vortex pinning, we would then conclude that if an upper critical field exists, then it should be clearly larger than 7 T in the temperature range of the figure. In spite of some interesting differences between the behavior obtained for bulk HOPG and WTGP, the similarities already suggest that the water treatment helps to produce a certain amount of interfaces between graphite grains, being the ori- gin for the whole hysteresis. Thermal annealing as well as pressing the WTGP are detrimental indi- cating that defects and/or hydrogen or oxygen at the interfaces could play an important role in the observed phenomena. VI. Discussion Superconductivity in carbon-based systems is a rather old, well recognized fact. This phenomenon was probably first observed in the potassium inter- calated graphite C8K [64] back in 1965. Since then, a considerable amount of studies reported this phe- nomenon in carbon-based systems, reaching critical temperatures Tc ∼ 10 K in intercalated graphite [65, 66] and above 30 K - though not percolative - in some HOPG samples [59] as well as in doped graphite and amorphous carbon systems [67–70]. Traces of superconductivity at Tc = 65 K have been recently reported in amorphous carbon powder that contained a small amount of sulfur [71]. Supercon- ductivity was found also in carbon nanotubes with Tc = 0.55 K [72] and 12 K [73] or possibly even higher critical temperatures [74, 75]. Superconduc- tivity with Tc ∼ 4 K in boron-doped diamond [76] and in diamond films with Tc ∼ 7 K [77] belong also to the recently published list of carbon-based superconductors. We should note, however, that superconductivity at room temperature in a disor- dered graphite powder has been already reported in 1974 [78], see also [79], a work that did not attract the necessary attention in the community. Whether quasi two dimensional interfaces play a role in the above mentioned carbon-based su- perconductors, one can probably rule out only for the intercalated graphite and doped diamond com- pounds, where the three dimensional superconduc- tivity is characterized by a relatively low critical temperature. We may speculate that the traces of superconductivity found in doped amorphous car- bon, disordered or ordered graphite powders may be related to some interfaces between well ordered graphite regions. The experience of the high tem- perature superconducting oxides already suggests that two dimensionality is advantageous to achieve higher critical temperatures. Apart from the usual transport and magneti- zation measurements used to characterize the su- perconducting state, there are scanning tunnel- ing spectroscopy (STS) results obtained on cer- tain disordered regions of a HOPG surface at T = 4.2 K that revealed an apparent energy gap ∼ 100 meV [80]. Although the overall curves resem- ble a superconducting-like density of states, the au- thors suggested that the gap originates from charg- ing effects. See further STS results and the discus- 050007-13 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi sion in [70]. Theoretical works that deal with superconduc- tivity in graphite as well as in graphene have been published in recent years. For example, p-type su- perconductivity has been predicted to occur in in- homogeneous regions of the graphite structure [50] or d−wave high-Tc superconductivity [81] based also on resonance valence bonds [82], or at the graphite surface region with rhombohedral stack- ing due to a topologically protected flat band [83]. For the graphite structure, the experimental evi- dence obtained in the last years suggests that high temperature superconductivity exists at certain in- terfaces or interface regions within the usual Bernal structure although the structure of the supercon- ducting regions remains unknown. One can further speculate that due to the high carrier concentra- tion that can be localized at those interfaces, they should be predestined to play a role in triggering superconductivity. Following a BCS approach in two dimensions (with anisotropy), for example, a critical temperature Tc ∼ 60 K has been estimated if the density of conduction electrons per graphene plane increases to n ∼ 1014 cm−2, a density that might be induced by defects and/or hydrogen ad- atoms [84] at the interfaces, or by Li deposition [85]. Further predictions for superconductivity in graphene support the premise that n > 1013 cm−2 in order to reach Tc > 1 K [86, 87]. On the other hand, the possibility to have high temperature su- perconductivity at the surface of or in the rhombo- hedral graphite phase [21, 83] – a phase that some- times is found in graphite samples [19, 20] – stim- ulates further careful studies of these hidden in- terfaces. In the last years, superconductivity has been found at the interfaces between oxide insula- tors [88] as well as between metallic and insulating copper oxides with Tc & 50 K[89]. Also, interfaces in different Bi bicrystals show superconductivity up to 21 K, although Bi bulk is not a superconductor [90, 91]. Finally, we think that some of the interfaces are also the origin for the metallic-like behavior of graphite samples as well as for the quantum Hall effect (QHE) found in some HOPG samples [48,92]. Because the existence, density as well as the intrin- sic properties of these interfaces depend on sample, we can now understand why the reproducibility of the QHE in bulk HOPG samples is rather poor. VII. Conclusion In this review, we have discussed the following ex- perimental evidence: Firstly, the temperature and magnetic field depen- dence of the electrical resistance of bulk and thin films of graphite samples and its relation with the existence of two dimensional interfaces. Secondly, the Josephson behavior of the current- voltage curves with an apparent zero resistance state at high temperatures in especially made TEM lamellae. Thirdly, the anomalous hysteresis in the magne- toresistance observed in graphite thin samples as well as its enhancement restricting the current path within the sample. Finally, the overall magnetization of bulk graphite samples, with and without interfaces, as well as wa- ter treated graphite powders. All this experimental evidence as a whole in- dicates the existence of superconductivity located at certain interfaces inside graphite samples. Al- though we cannot rule out other interpretations for some of the observations discussed in this work, the whole evidence suggests that superconductiv- ity should be the origin for all the phenomena dis- cussed here. Clearly, the situation is still highly un- satisfactory because several open questions remain, namely, the characteristics of the superconducting phase(s), from the structure to the main supercon- ducting parameters, as “simple” as the critical tem- perature and critical fields, the coherence and pen- etration lengths, etc. It is clear that further stud- ies are necessary in the future but the overall work done until now shows us the way to go. Acknowledgements - The author acknowledges the support provided by the Deutsche Forschungs- gemeinschaft under contract DFG ES 86/16-1 and the ESF-Nano under the Graduate School of Nat- ural Sciences “BuildMona”. The results presented in this review were part of the Ph.D. thesis of Heiko Kempa (section II), Srujana Dusari (section III) and Ana Ballestar (section IV) as well as the mas- ter thesis of Thomas Scheike (section V) done in the Division of Superconductivity and Magnetism of the Institute for Experimental Physics II of the University of Leipzig. The author thanks Dipl. 050007-14 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi Kris. Annette Setzer, Dr. José Barzola-Quiquia and Dr. Winfried Böhlmann for their experimen- tal assistance and support. The permanent sup- port as well as the discussions with Nicolás Garćıa are gratefully acknowledged. Special thanks go to Yakov Kopelevich with whom we started in the year 1999 and in a rather naive way the research of a new and unexpected world behind graphite. Note added in proof Since the submission of this manuscript, some new works related to the subject of this review were published. Tight-binding simulations done in Ref. [93] support the work done in Ref. [21] and found that surface superconductivity is robust for ABC stacked multilayer graphene, even at very low pair- ing potentials. Through the observation of persis- tent currents in a graphite filled ring-shaped con- tainer immersed in alkanes, the author in Ref. [94] claimed possible room temperature superconduc- tivity. For completeness, we include them in the reference list. [1] J W McClure, Energy band structure of graphite, IBM J. Res. Dev. 8, 255 (1964). [2] B T Kelly, Physics of graphite, Applied Sci- ence Publishers, London (1981). [3] A Grüneis, C Attaccalite, T Pich- ler, V Zabolotnyy, H Shiozawa, S L Molodtsov, D Inosov, A Koitzsch, M Knupfer, J Schiessling, R Follath, R Weber, P Rudolf, R Wirtz, A Rubio, Electron-electron correla- tion in graphite: A combined angle-resolved photoemission and first-principles study, Phys. Rev. Lett. 100, 037601 (2008). [4] I M Tsidilkovski, Electron spectrum of gapless semiconductors. Springer Series in Solid-State Sciences Vol. 116, Springer Verlag (1997). [5] R O Dillon, I L Spain, J W McClure, Elec- tronic energy band parameters of graphite and their dependence on pressure, temperature and acceptor concentration, J. Phys. Chem. Solids 38, 635 (1977). [6] J M Schneider, M Orlita, M Potemski, D K Maude, Consistent interpretation of the low- temperature magnetotransport in graphite us- ing the Slonczewski-Weiss-McClure 3D band- structure calculations, Phys. Rev. Lett. 102, 166403 (2009). [7] J Barzola-Quiquia, J L Yao, P Rödiger, K Schindler, P Esquinazi, Sample size ef- fects on the transport properties of mesoscopic graphite samples, Physica Status Solidi A 205, 2924 (2008). [8] M Inagaki, New carbons: Control of structure and functions, Elsevier (2000). [9] A Arndt, D Spoddig, P Esquinazi, J Barzola- Quiquia, S Dusari, T Butz, Electric carrier concentration in graphite: Dependence of elec- trical resistivity and magnetoresistance on de- fect concentration, Phys. Rev. B 80, 195402 (2009). [10] N Garćıa, P Esquinazi, J Barzola-Quiquia, S Dusari, Evidence for semiconducting behav- ior with a narrow band gap of Bernal graphite, New J. Phys. 14, 053015 (2012). [11] Y Ohashi, K Yamamoto, T Kubo, Shubnikov - de Haas effect of very thin graphite crystals, In: Carbon’01, An International Conference on Carbon, Pag. 568, The American Carbon Society, Lexington, KY, United States, (2001). [12] B C Camargo, Y Kopelevich, S B Hubbard, A Usher, W Böhlmann, P Esquinazi, Effect of structural disorder on the quantum oscillations in graphite (unpublished). In this work the authors show that in certain HOPG samples (SPI) of high grade, the density of interfaces is much lower than in, for example, Advanced Ceramics HOPG ZYA samples. In this new HOPG samples basically no SdH oscillations are found and the temperature dependence of the resistance shows a semiconducting behav- ior with saturation a low temperatures (2013). [13] M Orlita, C Faugeras, G Martinez, D K Maude, M L Sadowski, J M Schneider, M Potemski, Magneto-transmission as a probe of Dirac fermions in bulk graphite, J. Phys: Cond. Mat. 20, 454223 (2008). 050007-15 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi [14] N A Goncharuk, L Nádvorńık, C Faugeras, M Orlita, L Smrčka, Infrared magnetospec- troscopy of graphite in tilted fields, Phys. Rev. B 86, 155409 (2012). [15] J C González, M Muñoz, N Garćıa, J Barzola- Quiquia, D Spoddig, K Schindler, P Esquinazi, Sample-size effects in the magnetoresistance of graphite, Phys. Rev. Lett. 99, 216601 (2007). [16] N Garćıa, P Esquinazi, J Barzola-Quiquia, B Ming, D Spoddig, Transition from ohmic to ballistic transport in oriented graphite: Mea- surements and numerical simulations, Phys. Rev. B 78, 035413 (2008). [17] S Dusari, J Barzola-Quiquia, P Esquinazi, N Garćıa, Ballistic transport at room tempera- ture in micrometer-size graphite flakes, Phys. Rev. B 83, 125402 (2011). [18] P Esquinazi, J Barzola-Quiquia, S Dusari, N Garćıa, Length dependence of the resistance in graphite: Influence of ballistic transport, J. Appl. Phys. 111, 033709 (2012). [19] Q Lin, T Li, Z Liu, Y Song, L He, Z Hu, Q Guo, H Ye, High-resolution TEM obser- vations of isolated rhombohedral crystallites in graphite blocks, Carbon 50, 2369 (2012). [20] C H Lui, Z Li, Z Chen, P V Klimov, L E Brus, T F Heinz, Imaging stacking order in few-layer graphene, Nano Lett. 11, 164 (2011). [21] N B Kopnin, M Ijäs, A Harju, T T Heikkilä, High-temperature surface superconductivity in rhombohedral graphite, Phys. Rev. B 87, 140503 (2013). [22] P Esquinazi, J Barzola-Quiquia, D Spemann, M Rothermel, H Ohldag, N Garćıa, A Setzer, T Butz, Magnetic order in graphite: Exper- imental evidence, intrinsic and extrinsic dif- ficulties, J. Magn. Magn. Mat. 322, 1156 (2010). [23] H Kempa, Y Kopelevich, F Mrowka, A Setzer, J H S Torres, R Höhne, P Esquinazi, Mag- netic field driven superconductor-insulator- type transition in graphite, Solid State Com- mun. 115, 539 (2000). [24] Y Kopelevich, P Esquinazi, J H S Torres, R R da Silva, H Kempa, Graphite as a highly cor- related electron liquid, In: Advances in Solid State Physics, Vol. 43, Ed. B Kramer, Pag. 207, Springer-Verlag, Berlin (2003). [25] T Tokumoto, E Jobiliong, E Choi, Y Oshima, J Brooks, Electric and thermoelectric transport probes of metal-insulator and two-band mag- netotransport behavior in graphite, Solid State Commun. 129, 599 (2004). [26] X Du, S W Tsai, D L Maslov, A F Hebard, Metal-insulator-like behavior in semimetallic bismuth and graphite, Phys. Rev. Lett. 94, 166601 (2005). [27] Y Zhang, J P Small, W V Pontius, P Kim, Fabrication and electric-field-dependent trans- port measurements of mesoscopic graphite de- vices, Appl. Phys. Lett. 86, 073104 (2005). [28] See several reviews in: Graphite and precur- sors, World of carbon series, Vol. 1, Ed. P Del- haes, Gordon and Breach Science Publishers (2001). [29] T Scheike, P Esquinazi, A Setzer, W Böhlmann, Granular superconductiv- ity at room temperature in bulk highly oriented pyrolytic graphite samples, Carbon 59, 140 (2013). [30] Y Kopelevich, V Lemanov, S Moehlecke, J Torres, Landau level quantization and possi- ble superconducting instabilities in highly ori- ented pyrolitic graphite, Phys. Solid State 41, 1959 (1999). [31] H Kempa, H C Semmelhack, P Esquinazi, Y Kopelevich, Absence of metal-insulator transition and coherent interlayer transport in oriented graphite in parallel magnetic fields, Solid State Commun. 125, 1 (2003). [32] Y Kopelevich, P Esquinazi, J Torres, R da Silva, H Kempa, F Mrowka, R Ocaña, Metal-insulator-metal transitions, supercon- ductivity and magnetism in graphite, In: Stud- ies of High Temperature Superconductors, Vol. 45, Chap. 3, Pag. 59, Nova Science Publishers Inc. (2003). 050007-16 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi [33] M P A Fisher, Quantum phase transitions in disordered two-dimensional superconductors, Phys. Rev. Lett. 65, 923 (2000). [34] E Abrahams, S V Kravchenko, M P Sarachik, Metallic behavior and related phenomena in two dimensions, Rev. Mod. Phys. 73, 251 (2001). [35] Y Shapira, G Deutscher, Semiconductor- superconductor transition in granular Al-Ge, Phys. Rev. B 27, 4463 (1983). [36] J S Langer, V Ambegaokar, Intrinsic resistive transition in narrow superconducting channels, Phys. Rev. 164, 498 (1967). [37] D E McCumber, B I Halperin, Time scale of intrinsic resistive fluctuations in thin super- conducting wires, Phys. Rev. B 1, 1054 (1970). [38] P Esquinazi, N Garćıa, J Barzola-Quiquia, P Rödiger, K Schindler, J L Yao, M Ziese, Indications for intrinsic superconductivity in highly oriented pyrolytic graphite, Phys. Rev. B 78, 134516 (2008). [39] L Ji, M S Rzchowski, N Anand, M Thinkam, Magnetic-field-dependent surface resistance and two-level critical-state model for granular superconductors, Phys. Rev. B 47, 470 (1993). [40] Y Kopelevich, C dos Santos, S Moehlecke, A Machado, Current-induced superconductor- insulator transition in granular high-Tc super- conductors, arXiv:0108311 (2001). [41] I Felner, E Galstyan, B Lorenz, D Cao, Y S Wang, Y Y Xue, C W Chu, Magnetoresistance hysteresis and critical current density in gran- ular RuSr2Gd2−xCexCu2O10−δ, Phys. Rev. B 67, 134506 (2003). [42] S Dusari, J Barzola-Quiquia, P Esquinazi, Superconducting behavior of interfaces in graphite: Transport measurements of micro- constrictions, J. Supercond. Nov. Magn. 24, 401 (2011). [43] J Barzola-Quiquia, S Dusari, G Bridoux, F Bern, A Molle, P Esquinazi, The influence of Ga+ irradiation on the transport properties of mesoscopic conducting thin films, Nanotech. 21, 145306 (2010). [44] V Ambegaokar, B I Halperin, Voltage due to thermal noise in the DC Josephson effect, Phys. Rev. Lett. 22, 1364 (1969). [45] A Ballestar, J Barzola-Quiquia, T Scheike, P Esquinazi, Evidence of Josephson-coupled superconducting regions at the interfaces of highly oriented pyrolytic graphite, New J. Phys. 15, 023024 (2013). [46] J Barzola-Quiquia, P Esquinazi, Ferromagnetic- and superconducting-like behavior of the electrical resistance of an inhomogeneous graphite flake, J. Supercond. Nov. Magn. 23, 451 (2010). [47] A Ballestar, P Esquinazi, Highly oriented py- rolytic graphite TEM lamellae preparation to study transport properties of the internal in- terfaces, J. Visual. Exp. (in press). [48] Y Kopelevich, J H S Torres, R R da Silva, F Mrowka, H Kempa, P Esquinazi, Reentrant metallic behavior of graphite in the quantum limit, Phys. Rev. Lett. 90, 156402 (2003). [49] H J Gardner, A Kumar, L Yu, P Xiong, M P Warusawithana, L Wang, O Vafek, D G Schlom, Enhancement of superconductivity by a parallel magnetic field in two-dimensional superconductors, Nat. Phys. 7, 895 (2011). [50] J González, F Guinea, M A H Vozmedi- ano, Electron-electron interactions in graphene sheets, Phys. Rev. B 63, 134421 (2001). [51] K Scharnberg, R A Klemm, p-wave supercon- ductors in magnetic fields, Phys. Rev. B 22, 5233 (1980). [52] A Knigavko, B Rosenstein, Spontaneous vor- tex state and ferromagnetic behavior of type- II p-wave superconductors, Phys. Rev. B 58, 9354 (1998). [53] N Casan-Pastor, P Gomez-Romero, L C Baker, Magnetic measurements with a squid magnetometer: Possible artifacts induced by sample holder off centering, J. Appl. Phys. 69, 5088 (1991). [54] A Ney, T Kammermeier, V Ney, K Ollefs, S Ye, Limitations of measuring small magnetic 050007-17 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi signals of samples deposited on a diamagnetic substrate, J. Magn. Magn. Mater. 320, 3341 (2008). [55] M Sawicki, W Stefanowicz, A Ney, Sensi- tive SQUID magnetometry for studying nano- magnetism, Semicond. Sci. Tech. 26, 064006 (2011). [56] J Barzola-Quiquia, W Böhlmann, P Es- quinazi, A Schadewitz, A Ballestar, S Dusari, L Schultze-Nobre, B Kersting, Enhancement of the ferromagnetic order of graphite after sul- phuric acid treatment, Appl. Phys. Lett. 98, 192511 (2011). [57] M W McElfresh, Y Yeshurun, A P Maloze- moff, F Holtzberg, Remanent magnetization, lower critical fields and surface barriers in an YBa2Cu3O7 crystal, Physica A 168, 308 (1990). [58] T Scheike, W Böhlmann, P Esquinazi, J Barzola-Quiquia, A Ballestar, A Setzer, Can doping graphite trigger room temperature su- perconductivity? Evidence for granular high- temperature superconductivity in water-treated graphite powder, Adv. Mater. 24, 5826 (2012). [59] Y Kopelevich, P Esquinazi, J Tor- res, S Moehlecke, Ferromagnetic- and superconducting-like behavior of graphite, J. Low Temp. Phys. 119, 691 (2000). [60] S Senoussi, C Aguillon, S Hadjoudj, The con- tribution of the intergrain currents to the low field hysteresis cycle of granular superconduc- tors and the connection with the micro- and macrostructures, Physica C 175, 215 (1991). [61] M Borik, M Chernikov, V Veselago, V Stepankin, Anomalies of the magnetic properties of granular oxide superconductor BaPbl−xBixO3, J. Low Temp. Phys. 85, 283 (1991). [62] B Andrzejewski, E Guilmeau, C Simon, Mod- elling of the magnetic behaviour of random granular superconductors by the single junc- tion model, Supercond. Sci. Tech. 14, 904 (2001). [63] R Prozorov, Y Yeshurun, T Prozorov, A Gedanken, Magnetic irreversibility and relaxation in assembly of ferromagnetic nanoparticles, Phys. Rev. B 59, 6956 (1999). [64] N B Hannay, T H Geballe, B T Matthias, K Andres, P Schmidt, D MacNair, Supercon- ductivity in graphitic compounds, Phys. Rev. Lett. 14, 225 (1965). [65] T E Weller, M Ellerby, S S Siddharth, R P Smith, T Skippe, Superconductivity in the in- tercaled graphite compounds C6Yb and C6Ca, Nat. Phys. 1, 39 (2005). [66] N Emery, C Hérold, M D’Astuto, V Garcia, C Bellin, J F Marêché, P Lagrange, G Loupias, Superconductivity of bulk CaC6, Phys. Rev. Lett. 95, 035413 (2005). [67] R R da Silva, J H S Torres, Y Kopelevich, Indi- cation of superconductivity at 35 K in graphite- sulfur composites, Phys. Rev. Lett. 87, 147001 (2001). [68] Y Kopelevich, R R da Silva, J H S Torres, S Moehlecke, M B Maple, High-temperature local superconductivity in graphite-sulfur com- posites, Physica C 408, 77 (2004). [69] I Felner, Y Kopelevich, Magnetization mea- surement of a possible high-temperature super- conducting state in amorphous carbon doped with sulfur, Phys. Rev. B 79, 233409 (2009). [70] Y Kopelevich, P Esquinazi, Ferromagnetism and superconductivity in carbon-based systems, J. Low Temp. Phys. 146, 629 (2007). [71] I Felner, O Wolf, O Millo, High-temperature superconductivity in sulfur-doped amorphous carbon systems, J. Supercond. Nov. Magn. 25, 7 (2012). [72] M Kociak, A Y Kasumov, S Guéron, B Reulet, I I Khodos, Y B Gorbatov, V T Volkov, L Vac- carini, H Bouchiat, Superconductivity in ropes of single-walled carbon nanotubes, Phys. Rev. Lett. 86, 2416 (2001). [73] I Takesue, J Haruyama, N Kobayashi, S Chi- ashi, S Maruyama, T Sugai, H Shinohara, Su- perconductivity in entirely end-bonded multi- walled carbon nanotubes, Phys. Rev. Lett. 96, 057001 (2006). 050007-18 Papers in Physics, vol. 5, art. 050007 (2013) / P. Esquinazi [74] Z K Tang, L Zhang, N Wang, X X Zhang, G H Wen, G D Li, J N Wang, C T Chan, P Sheng, Superconductivity in 4 angstrom single-walled carbon nanotubes, Science 292, 2462 (2001). [75] G M Zhao, Y S Wang, Possible superconduc- tivity above 400 k in carbon-based multiwall nanotubes, arXiv:0111268 (2001). [76] E A Ekimov, V A Sidorov, E D Bauer, N N Mel’nik, N J Curro, J D Thompson, S M Stishov, Superconductivity in diamond, Na- ture 428, 542 (2004). [77] Y Takano, M Nagao, I Sakaguchi, M Tachiki, T Hatano, K Kobayashi, H Umezawa, H Kawarada, Superconductivity in diamond thin films well above liquid helium temperature, Appl. Phys. Lett. 85, 2851 (2004). [78] K Antonowicz, Possible superconductivity at room temperature, Nature 247, 358 (1974). [79] K Antonowicz, The effect of microwaves on dc current in an Al-Carbon-Al sandwich, Physica Status Solidi A 28, 497 (1975). [80] N Agrait, J Rodrigo, S Vieira, On the tran- sition from tunneling regime to point-contact: graphite, Ultramicroscopy 42–44, Part 1, 177 (1992). [81] R Nandkishore, L S Levitov, A V Chubukov, Chiral superconductivity from repulsive inter- actions in doped graphene., Nat. Phys. 8, 158 (2012). [82] A M Black-Schaffer, S Doniach, Resonating valence bonds and mean-field d-wave supercon- ductivity in graphite, Phys. Rev. B 75, 134512 (2007). [83] N B Kopnin, T T Heikkilä, G E Volovik, High- temperature surface superconductivity in topo- logical flat-band systems, Phys. Rev. B 83, 220503 (2011). [84] N Garćıa, P Esquinazi, Mean field supercon- ductivity approach in two dimensions, J. Su- percond. Nov. Magn. 22, 439 (2009). [85] G Profeta, M Calandra, F Mauri, Phonon- mediated superconductivity in graphene by lithium deposition, Nat. Phys. 8, 131 (2012). [86] B Uchoa, A H C Neto, Superconducting states of pure and doped graphene, Phys. Rev. Lett. 98, 146801 (2007). [87] N B Kopnin, E B Sonin, BCS superconductiv- ity of Dirac electrons in graphene layers, Phys. Rev. Lett. 100, 246808 (2008). [88] N Reyren, S Thiel, A D Caviglia, L F Kourk- outis, G Hammerl, C Richter, C W Schneider, T Kopp, A S Rüetschia, D Jaccard, M Gabay, D A Muller, J M Triscone, J Mannhart, Su- perconducting interfaces between insulating ox- ides, Science 317, 1196 (2007). [89] A Gozar, G Logvenov, L F Kourkoutis, A T Bollinger, L A Giannuzzi, L A Muller, I Bo- zovic, High-temperature interface supercon- ductivity between metallic and insulating cop- per oxides, Nature 455, 782 (2008). [90] F Muntyanua, A Gilewski, K Nenkov, J Warchulska, A Zaleski, Experimental mag- netization evidence for two superconducting phases in Bi bicrystals with large crystallite disorientation angle, Phys. Rev. B 73, 132507 (2006). [91] F Muntyanua, A Gilewski, K Nenkov, A Za- leski, V Chistol, Superconducting crystallite interfaces with Tc up to 21 K in Bi and Bi- Sb bicrystals of inclination type, Solid State Commun. 147, 183 (2008). [92] Y Kopelevich, P Esquinazi, Graphene physics in graphite, Adv. Mater. (Weinheim, Ger.) 19, 4559 (2007). [93] W A Muñoz, L Covaci, F Peeters, Tight- binding description of intrinsic superconduct- ing correlations in multilayer graphene, Phys. Rev. B 87, 134509 (2013). [94] Y Kawashima, Possible room temperature superconductivity in conductors obtained by bringing alkanes into contact with a graphite surface, AIP Advances 3, 052132 (2013). 050007-19