Evolution of shape and volume fraction of superconducting domains with temperature and anion disorder in (TMTSF)$_2$ClO$_4$

In highly anisotropic organic superconductor (TMTSF)$_2$ClO$_4$, superconducting (SC) phase coexists with metallic and spin density wave phases in the form of domains. Using the Maxwell-Garnett approximation (MGA), we calculate the volume ratio and estimate the shape of these embedded SC domains from resistivity data at various temperature and anion disorder, controlled by the cooling rate or annealing time of (TMTSF)$_{2}$ClO$_{4}$ samples. We found that the variation of cooling rate and of annealing time affect differently the shape of SC domains. In all cases the SC domains have oblate shape, being the shortest along the interlayer $z$-axis. This contradicts the widely assumed filamentary superconductivity along $z$-axis, used to explain the anisotropic superconductivity onset. We show that anisotropic resistivity drop at the SC transition can be described by the analytical MGA theory with anisotropic background resistance, while the anisotropic $T_c$ can be explained by considering a finite size and flat shape of the samples. Due to a flat/needle sample shape, the probability of percolation via SC domains is the highest along the shortest sample dimension ($z$-axis), and the lowest along the sample length ($x$-axis). Our theory can be applied to other heterogeneous superconductors, where the size $d$ of SC domains is much larger than the SC coherence length $\xi$, e.g. cuprates, iron based or organic superconductors. It is also applicable when the spin/charge-density wave domains are embedded inside a metallic background, or vice versa.


