The 3-Dimensional Fermi Liquid Description for the Iron-Based Superconductors

The quasiparticles in the normal state of iron-based superconductors have been shown to behave universally as a 3-dimensional Fermi liquid. Because of interactions and the presence of sharp Fermi surfaces, the quasiparticle energy contains, as a function of the momentum p\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{p}$$\end{document}, a term of the form (p-p0)3ln(|p-p0|/p0)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$( p - p_0)^3 \ln {( |p-p_0|/p_0)} $$\end{document}, where p=|p|\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p = | \varvec{p} |$$\end{document} and p0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_0$$\end{document} is the Fermi momentum. The electronic specific heat coefficient, magnetic susceptibility (Knight shift), electrical resistivity, Hall coefficient and thermoelectric power divided by temperature follow, as functions of temperature T, the logarithmic formula a-bT2ln(T/T∗)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a-b T^2 \ln {(T/T^*)}$$\end{document}, a,b\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a, \, b$$\end{document} and T∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T^*$$\end{document} being constant; these formulae have been shown to explain the observed data for all iron-based superconductors. It is shown that the concept of non-Fermi liquids or anomalous metals which appears in the literature is not needed for descriptions of the present systems. When the superconducting transition temperature TC\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{\mathrm {C}}$$\end{document} and the b / a value for the resistivity are plotted as functions of the doping content x, there appear various characteristic diagrams in which regions of positive correlation and those of negative correlation between TC\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{\mathrm {C}}$$\end{document} and b / a are interconnected; from these diagrams, we may make speculations about the types of superconductivity and the crossover between them.


Introduction
Studies of the superconductivity of iron-based superconductors have been explored to a certain extent. Concerning the nature of the normal state, however, the conclusions have not yet been settled. Many authors [1][2][3] in this field describe the system of quasiparticles as a non-Fermi liquid. The present author [4], however, states that the normal state of the substances can be described as a truly normal 3-dimensional (3D) Fermi liquid. Because of the Fermi liquid effect which is caused by the presence of interactions and the sharp Fermi surfaces, most of thermodynamic and transport properties, e.g. the electrical resistivity, behave as a − bT 2 ln (T /T * ) as a function of temperature T , where a, b and T * are constants.
According to the Green's function theory [5], the real part of the self-energy part, (ε), for the 3D Fermi liquid contain the ε 3 ln |ε| term, where ε is the quasiparticle energy measured from the Fermi energy μ; this contributes the ε 2 ln |ε| term to the density-of-states function and the T 2 ln (T /T * ) term to the electronic specific heat coefficient, T * being constant. The ε 3 ln |ε| dependence of has been really observed in Ba 1−x K x Fe 2 As 2 [6]; thus, the 3-dimensionality of the system has been confirmed.
As has been discussed in the preceding paper [4], the b/a values for the resistivity are free from the nature of impurities and represent solely the inherent characters of the electron system. If we plot b/a and the superconducting transition temperature T C as functions of the doping content x of the compounds, we obtain various kinds of characteristic diagrams for the arrangement of regions of positive correlation and those of negative correlation between b/a and T C . By observing the diagrams, we may speculate about the types of superconductivity and the crossover between those.
Since the basic substance of the iron-based superconductors shows generally antiferromagnetism, the doped materials are believed also to be antiferromagnetic. Concerning Ba(Fe 1−x Co x ) 2 As 2 compounds, however, by analysing the susceptibility [7,8] and specific heat data [9], we find the system to be nearly ferromagnetic. The general view that antiferromagnetic fluctuations dominate these systems should not be taken as universal.
Concerning BaFe 2 (As 1−x P x ) 2 compounds, Nakai et al. [10] have argued that this system exhibits the quantum critical point (QCP) near x = 0.3, since the effective mass determined by the de Haas-van Alphen effect [11] becomes very large. Shibauchi et al. [12] have also examined the QCP through measurements of the London penetration depth λ L near absolute zero; λ L is known as λ 2 L ∝ m s /n s , m s and n s being the mass and the number density of superconducting carriers. They have found that λ L shows a steep peak near x = 0.3 and concluded that, because of the quantum fluctuations at absolute zero, m s becomes very large and the system really undergoes the QCP. We have analysed the Knight shift experiment of the same system by Nakai et al. to find that the effective mass of carriers which is proportional to the density-of-states is not enhanced in the neighbourhood of x = 0.3; it is to be noted that how the interactions affect the effective mass varies with experiments.
We have already shown that the 3D Fermi liquid can possess the QCP in the presence of a magnetic field H . The free energy of the magnetized system has been shown to contain a term M 4 ln M, M being the magnetization [13,14]. Thus, in terms of dimensionless H , M and T , the magnetization relation is given by where λ is a positive constant representing the interaction strength. Because of this logarithmic form, the system is favourable for the appearance of a thermodynamic instability; the instability occurs when ∂H /∂M = 0 has two solutions; atT = 0 this leads to λ > λ C = e 5/3 /3 1.76 for the condition of metamagnetic transitions. We have already shown that Laves phase compound Y(Co 0.92 Al 0.08 ) 2 , λ = 1.88, [15] and MnSi under high pressures, λ = 1.85, [16] really undergo the transition, i.e. QCP. The temperature dependence of the Hall coefficient [1], magnetoresistance [1] and thermoelectric power [17] has also been shown to follow logarithmic formulae; these formulae explain clearly the experimental facts. When Wilson wrote the book, "The Theory of Metals" [18], the theory could not account for the experiment of the thermoelectric power. Now, the (inverse) v-shape curves for Pt and iron-pnictide compounds can be explained by αT + βT 3 ln (T /T * ) law, where α, β and T * are constants.
Among the literature related to the Fermi liquid description of the iron-based superconductors, we may list the works by Popovich et al. [19] and by Charnukha et al. [20]. These authors have shown that the specific heat and optical properties of high purity Ba 0.68 K 0.32 Fe 2 As 2 single crystal can be analysed without invoking the non-Fermi liquid concept. For description of the physics for highly correlated electron systems, we should list the work by Yamada and his colleagues [21,22]. They have shown that the normal state of high-T C -cupurates can be thoroughly described and clarified on the basis of the Fermi liquid theory. Matsuda [23] has pointed out that "the breakdown of the Fermi liquid theory" which appears frequently in the literature can be avoided through proper recognition to the differences of effective masses.

