Charge transport in a polar metal

The fate of electric dipoles inside a Fermi sea is an old issue, yet poorly explored. Sr1−x\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{1-x}$$\end{document}Cax\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{x}$$\end{document}TiO3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{3}$$\end{document} hosts a robust but dilute ferroelectricity in a narrow (0.0018<x<0.02\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.0018\ <\ x\ <\ 0.02$$\end{document}) window of substitution. This insulator becomes metallic by removal of a tiny fraction of its oxygen atoms. Here, we present a detailed study of low-temperature charge transport in Sr1−x\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{1-x}$$\end{document}Cax\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{x}$$\end{document}TiO3−δ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{3-\delta }$$\end{document}, documenting the evolution of resistivity with increasing carrier concentration (n\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n$$\end{document}). Below a threshold carrier concentration, n*(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}^{* }(x)$$\end{document}, the polar structural-phase transition has a clear signature in resistivity and Ca substitution significantly reduces the 2 K mobility at a given carrier density. For three different Ca concentrations, we find that the phase transition fades away when one mobile electron is introduced for about 7.9±0.6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$7.9\pm 0.6$$\end{document} dipoles. This threshold corresponds to the expected peak in anti-ferroelectric coupling mediated by a diplolar counterpart of Ruderman–Kittel–Kasuya–Yosida (RKKY) interaction. Our results imply that the transition is driven by dipole–dipole interaction, even in presence of a dilute Fermi sea. Charge transport for n<n*(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\ <\ {n}^{* }(x)$$\end{document} shows a non-monotonic temperature dependence, most probably caused by scattering off the transverse optical phonon mode. A quantitative explanation of charge transport in this polar metal remains a challenge to theory. For n≥n*(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\ge {n}^{* }(x)$$\end{document}, resistivity follows a T-square behavior together with slight upturns (in both Ca-free and Ca-substituted samples). The latter are reminiscent of Kondo effect and most probably due to oxygen vacancies.


INTRODUCTION
The concept of a polar or "ferroelectric" metal was first proposed by Anderson and Blount in 1960s. 1 They considered a continuous structural-phase transition breaking the inversion symmetry and leading to the appearance of a polar axis in a metal. This appears counter-intuitive since one expects mobile electrons to strongly screen electric field. Recently, however, "ferroelectric" metallicity was reported in LiOsO 3 . 2 It was found that this stoichiometric metal shows a structural-phase transition to non-centrosymmetric rhombohedral phase (R3c) below 140 K. The polar structural transition, which manifests itself as a kink in resistivity of this metal, is analogue to what occurs in its insulating ferroelectric cousins LiNbO 3 and LiTaO 3 . 2 Puggioni and Rondinelli 3 have argued that such polar metals can be found when there is an unusually weak coupling between mobile electrons and transverse optical phonons, which drive ferroelectricity. This scenario has been documented in LiOsO 3 by recent ultrafast spectroscopy measurements. 4 Paraelectric solids close to a ferroelectric transition 5 provide an alternative platform for a meeting between metallicity and ferroelectricity. One example is PbTe, a narrow-gap semiconductor close to a ferroelectric instability. The unavoidable presence of doping defects makes available PbTe samples dilute metals. Isovalent substitution of Pb by Ge leads to a structural-phase transition from cubic to a non-centrosymmetric rhombohedral phase, 6 which would have been ferroelectric in absence of mobile electrons. The ferroelectric-like transition was revealed by X-ray diffraction, inelastic neutron, and Raman scattering. [6][7][8] Charge transport in presence of local dipoles was studied several decades ago. [9][10][11][12] SrTiO 3 single crystals, in contrast to PbTe, can be made stoichiometric enough to be insulating. Proximity to a ferroelectric quantum critical point 13,14 is manifested by a large electric permittivity of SrTiO 3 (ε r > 20000). 15 As a consequence, this insulator displays a number of intriguing properties. 16 A FE state emerges upon substitution of a tiny fraction of Sr ions with Ca. 17 Moreover, this quantum paraelectric can become a dilute metal (with a carrier concentration as small as %10 16 cm −3 ) upon oxygen reduction. 18,19 The dilute metal undergoes a superconducting transition below 0.3 K. [20][21][22] Rischau et al. have recently found that a superconducting phase coexists with a FE-like instability in n-doped Sr 1Àx Ca x TiO 3Àδ , and superconductivity and ferroelectricity (FE) may be intimately linked. 23 The FE transition of insulating Sr 1Àx Ca x TiO 3 was found to survive in metallic Sr 1Àx Ca x TiO 3Àδ . Being a metal, the latter does not show a bulk reversible electric polarization and cannot be a true ferroelectric. Nevertheless, it shows anomalies in various physical properties at the Curie temperature of the insulator. For example, Raman scattering found that the hardening of the FE soft mode in the dilute metal is indistinguishably similar to what is seen in the insulator. 23 The anomaly in resistivity was found to terminate at a threshold carrier density (n Ã ), near which the superconducting transition temperature was enhanced 23 providing evidence for a link between superconducting pairing and ferroelectricity, a subject of present attention. [24][25][26][27][28][29][30][31][32][33][34] In this paper, we present a study of low-temperature electrical resistivity in dozens of Sr 1Àx Ca x TiO 3Àδ single-crystals with x ¼ 0, 0:22%, 0:45%, 0:9% documenting in detail the evolution of charge transport with increasing concentrations of electric dipoles and charge carriers, a question, which was not addressed in depth by a previous study. 23 In Sr 1Àx Ca x TiO 3Àδ , metallicity and ferroelectricity are both dilute. Therefore, the distance between dipoles and mobile electrons can be separately tuned, but kept much longer than the interatomic distance. We find that the magnitude of low-temperature mobility is significantly reduced below n Ã ðxÞ and the mean-free-path gently peaks near n Ã ðxÞ, where the quadratic temperature dependence of resistivity is restored. Moreover, we find that n Ã is proportional to x, implying the threshold density occurs at a fixed ratio between the inter-carrier and the interdipole distance. We will argue that these features are all consistent with the hypothesis of a dipolar RKKY interaction, which was theoretically proposed a quarter-century ago. 35 At much higher carrier densities (n > 5 10 19 cm −3 ), we observe a slight upturn in low-temperature resistivity of Ca-free and Ca-substituted samples and attribute it to a Kondo effect associated with oxygen vacancies.

RESULTS
The upper panel (a-l) of Fig. 1 plots, the low-temperature dependence of resistivity for Sr 1Àx Ca x TiO 3Àδ at low n and x ¼ 0:22%, 0:45%, 0:9%. In Fig. 1a, e, i, resistivity shows an anomaly at lowest n for all three groups of samples. Increasing ndoping, the anomaly shifts to lower temperatures, evolves into a minimum and finally disappears at a threshold doping (n Ã ðxÞ), seen in Fig. 1b-d, f-h, j-l. Rischau et al. found that in selected samples at x ¼ 0:2% and 0:9%, this anomaly occurs close to where the structural-phase transition was detected by Raman spectroscopy, sound velocity, and thermal-expansion measurements. 23 Most recently, thermal-expansion measurements documented the evolution of the structural-phase transition in Sr 1Àx Ca x TiO 3Àδ (x ¼ 0:9%), starting from the insulating phase and extending deep into the metallic phase. 36 The study found that below n Ã , the transition temperature and the magnitude of the transition-induced anomaly continuously decreased with increasing carrier concentration. Above n Ã , a small residual anomaly with a concentration-independent temperature scale was observed to survive. No sign change in the thermal-expansion coefficient (α) was observed at n Ã . Thus, the corresponding Grüneisen ratio α=C p (C p is the specific heat) does not change sign. This implies either  17 The green, blue, red balls mark the T c of polar structural-phase transitions in n-doped Sr 1Àx Ca x TiO 3Àδ with x ¼ 0:22%, 0:45% and 0:9%, respectively. the absence of a quantum critical point at n Ã or its insensitivity to uniaxial pressure. 36 Figure 1m shows a 3D plot of the phase diagram of Sr 1Àx Ca x TiO 3Àδ with two tunable parameters (n and x). In the n ¼ 0 plane, Sr 1Àx Ca x TiO 3 becomes ferroelectric when x c > 0:18%. Above this critical concentration of Ca substituents, T c scales with x following T c jx À x c j 1=2 . 17 The anomaly in resistivity of n-doped samples is shown with green, blue, and red symbols for x ¼ 0:22%, 0:45%, 0:9%, respectively. The anomaly caused by the structural-phase transition (strictly speaking, only a true FE state at n ¼ 0) shifts to lower temperatures with increasing n. At a threshold doping (n Ã ð0:22%Þ % 5 10 18 cm −3 , n Ã ð0:45%Þ % 1 10 19 cm −3 and n Ã ð0:9%Þ % 2 10 19 cm −3 ), resistivity becomes metallic down to lowest temperatures (see panels 1d, h, and l).
The low-temperature resistivity of strontium titanate follows a simple quadratic temperature dependence: ρ ¼ ρ 0 þ AT 2 . 37,38 Elastic scattering is represented by residual resistivity, ρ 0 , intimately linked to the asymptotic low-temperature mobility set by disorder. Inelastic scattering among electrons is expected to give rise to the AT 2 term. However, the quadratic temperature dependence of resistivity has several puzzling features. One is that it persists even in the dilute regime where the Fermi surface is too small to allow Umklapp scattering among electrons. 37 The other is that this T-square resistivity extends beyond the degeneracy temperature and the exponent gradually enhances with warming. 39 Nevertheless, the magnitude of T-square resistivity prefactor corresponds to what is found empirically in other metals (dilute or correlated) with a similar Fermi energy. 22 Moreover, this prefactor evolves smoothly and universally with carrier concentration and has a similar amplitude in other quantum paralectric perovskite oxides, 40 such as EuTiO 3 and KTaO 3 in spite of their different phonon spectra. The origin of this T-square resistivity is still debated. 22,41 Both the elastic and inelastic terms of resistivity are affected by the passage across n Ã ðxÞ. The quadratic temperature dependence is drastically altered upon the introduction of Ca atoms. In Fig. 2, the resistivity of Sr 1Àx Ca x TiO 3Àδ with x ¼ 0:22%, 0:45%, and 0:9%, is plotted vs. T 2 and compared with Ca-free samples of similar n. As seen in the upper panel, which shows typical data for n < n Ã ðxÞ, the resistivity shows a minimum followed by an upturn, in presence of the structural-phase transition. Note that the non-monotonic temperature dependence of mobility cannot be a simple consequence of the temperature dependence of the electric permittivity (see Supplementary Fig. 1). In contrast, in SrTiO 3Àδ with similar n, T 2 dependence of resistivity persists down to lowest temperatures. 37,38 The lower panel of Fig. 2 shows the behavior at n ' n Ã ðxÞ. One can see that the T 2 resistivity of Cadoped samples is restored. For more data at n > n Ã ðxÞ, please refer to Supplementary Fig. 2. Moreover, the slope is similar in Casubstituted and Ca-free samples. In other words, the prefactor of quadratic resistivity, which strongly depends on carrier concentration but not on residual resistivity, is similar in Ca-substituted and Ca-free samples at the same carrier density. Let us note that this strict Fermi-liquid behavior at n Ã ðxÞ indicates that it is not a "quantum critical point" where a "non-Fermi-liquid behavior" is commonly sought and often found. This is in agreement with the absence of a sign change in thermal expansion. 36 The other implication of this observation is that below n Ã ðxÞ, electrons suffer additional inelastic scattering in presence of aligned dipoles. The destruction of the structural-phase transition at n Ã ðxÞ suppresses this additional mechanism and restores the T 2 resistivity associated with electron-electron scattering. 37 This implies the existence of additional scattering in the polar metallic state, totally smearing T-square resistivity.
The presence of FE-like order affects the elastic scattering of the carriers too. The Hall mobility (μ) at 2 K for Sr 1Àx Ca x TiO 3Àδ with x ¼ 0, 0:22%, 0:45%, 0:9% is presented in Fig. 3. The mobility μ is extracted from Hall resistivity ρ yx and longitudinal resistivity ρ xx using μ ¼ ρ yx ρ xx B . As seen in the figure, μ increases with the decrease in carrier density (n) following an approximate power law μ $ n Àα shown by the dashed lines. For SrTiO 3Àδ , the power law behavior persists down to n $ 2 10 18 cm −3 . Below this concentration, it begins to saturate and then drops at even lower carrier concentrations with the approach of the metal-insulator transition. 18,21 In many doped semiconductors, the mobility decreases with increasing carrier concentration. 42,43 This has been often discussed in the framework of ionized impurity scattering. 44 Several features distinguish the 2 K mobility of metallic strontium titanate from ordinary doped semiconductors. First of all, the dependence of mobility with carrier concentration is very steep (i.e., α in μ $ n Àα is close to unity). Second, this mobility is strongly temperature dependent and passes from a room-temperature value of 5 cm 2 V −1 s −1 to 20,000 cm 2 V −1 s −1 at liquid He temperature. 39,45 This latter value implies that the low-temperature carrier mean-freepath (') becomes much longer than the interdopant distance (l ee ¼ n À1=3 ).
Dingle provided an expression for low-temperature mobility in a degenerate-doped semiconductor, based on the Born approximation for scattering. 46 The formula can be written as where n is the carrier concentration and a Ã B is the effective Bohr radius. This expression does not seem to provide a good description of our data, as seen in Fig. 3a.
Recently, it has been argued 47 that these peculiar features of mobility in dilute metallic strontium titanate can be traced back to the long effective Bohr radius, a Ã B , of the parent insulator. The large electric permittivity 15 elongates the Bohr radius to 600 nm, which is to be compared to 1.5 nm in silicon. This in turn affects the Thomas-Fermi screening length of the metal, which depends on it: The combination of a large a Ã B and a small k F elongates the screening length and therefore short-distant irregularities in the dopant distribution are smoothed out. This simple approach yields Fig. 2 Resistivity plotted as a function of T 2 at low temperatures for Sr 1Àx Ca x TiO 3Àδ with x ¼ 0:22%, 0:45%, and 0:9%, compared with x ¼ 0 at similar carrier densities. When n < n Ã ðxÞ (a-c), the Casubstituted samples show a strong deviation from the T 2 temperature dependence seen in Ca-free samples. On the other hand, when n % n Ã ðxÞ (d-f), the T 2 temperature dependence is maintained in Ca-substituted samples and the slope remains the same.
this expression for mobility: 47 (2) Figure 3a confirms that this expression gives a surprisingly good account of the variation of 2 K mobility with n in SrTiO 3Àδ . 47 As seen in Fig. 3b-d, Ca substitution leads to a slight decrease in the exponent of the power law exponent (α). More importantly, a clear deviation from the power law behavior occurs below n Ã ðxÞ marked by vertical arrows. In other words, the carrier mobility is significantly reduced when the system orders. Figure 3e shows the mean-free-path (') calculated through ' ¼ μ_k F =e, in which μ is taken from the Hall mobility, k F is the Fermi wave vector deduced from k F ¼ ð3π 2 nÞ 1=3 by assuming a single isotropic Fermi surface. It presents a mild maximum marked by three arrows at n Ã ðxÞ. Well above this threshold density of n Ã ðxÞ, the mean-free-path of the Ca-doped and Ca-free samples gradually merges. Thus, we conclude that the presence of the FE-like order inside the metal has drastic consequences for elastic scattering of electrons too. Figure 3f shows n Ã extracted from our data as a function of n Ca . One can see that the two concentrations are proportional to each other. In other words, the threshold inter-electron distance (l Ã ee ) linearly scales with the average distance between Ca ions (n À1=3 Ca ), the slope of which indicates that the destruction of the FE-like order happens when there is one mobile electron per 7.9 ± 0.6 Ca ions.

DISCUSSION
Thus, we find that below a threshold concentration: (i) an additional mechanism for inelastic scattering sets in; (ii) the lowtemperature mobility is significantly reduced. Moreover, (iii) this threshold density for the destruction of the polar metal is proportional to Ca concentration; and (iv) at this density carrier mean-free-path gently peaks. We are now going to argue that these observations support a picture in which off-center Ca sites generate electric dipoles interacting with each other inside a Fermi sea.
Let us begin with a fundamental question: what drives the FE transition in the insulating Sr 1Àx Ca x TiO 3 for 0:0018 < x < 0:02? According to the scenario invoked by previous authors, 17,48 an offsite Ca atom generates an electric dipole, which couples to the soft transverse optical phonon of the host lattice. Theory 49 had predicted that cooperative phenomena between paraelectric defects in an easily polarizable crystal lead to polar clusters whose size is set by the polar correlation length of the host lattice (r c ). The latter depends on the velocity and the frequency of the soft mode of the host polarizable lattice (r c ¼ v s =ω 0 ). 50,51 With cooling r c increases, the interaction between clusters becomes stronger and at a sufficiently low temperature, cluster percolation leads to a ferroelectric order. A transverse Ising model, treating dipoles as pseudospins, has been successful in describing the Sr 1Àx Ca x TiO 3 phase diagram. 52,53 In an ordinary ionic crystal, there is no mean-field basis for favoring ferroelectricity. This is because the electric field generated by a dipole does favor opposite alignments for a neighboring dipole along perpendicular orientations (Fig. 4a). In an easily polarizable crystal, on the other hand, the interaction becomes ferroelectric over a polarization length (r c ), thanks to the presence of soft transverse optic phonon modes. Since r c is significantly longer than the typical distance between dipolar impurities, sizable ferroelectric clusters [49][50][51] (Fig. 4b) become energetically stable. This paves the way for the emergence of long-range order.
This approach provides a satisfactory explanation for two experimental observations. The first is the presence of a ferroelectric order triggered by off-center impurities in quantum paralectrics such as SrTiO 3 and KTaO 3 54 , and its absence in ordinary ionic crystals with no soft mode such as alkali halide crystals. 50 The second observation is the existence of a critical concentration for the emergence of ferroelectricity in Casubstituted SrTiO 3 . 17 Let us note, however, that extrinsic Ca ions in strontium titanate lattice can also interact with each other via the strain field, they create by reducing the size of the lattice cell which they occupy. The latter feature has been invoked to explain the weakening of FE order for x < 0:02. 55 In this context, our observations raise two questions specific to metallic Sr 1Àx Ca x TiO 3Àδ . First of all, at small n, how can a phase transition driven by the dipolar interaction among Ca sites hardly be affected by the capacity of mobile electrons supposed to screen any electric field? A second question arises by the existence of n Ã ðxÞ. Why does the order eventually fade away and what does this threshold density correspond to?
The first question brings to mind a remark by Landauer 56 according to which, neither electric field nor carrier concentration can be strictly homogeneous in a metal containing defects. This is particularly true in our case, where the Thomas-Fermi screening length is very long by a combination of long Bohr radius and small Fermi wave vector. When the carrier density is as low as 10 18 cm −3 , Eq. (1) yields r TF ' 80 nm, much longer than the distance between the dipoles or the r c of the insulator. The Fermi sea is too dilute to impede dipolar interaction. The attenuation of mobility implies enhanced electric-field inhomogeneity brought by the alignment of dipoles, providing additional support for picturing persistent electric interaction in spite of a Fermi sea.
The second question leads us to an additional interaction mechanism offered to electric dipoles in presence of a Fermi sea (see Fig. 4c). This was proposed first by Glinchuk and Kondakova in 1992. 35 The opposite charges of an electric dipole are expected each to generate Friedel oscillations 57 inside a Fermi sea, which is a distinct source of dipolar interaction. Indeed, it is established that localized magnetic spins inside a Fermi sea generate similar oscillations and interact through what is known as the Ruderman-Kittel-Kasuya-Yosida (RKKY) mechanism. 58 Glinchuk and Kondakova proposed a dipolar analogue of RKKY interaction. This dipole-dipole interaction (V dd R ) depends on the magnitude of the electric dipole moment and has an alternating sign as a function of the inter-dipole distance (l dd ) and Fermi wave vector (k F ): 35 This expression for interaction between two dipoles is to be compared with what is expected in vacuum: (4) or in a highly polarizable insulator with a polarization correlation radius of r c : 50 According to Eq. (3), the parallel alignment of nearest dipoles becomes energetically unfavorable when cosð2k F l dd Þ ¼ À1 or 2k F l dd ¼ π. Assuming an isotropic Fermi surface (k F ¼ ð3π 2 nÞ 1=3 ), one finds that this will happen when n d n ¼ 24 π . In other words, in presence of a single-band isotropic metal, the anti-ferroelectric coupling between neighboring dipoles is expected to peak when there is 24 π ' 7:6 dipoles per mobile electrons. Thus, the simplest picture of RKKY-like interaction between electric dipoles expects a proportionality between dipole density and the charge carrier density required to achieve maximum antiferroelectric coupling between neighboring dipoles and this proportionality is simply close to 8 (which is because we are in three dimensions, and for each dimension one needs two dipoles per one electron for destructive interference). It is worth noting the following two facts. First of all, Eq. (3) does not include the angular dependence of the interaction between neighboring dipoles, which can have different orientations. Second: the dopants are randomly distributed in the real system, which should be taken account in a more sophisticated theoretical model. Thus, one may conjecture that the n d n ratio remains the same after averaging the space and angle dependence of interactions. If not, it would be an amazing coincidence that for three different Ca contents, the experimentally observed n d n corresponds to what is expected in the simple approach neglecting these details. Further theoretical efforts are required to clarify this issue.
One may be tempted by an alternative picture in which the FElike transition is destroyed, because the screening length shrinks with increasing carrier concentration. However, since this evolution is very slow (r TF / n À1=6 ), it would be very hard to explain the linear proportionality between n Ã and x. . a In an ordinary ionic crystal, the interaction is ferroelectric along one orientation and anti-ferroelectric along the perpendicular one, impeding the emergence of long-range ferroelectricity. b In a highly polarizable lattice, the interaction is ferroelectric over a length scale, r c , longer than the interdopant distance, l dd . This allows the formation of large ferroelectric cluster composed of many dipoles. The typical size of these clusters is set by r c . [49][50][51] c In the presence of a Fermi sea, ferroelectric and anti-ferroelectric coupling between dipoles alternate radially. According to our observations, anti-ferroelectric coupling between dipoles destroys the ferroelectric order in Sr 1Àx Ca x TiO 3Àδ at a carrier density of n Ã ðxÞ, which corresponds to 2k F l dd ¼ π.
As mentioned above, when carrier density is below n Ã ðxÞ, we are in presence of a polar metal where mobile electrons and aligned dipoles co-exist. Explaining the transport properties of the metal in this regime emerges as a challenge to future theoretical works. The complex temperature dependence of electrical resistivity (markedly different from quadratic) is to be understood in a picture where electrons are scattered off a hardened ferroelectric mode.
What happens to the electric dipoles when the carrier density exceeds n Ã ðxÞ? According to the thermal-expansion measurements, 36 a smeared anomaly continues to survive. One may speculate that a non-percolative dipole glass 50 persists at higher densities. In our study of resistivity, we cannot see any significant difference between Ca-free and Ca-substituted strontium titanate when n > n Ã ðxÞ. On the other hand, we do see small upturns in resistivity, reminiscent of the Kondo effect, which constitutes yet another source of information.
Low-temperature resistivity of Sr 1Àx Ca x TiO 3Àδ (with x ¼ 0:9%) above n Ã ðxÞ is shown in Fig. 5a-d. One can see that a less pronounced upturn re-appears above n Ã ðx ¼ 0:9%Þ. Similar data for Ca-free samples are shown in Fig. 5e-h for comparison. No clear upturn can be seen in Fig. 5e, f with n below 5:2 10 19 cm −3 . With further increase of carrier concentration, a small upturn appears. No difference between Ca-doped and Ca-free samples can be detected. We note, however, that these upturns in resistivity correspond to a relative change of several 10 À3 . This is orders of magnitudes smaller than what can be seen in Fig. 2 for n < n Ã . Their small magnitude as well as their presence in Ca-free samples point to an origin which is not the one seen in the presence of the phase transition and discussed above.
As seen in Fig. 5c, g, the temperature dependence of these upturns is roughly logarithmic. The logarithmic temperature dependence of resistivity is reminiscent of the Kondo effect arising when a localized spin couples to a Fermi sea. 59 It is now known that the Kondo effect can be present whenever a Two-Level-System (TLS) is embedded in a Fermi sea 60 and experiments 61 63 If this is a Kondo effect, what can be its origin? Oxygen vacancies are the principal suspect. By distorting the TiO 6 octahedra, they can generate degeneracies in different degrees of freedom. Moreover, they are suspected to host d 0 magnetism. 64 One cannot forget, however, the unavoidable presence of magnetic impurities at the ppm level in SrTiO 3 . 65,66 The absence of this upturn below 5 10 19 cm −3 and its gradual appearance at higher densities begs explanation. The evolution of the T-square resistivity prefactor may provide one. In presence of a large temperature-dependent resistivity, a small sub-percent upturn may be undetectable. However, this hypothesis fails to explain the gradual fading of the effect at even higher carrier concentrations. Why the upturn is most prominent when n ' 2 10 20 cm −3 , which roughly corresponds to one vacancy out of 100 oxygen atoms? A possible clue is provided by the wellestablished fact that in general, the Kondo physics of magnetic impurities above a threshold concentration is replaced by a spinglass state. In the case of iron impurities in silver, 67 this threshold is about a percent. The same may happen to the Kondo effect of oxygen vacancies. Note that this is yet another argument against extrinsic magnetic impurities as the source of the Kondo effect.
In summary, we found that the existence of the FE order in insulating Sr 1Àx Ca x TiO 3 deeply affects charge transport in dilute metallic Sr 1Àx Ca x TiO 3Àδ . Below a threshold of carrier concentration n Ã ðxÞ, inelastic and elastic scattering of carriers are both amplified. This threshold density scales with Ca content. Both these features can be explained by invoking the survival of electric dipole interaction in presence of mobile electrons.
Inside the Fermi sea, the dipolar counterpart of RKKY interaction can give rise to a strong anti-ferroelectric coupling at a carrier density proportional to the dipole density. There is striking match between the experimentally resolved proportionality (7.9 dipoles per electron) and what is expected in the simplest RKKY dipolar picture.
At higher carrier concentrations, deep inside the metallic state, clear but tiny upturns of resistivity are detectable, which are reminiscent of the Kondo effect and most probably associated with oxygen vacancies.

METHODS
SrTiO 3 and Sr 1Àx Ca x TiO 3 single-crystals with x ¼ 0:22%; 0:45%, and 0:9% were commercially obtained. Mobile electrons were introduced by generating oxygen vacancies. In order to remove oxygen atoms, samples were sealed with high-purity Ti powder (99:99%) in an evacuated quartz tube and annealed in an oven. Samples with various electron carrier concentrations (from 10 18 to 10 21 cm −3 ) were obtained by varying the annealing temperature between 700 and 900°C. Resistivity was measured through the four-terminal method down to 1.8 K in Quantum Design PPMS. Hall effect is measured concomitantly under magnetic field to determine the carrier density (n).

DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request. Fig. 5 Temperature dependence of resistivity at low temperatures above n Ã . Data for Sr 1Àx Ca x TiO 3Àδ with x ¼ 0:9% (a-d) is compared with x ¼ 0 (e-h). The vertical blue bars correspond to a value of 0.2 μΩ.cm. Note the presence of small upturns in both cases corresponds to a sub-percent change in resistivity.