Introduction
Raising the superconducting transition temperature (T c ) has been the goal of active research for a century. Compounds like cuprates [2,3], iron based superconductors [4], organic superconductors (hereafter denoted as OrS) [5] are some major high-T c superconductors at ambient pressure. These materials have several common properties: (i) layered crystal structure and, hence, high conductivity anisotropy; (ii) interplay between various types of electron ordering, i.e. between spin/chargedensity wave and superconductivity; (iii) spatial inhomogeneity. Therefore, many effects and methods, both experimental and theoretical, are common for these materials. A recent and good review on all these systems is given by Stewart [6]. Our paper concerns the OrS and is devoted to two problems: (i) show that Maxwell-Garnett approximation can be used to estimate superconducting volume fraction in coexistence regime of superconducting, metallic and spin/charge density wave phases; (ii) analysis of corresponding experimental data [7,8] in the organic superconductor (TMTSF) 2 ClO 4 to study the effect of cooling rate and of disorder on the volume fraction, shape and size of SC domains.
(TMTSF) 2 X series belongs to quasi-1D OrS and has been widely studied for 40 years [5,[9][10][11]. Many important effects have been discovered and investigated on these compounds, e.g. angular magnetoresistance oscillations (AMRO) in quasi-1D metals [10][11][12], field-induced spin-density waves (FISDW) [13,14], etc [10,11]. An interesting and puzzling property of these materials related to our subject is that, with the increase of pressure for (TMTSF) 2 PF 6 or of anion ordering for (TMTSF) 2 ClO 4 , the superconductivity first appears along the least conducting z-axis, while along the most conducting x-direction only in the last turn [7,8,15,16]. Recent experimental study on (TMTSF) 2 PF 6 [15][16][17][18] and (TMTSF) 2 ClO 4 [7,8] has shown that spin density wave (SDW), superconducting (SC) and metallic phase coexist in form of segregated domains. With the increase in pressure (for (TMTSF) 2 PF 6 ) [15,16,19] or in anion ordering (for (TMTSF) 2 ClO 4 ) [7,8] the volume fraction of SC or metallic phase increases. arXiv:2012.03231v2 [cond-mat.supr-con] 1 Jul 2023 Few theories has been suggested to describe the coexistence regime in these materials, especially in (TMTSF) 2 PF 6 . One of them is the application of SO(4) symmetry that explains SDW and SC coexistence [20] but does not account for the observed hysteresis [17], for the strong enhancement of the upper critical field H c2 [21] and for the anisotropic SC onset [15,16,19] in the coexistence phase. Less exotic theories suggest a separation of SDW and SC in a coordinate [15][16][17][21][22][23][24][25] or momentum space [24]. The momentum SC-SDW separation assumes a semi-metallic state in a SDW phase, where, due to the imperfect nesting, small ungapped Fermi-surface pockets appear and become superconducting [24]. The H c2 enhancement can be explained in both scenarios [24,25], but the observed hysteresis suggests a spatial SDW/SC separation [17]. To explain the H c2 enhancement [21] the SC domain width must not exceed the in-plane SC penetration depth λ ab ≡ λ. In (TMTSF) 2 ClO 4 the in-plane penetration depth is [26][27][28] λ ab (T = 0) ≈ 1 µm, and the out-of-plane penetration depth is [29] λ bc (T = 0.19 K) ≈ 40 µm. One possible mechanism of the formation of such narrow domains in the SDW state could be the soliton phase [15,22,23,25,30]. It suggests that SDW order parameter becomes non-uniform with metallic domain appearing perpendicular to the highest conducting x-axis. However, the width of these soliton-wall domains, being of the order of SDW coherence length ξ SDW ∼ 30 nm, is too small to be consistent with the recent observation of AMRO and FISDW in (TMTSF) 2 PF 6 [16] and in (TMTSF) 2 ClO 4 [31], suggesting the domain size d > 1 µm. Moreover, the soliton-phase scenario accounts for the SC suppression along the most conducting x-axis only, but it could not explain why SC first appears along the least conducting z-axis, because in this scenario the soliton walls are extended along both y and z-axes, which should result to SC along both these directions.
A probable explanation of this anisotropic SC onset was proposed recently [1]. It is based on two ideas. First, in anisotropic media the isolated SC islands increase conductivity much stronger along the least conducting direction than along the others, as observed in FeSe [32,33] and described [32][33][34] using the Maxwell-Garnett approximation (MGA) [35] for small volume fraction φ of SC phase. However, this MGA theory cannot explain the anisotropic zero-resistance onset. For this we need the second idea, which takes into account the finite sample size L as compared to the size d of SC grains. If the sample shape is very anisotropic, e.g. a thin plate, then the current percolation via the SC grains, responsible for the zero-resistance onset, is most probable along the shortest sample dimension, [1] i.e. along the sample thickness, and least probable along the sample length. These two ideas were applied to explain [1] the experimental data [15,16,19] in PF6, taken on thin elongated samples with dimensions [15,19] 3 × 0.2 × 0.1 mm 3 . As (TMTSF) 2 ClO 4 samples are usually flat shaped either [7,8], similar ideas can be applied to analyze the resistivity experimental data there too. In this paper using the MGA we find superconducting volume ratio φ and shape of inclusions from available experimental data. We also investigate how disorder and cooling rate affect the inclusions' shape and size. This knowledge may help to better understand the microscopic structure and electronic properties of the SDW/SC coexistence phase in organic superconductors.
In Sec. 2 we briefly describe the important properties of (TMTSF) 2 ClO 4 and argue that MGA can be applied to estimate temperature dependence of SC volume ratio (φ) in the presence of all 3 phases, i.e. metallic, SDW and SC. In Sec. 3 we describe the theoretical model in MGA. Further, in Sec. 4 we apply our theoretical model to analyze the experimental data on (TMTSF) 2 ClO 4 . Using the experimental data from Ref. [8] we find out how cooling rate effects φ. Similarly, using the experiments of Ref. [7] we analyze the evolution of aspect ratio of SC inclusions with sample disorder. Finally, in Sec. 5 we discuss the main results of our investigation and their consequences.

Material
(TMTSF) 2 ClO 4 is the only member of Beechgard salts which becomes superconducting at ambient pressure [36,37]. It is a quasi-1D superconductor in which cooper-pair ordering (SC) coexists with insulating Pierels ordering (SDW) 1 . This very competition between SC and SDW is the key to understand the unconventional superconductivity [38][39][40][41]. Among quasi-1D superconductors (TMTSF) 2 ClO 4 is the only compound in which superconductivity can be controlled by the cooling rate of samples, which affects the disorderliness of ClO 4 anions. Slowly cooled (TMTSF) 2 ClO 4 samples (relaxed state) undergo an SC transition at T c ≈ 1.3 K [42,43]. However, when cooled very fast (quenched state) (TMTSF) 2 ClO 4 has an insulating SDW transition at T SDW ≈ 4 − 5 K [44][45][46]. This behavior can be ascribed to structural change that occurs at the anion ordering temperature [47] T AO ≈ 24.5 K. For T > T AO the noncentrosymmetric tetrahedral ClO 4 anions, which are located at the inversion centers [48,49], preserve the inversion symmetry due to thermal motion of ClO 4 anions. ClO 4 anions randomly occupies one or other orientations, hence, on average the inversion symmetry is preserved [48]. For T < T AO , if the sample is cooled fast enough then the randomness of orientation of ClO 4 anions is preserved [45]. However, if the sample is cooled slowly through T AO , ClO 4 anions along a,c-axes are ordered uniformly; and along b-axis ordered alternatively [47,50,51]. The anion ordering introduces the new wave vector and the Fermi-surface folding. It disturbs the Fermi-surface nesting, preventing the SDW and favoring SC.
Recent experiments strongly support the presence of SC inclusions embedded in the background of metallic/SDW phase when (TMTSF) 2 ClO 4 is cooled at intermediate rate [7,8]. This means granular superconductivity for partially ordered samples. From the crystallographic point of view, in this coexistence phase the domains of ordered ClO 4 anions are embedded inside disordered background. In these domains the anions have alternating ordering pattern, being disordered outside the domains [7,8,47]. The electronic state inside these domains is assumed to remain metallic even at T < T AO , and at T < T * these ClO 4 -ordered grains becomes superconducting [8]. Here T * is the superconducting onset temperature. The volume fraction φ of superconducting phase increases with the decrease of temperature. For slow or intermediate cooling rate at T = T c the phase coherence of these SC islands establishes in the entire sample, leading to its zero resistance. In the temperature interval T c < T < T * the effective medium model for highly anisotropic heterogeneous layered compounds [35], developed and applied by the authors to various anisotropic compounds in Refs. [32][33][34]52], can be used to estimate the volume fraction and the shape of SC inclusions inside the samples.