The Logarithmic Quasiparticle Energy in 3D Systems
We shall describe here the self-energy part of the quasiparticles in the normal state, ( p, ε), as required for the 3D Fermi liquid, where p is the momentum.
In 1958, Galitsky [24] showed that the quasiparticle energy of interacting fermions contains a ( p − p 0 ) 3 ln (| p − p 0 |/ p 0 ) term, where p = | p|, p 0 being the Fermi momentum. Subsequently, Migdal [5] discussed comprehensively the analytic form of the self-energy part in the Green function G( p, ε). Examining the lowest order diagrams for , he showed that Im takes the form where γ, α and α are constants. Here (ε) is not an analytic function of ε; the analytic function in the upper half plane, σ (ε), is given by which states that, because of the presence of the Fermi surface, Im σ has a discontinuity equal to 2αε 3 when the sign of ε changes. Using the dispersion relation for σ (ε), he derived where μ is the Fermi energy (chemical potential), Re 0 (ε) is expressed in a power series of ε, and Im 0 (ε) is given by Eq. (2.1). The ε 3 ln |ε| term in Eq. (2.3) confirms the result of Galitsky. It should be noted that the logarithmic αε 3 ln (|ε|/μ) term arises exclusively from the presence of the sharp Fermi surfaces in 3D systems. From the ε 3 ln |ε| dependence of the quasiparticle energy, we have ε 2 ln |ε| dependence of the density-of-states function, and the T 2 ln T dependence of the electronic specific heat coefficient γ (T ). The T 2 ln T dependence of γ (T ) has been confirmed in liquid He 3 [25,26], heavy fermion compounds [27,28], and iron-based superconductors [4]. Here it should be remarked that, for heavy fermion CeAl 3 [27], an extremely clear T 2 ln T term of γ (T ) has been observed.

