Pergamon
Solid State Communications, Vol. 92, Nos i-2, pp. 89-100, 1994 Elsevier Science Ltd Printed in Great Britain 0038-1098(94)00492-7 0038-1098/94 $7.00+.00 ORGANIC SUPERCONDUCTORS DJ~rome
Laboratoire de Physique des Solides (associ6 au CNRS) Universit6 Paris-Sud, 91405 Orsay, France
The discovery of organic superconductors in 1980 was an happy conclusion of a search for high Tc superconductors (although Tc was modestly in the Kelvin range). The first generation of organic superconductors belonged to a large family of quasi-one-dimensional conducting cation radical salts. Most of them exhibit a wide range of new phenomena, including the competition between ground states, the influence of a magnetic field on the non ordered phase of a quasi-one-dimensionalelectron gas and the onset of spin density wave phases at high magnetic field with quantized Hall effect. The extensive study of the (TM)2X series has shown that electron interactions are repulsive and of the order of the electron bandwidth. However, the nature of the ground states relies essentially on the balance between charge localisation boosted by the Umklapp scattering and the interchain hopping inte~al. Second and third generation organic superconductors axe two and three dimensional molecular conductors respectively with maximum Tc of 12K (ET)2X and 33 K (fullerides). Keywords:A.organic crystals,superconductors
I.
Introduction
fulvalene molecule [21.These molecules tend to form stacks in the crystalline state. Such a packing optimizes the overlap between molecular orbitals along the stacks. When the Coulomb repulsion (U) does not overcome the energy gained by the band formation, conducting properties can be observed. However, due to the strong directionality of the intermolecular couplings a pronounced one dimensional character is obtained in most organic conductors (they are called one dimensional,(1-D), or quasi one dimensional, (Ql-D), conductors when the small interchain coupling is taken into account.
Crystalline conductors of molecular origin are basically different from the regular metallic crystals. They are solids made of building blocks (the molecules) reflecting on the one hand the properties of the molecular species (ionization energy or electron affinity, characteristic infrared or Raman active vibration modes, NMR shifts related to individual atomic sites,...) and on the other hand the physical properties of a metal (dp/dT > 0, small Hall coefficient, plasma frequency in the optical reflectance,...).
Salts of the tetramethyltetrathia(selena)fulvalene, i.e. TMTT(S)F molecule [31 are prototype materials for molecular conductors [41 figure(l). They reveal some of the most typical behavior for 1-D physics which will be reviewed in this article 15.61.
Contrasting with conventional molecular crystals made of neutral organic molecules held together by weak Van der Waals forces, organic conductors contain molecules with unpaired carriers in x-molecular orbitals presenting an open shell configuration. Such a situation comes generally from a partial oxidation (reduction) of donor (acceptor) molecule in the formation of a salt with an inorganic anion (cation). In addition, a strong intermolecular overlap of ~orbitals makes the electron delocalization possible over all molecular sites in the crystals (although in some cases the delocalization occurs preferentially along selected cristallographic directions.
These molecules also give rise to higher dimensionality conductors. In that case, the planar molecules form dimers with a strong overlap between heteroatoms both within each dimer and between neighbouring dimers. The dimers are packed in planes according to a chess-board pattern. This crystal structure gives rise to a large metal-like conductivity within 2-D conducting molecular layers and about 104 less conduction between the sheets 171.
Molecular conductors which are the subject of this article should not be confused with extended conjugated polymers and even graphite, where the x-electron system of the extended molecules provides the conducting path-way. The development of molecular conductors as well as the subsequent development of conducting polymers has been stimulated to a large extent by the suggestion made by Little in 1964 that the arrangement of chain conductors in a polarizable medium could promote superconductivity at high temperature [11.
Tile recently discovered fullerene conductors bear much similarity with organic conducting salts. Neutral C60 molecule possesses a closed shell electronic structure giving rise to an insulating solid. But, very much like organic molecules, when carriers are added to the neutral molecule through its reduction by an alkali cation the salt thus formed may become conducting 181.
Molecular conductors which are surveyed in this article usually derive from molecules which are small modifications of the flat tetrathiafulvalene (TI'F) prototype
Similarly to 1 and 2-D molecular conductors the conductivity of fullerenes is due to the intermolecular overlap of x-molecular orbitais. For all these molecular 89
90
ORGANIC SUPERCONDUCTORS
Cl~j~Se
Se_/CHa
c, J's. TMTSF
%
Vol.92, Nos 1-2
Table I. Parameters of the calculated dispersion relations for (TM'I'rF)2Br and (TMTSF)2PF6, figure (2) for the labels) using their room temperature and ambient pressure structures (private communication from E.Canadell and to be published in J.Physique I). The calculations were carried out using a basis set of double-~ Slater type orbitals.Numbers in brackets are the band parameters obtained from the two dimensional multi-~ calculation in ref 11olAll values are in eV.
(TM'I~F)2Br
(TMTSF)2PF6
......................................................................
WI Wn Wa Wb Wb' Wc Wc'
1) TMTSF donor molecule and a side view of the prototype (TMTSF)2X structure.
0.316 (0.372) 0.417 0.797 (0.859) 0.052 (0.044) 0.257 0.002 0.004
0.518(0.595) 0.573 1.218(1.366) 0.140 (0.114) 0.248 0.003 0.002
We shall see that the half filled band character together with the amplitude of the interchain overlap are both important parameters which govern all physical properties in the (TM)2X family.
conductors counter ions do not play any significant role at first sight in the intermolecular conduction process. This article intends to present a very short review of some physical aspects of organic conductors in which superconductivity has been stabilized 19]. The low dimensionality of these Q-1-D conductors makes the physical properties of these materials extremely diversified and fascinating. Therefore, besides superconductivity, the role of dimensionality,correlations and band filling will also be discussed.
2.
Materials and Electronic Structure
The salts of the (TM)2X family where TM is the molecule TM'Iq'F, TMTSF or any mixed (S-Se) molecule represent the archetype of I-D conductors, figure (1). Superconductivity, when it exists in these materials is restricted to the low temperature domain ( T < 2K) but a wealth of other instabilities are also observed 161. All isostructural (TM)2X compounds exhibit a rather simple band structure with a conduction band built from the highest occupied molecular orbital (HOMO) of the isolated molecule. The conduction band has a width ranging between 0.4 and 1.2 eV [lOl figure (2) and its filling is formally 1/4 since one carrier (hole) is provided to each one dimensional unit cell containing two TM entities by the ionization of the monoanton X. However, the three dimensional cristallographic structure of the anions induces a slight dimerization of the bond lengths (and of the overlaps) between molecules along the stacking direction giving rise to a splitting of the conduction band into an half-filled upper band and a filled lower band. The half-filled band character is related to the amount of band dimerization [111; it is generally more pronounced for sulfur than for selenium compounds and it is known to be reduced under pressure, figure (3). Some band parameters of these salts are listed in Table I. Typical band widths are in'the range of 0.4 and 1.2 eV for sulfur and selenium salts respectively.
3. I n s t a b i l i t i e s in Theoretical Features.
the
I-D
Electron
Gas,
To create a superconducting phase that is stable below a certain critical temperature, one must overcome the su-ong tendency shown by the low dimensional electron gas (in particular those with a pronounced 1-D character) to undergo a transition toward an insulating ground state exhibiting a density wave periodic potential with wave vector Q The tendency towards instability in a I-D electron gas is already visible without interactions as shown by the temperature dependence of the non-interacting electron susceptibility I131: tl~- tO" , ) - , n ( k ) - n ( k + q ) Z,, ~q,
(1)
) = "7" e ( k ) - e(k + q)
where e(~.) and n(k) are the excitation energy counted
-i0.0
t~
-I0.2
c LIJ
-10.4-
/
f ,t
\
-i0.6
X
F
¥
2) Energy dispersion of (TMTSF)2PFb and Fermi surface.
interaction is usually approximated by two coupling constants gl and g2 which are respectively the backward (q = 2kF) and forward (q = 0) scattering amplitudes It61.
X
:.':h ×
In addition, there are "also g3 terms involving the transfer of two particles from one side of the Fermi surface to the other. Since the total momentum transfer is 4kF, this process is allowed only if 4kF is a reciprocal lattice vector (Umklapp scattering) 1171.
v
-m
0
-~
~ p='/2,2k1=: i ~, i~: ]J~ k za[mod2~a
3) Illustration (schematic) of the structural dimerization giving rise to an upper half-filled band in (TM)2X compounds.
from the Fermi level and the Fermi occupation factor respectively. At low temperature eq.(l) leads to logarithmic divergences, Z~o°(2kr,T)
=
91
ORGANIC SUPERCONDUCTORS
Vol.92, Nos 1-2
Iog(E--T~)
and
Z'on ( q - 2k~. , T = 0,o~ = 0 ) - log(,,,
(q - 2k F ) EF
where ke is the wave vector which nests the Fermi surface of a 1-D electron gas. In eq. (1) the electron dispersion has been linearized in the vicinity of the Fermi energy, namely
The mean-field treatment of the 1-D electron gas amounts to neglecting the mixture of Peierls and Cooper channels and evaluating correlation functions via the perturbation expansion of ladder diagrams. This summation procedure is not appropriate for the I-D electron gas as it leads to finite transition temperatures (against the Landau criterion). However, when superconducging and electron-hole instabilities are treated on equal footing the famous (gl, g2) phase diagram for the most stable ground state of a 1D electron gas reveals a common border (gl = 2g2) between Peierls and Cooper instabilities (g3 = 0 so far). For gl < 0, the only stable modulated state in the one exhibiting a charge density wave oscillation (CDW) whereas there is a gap in the spin excitations. On the other hand, the spin modulated ground state is stable only for gl > 0. There are low-lying spin excitations and the spin susceptibility is finite at T = 0.
(2)
Response functions of all instabilities contain at all orders in perturbation logarithnaic terms such as log (OWEF), as well as a mixture of Peierls and Cooper divergences.
Because the FS of TM2X compounds is open, the major concepts of 1-D physics should apply to these materials, namely, (i) the absence of long range order at finite temperature 1131 (ii) the existence of equally divergent correlations between electron states at (+k,l")and (-k,,l,) 1141 leading to superconducting pairing (Cooper channel) and electron and hole states at k and k + O leading to density wave instabilities (Peierls channel) 1141.These two pairing channels coexist in the I-D electron gas but the divergence of the electron-hole (density wave) channel is, however, limited at low temperature by the existence of a finite coupling tz between near neighbouring chains, which gives rise to a cross over temperature below which 1-D physics no longer applies [ 151 namely
The only relevant treatment of the 1-D electron gas is the one which retains the mixture of channels as the temperature is lowered. Such a situation is achieved when the band width EF is renormalized towards a smaller value at low temperature so that the various log (T/EF) terms of the perturbative expansion become less divergent. The renormalization of the bandwidth is also accompanied by a renormalization of the various coupling constants gi 1161. Spin excitations (the static su~eptibility) rely only on the backward scattering gl whereas the charge excitations depend on the set of parameters g3, gl-2g2. Renormalization procedures lead to hyperbolic trajectories 1171 as shown in fig. (4).
f(k) = fo + V~.(Ik,l-k~-~
T~ - t>/~
g2 _ (gl - 2g2)2 = const (3)
and the I-D bare susceptibility saturates at low temperature.
If g 3 l < g l - 2 g 2, then g3 scales to zero and the system evolves towards an attractive pairing instability as T ~ 0.
Z:,° (T < T ,..) - I o g ( ~ E F )
(4)
On the other hand, for g3 > gl - 2g2 i.e. strong Umklapp situation for an half-filled band, g3 scales towards strong
The divergence of the Cooper channel is not affected by the interchain coupling since the equivalence between + k a n d - k states is preserved even in the presence of a finite transverse coupling. Going to the situation of the interacting 1-D electron gas, and focussing on low energy excitations, the only pertinent wave vectors to be taken into account for the Coulomb interaction are those which are close to q = 0 and q = 2 k F since the electron-hole pair excitation energy is zero for these vectors. Consequently, the Coulomb
coupling with the concomitant opening of a gap Ap in the charge excitation spectrum below the temperature Tp. However, gl(T) renormalizes to zero when gl > 0, following the equation, g| gl(T) =
gl 1- ~
T
(5)
In
F EF The spin susceptibility is not affected by the charge
92
ORGANIC SUPERCONDUCTORS
localization below Tp. For the strong or weak g3 situation )~s(T) then reads:
Vol.92, Nos 1-2
and the single particle density of states vanishes at EF according to a power law singularity N(to)-Jto[ a
(6)
Zs(T) = ZS(0)(I + gl (T))
where XS(0) is the zero temperature susceptibility enhanced above the bare band value by a Stoner factor. This phenomenon is known as the spin-charge separation of the 1-D electron gas. Experimentally, one should not be able to tell from the only knowledge of the temperature dependent susceptibility whether the chain is a metallic conductor or a linear chain of localized spins. It can be shown that the long distance decay of the correlations functions of the interacting 1-D electron gas is determined by a unique exponent which also governs their temperature dependence HsI namely
ZDw(2kF,T)aT(Kp-1)
and
ZscaT (1/ Kp-I) (7)
Susceptibilities follow power laws dependences at low temperature. Unlike the mean field treatment, no phase transition is observed at finite temperatm'e. Density waves or pairing correlations are divergent at low temperature for Kp < 1 and > l, respectively. The exponent Kp depends on the bare values of the coupling constants, namely, Kp =
1+g'-2g~-
Therefore, Kp = l con'esponds to the non2rtvF interacting situation and Kp = 0 represents the strong coupling limit with a large gap in the charge degrees of freedom. In addition, the momentum distribution no longer exhibits a step at the Fermi wave vector; this a a Luttinger liquid.
nk = nkF -flsign(k - kF)k - kF a
A(q,~)
(8)
What is expected from an experimental study of I-D conductors is the search for the most characteristic features of 1-D physics, namely, the existence of a I-D regime in temperature, the spin-charge separation, evidences for a Luttinger liquid and the determination of power laws for divergent susceptibilities. 4.
I-D Physics : Experimental
A direct consequence of the pronounced I-D character of the TM2X series is the existence of a wide variety of behaviours that can be observed for the transport properties, depending on parameters such as the chemical composition of the organic molecule or the inorganic anion and the hydrostatic pressure in various members of the family. The incipient instabilities that develop at high temperature in the I-D regime depend on the intra-chain interactions and on the strength of the Umklapp scattering (g3) whereas the nature of the long range that is stabilized at low temperature depends also on the coupling between the chains. The generic diagram in figure (5) displays the v~uiety of 1-D regimes and ground states that can be observed in the TM2X series [61. Special attention has been paid to some key compounds labelled by letters in figure (5). At the left side of the diagram, (TMTI'F)2PF6 (a) presents a poor conductivity o'. = 30f)-~cm-~ but a metal like conductivity at high temperature and a charge localized (CL) behaviour below Tp = 250K. It is the compound in which the halffilling character is most pronounced within the TM2X series since the overlap alternation amounts to A t . / t . = 0.38. The localization which takes place below To is due to repulsive interactions between 1-D electrons belonging to the half-filled band [191 figure (6). The average distance between electrons is a (the periodicity along the stacking direction). In the same temperature regime the localized electrons of the Heisenberg chain couple to phonons and give rise to 1-D lattice fluctuations at wave vector 2kF (periodicity 2a) detected by X-ray experiments.
~100 (b)
(a)
where c t = l ( K + - L - 2 ) 4 ¢ K/,
CL
v
#
-..
CONDUCTOR
IQ, i
I
Kp
4) Diagramatic representation of the pair susceptibilities. Cooper A(tI, to) and Peierls x (2kF + q, to) at lowest order (a) and second order (b) showing the mixture of channels. Renormalization paths for g3 and the charge degree of freedom exponent Kp.
t o
t~ bc
PRESSURE
t d
t 5 kbar e
5) Generalized phase diagram for the (TM)2X series. Spin-Peierls (SP), spin density wave (SDW) and superconductivity (SC) are indicated together with the zero pressure location of some prototypical compounds, see text.
93
ORGANIC SUPERCONDUCTORS
Vol.92, Nos 1-2
Below about 20 K these 1-D lattice fluctuations order three-dimensionally and the system undergoes a phase transition towards a spin-Peierls phase with the loss of the spin degrees of freedom. NMR experiments have shown an activated behaviour for the spin-lattice relaxation time related to the spin-Peierls gap in the spin degrees of freedom 12Ol. As one moves to the right in the diagram not much is noticed for the temperature dependent susceptibility but the g3 term (proportional to the band dimerization, g3=gl(Ao/EF) becomes less important and the charge localization is sllifted to lower temperatures, To = 100K for (TMTTF)2Br, (b) 1211. Coupling of weakly localized electrons to phonons looses in strength, although some 2kF I-D lattice fluctuations are still observed below Tp 1221. The weakly localized electron gas undergoes a transition towards a ground state which is characterized by a modulation of the spin density at wave vector 2kF. NMR experiments have demonstrated both the magnetic nature of the ground state in (TMTTF)2Br and its commensurate character 123.241.The localized moments order antiferrromagnetically along transverse directions and the periodicity is doubled accordingly. The commensurate character of the SDW ground state is also attested by the activation of the relaxation rate since no low lying excitation modes (phasons) are present in this x 2 connmensurate ground state to couple to nuclear spins. The behaviour of the Bechgaard salt (TMTSF)2PF6 (c) 131is interesting to compare with that of (TMTIT)2Br as its conductivity is metal-like down to a phase transition towards a SDW ground state, figure (7). No 1-D localization is observed before the transition at TSDW (= 12 K). The additional periodicity in the exchange potential in the SDW phase is responsible for the opening of a gap at the Fermi level accompanied by a metal-insulator transition 1251. No doubts the concept of Fermi surface is valid for the compound (TMTSF)2PF6. Furthermore, it is the optimum nesting condition of the Q-I-D Fermi surface which decides for the actual value of the magnetic modulation wave vector Q. 13c and IH-NMR spectra have unambiguously shown that O is incommensurate with the underlying lattice 12t~.27.2~1. Furthermore, the existence of low lying phason
modes (although pinned by impurities and defects) makes the NMR relaxation rate temperature independent in the SDW state. The main difference between (TMTSF)2PF6 and (TM'I~F)2Br is the weak localization which takes place at high temperature for the latter compound. The interchain interaction which is required for the onset of long range order in ( T M T T F ) 2 B r can be ascribed to an antiferromagnetic coupling interaction between weakly localized spins provided by the exchange of a 2kF spin fluctuation between chains 1291.The existence of such an interchain coupling suggests that the correlated jump of an electron-hole pair (a 2kF spin fluctuation) between neighbouring chains is easier than the jump of an individual particle : the electron-hole cross-over temperature Tx2 is thus higher than the single particle cross-over at Txi. Actually, it has been suggested that the single particle tunnelling can be impeded by I-D correlations. The crossover temperature thus becomes instead of eq.(3) [291: 1-Kp
T ~ T:(~,)
%
(9)
The reduction of Txl in eq.(9) can be quite significant even for relatively small repulsive interactions ( Kp <; 1) since t.I./EF = 0.1 and may lead to a suppression of the single particle cross-over in agreement with the NMR determination for a lower limit value in the vicinity of 10 K I301. Direct evidences for the existence of a 1-D regime and the role of electron-electron interactions in the (TM)2X series are given by the extensive NMR studies which have been conducted in various compounds. NMR spin lattice relaxation probes the low energy excitation modes of the I-D electron gas at q = 0 and q = 2kF. Therefore, instead of a summation over all q-vectors in the Moriya's formulation 1321for T 1 the relaxation rate in a I-D electron gas simplifies 1331
l T, TI-' = C,,x~ (T)
t.2
(lO)
C,
+
(TMTTF)~Br
X, ( a . u . ) 1.0
I0 I'
2.5
0.8 tee
•
o~
•
0.8 0.4.
10 "d
10 ~ o
50
I! :
100_ 150 ' Teaml~nnztax~
ee °
° e °
~U
',,pit. Br lO0 ' (~O
eg o
....... 200
m X
1.5
30O 250
,.-, ,"
."
3O0
1°o 6) Temperature dependence of the ESR susceptibility of
( T M T T F ) 2 B r and p(T) (TMTFF)2Br.
,
E
\,
10 "~
0.2
0.0
2.0
i
1
"
• • • •
•r (K)
|0 #
~
(TMTSF)zPF,
'
16o 26o TEMPERATURE (K) '
'
300
for (TMTTF)2PF6 and 7) Resistivity and susceptibility of (TMTSF)2PF6.
94
Vol.92, Nos 1-2
ORGANIC SUPERCONDUCTORS
The first term it) eq.(9) is linked to the uniform (q = 0) spin excitations of the I-D electron gas (Co is a constant related to the hyperfine interaction), whereas the second term probes the q = 2kF spin fluctuations. CI (T) is related to the 2kF spin susceptibility, namely CI (T) ~ XSDW (2kF, T) Qt T ~ ' Hence, the contribution to T I ' I coming from 2kF spin fluctuations becomes (-~1)2,F ,,T ~
(1 l)
Figure (8) shows that the temperature independent contribution to l/T1 (extreme situation of strong coupling with Kp= 0 for a 1-D quantum antiferromagnet ) prevails in the whole temperature domain below T o for (TM'i~F)2PF6 and (TM'I'TF)2Br salts. The temperature dependence of XSDW (2kF, T) has been approached with much accuracy in the conductor (TMDTDSF)2PF6 (c in fig. 5) built with a mixed sulfurselenium TM molecule I341 figure (9). There, the strong coupling limit, Kp ---->0,is reached below Tp/2. Some temperature dependence of Kp is expected at higher temperature according to the renormalization trajectories displayed in fig. (4). There is also evidence from NMR data for significant 1-D 2kF spin fluctuations even in compounds exhibiting a metal-like conduction down to the onset of long range ordering at low temperature. For example 2kF spin fluctuations contribute to T~-~ below 100 K and 30 K in (TMTSF)2PF6 and (TMTSF)2CIO4 respectively [311. The data in fig. (8) show that a renormalized "Korringa" behaviour is recovered below 10 K. Therefore, that temperature can be taken as the lower limit value Tx] for the I-D to 3-D single particle cross over regime. The cross-over regime extends up to about 30 K in (TMTSF)2CIO4. Thus the enhancement of T 1-1 observed between 30 and 10 K is due to spin fluctuations which evolve from I-D to 2 or 3-D on cooling towards low temperatures.
Recent photoemission experiments performed on (TMTSF)2PF6 at 50 K have failed to show a Fermi step and a finite single particle density of states at the Fermi level [351. Instead, N(E) is proportional to (E-EF) 1.2 which leads, according to eq.(8), to K 9 = 0.2. Such a small value for Kp represents admittedly a rather strong repulsive coupling situation 1181which would be in qualitative agreement with the low tempet-ature TI "1 data in the same compound but yet hard to reconcile with the metal-like character of the conductivity which shows no sign of localization below 100 K [361. Power law divergences of XSDW have also been measured in the 1-D conductor Tl'F-[Ni(dmit)2] where the band filling is different from 1/2 [371.There again Kp = 0.4 although the TTF chain contributing to the conduction reveals no sign of localization 1381. NMR and susceptibilities have strongly contributed to the understanding of the physical properties of (TM)2X conductors. Thus at ambient pressure the ratio gl/ltVF is found to be close to unity and it is essentially the same for all compounds in the (TM)2X family [311. Furthermore, the amplitude of the AF spin fluctuations and the charge degrees of freedom are governed by the Umklapp coupling term g3 and the combination g 1-2g2. The Umklapp term is dominant for sulfur compounds leading to a localized 1-D antiferromagnet at low temperature but weakens gradually for selenium system (or under pressure). NMR has also provided indications for the dimensionality cross-over regime extanding between 30 and 10 K in selenium compounds. The onset of the SDW ground state of (TMTSF)2PF6 is linked to an incommensurate vector Q nesting optimally the Q-I-D Fermi surface. Inter-chain interactions between first neighbour chains amounts to adding a term such as coskjb in eq.(2). Deviations to Q-1D perfect nesting (occuring under pressure) can be described including also next ,)ear neighbour interactions, namely i391:
E(k)=E,+V/.(l~,l-~.)+2ticosk,b+2t~lcos2kb
(12)
with t t
~ I '''II'
('rllfr~:~l~'! 'IC'
I'TID
' /(
The stability of the SDW ground state is suppressed as soot) as the deviation to perfect nesting becomes
,,/,/,"'I "T
t
I--
i
.--.-
• 2
6
.
I
.
.CtN v
o~
/
2,
.
I
"k:2
x,'T(a.u.)
3
8) 13C spin lattice relaxation rate versus TZs(T) for (TMTTF)2PF6 and (TM'ITF)2Br undergoing a charge localization below Tp. No localization is observed in (TMTSF)2X compounds (dashed line). T~-tin ( T M T S F ) 2 C I O 4 reveals a strong deviation from a "Korringa" behaviour below 30 K 0n.sen).
|
200
100
T/K 9) Power law temperature dependence of XSDW(2kF,T) derived from NMR in (TMDTDSF)2PF6.
Vol.92, Nos 1-2
ORGANIC SUPERCONDUCTORS 1
......
,
. . . . . . . .
,
.
.
.
.
.
.
.
.
"
'
"
95
"
(TMTSFI,CIO~
10
/
,...,'10
~o
"3
hi
"10 i-
10
10
I
1
TEmPERAtURE(K)
1
T~(KI ]
oo
10) Superconducting transition of (TMTSF)2PF6 under pressure and (ET)2Cu(SCN)2 at ambient pressure. Also shown, the specific heat anomaly and the Meissner expulsion in (TMTSF)2CIO4.
o significant, i.e. t±1 > Tsn w where Ts~w is the 1-D mean-field temperature for the onset of the SDW state 14Ol.
Once the SDW ground state is removed increasing the deviation to perfect nesting in eq.(1 1) either under pressure or chemically, superconductivity becomes stable at low temperature with a maximum critical temperature smaller than 2K 19.4tl figure (10). As shown by NMR experiments, 2kF spin fluctuations are still very active when superconductivity sets in. The intra-chain repulsive interactions should prevent s-wave superconducting order. However, the exchange of 2kF spin fluctuations between chains could protnote an attractive pairing for electrons belonging to neighbouring stacks 142]. Such a pairing mechanism would give rise to an anisotropic superconducting gap za(k)= ~o cosk~b with lines of zeros on
results, to be published in Advanced Materials, 1994). It is still premature to state that organic superconductivity of I-D conductors is non-conventional and mediated by a magnetic exchange although several experimental results support this possibility : i) NMR in the superconducting phase revealed a power law temperature dependence of T1-1 and has failed to show a coherence peak at T c 1451. ii) Superconductivity is present in all members of the (TM)2X family but seems to ve very sensitive to small concentration of impurities, cristalline disorder or packing defects 1461. iii) The stability of the superconducting ground state takes advantage of the proximity with the antiferromagnetic ground state as attested by the large negative pressure coefficient 8 1 n T / Sp = -25%kbar -t 1471.
the Fen'hi surface. Figure (10) displays the superconducting transition of (TMTSF)2PF 6 salts observed by resistivity data. What is remarkable in fig.(10) is the sta'ong temperature dependence of p(T) above To, unlike the behaviour of regular metals for which the resistivity is limited at low temperature by temperature independent elastic scattering. The observation of a Meissner flux expulsion 1431 and of a specific heat anomaly 1441support the picture of bulk superconductivity, figure (10). On a C/T vs T 2 plot, the specific heat of (TMTSF)2CIO4 displays a very large anomaly at 1.2K and follows the law: C/T = y + ~T, with -/ = 10.5 mJ. tool "1 K "2 corresponding to N(EF) = 2.1 states.eV -1 tool-1 {N(EF) = 3~'/2/r2k~l and 13 = 11.4 mJ. tool- I K-4. The existence of superconductivity in 'all members of the (TM)2X series has recently been confirmed, (Orsay
A t t e m p t s to force the s t a b i l i z a t i o n of superconductivity at lower pressures hoping to increase Tc in the compound (TMTSF)2ReO4 (where a pressure of 10kbar is necessary to prevent the Red4- anion ordering to open a gap at the Fermi level) have failed to raise Tc above 1.2 K 1481.This behaviour illustrates the close interplay that exists between magnetism and superconductivity in the (TM)2X series.
5.
High Magnetic Fields Properties
Starting fi'om the superconducting state, increasing magnetic field perpendicular to the stacking axis gives rise to spectacular effects. First, the superconducting state is destroyed above a critical field which is very anisotropic (H'~~O.IT,H~-2T) 1491. At higher fields, an antiferomagnetic ground state is recovered above a threshold field ~ 3T for field//c* axis 15ol. Increasing the field further induces a sequence of phase transitions between SDW ground states. The major peculiarity of these phases is
96
ORGANIC SUPERCONDUCTORS -iO0
.....
. . . . , ....
, ....
Vol.92, Nos 1-2
, .... e/
(TNTSF~PF
('rxTsFI~r~
e
0
~
/.
S
-50
k.
.J
xO.
e/
/
:
..//
e[
"~i
/
/
1
!
o
e(
5C
5
10 L5 :~0 ~GNETZC FZELO (T)
fi'F: u ,
25
0
tO 14 ( T e s l s I
20
1 1) Quantized Hall voltage of (TMTSF)2PF6 (under 9 kbar) and phase diagram. Integers show the number of filled Landau levels in each subphase.
provided by the field-independent Hall resistance within each subphase 1511 which is quantized at the value Pxy = h/2Ne 2, where N is a small integer, N = 1, 2,...5.. ( T M T S F ) 2 P F 6 ) gives a remarkable illustration for the behaviour of the quantized Hall resistance 152,531, figure (11). Field-induced spin density wave (FISDW)phases correspond to semimetallic phases with a density of carriers which is about 102 times smaller that that of the Q-I-D metallic phase at H < 3T. Above 18T, the resistance jumps to a very large value and the Hall resistance becomes large and no longer follows the quantized law.The phase above 18T is thus ascribed to a semiconducting SDW phase with no free carriers at zero temperature i.e. N = 0 (unlike subphases with N > 0 which do contain free carriers). Similar restorations of SDW phases under magnetic fields have been achieved with other members of the (TMTSF)2X series, X = CIO4- I5.1.551and ReO4- [561. The C104- compound has raised much interest since the Hall plateaus can be indexed by the quantum numbers N = 1, 3, 6 and so forth. The plateau N = 2 is clearly missing and a sign reversal of the Hall resistance is observed between N = 3 and N = 6 phases. The N = 1 phase of this compound extends from 8 T up to 27 T [571 and the N = (I phase has not yet been clearly identified but is probably shifted above 27 T [58.591 see figure (12). The giant reentrance of TFISDW in the CIO4- phase diagram ,seems to be specific to that compound. No such anomaly for the stability of the FISDW phases has been detected in the prototype system (TMTSF)2PF6. The particular ordering of the CIO4 anions is very likely responsible for this phenomenon. Secondly. as soon as a moderate magnetic field is applied along the c* direction the temperature dependence of the longitudinal resistance no longer remains metal-like. A shallow minimum is observed instead at low temperature and rises continuously with magnetic field 16Ol.Clearly, the magnetic field induces a weak electron localization which bears some similarities with the charge localization observed below Tp in a sulfur compounds such as (TMTTF)2Br at zero field, see fig (5). Furthermore, the localization is most efficient when H//c* and no localization is observed up to 27T as H//b'. Figure (12) displays the phase diagram of (TMTSF)2CIO4 showing both the various ground states corresponding to FISDW phases and the field-induced
localization occuring in the undistorted phase at high temperature, comparing the development of the localization under magnetic field and the diagram of FISDW phases it is clear that a feature such as the reentrance of the FISDW states above 15 T is not observed for the localization. This
va (.V) I
i
I
i
I
(~F):C'104 -20
'(2)
(I)
kj i
20
I
i
I
I
lO '
,
,
,
,
,
,
,
,
,
,
,
,
(13fTSIr)=C10,
20
"
,
,
|If
- st.ti~-"
il
°-1o ,
r,
z-]
~
CL
-
State
I"
o;
,i 2
4
•
•
tO
12
14
IS
18
!0
112
2t
~
~
3O
32
H (TZS~)
12) Phase diagram of (TMTSF)2CIO4 showing the FISDW ground states and the restoration of the localization under magnetic field (H//c*). Quantized Hall voltage of (TMTSF)2CIO4.
97
ORGANIC SUPERCONDUCTORS
Vol.92, Nos 1-2
feature enables us to rule out (at least at high fields) 3-D critical fluctuations in the vicinity of the AF phases as the source of localization. Instead, we have proposed that the localization can be explained by the magnetic field limiting the amplitude of the ta'ansverse quasi-classical motion of Q1-D electrons moving along the stacks. The dimensionality cross-over (or the effective interchain overlap) is thus decreased under magnetic field. Hence, the localization at Tp related to a charge gap in a I-D electron gas is restored by the magnetic field. At zero magnetic field, the cross-over regime (20-30K) is reached at a temperature which is higher than Tp and the charge localization of the I-D conductor looses its significance. FISDW phases have been given a very elegant theoretical treatment [6[I called the "standard model" which also relies on the field induced one dimensionalization of Q1-D electron trajectories under a transverse magnetic field 1621.The divergence of X~o~ (Q)at low temperature which is suppressed by pressure together with the SDW ground state in (TMTSF)2PF6, can be reactivated as soon as a magnetic field is applied to Q-I-D electrons along the direction c* of lowest conductivity 1391. The spin susceptibility X.w,.(q,H) thus displays logarithmically divergent peaks as T approaches zero for the set of different wave vectors defined by qx = 2kF + NeHb/h, where N is an integer. Crudely speaking, the field restores enough 1-D character in the energy dispersion for the SDW distortion to become stable. A sequence of phases must be crossed before the N = 0 phase (which is believed to correspond to the insulating SDW ground state which is stable at anabient pressure) can be reached above 18T. The integer variable N which is related to the deviation of qx fi'om 2kF labels a semimetallic subphase containing N fully occupied Landau levels below the Fermi level. As the degeneracy of each Landau level increases linearly with the magnetic field, the density of carriers in the semimetallic pockets increases accordingly and the Hall resistance remains fixed at the quantized value. The situation where the Fermi level falls between N filled levels and upper empty levels minimizes the diamagnetic energy of the 2D carriers. This situation prevails in a finite range of magnetic fields as long as the wave vector of the magnetic modulation of a given
subphase can vary linearly with the field. First order phase transitions, with a finite jump of the modulation vector, are expected between various subphases, these transitions are in agreement with the hysteresis shown in fig. (12).
6.
Spin Density Wave Collective Conduction
The broken-symmetry ground state characterized by the magnetic modulation (SDW) of wave vector Q is responsible for the activated behaviour of the conductivity in the semiconducting phase of (TMTSF)2X compounds through the single particle gap 2A. Moreover, two kinds of collective excitations arise in the SDW condensate [631. They are (i) magnon excitations that manifest themselves in the antiferromagnetic resonance and (ii) the phase mode of the condensate, which is gapless provided the spin rnodulation is non-commensurate with the underlying lattice. In (TMTSF)2PF 6, the wave-vector corresponding to the optimum nesting of the Fermi surface is shown by NMR to be incommensurate although it is very close to the commensurate value (0.5, 0.25, 0). As long as Q is incommensurate with the underlying lattice a rigid motion of the density wave can be achieved without any cost in energy. This translational invariance gives rise to a collective conu'ibution to the conductivity: o(o~) = (ne2/m * )6 (o~C0o) where n is the density of condensed carriers in the SDW state, m* is the effective mass of the condensate and coo is the oscillation frequency of the collective mode, either zero for an ideally pure material, or a finite frequency when the collective mode is pinned by commensurability or impurities 1641. The interaction between the SDW and non magnetic impurities provides a finite pinning energy. A threshold electric field ET must, therefore, be reached before the condensate can contribute to the cot~duction collectively. ET is given by the balance between the energy provided by the electric field when the condensate is moved by the wavelength 27t/kF and the pinning energy, namely ETO~C0o2. In all SDW phases of the (TMTSF)2X series ET lies in the millivolt per centimeter range 1651. As shown in figure (13), the conductivity of the
i
600f
(TMTSF)2AsF6
211
I11
I,
/~
•
/
|
112 .
£
•
t
'.
h.__
s12
s/t
I
Electricfield(mV/cm)
i
113
,
,
,oT
500 1000 1500 2000 rxtemolfrequencyvl (kHz) I
l
13) Typical curJ'ent-vohage characteristics showing the onset of SDW sliding above the threshold field and the beating between an external AC current and the current oscillation of the sliding SDW state (Shapiro steps).
I
i
I
I
I
I
98
ORGANIC SUPERCONDUCTORS
Vol.92, Nos 1-2
SDW phase becomes field dependent above ET. Then the total current comprises a regular component excited through the single particle SDW gap and a collective term 1SDW. A direct evidence for a spin density wave transport is the detection of a current oscillating at a frequency which is proportional to the DC current catTied collectively 1661.The recent observation of such oscillations, the hannonic and sub-harmonic locking of this oscillation to an external ac source and a motional narrowing of the NMR spectrum in the sliding SDW state, figure (14) has provided a finn evidence for the existence of a novel collective transport in an SDW condensate I261.
(commensurate for (TMTTF)2Br). The ground state becomes an inconlmensurate SDW state due to a nested Fermi surface if the system becomes 2 or 3-D before the I-D charge localization can be active (prototype (TMTSF)2PF6). Small deviations to perfect nesting at T
7.
Charge localization in 1-D conductors remains always behind the curtain since superconductivity turns into a localized state when magnetic fields in the ten Tesla range are applied perpendicular to the conducting direction.
Conclusion and other Remarks
The study of new classes of conductors in which organic superconductivity has been discovered has contibuted to a better experimental and theoretical understanding of the physics of the I-D electron gas. Superconductivity in the (TM)2X series is only one instability among many possible ground states whose stability is governed by the competition between the charge gap (enhanced by Umklapp scattering terms g3) and the amplitude of the trzmsversecoupling t~. lntrachain repulsive interactions are important for all compounds (g l//tVF= 1) but they give rise to a clear cut spin-charge separation only when the Umklapp contribution is dominant. For strong charge localization at wave vector 4kF the intrachain exchange coupling leads to strong spin-Peierls fluctuations and concomitant I-D 2kF lattice softening which promote a spin-Peierls lattice instability (prototype example (TMTSF)2PF6). When the charge localization is weaker, the interchain AF exchange favors the propagation of 2kF spin correlations and the ground state becomes SDW
The angular dependent magnetoresistance observed in 1 and 2D conductors were not discussed in that paper. Very briefly.the phenomenon first discovered in I-D (TMTSF)2CIO4 1681has been studied in more details with (TMTSF)2PF6 [(~9.7o1.The sharp dips of magnetoresistance which are observed at angles which fullfill "magic angles" conditions, figure (15) have been explained in terms of particular electron paths which do not sweep the entire Fermi surface for these particular orientations of the magentic field but are restricted to trajectories going through the nodes of the reciprocal lattice [711. Similarly, we did not discuss the properties of 2-D organic conductors in which superconductivity can be stabilized up to 12K. Giant magnetoresistance oscillations have been detected ( [3H-(ET)213 1731 [3-(ET)21Br2 1741,K(ET)2Cu(SCN)2
1751 and ¢ x - ( E T ) 2 X , with X =
TIHg(SeCN)4) ITMand.peaks of magnetoresistance are "also visible at some peculiar orientations of the magnetic field [771.They have ben ascribed to the optimization of the Fermi surface cross section stationarity achieved at these angles 1781.Electron COtTelationsare still remaining rather strong in these materials. They are detected through the blue shift of the o(o)) oscillator strength 1791, a non-Korringa behaviour of the spin-lattice relaxation rate, a large pressure coefficient I=l.O mA ]''
D
I'"
I''
I''
I"'
[ ' ' I
121 e.,
I = 0 . 5 mA
,=
8
E 4
I=0 m A
--500
0
T=41 K
500
I000
frequency (ppm)
14) 13C-NMR spectrum in the SDW phase of (TMTSF)2PF6 showing the motional narrowing arising when a current is passed through the .sample.
---0
1
e I.,
I i i I,
-90 -60 -30
b
, I..Ii,
0
30
I.,
60
I
90
e(degree) 15) Angular oscillation of the magnetoresistance (Lebed's oscillation) in (TMTSF)2PF6, art,, C.°r [691•
Voi.92, Nos 1-2
ORGANIC SUPERCONDUCTORS
of the susceptibility 18Oland in some cases the onset of SDW ground states 18tl. Superconductivity has been raised up to 33 K in fullerene compounds 1821. The response of the magnetic properties of these 3-D molecular conductors to a pressureinduced modification of the bandwidth has revealed a very minor role of electron correlations which can be related to the remarkably large area available on each organic molecule for electron delocalization 1831. The pairing mechanism in
99
these 3-D superconductors is probably mediated by high energy instraball vibrations Thus Tc is governed by the intermolecular distance which depends on the volume of the alkali atom or on pressure [841. A c k n o w l e d g m e n t - I gratefully acknowledge the contribution of a large number of colleagues at Orsay and in several laboratories all over the world. The research activity has been partly supported by DG XII and DG XIII European Actions and DRET and CNRS contracts.
REFERENCES |1 l W.A. Little, Phys. Rev.,A 134, 1416(1964). [2 1 F. Wudl et al, J. Chem. Soc. Chem. Comm., 1453 (1970). [3 1 K. Bechgaard et al, Solid State Comm., 33, 1119 (1980). [4 I N. Thorupetal. ActaCryst. B 371236 (1981) [5 1 D. J~rome and H.J. Schulz, Adv in Physics, 31, 299 (1982). [6 ] D. J6rome, Science, 252, 1509 (1991). [7 ] H. Urayama et al, Chem. Lett., 1988.55 (1988). [8 I A.F. Hebbard et al, Nature, 350, 600 (1991). 19 ] D. J~rome et al, J. Phys. Lett.,. 41, L-95 (1980). [10 l P.M. Grant, J. Physique 44, C3-847 (1983) [ 11 ] J-P. Pouget in Low Dimensional Conductors and Superconductors p. 17,D.J6rome and L.G. Caron editors, Plenum Press New-York (1987). [12 ] J. Friedel and D. J6rome, Contemp Physics, 23, 583 (1982). [13 ] L.D. Landau and E.M. Lifshitz, Statistical Physics, p. 482, Pergamon Press (1959). [14 ] Y.A. Byschkov et al, Sov. Plays. JETP, 23, 489 (1966). [15 l V.J. Emery, J. Physique, C3, 44. 977 (1983). [16 ] J. Solyom, Adv. in Physics, 28,201 (1979). [ 17 ] V. Emery in "Highly Conducting One-Dimensional Solids", p. 247, J.T. Devreese, R.P. Evrard and V.C.E. van Doren editors, Plenum Press (1979). [18 ] H.J. Schulz, Int. J. Mod. Physics, B 5. 57 (1991). [19 ] V. Emery et al, Phys. Rev. Lett., 48, 1039 (1982). [20 ] F. Creuzet et al, J. Phys. Lett., 43, L-755 (1982). [21 1 C. Coulon. Thesis Universit6 de Bordeaux 1981 and J. Physique, 44. C-3, 885 (1982). [22 1 J.P. Pouget, private communication [23 l F. Creuzet, Thesis Universit6 d'Orsay, (1987). [24 ] F. Barthel et al, Europhys. Lett., 21, 87 (1993). 125 ] A.W. Overhauser, Plays. Rev. Lett. 4_, 462 (1960) [26 ] E. Barthel et al Plays. Rev. Lett. 71, 2825 (1993). [27 ] J.M. Delrieu, M. Roger, Z. Toffano, A. Moradpour and K. Bechgaard, J. Physique, 47,839 (1986) [28 ] T. Takahashi, Y. Maniwa, H. Kawamura and G. Saito, J. Plays. Soc. Japan 55, 1364 (1986) [29 ] C. Bourbonnais in "Low Dimensional Conductors and Superconductors", p155,D. J6rome and L.G. Caron eds NATO-ASI Series, Plenum l:h'ess NewYork (1987). [30 1 C. Bourbonnais et al. Phys. Rev. Lett. 62, 1532 (1989) [31 ] P. Wzietek et al, J. Physique I, 3, 171 (1993). [32 ] T. Moriya, J. Phys. Sot.Japan, 18, 516 (1963). [33 ] C. Bourbonnais J. Physique I, 3, 143 (1993) [34 l B. Gotschy et al, J. Physique I, 2, 677 (1992). [35 ] B. Dardel, Europhys. Lett. 24 687 (1993) [36 l T. Giamarchi, Phys. Rev., B44, 2905 (1991) and Phys. Rev. B 46, 342 (1992) [37 ] C. Bourbonnais et al, Europhys. Lett. 6 177 (1988) 138 ] M. Bousseau et al, Nouv. Jour. de Chimie 8, 3 (1984)
[39 ]
140 I [41 ] [42 ] [43 I [44 1 [45 1 [46 1 [47 ] [48 ] J49 1 [50 I 151 ] [52 1 [53 l [54 ] 155 l [56] [57 ] J58 ] [59 [60 [61 [62
1 ] ] ]
[63 ]
164 ] 165 ] [66 1 [67 ]
168 l 169 l
G. Montambaux in Low Dimensional Conductors and Superconductors, p. 233, D. J~rome and L.G. Caron editors, Plenum Press New-York (1987). K. Yamaji, in Organic Superconductors, T.lshiguro and K. Yamagi, Springer Verlag (1990). K. Yamaji, J. Phys. Soc. Japan, 51, 2787 (1982). S.S.P. Parkin et al, J.Phys.C,14, 5305, (1981) C. Bourbonnais and L.G. Caron, Europhysics, 5, 209 (1988). K. Andres et a l Phys. Rev. Lett., 45, 1449 (1980) and D. Mailly et al, J. Physique, 44, C-3, 1037 (1983). P. Garoche et al, J. Phys. Lett., 43, L-147 (1982). M. Takigawa et al, J. Phys. Soc. Jap., ,56, 873 (1987). S. Bouffard et al, J. Phys., C15, 2951 (1981). H.J. Schulz et al, J. Physique Lett.42, L-51 (1981) S. Tomic and D. J6rome, J. Plays. Cond. Matt., 1, 4451 (1989). J.F. Kwak et al, Phys. Rev. Lett 46 1296 (1981) T. Takahashi et al, J. Physique Lett. 43, L-565 (1982) M. Ribault et al, J. Physique. Lett. 44, L-953 (1983) J.R. Cooper et al, Plays. Rev. Lett.. 63, 1988 (1989). S.T. Hannahs et al, Phys. Rev. Lett., 63, 1988 (1989). P.M. Chaikin, Phys. Rev. Lett. 51, 2333 (1983) M. Ribault et al, J. Physique Lett 45, L-935 (1984) W. Kangetal,Phys. Rev.,B43,, 11467 (1991). W. Kang and D. J6rome, J. Physique I, 1, 449 (1991). M.J. Naughton et al, Plays. Rev. Lett., 61,2276 (1988). R.C. Yu et al, Phys. Rev. Lett. 65, 2458 (1990) Orsay data (to be published, 1994) M. H6ritier et al, J. Phys. Lett., 45, L-943 (1984). L.P. Gork'ov and Lebed, J. Phys.Lett., 45, L-433 (1984). For a review see "Charge density waves in Solids", L.P. Gork'ov and G. Grtinel' editors, North Holland (1989) and G. Grfiner, Synth. Metals, 43, 3767 (1991). P.A. Lee et al, Solid State Com., 14, 703 (1974). S.Tomic et al, Phys. Rev. Lett., 62, 462 (1989). G. Kriza et al, Phys. Rev. Lett., 66, 1922 (1991). A recent review of I-D Organic Conductors can be found in D. J6rome, in Organic conductors, p. 405 J.P, Farges editor, M. Dekker, New-York (1994) and C. Bourbonnais, in Highly Correlated Fermion Systems, B. Douqot and R. Rammal editors, Elsevier (to be published 1994) G.S. Boebinger et al, Phys. Rev. Lett., 64, 591 (1990). and M.J. Naughton et al, MRS-Symposia Proceedings, N ° 173, p. 257, (1990). W. Kang et al, Plays. Rev. Lett., (22, 2827 (1992).
I00 [7(1 I I71 1 I721 173 1 174 ]
ORGANIC SUPERCONDUCTORS
K. Belmia et al, Europhysics Lett. 25, 285 (1993). A.G. Lebed, JETP-Lett., 43, 174 (1986) J.M. Williamsetal, Science, 252, 1501 (1991). W. Kang et al, Phys. Rev. Lett. 62, 2559 (1989). M.V. Kartsovnik et al, JETP Lett., 48, 541 (1988) and Soy. Physics JETP 70, 735 (1990) 175 ] T. Sasaki et al Solid State Comm. 76, 507 (1990), C.P. Heidman et al, synth-Metals 41,2029 (1991) 176 l L.I. Burarov et al, J. Physique I, 4 441 (1994) [77 I M.V. Kartsovnik et al, J. Physique !, 2, 89 (1992) 178 I K. Yamaji, J. Phys. Soc. Japan, 58, 1520 (1989).
179 ]
Vol.92, Nos 1-2
C.S. Jacobsen, Thesis, Technical Univel,'sity of Denmark (1986) and C.S. Jacobsen et al, Phys. Rev. Left., 53, 194 (1984). 180 I H. Mayaffre et al, (to be published, 1994) 181 l See T. lshiguro and K. Yamaji, Organic Superconductors Springer Verlag (1990) 182 I R. Haddon et al, R.C. Haddon, Accounts of Chem. Research, 2,.~., 127 (1992). 183 1 R. Kerkoud et al. Europhysics Lett., 25, 379 (1994) and G. Quirion et al, Europhysics Lett., 21, 233 (1993). 184 ] R.M. Fleming et al, Nature, 352, 787 (1991)