Method
To analyze the temperature dependence of resistivity at T c < T < T * we use the Maxwell-Garnett approximation (MGA) [35], valid when the volume fraction of SC phase φ 1. Note that the condition φ 1 is fulfilled in a wide range of parameters, because in 3D anisotropic samples even the SC percolation threshold φ c , corresponding to the onset of nearly zero resistance, is considerably smaller than unity. This is illustrated in our recent work [1], where φ c in needle-shaped flat (TMTSF) 2 PF 6 samples from Refs. [15,19], typical for organic superconductors, was calculated numerically, assuming a rectangular sample shape and ellipsoidal SC inclusions of varying size at randomly distributed but fixed positions. This calculation has shown the strong anisotropy and small value of these φ c (see Figs. 3 and 4 of Ref. [1]). Thus, for SC domain size d = 40 µm the percolation threshold along the x-axis φ c ≈ 0.3; for d = 15 µm φ c ≈ 0.15 (see Fig. 3 of Ref. [1]). For y and z axes φ c is even smaller. In all these cases φ c 1, suggesting a wide range where MGA can be applied. The (TMTSF) 2 ClO 4 samples, studied in the present paper and in Refs. [7,8], are also flat and needle-shaped, similar to (TMTSF) 2 PF 6 samples from Refs. [15,19].
In recent works [32][33][34]52] using MGA the superconducting volume ratio φ was found when SC inclusions were embedded inside metallic background. Similarly, here we also assume SC inclusions are embedded inside a background phase. However, here the background phase consists of SDW and metallic phases. This is permissible in MGA approximation as long as we know the effective conductivity of this mixed background phase. We make two more assumptions: (i) the SC inclusions are of ellipsoidal shape, which simplifies the calculations and allows deriving analytical formulas for conductivity; (ii) the size d and distance between SC inclusions l are much greater than the SC coherence length ξ, so that the SC proximity effects and the Josephson coupling between the SC grains do not change the results considerably. Both in (TMTSF) 2 PF 6 and (TMTSF) 2 ClO 4 the size d 1 µm of metal/SC inclusions is indeed much larger than the SC coherence length [8,53] ξ = 70, 30, and 2 nm along the a, b, and c axes, respectively. Such a large metal/SC grain size is evidenced by the observation [16,31] of angular magnetoresistance oscillations (AMRO) in the mixed phase of these compounds. If φ 1, one can also take l ξ. Due to the proximity effect a shell of SC condensates around SC inclusions with thickness ∼ ξ gets created, which changes the effective size and shape of SC islands. Another quantum mechanical effect is the Josephson coupling between SC grains. The Josephson coupling energy E J depends directly on Josephson junction current I c , E J ≡hI c /2e, which exponentially decreases with the increase of distance l between neighbouring SC grains: I c = I 0 exp(−l/ξ). If E J T, the Josephson coupling can be disregarded [54]. As we assume l, d ξ, both these effects can be neglected. The conductance through N-S boundaries of a normal metal and SC may increase, maximum two times, due to the Andreev reflection 2 . However, this increase of the interface conductance should not considerably change the effective conductivity of the whole sample because the main voltage drop (at a given current) comes not from the N-S interfaces, but from the large non-SC parts of the sample. Even if we take an infinite N-S interface conductance, this does not much affect the total sample resistance.
Considering the above facts we neglect these small quantum mechanical effects, and treat conductivity due to embedded SC inclusions using the classical MGA approximation and percolation.