The Kink Phenomenon of the 3D Quasiparticle Self-Energy
The quasiparticle behaviours, energies and states around Fermi surfaces of (Ba, K)Fe 2 As 2 and (Sr, K)Fe 2 As 2 were thoroughly examined by Wray et al. [6] through the angle-resolved-photoemission spectroscopic (ARPES) experiments. They measured the self-energy of quasiparticles in the superconducting and normal states and found characteristic kinks in Re (ε) around ε = 40 meV. They expected that the kinks might be associated with some collective modes such as phonons and spin fluctuations. They have found, however, that phonons are inadequate and the relation to the collective spin modes cannot be confirmed.
Here we analyse their ARPES results. In Fig. 4 of Ref. [6], the quasiparticle dispersion curves, Re (ε) and Im (ε)/2v F , v F being the Fermi velocity, for normal and superconducting states of (Ba, K)Fe 2 As 2 are shown. We carefully read the experimental curves at 11K of this figure, on the basis of Eq. (2.3), to obtain where we have used v F = 0.7 eV· Å from Ref. [6]. In deriving Eq. (3.1), the zero point, Re (ε) = 0, is shifted to ε = −10.4 meV; this is caused possibly by experimental errors and an additional energy of breaking the Cooper pair. Re (ε), Eq. (3.1), is drawn in Fig. 1. We should comment the accuracy of Eq. (3.1); the root-mean-square (rms) error of Eq. (3.1) to the original curve is found to be 0.39%. If the curve is expressed, without the logarithmic ε 3 ln |ε| term, by a power series up to ε 3 , ε 4 or ε 5 , the corresponding rms error is 1.93, 0.59 or 0.43%; this demonstrates the existence of the logarithmic term to be essential. The energy of the kink point is determined as 33.9 meV from Eq.  Real part of the self-energy, Re (ε), for quasiparticles in (Ba, K)Fe 2 As 2 at 11K is plotted as a function of ε. Circles show the experimental data [6], and the solid line is the theoretical curve, Eq. We have to note that, contrary to the general belief, the energy of the kink (34 meV) has no special meaning; if we assume that the ε 3 ln |ε| term is related to a collective or coherent motion with the characteristic energy σ , then we may have a term aε 3 ln (|ε|/σ ), to which we have to add ordinary bε 3 term from Eq. (2.3), where a and b are constants; here a arises exclusively from the interactions, while b is the sum of an interaction-free (band) term and an interaction-dependent term. Since the kink energy is a fraction of σ e −b/a , which is usually quite different from σ , since e −b/a ranges from 0 to ∞, depending on the interaction strength.