Theory
We use the Maxwell-Garnett approximation (MGA) [55,56] to find the SC volume ratio φ. We derive a general formula for effective conductivity in anisotropic heterogeneous media, where unidirectional ellipsoidal SC (SDW) inclusions are embedded inside the background SDW/metallic (metallic) phase with anisotropic resistivity. We denote the diagonal effective conductivity tensor of this materials as diag(σ xx , σ yy , σ zz ). Similarly, the background conductivity tensor of this material is denoted as diag(σ b xx , σ b yy , σ b zz ). Background phase consists of both metallic and SDW phase, or only of the metallic phase. Since the MGA in its standard form [35] can be applied only to an isotropic medium, we use the coordinate-space dilation to transform the initial problem with anisotropic conductivity into isotropic one [32][33][34] (see appendix A for the details of this mapping and Eqs. (A4) and (A5) for the definition of dilation coefficients µ and η).
Before explaining the idea of MGA [35,55,56] we note that the stationary-current equation for electrostatic potential V(r), coming from the continuity equation for the electric current J i in heterogeneous media with coordinate-dependent diagonal conductivity σ(r) 3 is completely equivalent to the electrostatic equation in heterogeneous media with coordinate-dependent dielectric constant ε(r), as comes from the Maxwell's equations: where E i and D i are the electric and displacement fields, and ∇ i ≡ ∂/∂x i . Hence, the problem of an effective dielectric constant of such a medium with heterogeneous dielectric function ε(r) is equivalent to the problem of effective conductivity of a heterogeneous medium with the same coordinate dependence of σ(r).  To explain the idea of MGA 4 [35,55,56] of calculating the effective dielectric constant of heterogeneous isotropic media with a small volume fraction φ of inclusions of the second phase, we refer to Fig. 1. The background phase is represented as "2", and the inclusion phase as "1" in the figure. We take a sphere of our heterogeneous media. Inside this sphere the inclusions are embedded as shown in Fig. 1. Let this sphere be embedded in an infinite medium of background phase. An electric field E is applied to the infinite medium. Let E be the electric field at a distant point P. E includes polarization due to individual inclusions. We denote the inclusion volume ratio as φ and the background phase volume ratio as (1 − φ). We assume that the heterogeneous sphere with both inclusion and background phase can be substituted by a sphere of homogeneous phase with an effective dielectric constant ε e f f . Then the electric field at point P can be found in two ways: (i) By summing the polarization effect due to individual inclusions; the corresponding electric field at point P we denote as E . (ii) By taking the polarization effect at point P of homogeneous sphere with effective dielectric constant ε e f f ; the corresponding electric field at point P we denote as E . The MGA assumes these fields to be the same, i.e. E = E , which gives an equation on ε e f f 5 . Using the well-known formula for the polarization of dielectric ellipsoid in isotropic medium, and replacing ε ii by σ ii , we find the following equation (see Eqs. 18.9 and 18.10 of Ref. [35]) for the effective conductivity σ * i along main axis i in the mapped space 6 : Here A i are the depolarization factors of a dielectric ellipsoid with semiaxes a i in the mapped space: The analytical solution of this integral can be found in terms of incomplete elliptic integrals of first and second kind 7 . If the inclusions have finite conductivity, then solving Eq. (3) for σ * i we obtain If the inclusions are superconducting, we take σ isl → ∞. Eq. (3) in this case simplifies to Solving Eq. (6) for σ * i we obtain Both Eqs. (5) and (7) gives the effective conductivity σ * in the mapped space. Although the background conductivity in the mapped space is isotropic, the effective conductivity σ * is anisotropic because of the anisotropic ellipsoidal shape of inclusions. The effective conductivity of original anisotropic material in real space is found by the reverse mapping. It is 5 See Sec. 18.1.1 of Ref. [35] for complete discussion on MGA. 6 Eqs. (3),(5) and (8) assume that in the mapped space the conductivities of both phases, σ isl and σ b , are isotropic.
For a generalization to anisotropic σ isl or σ b in the mapped space see Refs. [58,59]. In the case of superconducting inclusions, applied in Sec. IV to analyze the experimental data, σ isl = ∞ is naturally isotropic. 7 See appendix B of Ref. [34]. done via multiplying the effective conductivity matrix in the mapped space σ * = diag(σ * xx , σ * yy , σ * zz ) by the inverse mapping coefficient matrix diag (1, µ, η). Hence, multiplying Eq. (5) by diag(1, µ, η), we find the effective conductivity in real space: Here σ b ii is the background conductivity along the axis i in real space. Similarly, multiplying Eq. (7) by diag(1, µ, η), we find the effective conductivity of original inhomogeneous material with SC inclusions of volume fraction φ: From Eq. (9) one can express the volume fraction φ via the effective σ ii and background σ b ii conductivities and depolarization factor A i along the same axis: Eq. (8) is helpful when the conductivities of background and inclusion phases are both finite. Hence, it can be used to find the effective conductivity of heterogeneous material, when, e.g., SDW domains are embedded inside a metallic background, or vice versa. Eq. (9) can be used when superconducting inclusions are embedded inside a background phase of finite conductivity. Below we use Eqs. (9) and (10) to analyze experimental resistivity data in (TMTSF) 2 ClO 4 in the mixed SC/SDW state. Here the background phase is made up of metallic as well as SDW phases. Resistivity along the corresponding axis is found by taking the inverse of Eq. (8) and Eq. (9).

Analysis of experimental data in (TMTSF) 2 ClO 4
We consider partially ordered (TMTSF) 2 ClO 4 samples. We denote T * as superconducting onset temperature. In these compounds for T > T * there is no superconductivity. However, for T < T * the domains containing ordered ClO 4 anions partially transform to superconducting inclusions [8]. These ordered domains are embedded inside the phase of unordered ClO 4 anions, where SDW prevails but may coexist with metallic phase. Further cooling results in the increase of SC volume fraction φ and in the formation of coherent clusters of SC inclusions. At T = T c < T * a complete SC channel gets opened [7,8], i.e. the SC phase coherence establishes in the whole sample. In Sec. 4.1, we use our theory to calculate the SC volume ratio φ. In Sec. 4.2, we study the influence of cooling rate on SC volume ratio. In the end, in Sec. 4.3, we find the approximate shape of SC inclusions in various disordered samples.

Application of MGA theory to describe resistivity and to find superconducting volume ratio φ
To find a typical temperature dependence of SC volume ratio φ(T) in (TMTSF) 2 ClO 4 at cooling rate − dT/ dt ≤ 100 K/min we choose the sample #4 in Fig. 2 of Ref. [7] 8 due to the availability of experimental resistivity data on ρ xx (T, H = 2 T) in a magnetic field H = H z , shown in Fig. 2(c) of Ref. [7]. The magnetic field destroys superconductivity, and we can use these data to find the conductivity of background phase σ b xx (T) = 1/(ρ xx (T, H = 2 T) − ∆ρ xx ) (see the inset in Fig. 3), where the offset ∆ρ xx = ρ xx (T * , H = 2 T) − ρ xx (T * , H = 0) accounts for magnetoresistance of metallic phase at H = 2 T. A magnetic field H 500 Oe is usually enough to destroy superconductivity in (TMTSF) 2 ClO 4 [60,61], but we take the data at H = 2 T where the SC effects can be safely ignored. Since the experimental data on ρ zz (T) under magnetic field for the same samples are absent, the background-phase conductivity σ b zz along the z-axis is found by extrapolating the The effective medium resistivity ρ xx containing SC, metallic and SDW phases (green squares), and the background resistivity ρ b xx in the absence of SC islands (blue triangles), taken from the data in Fig. 2c of Ref. [7] in magnetic field. metallic ρ zz (T) resistivity to low temperature by a second-order polynomial, similar to Ref. [61]. Here the second-order term comes from the electron-electron scattering at low temperature [61,62]. We take the x-axis as the reference axis for mapping to isotropic medium, i.e. σ b xx is taken as the background isotropic conductivity in the mapped space: σ b = σ b xx . According to Eq. (A4), the mapping coefficient along the z-axis is defined as η = σ b zz /σ b xx . The SC volume ratio φ is found from Eq. (10) for i = x. The calculated volume ratio as a function of temperature is plotted in Fig. 2. Substituting φ found from Eq. (10) for i = x to Eq. (9) for i = z, we predict the effective resistivity along the z-axis, σ zz (T). Its comparison with the experimental data from Fig. 2b of Ref. [7] is shown in Fig. 3. We see a rather good agreement. In this calculation we have one fitting parameter -the unknown ratio of the semiaxes a x and a z of ellipsoidal SC inclusions. We found that at high temperature T ≈ 1.2 K, typical inclusions have the aspect ratio a z /a x ≈ 0.16. At low temperature T ≈ 0.5 K, the inclusions have aspect ratio a z /a x ≈ 0.85. It means that with a decrease in temperature the SC inclusion becomes more isotropic along x and z-axes, i.e. a z /a x → 1. It may indicate the formation of large and almost isotropic clusters of small SC inclusions. Note that at any temperature the found aspect ratio a z /a x is much larger than the ratio of coherence lengthes in (TMTSF) 2 ClO 4 , ξ z /ξ x ≈ 0.03. This supports the fact that in (TMTSF) 2 ClO 4 the heterogeneity and SC islands originate from disorder and anion ordering rather than from usual SC fluctuations, because for SC fluctuations a z /a x ∼ ξ z /ξ x . Fig. 3 and the inset in Fig. 2 show the temperature dependence of resistivity for the same sample along the z and x axes correspondingly. From their comparison one observes that the resistivity drop near T c for ρ zz is much stronger than along for ρ xx . This feature originates from the strong anisotropy of background-phase resistivity ρ b ii and is naturally described in generalized MGA theory [32][33][34]. The qualitative interpretation of this anisotropic resistivity drop due to SC onset is illustrated in Fig. 1 of Refs. [32] or [33]. Because of high resistivity ρ b zz , the interlayer current mainly flows via the SC islands serving as shortcuts for this current direction. The effective resistivity ρ zz is then determined by the much smaller intralayer resistivity and by the typical length of in-plane path between two close SC domains, which is inversely proportional to the SC volume fraction φ.

Effect of cooling rate on superconducting volume ratio
The cooling rate of (TMTSF) 2 ClO 4 samples controls the fraction of ClO 4 -ordered domains. At slow cooling the ClO 4 anions have enough time to relax into ordered state. At fast cooling the thermal disorder remains in the samples, so that both ordered and disordered domains coexist. It was corroborated by resistivity [63], specific heat [64] and x-ray scattering [65] experiments.
The volume fraction φ o of anion-ordered domains as a function of cooling rate has been studied using the x-ray scattering [65,66] and, recently, by resistivity and magnetic susceptibility [8] measurements. The corresponding results are compared in   [8] in the range of cooling rate 0 < − dT/ dt < 20 K/min. Several assumptions are made in extracting φ o from these experiments [8,65]. First, (i) all these results assume that at the slowest cooling rate − dT/ dt| min ∼ 1 K/min all ClO 4 are ordered at low temperature. It is not evident, as some degree of anion disorder may remain. Second, (ii) in the estimate of volume fraction φ o from the resistivity measurements in the mixed SDW/metal phase in Ref. [8] the following equation (see Eq. 1 of Ref. [8]) for the effective conductivity σ zz along z-axis has been used: where ρ min = 1/σ min = 0.03 Ω·cm is taken as a residual resistance of the sample with lowest cooling rate and ρ max = 1/σ max = ρ min + ∆ρ c * = 0.26 Ω·cm is determined assuming that the difference ∆ρ c * = ρ max − ρ min = 0.23 Ω·cm is equal to the jump of resistivity at anion-ordering temperature T AO = 24.5 K due to the scattering by anion disorder. In fact, at low temperature the anion disorder has much stronger effect on conductivity than just the electron scattering by this disorder itself, because it also favors the formation of insulating SDW state. Even if the fraction of insulating SDW domains is about one half, as in Fig. 4 of Ref. [8], it may considerably affect the electron conductivity. In addition, (iii) Eq. (11) does not take into account the conductivity anisotropy of (TMTSF) 2 ClO 4 , which strongly enhances the effect of metal/SC domains on resistivity along the least-conducting axis [32][33][34], as given by Eq. (8). (iv) The extraction of φ o from the magnetic susceptibility χ(T) data, especially at rapid cooling rate when the SC volume fraction φ 1, depends strongly on the size and shape of SC domains [34,54]; therefore the assumption [8] .02 K/min is not valid when the size of SC domains is smaller than the London penetration depth.
In this subsection we estimate the volume fraction φ of SC phase in (TMTSF) 2 ClO 4 using the resistivity data from Ref. [8] and applying Eqs. (9) and (10) derived in MGA. Note that the volume fractions φ of SC phase and φ o of anion-ordered phase may differ, e.g., because the former depends on temperature. Since the MGA approximation is valid only at φ 1, it can only give φ(T) at T > T c . Thus, our method of estimating φ works better for higher cooling rate when T c is lower, therefore our results are rather complimentary to those in Ref. [8]. However, the cooling-rate dependence of φ(T > T c ) also gives the general tendency. Note that at cooling rate dT/ dt = 100 K/min and some annealing, by extrapolating the φ(T) curve in Fig. 3 to T = 0 we obtain φ(T → 0) ≈ 0.3, which is in a good agreement with other data in Fig. 4 of Ref. [8].
To observe the influence of cooling rate on SC volume ratio φ(T), we use resistivity data from the inset of Fig. 1(d) of Ref. [8]. These experimental values are taken as the effective resistivity ρ zz along z-axis for different cooling rates. For background-phase resistivity ρ b zz along the z-axis we use the 2nd order polynomial fit of metallic resistivity at T > T * . Substituting these ρ zz (T) and ρ b zz in Eq. (10) we obtain φ(T) for various cooling rates. Unfortunately, in Ref. [8] there is no resistivity data along other two axes which would allow us to find the ellipsoid aspect ratio. Therefore, in Fig. 4 we take the depolarization factor A i = 1/3, i.e. the spherical inclusions in the mapped space instead of ellipsoidal. This choice corresponds to ellipsoid semiaxes a i ∝ σ ii (T * ) ∝ ξ i in real space, as one expects for SC fluctuations. 9 The obtained φ(T) at cooling rates 0.02 K/min, 0.052 K/min, 2.5 K/min, 7.6 K/min and 18 K/min are shown in Fig. 4. The curves in Fig. 4 are similar to those in Fig. 4 of Ref. [8], but the values of SC volume fraction φ are expectedly smaller than φ o in Ref. [8], because at T > T c only a fraction of ClO 4 -ordered domains becomes superconducting. But φ increases with decreasing temperature and, probably, reaches φ o at T → 0.  Figure 4. Dependence of SC volume ratio φ on cooling rate calculated using Eq. (10) for different temperatures. For this calculation the experimental data from the inset of Fig. 1(d) in Ref. [8] are used.