The Electrical Resistivity
We have already stated the temperature dependence of the electrical resistivity, ρ(T ), for highly correlated electron system in Ref. [4]. Because of the importance of the problem, here we restate the essence of the problem. For evaluating ρ(T ), we need the knowledge of the density-of-states function for 3D quasiparticles, ν(ε). The quasiparticle energy with momentum p for interacting 3D fermion systems is given, by Galitsky [24], in terms of the scattering length for 2 quasiparticles, a, where s = ( p − p 0 )/ p 0 , and O(s, s 2 , s 3 ) denotes a sum of terms of order s, s 2 and s 3 ; 2 represents the interaction strength of the system, where m is the mass of a particle and k F = p 0 /h is the Fermi wave number. From this, ν(ε) for both spins is derived as where ε * 0 = p 2 0 /2m * , m * = p 0 /(∂ε/∂ p) p 0 is the effective mass of a quasiparticle. The electrical resistivity arises mainly by the Coulomb scattering and impurity scattering. In terms of the corresponding conductivities σ C and σ imp , we have These conductivities are given, in terms of the corresponding relaxation time τ (ε), as where e is the electronic charge, v 2 = (∂ε/∂ p) 2 is the velocity squared of the quasiparticle, and f is the Fermi distribution function, f = (e βε + 1) −1 with β = 1/k B T , k B being the Boltzmann constant.
Here v 2 in Eq. (4.4) contains, according to Eq. (4.1), a logarithmic term where v 0 = p 0 /m, and s 1 is a constant. For more general Fermi liquids, since where η 1 , η 2 and η * are constants. Concerning ν(ε), we can write, as a generalization of Eq. (4.2), where ν 1 , ν 2 and ν * are constants. Concerning the Coulomb scattering, because of the Umklapp processes, τ (ε) is given by τ 0 (ε)/T 2 , and the conductivity is given by σ C = (r T 2 ) −1 , where τ 0 (ε) is a function of ε and r is a constant [29]. If we consider the effect of logarithmic terms in v 2 ν(ε) of Eq. (4.4), there appears a term of the order T 4 ln T for σ C , which will be ignored in this paper. Thus we have, for the Coulomb resistivity, (4.8) Concerning the impurity scattering, the form of the relaxation time τ (ε) varies with assumptions; whether the scattering of a quasiparticle with impurities is treated as the single particle scattering, the Born approximation or unitary approximation, τ (ε) becomes independent of, inversely proportional to, or proportional to ν(ε). Thus, we have where φ 1 and φ * are constants, and φ 2 takes −η 2 + ν 2 , −η 2 , or −η 2 + 2ν 2 for above three assumptions. From Eqs. (4.4) and (4.9), σ imp is obtained as where the order T 2 term, ( The total resistivity, Eq. (4.3), is given by It should be noted that the r T 2 Coulomb resistivity is included in the second term of Eq. (4.11).
If we make the ratio ρ 1 /ρ 0 , since ρ 1 /ρ 0 = σ 1 /σ 0 = (π 2 /3)k 2 B φ 2 , it does not depend on the character of impurities, but reflects solely the nature of the electron system; as seen from Eq. (4.10), ρ 1 /ρ 0 ∝ φ 2 represents the relative magnitude of the The observed resistivity of the normal state of LSCO cupurates [30], Ce-based oxysulfide [31] and iron-based superconductors [4] has been confirmed to follow clearly the Δρ ∝ T 2 ln T law. It should be remarked that extremely clear agreement between theory and experiment has been attained in a heavy fermion compound, CeCu 4 Ga [28].

Nearly Ferromagnetic Ba(Fe 1−x Co x ) 2 As 2
It is generally believed that, in the iron-pnictide materials, the antiferromagnetic fluctuations play an important role for the behaviours of this system. We here examine the type and degree of the fluctuations by taking the Ba(Fe 1−x Co x ) 2 As 2 system as an example.
In order to define the magnetic properties in the normal state, we have to find the value of Landau parameter F a 0 given by where χ(0) and γ (0) = (2π 2 /3)k 2 B ν(0) are the magnetic susceptibility and electronic specific heat coefficient at 0 K, and μ B is the Bohr magneton. In Eq. (5.1) F a 0 , which ranges −1 < F a 0 < ∞, specifies whether the system is nearly ferromagnetic, paramagnetic or nearly antiferromagnetic.
Here we shall find the F a 0 value for x = 0.085 specimen whose T C , 23 K, is near the highest in this series. The magnetic susceptibility for this specimen has been determined, on the basis of Wang et al.'s data [7], in our preceding paper [4]; This T 2 ln (T /T * ) variation can be confirmed by comparing with the Knight shift, K (T ), experiment by Ning et al. [8]. From their data for x = 0.09 in Fig. 3, we find It is remarkable that the relative temperature variations for χ(T ) and K (T ) are almost the same in spite of data by different groups.
The corresponding γ (0) value can be obtained from the electronic specific experiments by Yoshizawa et al. [9,32]. Their data have shown that, for the range 50 K T 100 K, γ (T ) values do not change for changes of doping x, 0 ≤ x ≤ 0.245. Thus, we may assume that, if the system stays non-superconducting, the γ (T ) values for 0 K ≤ T 50 K also should not change with x; γ (0) value for x = 0.085 is the same as that for the non-superconducting sample x = 0.245. For x = 0.245, their data of γ (T ) can be precisely fitted by Thus, by putting χ(0) = 5.88 × 10 −4 emu/Oe·mol and γ (0) = 12.5 mJ/mol·K 2 in Eq. (5.1), F a 0 for x = 0.085 is obtained as It is remarkable that the system is not nearly antiferromagnetic, but nearly ferromagnetic. The spin susceptibility is enhanced by (1 + F a 0 ) −1 = 3.4 times compared with the value evaluated by the real (exact) density-of-states. The frequently mentioned statement that, because of the antiferromagnetic fluctuations, the system tends to be a non-Fermi liquid is not applicable to this system.

The 3D Fermi Liquid Description for the BaFe 2 (As 1−x P x ) 2 System
Whether the system of quasiparticles in the normal state of iron-pnictide superconductors is a Fermi liquid or non-Fermi liquid has been a subject of principal concern for a long time. One of the most notorious iron-pnictides is the BaFe 2 (As 1−x P x ) 2 system in the neighbourhood of x = 0.3.
Kasahara et al. [1], Nakai et al. [10] and Shishido et al. [11] have performed extensive experimental studies using high-quality single crystals of this substance; measurements of the resistivity, Knight shift, spin-lattice relaxation time, 1/T 1 T , and de Haas-van Alphen effect have been taken. They have asserted that the 2-dimensional antiferromagnetic fluctuations play a central role for the non-Fermi liquid behaviour of this system; in the superconducting dome, the system has an antiferromagnetic critical point at T = 0 and x = 0.33; this is also a quantum critical point (QCP). For reasons of the non-Fermi liquid nature of the system, they have presented 3 points; outside and inside of the dome the effective mass of the quasiparticle is greatly enhanced near x = 0.33; the resistivity obeys the T -linear law and the Kohler's rule is violated.
On the contrary, we shall show that the system can be described in terms of the 3D normal Fermi liquid; outside of the superconducting dome in the T vs x phase diagram, the T -dependence of the resistivity, Hall coefficient and Knight shift follows the T 2 ln (T /T * ) variation. The reason why the Kohler's rule is not valid in this system is presented. Even inside of the dome, we can show that, from the analysis of the Knight shift experiment, the itinerant carriers behave as a 3D normal Fermi liquid. Fig. 2 Temperature dependence of the electrical resistivity of BaFe 2 (As 0.67 P 0.33 ) 2 from Ref. [1]. As described in the text, the T -dependence is analysed on the basis of 10 experimental points of the curve First we shall analyse their resistivity data for BaFe 2 (As 0.67 P 0.33 ) 2 in Fig. 3 of Ref. [1]. The authors of this Ref. have reported the results on the basis of the empirical law; n = 1.0 or 1.05 in the formula ρ(T ) = ρ 0 + AT n , where ρ 0 , A and n are constants. They have claimed that this system is a non-Fermi liquid, since n is highly deviated from 2. We reexamine their results; as shown in Fig. 2, we choose ten representative points from 40 K to 130 K on the curve, carefully read the points, and try to find the best fitting value of n. We have selected six n-values from n = 1.00 to 1.15, and evaluate the rms error in %, δ, as a function of n to find, δ = 0.493 + 40.75(n − 1.069) 2 ; the best fitting value of n is 1.07, the error being 0.49%, and the resistivity is given by Correspondingly, the resistivity formula for the T 2 ln T variation is obtained as and the error is 0.164% which is just 1/3 of that of the T 1.07 law. The T 2 ln T law reproduces the experimental data the most accurately. Generally stating, if one adopts the ρ(T ) = ρ 0 + AT n law, ∂ 2 ρ/∂T 2 ∝ n(n − 1)T n−2 should be positive for n > 1. If one observes the experimental data for Ba 2 Fe 2 (As 1−x P x ) compounds, Fig. 1a of Ref. [1], 6 curves from x = 0.14 to 0.41 are convex for the temperature range, T C T ≤ 300 K, and ∂ 2 ρ/∂T 2 < 0. Thus, the ρ ∝ T n law cannot describe the experimental data. Concerning a difficulty concerning the Hall coefficient, R H , which exceeds 1/ne, n being the number of electrons (or holes), they have argued that quasiparticle scattering times are different on different parts of Fermi surfaces. It should be remarked that the simple two-band model is not applicable to this system. The difficulty has arisen from

ρ(T, 0)}/ρ(T ) is given by a function of μ B H/ρ(T ) alone. Hereafter we define ξ = (μ B H/ρ(T )) 2 , η = Δρ(H )/ρ(T ) and the Bohr magneton μ B ≡ 1. If the Kohler's rule is valid, η(ξ )
should be a universal function of ξ ; their experimental data, Fig. 3 of Ref. [1], however, have been shown not to satisfy this requirement; η depends on both ξ and T . Now we shall prove that the Kohler's rule cannot be applied to the 3D Fermi liquid. Since, in the theory of magnetoresistance [33], the effect of melting of the Fermi surface by H 2 is the same as that by T 2 , ρ(T, H ) is given by where ρ's, a's and b's are constants. If we expand this in powers of H 2 up to H 4 , we have where and Here ρ(T ) = ρ(T, H = 0). The validity of Eq. (6.6) can be confirmed by examining whether or not A 1 /ρ(T ) is a linear function of ln T , and A 2 T 2 /ρ(T ) 3 is a constant or not. By analysing their experimental data of x = 0.33 sample for the range 40 K ≤ T ≤ 100 K, we find that A 1 /ρ(T ) is precisely fitted by and A 2 T 2 /ρ(T ) 3 is constant, 0.0018 μ ·cm·K −2 ·T −4 within 8% scattering. Thus, the violation of the Kohler's rule is not due to the non-Fermi liquid effect but represents the intrinsic nature of the 3D Fermi liquid. Nakai et al. [10] have discussed the appearance of QCP in this compounds in the vicinity of x = 0.3. They have found experimentally that the spin-lattice relaxation, T 1 , follows 1/T 1 T = A + C/(T + θ), A, C and θ being constants; θ changes sign at x ∼ 0.3 where the antiferromagnetic transition occurs. This transition persists to absolute zero, thus the appearance of QCP. Shishido et al. [11] have determined the effective mass m * through measuring the level spacing of the de Haas-van Alphen experiment, eh H/m * c, where e and c are the electronic charge and light velocity. Their results are reproduced in Fig. 4; m * /m increases steeply from 1.9 to 3.2 for decreasing x from 0.64 to 0.41. Fig. 4 Effective mass ratio m * /m and the Knight shift ratio κ(x) ≡ K spin (0, x)/K spin (0, x = 0.64) for BaFe 2 (As 1−x P x ) 2 are plotted as functions of x. Solid circles show the experimental data, reproduction of Fig. 3b of Ref. [11], and open circles are deduced from the experimental data of Ref. [10]; solid lines are smoothed curves Shibauchi et al. [12] have also examined the presence of QCP for the same substance. They have measured the London penetration depth λ L near absolute zero by three different methods to find that λ L shows a steep peak near x = 0.3; λ 2 L is proportional to the mass of superconducting carriers. They have concluded that this effective mass enhancement is caused by the quantum fluctuations arising from the QCP transition at absolute zero. In relation to this, Yanase et al, [21] have shown theoretically that λ L becomes longer if the vertex correction for the transport equation is rightly considered.
We reexamine the effective mass enhancement on the basis of Nakai et al.'s experiment for the Knight shift K (T ). K (T ) consists of the T -independent chemical shift, K chem , and T -dependent K spin (T ) which is proportional to the spin susceptibility χ(T ); K chem is estimated as 0.018 ± 0.019 %. We carefully read the experimental data and find that the data follow the T 2 ln T dependence; for instance, x = 0.33 specimen above T C , we obtain where 0.128 % is the extrapolated value to 0 K from the curve above T C ; K spin (0) = 0.110 %. It is known that the susceptibility χ(0) and hence K spin (0) is proportional to (m * /m)/(1 + F a 0 ). In Fig. 4 we plot κ(x) ≡ K spin (0, x)/K spin (0, x = 0.64) on the basis of 5 experimental points from x = 0.27 to 0.64. It is seen that κ(x) is almost constant. Thus, we conclude either that m * /m and F a 0 simultaneously tend to infinity near x = 0.33 because of the antiferromagnetic transition, or that the both quantities hardly change with x. Our conclusion is that the effective mass enhancement is not confirmed here. Circles are experimental data [34], and solid lines show theoretical fit to the data from Eqs. (6.11) and (6.12) in the text It should be remarked that how the interactions affect the effective mass varies with experiments. In the de Haas-van Alphen experiment, m * is determined by measuring the level spacing of Landau levels, while m * obtained from the susceptibility is a quantity proportional to the density-of-states on the Fermi surface. These two effective masses are irrelevant.
It is important to clarify the nature of the superconducting state of this system. We pay attention to the Knight shift experiment shown in Fig. 5, the reproduction of Fig.  1 of Ref. [34]. It should be noted that, for x = 0.33 and 0.41, T < T C , the system does not show perfect diamagnetism but show T -dependent paramagnetism; the superconductivity is nodal, and appreciable numbers of itinerant carriers are present. For x = 0.33, K (T ) decreases almost linearly with T from 0.13% to nearly zero within the uncertainty of K chem ; the system seems to be nearly antiferromagnetic. By carefully observing the experimental data for x = 0.41, we can find that the data follow the T 2 ln T variation. It is surely shown that K (T ) for x = 0.33 and 0.41 follows T 87.5 %, (6.11) and Thus, the system of itinerant carriers in the superconducting dome is truly a 3D normal Fermi liquid. From the above Knight shift data, we can estimate the number density of nonsuperconducting carriers, n * , in the superconducting state. Since χ(0), and hence K spin (0), is proportional to n 1/3 , on the crude assumption that (m * /m)/(1 + F a 0 ) is near to 1, from Eqs. (6.10), (6.11) and (6.12), we obtain n * /n = (0.0086/0.128) 3 [17], and solid lines are the theoretical fit to the data from Eq. (7.5) in the text. The inset shows the Q(T ) curve for Pt which appeared in the Wilson's book [18]; the theory, Eq. (7.4), fits the curve precisely up to about 150 K carriers; Nakai et al. [10] have reported that the analysis of the nuclear spin-lattice relaxation time leads to n * /n = 0.34 for x = 0.33. In order to confirm these values, we should consult other experiments such as the specific heat measurement.

The Thermoelectric Power of LaFe 1−x Co x As O 0.89 F 0.11
In the classic weighty book by A. H. Wilson, "The Theory of Metals", published in 1953, he described precisely the thermoelectric power of metals [18]. He showed there experimental data on the thermoelectric power Q of three metals; the Q vs T curve of Pt is reproduced in the inset of Fig. 6. Being faced with these data, he wrote, "The behaviour at very low temperatures is still more complicated and cannot at the moment be reconciled with the theory". The inverse-v shape, or peak, structure of the Q vs T curve at low temperatures was not able to be explained at that time.
The thermoelectric power is defined as an average of ε given by Here we need to consider higher order term in Eq. (4.9); where φ 3 and φ * * are constants. After taking the average, Q is obtained as up to order T 3 , where T * Q is a constant which depends on φ * , φ * * , φ 2 and φ 3 . We apply Eq. (7.3) to reproduce Q of Pt, and obtain This formula reproduces precisely the experimental data in the inset of Fig. 6 up to 150 K; thus, the appearance of the peak is explained by the Fermi liquid model. In order to get agreement at higher temperatures, we need to add the higher order term, q 2 T 5 ln (T /T * * Q ), where q 2 and T * * Q are constants. We have already shown that the thermoelectric power of a prototype heavy fermion compound, CeAl 3 , can be explained by this law [27].
Concerning the experimental study of iron-pnictide-compounds, we here take Q of La Fe 1−x Co x AsO 0.89 F 0.11 examined by Sato et al. [17]. On the basis of extensive experimental data including Q and R H , they have described these substances as an "anomalous metal". Here we shall show that the substances are not anomalous; the T -dependence of Q for this system is analysed by the form In Fig. 6, their data for x = 0.0125, 0.05 and 0.1 in Fig. 4a [17] and our analysis are presented. As shown in the figure, the data for the temperature range above the crystalline structural transition temperature, 75 K T 250 K, can be precisely fitted by Eq. (7.5) with q 1 = 1.28, 1.03 and 0.631 μV/K 2 , T 1 = 223, 233 and 263 K and T * Q = 479, 491 and 504 K; these q 1 values (the absolute value) are from 4 to 9 times larger than that of Pt, but about 1/30 of that of CeAl 3 .
We have also analysed the experimental data of the Hall coefficient of this system and find that the data follow the a − bT 2 ln (T /T * ) law. Thus, our conclusion is that the quasiparticles of this system are not "anomalous" but form a real 3D Fermi liquid.