Effect of disorder on the shape of superconducting inclusions
In Fig. 4 we investigated the effect of cooling rate on φ. However, along with φ the shape of SC domains also plays an important role. Recent work [1] has shown that the probability of percolation along the shortest sample dimension, i.e. along sample thickness, is higher than along other directions. It was corroborated by the experiment on FeSe [33,52], where by reducing the z-axis thickness of sample from 300 nm to ∼ 50 nm, one raised T c from 8 K to 12 K [52]. The (TMTSF) 2 ClO 4 samples are also usually flat. This effect of anisotropic superconductivity onset also depends on the shape of SC inclusions [1]. The knowledge of the shape of SC domains in (TMTSF) 2 ClO 4 is also helpful to better understand the mechanism of their formation. In Sec. 4.1 we found a z /a x for sample # 4 in Ref. [7] at cooling rate 100 K/min. Below we find a z /a y and a y /a x for the samples cooled at rate 600 K/min with various times of subsequent annealing. This gives the effect of disorder on the shape of SC inclusions.
To study the evolution of aspect ratios a z : a y : a x with disorder we use the experimental data from Figs. 3 and 4 of Ref. [7]. Unfortunately, the curves with equal numbers in these two figures correspond to different samples. Thus, we do not have the data on resistivity along all three axes for the same sample and parameters, required to determine the full shape of ellipsoidal inclusions. However, we use the fact that the depolarization factors A i in Eq. (4) depend most strongly on the semiaxis a i along the same direction, which allows us to vary only one parameter for each fit.
First we find the evolution of a z /a y with disorder. From the resistivity data along the z-axis, given in Fig. 4a of Ref. [7], using Eq. (10) we find φ(T) for different degrees of disorder. Using these φ(T) in Eq. (9) we predicted resistivity along y-axis. From best fit values of predicted and experimental resistivity along y-axis, given in Fig. 4b of Ref. [7], we find the ratio a z /a y for various degrees of disorder. The results are shown in Fig. 5a and reveal that at low disorder, up to sample #9, the ratio a z /a y almost remains same. However, at higher disorder the ratio a z /a y increases.
Similarly, to find the evolution of a y /a x with disorder we use the data from Fig. 3 of Ref. [7]. As before, from the resistivity data along y-axis using Eq. (10) we find φ(T) for different degrees of disorder. Using this φ(T) we predict resistivity along x-axis. We change the semiaxes of ellipsoids along yand x-direction, so that the theoretical and experimental values of resistivity agree. Thus obtained ratio of a y /a x is shown in Fig. 5b. . The dependence of aspect ratios a z /a y (Fig. a) and a y /a x (Fig. b) of superconducting domains in (TMTSF) 2 ClO 4 on disorder at two temperatures T = 0.8 K and 0.6 K, calculated using Eqs. (4,9,10) and resistivity data from Figs. 4 and 3 of Ref. [7], taken at cooling rate 600 K/min and various annealing times. At longer annealing time, i.e. at weaker disorder, the shape of SC domains is more anisotropic.