The ρ 1 /ρ 0 and T C vs x Diagrams and Superconductivity Types
As mentioned in Sec. 4, the ρ 1 /ρ 0 values for ρ(T ) represent the inherent nature of electrons; from the ρ 1 /ρ 0 values we may obtain precious information concerning the electron system. As described in the preceding paper [4], we have found that, for Ba(Fe 1−x Co x ) 2 As 2 , the ρ 1 /ρ 0 values have a close correlation with the superconducting transition temperature T C as a function of the doping content x. In the original BCS theory, the Cooper pair is formed between two quasiparticles whose momenta differ by 2 p 0 , while the logarithmic quasiparticle energy, 3 ln , arises mainly from the interactions of two quasiparticles whose momenta are almost equal or differ nearly by 2 p 0 ; thus, the strength and the nature of this interaction is closely connected with T Cvalues. It is frequently stated in the literature that the spin density fluctuations produce s +− -wave superconductivity [35,36], while the charge density (or orbital) fluctuations produce s ++ -wave superconductivity [37,38]. It is also stated that nodeless s ++ superconductors have relatively high T C values while nodal s ++ superconductors have low T C values.
When our analyses proceed to other compounds, however, we shall find that the correlation between ρ 1 /ρ 0 and T C is not always positive; each compound has a proper characteristic diagram in which the regions of positive, negative and no correlation are variously intermingled. The variety of the ρ 1 /ρ 0 and T C vs x diagrams reflects the variety of the electron correlation; the point of correlation change may suggest the crossover of the superconductivity types.
In order to confirm above views, we shall present to analyse the ρ 1 /ρ 0 and T C vs x diagrams of various superconductors.
First, in Fig. 7a, we reproduce the T C and ρ 1 /ρ 0 vs x diagram of Ba(Fe 1−x Co x ) 2 As 2 for the range 0 ≤ x 0.18 deduced from Chu et al.'s experiment [39]; Wang et al.'s experiment [7] yields essentially the same figure. Here we see a clear positive correlation between ρ 1 /ρ 0 and T C ; the superconductivity seems to arise mainly via intra-band parings. Since the single particle excitations are not observed in the specific heat experiment, the superconductivity seems to be nodeless s +− -wave state. As has been discussed in 5., however, because this system is nearly ferromagnetic, we should carefully examine the possibility of the p-wave superconductivity similar to liquid He 3 and UPt 3 .
Next we take up PrFeP 1−x As x O 0.9 F 0.1 system for which Takemori et al. [40] and Miyasaka et al. [3] have measured ρ(T ); in Fig. 7b, the ρ 1 /ρ 0 and T C vs x diagram is shown from the ρ(T ) data by Takemori et al. for the range 100 K ≤ T ≤ 300 K; Miyasaka et al.'s data give almost the same diagram. The figure shows that the system can be classified to 2 regions of x; for 0 ≤ x 0.6, T C is relatively low and ρ 1 /ρ 0 and T C show a negative correlation, while for 0.8 x ≤ 1.0, T C is high and they show a positive correlation. Miyasaka et al. have mentioned anomalies of this system for the range 0.6 x ≤ 0.8; they have stated that, although the antiferromagnetic fluctuations seem to play the principal role for this superconductivity, the spin fluctuation theory cannot explain observed anomalies. By observing Fig. 6b, we see that, just in this range, the system seems to exhibit the transition; the crossover of states from s ++ to s +− -wave occurs.
In Fig. 7c, the ρ 1 /ρ 0 and T C vs x diagram for BaFe 2 (As 1−x P x ) 2 system from Kasahara et al.'s data [1] is presented. At a glance, T C and ρ 1 /ρ 0 have an intimate correlation as functions of x. By a close look, however, although ρ 1 /ρ 0 increases by 5 times for the change of x from 0.27 to 0.33, T C keeps constant at about 30 K; there is no correlation for the range of x up to 0.40. Kontani and Onari [37] have described, with consideration of the electron-phonon interaction, the superconducting state showing a crossover from s +− -wave to s ++ -wave state as the impurity concentration increases.
Here it should be remarked that, as shown in 6., the Knight shift inside the super-  [41], and f LaFe 1−x Co x As 0.89 Fe 0.11 compounds from Sato et al.'s experiment [17]; solid and open circles show T C and ρ 1 /ρ 0 values, and solid lines are smoothed curves conducting dome exhibits partly the normal state values; this may be the proof of the p-wave superconductivity for which we need more careful observation.
In Fig. 7d, the ρ 1 /ρ 0 and T C vs x diagram for LaFeP 1−x As x O 0.9 F 0.1 , from Miyasaka et al.'s experiment [3], is shown. Here clear correlation between T C and ρ 1 /ρ 0 is not Fig. 8 A schematic diagram for obtaining higher-T C superconductors on the basis of the ρ 1 /ρ 0 and T C vs x diagrams seen; low T C and high T C regions are clearly separated in the vicinity of x = 0.3 where the crossover from s ++ to s +− -wave state seems to occur. In view of high-T C and low ρ 1 /ρ 0 values for x 0.4, here the interband pairing seems to be dominant.
The diagram for NdFe 1−x Ru x AsO 0.89 F 0.11 is presented in Fig. 7e from Lee et al.'s experiment [41]. Remarkably, here is found no correlation between T C and ρ 1 /ρ 0 ; although, for a change of doping x from 0 to 0.1, ρ 1 /ρ 0 decreases to 1/60, T C decreases only by 10 K(2%); in spite of that the small Ru-doping affects tremendously the electronic states, it hardly affects the electrons (or holes) which are associated with superconductivity. Lee et al. have stated that, in the small x-region, the lowering of T C is too small to be explained by the s +− -wave state. We cannot judge, by simply observing Fig. 7e, whether, in the vicinity of x = 0.3, the crossover from s +− to s ++ -wave state really occurs or not; Eliashberg equation for 3D system may give the definite answer.
Finally, the diagram for LaFe 1−x Co x AsO 0.89 F 0.11 is presented in Fig. 7f from Sato et al.'s experiment [17]. Here we find very clear correlation between ρ 1 /ρ 0 and T C ; when the ρ 1 /ρ 0 value becomes negative, the superconductivity vanishes. From the diagram, it is clearly seen that the interaction of two quasiparticles whose momenta differ by 2 p 0 is attractive for 0 ≤ x ≤ 0.075.