Discussion
In this paper we propose a method based on MGA to investigate the microscopic parameters of heterogeneous superconductors from resistivity data. We apply this method to the organic superconductor (TMTSF) 2 ClO 4 , where SC coexists with SDW in the form of isolated domains. Using our method we study the SC volume fraction φ and the shape of SC islands as a function of external parameters, such as temperature and ClO 4 anion disorder, which can be experimentally controlled by the cooling rate through the anion-ordering transition at T AO ≈ 24.5 K [8,65] or by the annealing of (TMTSF) 2 ClO 4 samples [7].
For the best use of proposed method, one needs the following experimental data: (i) Temperature dependence of resistivity ρ ii (T) along each of non-equivalent main crystal axes 10 ; (ii) ρ ii (T, H 0 ) in a magnetic field H 0 > H c destroying SC, to get the resistivity ρ b ii (T) of the background homogeneous phase. In the absence of ρ ii (T, H 0 ) or in the case of non-SC inclusions, one needs to make an extrapolation of ρ ii (T) from T > T * to lower T in order to get ρ b ii (T), which is less accurate. In the case of SDW/CDW inclusions one can also apply an external pressure destroying SDW/CDW to get ρ b ii (T). If also (iii) magnetic susceptibility data are present, especially for all non-equivalent magnetic-field orientations, they help to independently check the obtained microscopic parameters and allow the estimate of the average size of SC inclusions as compared to the SC penetration depth [33,34]. Ideally, all these data are available for several values of external parameters that one is interested in, for example, at each studied cooling rate of (TMTSF) 2 ClO 4 . Unfortunately, in spite of an active experimental investigation of (TMTSF) 2 ClO 4 by resistivity measurements, e.g. performed recently in Refs. [7,8], this full set of data is absent. Nevertheless, we have analyzed the available data from Refs. [7,8] to make some physical predictions concerning the mixed SC/SDW phase in (TMTSF) 2 ClO 4 .
In Figs. 2 and 4 we show the obtained SC volume fraction φ for various cooling rates and temperatures T c < T < T * . Fig. 3 illustrates how well the MGA model typically fits the experimental data. Figs. 5 and 6 show the evolution of the shape of SC grains with the change of disorder by the annealing of rapidly cooled samples.
From Fig. 5a we observe that, for partially ordered samples #2, #7, #9 the aspect ratio a z /a y ≈ 0.05 of SC domains at temperature T ≈ 0.6 − 0.8 K depends weakly on disorder. The corresponding SC volume ratio φ ≈ 0.1 also weakly depends on disorder. At shorter annealing time, i.e. at larger disorder, the SC volume ratio decreases to φ ≈ 0.03, while the aspect ratio increases to a z /a y ≈ 0.13. Another aspect ratio a y /a x ≈ 0.4 depends much weaker on disorder, as shown in Fig. 5b. Note that the obtained ratios a x : a y : a z at cooling rate 600 K/min are close to those of SC coherence lengthes. In Fig. 2 we found that for partially ordered sample #4 at cooling rate 100 K/min the aspect ratio of SC domains a z /a x ≈ 0.85 at temperature T ≈ 0.5 K, which differs considerably from what we get at 600 K/min even for long annealing time. This means that the decrease of cooling rate, as in Refs. [8,65] is not completely equivalent to the increase of annealing time at T < T AO used in Ref. [7]: they have similar effect on SC volume fraction φ but different effect on the shape of SC domains. Therefore, it would be very interesting to study their effect on the SC domain size, which can be extracted from the simultaneous magnetic susceptibility measurements as done for other compounds [33,34].
In organic superconductors (TMTSF) 2 ClO 4 [7] and (TMTSF) 2 PF 6 [15,16] superconductivity onsets anisotropically, i.e. first along the highest conducting z-axis and only in the end along the lowest-conducting x-axis. This behavior was first explained by assuming filamentary SC inclusions elongated along z-axis [15], but this hypothesis received neither theoretical nor experimental proof till now. Our analysis also shows that the SC domains are not elongated along z-axis but, on contrary, are oblate. Nevertheless, we predict much stronger decrease of resistivity along the least conducting direction (compare Fig. 3 and the inset in Fig. 2), similar to experimental observations in (TMTSF) 2 ClO 4 [7] and in many other heterogeneous anisotropic superconductors [32][33][34]. In our recent work [1] we proposed a simple model to explain the anisotropic zeroresistance onset also. We have shown [1] that the percolation probability along SC islands in needle or flat shaped samples is the highest along the shortest direction. A schematic illustration of this idea is given in Fig. 6. The same idea can be applied to (TMTSF) 2 ClO 4 to explain the anisotropic onset of superconductivity. Usually, the (TMTSF) 2 ClO 4 samples are much shorter along the interlayer z-axis than along other two, e.g., the dimensions of samples in Ref. [7] are 3 × 0.1 × 0.03 mm 3 . Hence, the probability of percolation for the ellipsoidal inclusion, even with the obtained anisotropic aspect ratios a z /a x and a z /a y , will be the highest along the z-axis. Our preliminary calculations of percolation threshold for such flat elongated samples confirm this statement and also show that the zero resistance, i.e. the percolation threshold, can be achieved even when the SC volume ratio φ = φ c 1. Hence, without invoking a filamentary SC one can easily explain the anisotropic onset of superconductivity in organic metals.
x z Figure 6. Schematic illustration of superconducting ellipsoid inclusions embedded inside a long thin conductor with dimensions L z L x . It shows that the percolation probability along z is higher than along x-axis if a z /a x > L z /L x .
In our semi-classical theory based on MGA we neglect the dynamic fluctuations. Probably, the position and size of domains may fluctuate. But in the studied organic metal the competition between SDW and superconductivity, leading to the domain formation, is governed by the disorder of ClO 4 anion orientation, which is frozen well below the anion-ordering transition temperature T AO ≈ 24.5 K. Hence, below T SDW ≈ 4 − 5 K we may consider the metal/SC and SDW domains as time-independent.
For all our calculation in Sec. 4 we have used Eqs. (9) and (10) because the inclusions are superconducting. However, the similar approach and Eq. (8) can be used when the conductivity of inclusions is finite, e.g., when SDW inclusions are embedded inside a metallic background or vice versa. This is a usual occurrence in organic superconductors [5]. Thus, instead of using Eq. (11) or similar phenomenological formulas to analyze the resistivity data in mixed metal/SDW or metal/CDW phases in organic metals, we recommend to use MGA formulas (8) and (4), which take into account the strong anisotropy of layered organic metals and the actual shape of domains. Apart from that, in several high-T c superconductors as Bi 2 Sr x Ca 3−x Cu 2 O 8 [67], La 2−x Sr x CuO 4 [68], Ba(Fe 1−x Co x ) 2 As 2 [69][70][71], and in many other materials, e.g., HfTe 3 [72], there is also a spatially-separated coexistence of metallic and SC or SDW/CDW phases. Hence in all these heterogeneous materials, provided the conductivity of both metallic and SDW/CDW phases is known and the domain size exceeds the coherence length, one can estimate the second-phase volume ratio and the domain shape from resistivity data using the above method.