A Program for Obtaining Higher-T C Superconductors
By observing diagrams presented in the preceding section, we can pursue the possibility of obtaining higher-T C superconductors. In Fig. 7b for PrFeP 1−x As x O 0.9 F 0.1 , if we can find a substance whose electronic structure (Fermi surfaces) is situated on a virtual system x A (> 1) as the continuation of the ρ 1 /ρ 0 curve for 0.8 x ≤ 1, as shown in Fig. 8, then we expect the corresponding T C value, T A , may be much larger than T C at x = 1, T M ; thus, the problem is whether or not we may find compounds situating on this line by using As, Sb or/and Bi. The same argument can be applied to the compound NdFeP 1−x As x O 0.9 F 0.1 .
Furthermore, on the basis of Fig. 7e for NdFe 1−x Ru x AsO 0.89 F 0.11 , if we can find the virtual system x < 0, as the smooth continuation of the ρ 1 /ρ 0 curve to x < 0, then we may get higher T C substances; the problem is whether or not we can find such transition elements in place of Ru in NdFe 1−x Ru x AsO 0.89 F 0.11 .

Conclusions
We have shown that the itinerant carriers (electrons or holes) in the normal state of all iron-based superconductors form strictly a 3D Fermi liquid. Because of interactions and the presence of sharp Fermi surfaces, the thermodynamic and transport properties obey the logarithmic formula with respect to T 2 , H 2 or M 2 , M being the magnetization; the electronic specific heat coefficient, magnetic susceptibility, electrical resistivity, Hall coefficient and thermoelectric power divided by temperature follow the formula a − bT 2 ln (T /T * ). This formula explains observations for the normal state of all iron-based superconductors. Remarkably, in certain cases, even the itinerant carriers inside the superconducting dome behave as a 3D normal Fermi liquid.
In most papers, in this field the authors describe the resistivity through the formula ρ = ρ 0 + AT n , ρ 0 , A and n being constant; when n deviates from 2, they call the system a non-Fermi liquid. In view of the present paper, the notion of non-Fermi liquids or anomalous metals is not needed.
Finally, in order to determine theoretically the superconductivity type and transition temperature, we have to solve the Eliashberg equation on the basis of the 3D system.