Conclusions
Using the Maxwell-Garnett approximation we estimate the volume fraction and the shape of superconducting domains from resistivity measurements. We apply this method to investigate the heterogeneous electronic structure in organic superconductor (TMTSF) 2 ClO 4 , where superconductivity coexists with the spin-density wave in the form of isolated domains. This material is especially important to study such coexistence because it appears even at ambient pressure and can be easily controlled by changing the cooling rate or annealing time of the samples. From available resistivity data we study the evolution of the volume fraction and of the shape of superconducting domains in (TMTSF) 2 ClO 4 with disorder and temperature. Our method applies not only to superconductors, but also when the density-wave or other-type domains are embedded inside a metallic background, or vice versa.
Author Contributions: P.G. conceptualized and developed the Methodology. Formal analysis, validation, manuscript writing and review was done by K.K, P.G and V.K. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A Anisotropic dilation of the problem of static current distribution
Anisotropic medium is converted to isotropic one by mapping the real space to a mapped space, where the solution is simpler. The mapping should satisfy the following conditions: (i) The conductivity of background phase in the mapped space should be isotropic, and (ii) the electrostatic continuity equation should be satisfied in the mapped space with the same solution.
In our notations, σ b xx , σ b yy , σ b zz are the constant conductivity components of background phase in the original heterogeneous medium. Let J and V be the current density and the applied potential respectively. The electrostatic continuity equation in the background phase is then written as Heterogeneity is hidden in the boundary conditions on the surface of each grain and of the sample. Let x , y and z be the axes in mapped space, where conductivity should be isotropic: σ b x x = σ b y y = σ b z z = σ b . If J and V are the current density and electrostatic potential respectively in the mapped space, the continuity equation in the mapped space is written as The condition (ii) for our mapping means that the solution V(x, y, z) = V (x , y , z ) (A3) of Eqs. (A1),(A2) is the same. This solution completely determines the effective conductivity of heterogeneous medium, as it also gives the current density via Eqs. (A1),(A2) and (A3) are consistent if the mapping is the anisotropic scaling, e.g., with constant mapping coefficients µ and η determined by conductivity anisotropy in real space: Instead of (A4) one could choose a product of the mapping (A4) and of any isotropic scaling by a factor α with the simultaneous change of σ b → α 2 σ b . We have chosen σ b = σ b xx , so that µ, ν < 1. The current components in the real space r = (x, y, z) and in the mapped space r = (x , y , z ) are related as J x (r) = J x (r ), J y (r) = √ µJ y (r ), J z (r) = √ η J z (r ).
The shapes of inclusions are not preserved during this mapping procedure. For example, a sphere with radius a x described by the equation x 2 /a 2 x + y 2 /a 2 x + z 2 /a 2 x = 1 in non-homogeneous medium will transform to an ellipsoid described by the equation x 2 /a 2 x + y 2 /a 2 y + z 2 /a 2 z = 1 in mapped space, where the semiaxes are given by a x = a x , a y = a x / √ µ, a z = a x / √ η.
If z is the lowest conducting axis, and the highest conducting axis is x, then a sphere in real space transforms to an ellipsoid elongated along z-axis. Due to the temperature dependence of conductivity anisotropy, the coefficients µ, η and the shape of inclusions change with temperature either. If in real space the inclusion is ellipsoid with semiaxes a = a x , b = βa x and c = γa x , then it transforms to ellipsoid in mapped space with semiaxes a x = a x , a y = a x β/ √ µ, a z = a x γ/ √ η.
In MGA we take ellipsoidal inclusions with fixed aspect ratios β and γ, but varing size.