# Anomalous Transport in a Superfluid Fluctuation Regime

###### Abstract

Motivated by a recent experiment in ultracold atoms [ S. Krinner et al., Proc. Natl. Acad. Sci. U.S.A 113, 8144 (2016)], we analyze transport of attractively interacting fermions through a one-dimensional wire near the superfluid transition. We show that in a ballistic regime where the conductance is quantized in the absence of interaction, the conductance is renormalized by superfluid fluctuations in reservoirs. In particular, the particle conductance is strongly enhanced and the plateau is blurred by emergent bosonic pair transport. For spin transport, in addition to the contact resistance the wire itself is resistive, leading to a suppression of the measured spin conductance. Our results are qualitatively consistent with the experimental observations.

###### pacs:

05.60.Gg, 67.85.-dTransport measurements often play crucial roles in revealing the fundamental nature of matter. In condensed matter physics, superconductivity, the Kondo effect, and the quantum Hall effect were all discovered with transport measurements. A two-terminal setup realized in ultracold atoms has opened up yet another avenue to explore strongly correlated systems through transport Brantut et al. (2012); Stadler et al. (2012); Krinner et al. (2013); Brantut et al. (2013); Krinner et al. (2015); Husmann et al. (2015); Krinner et al. (2016). In ultracold atoms, quantum transport typically occurs in the clean limit. The bulk conductivity cannot distinguish between different quantum states, since the -sum rule and the momentum conservation dictate that the conductivity involves the delta-function singularity at zero frequency whose weight does not depend on detailed states of matter Enss and Haussmann (2012). However, transport through a constriction such as a quantum point contact allows one to distinguish between different states due to the breakdown of the momentum conservation at the constriction. Indeed, such a setup has unveiled different transport properties for non-interacting Krinner et al. (2015) and superfluid fermions Husmann et al. (2015).

Recently, particle conductance and spin conductance have been measured with a quantum point contact in ultracold fermions Krinner et al. (2016). There, with increasing attractive interaction, both of them deviate significantly from quantized values Imry (2002) just above the superfluid critical temperature. More specifically, compared with the non-interacting limit, the particle conductance is enhanced, whereas the spin conductance is suppressed. These remarkable features stand in sharp contrast with the conventional wisdom that a conductance is not renormalized by an interaction in a one-dimensional wire Tarucha et al. (1995); Maslov and Stone (1995); Ponomarenko (1995); Safi and Schulz (1995); Kawabata (1996); Shimizu (1996). Meanwhile, different from condensed matter situations, fluctuation and interaction effects in reservoirs may be significant in cold atom experiments.

On another front, to cope with effects of an interaction at reservoirs in ballistic transport presents a theoretical challenge, since Landauer’s approach does not operate with an interaction, and phenomenological approaches Maslov and Stone (1995); Ponomarenko (1995); Safi and Schulz (1995); Kawabata (1996); Shimizu (1996) used to explain the interaction effect in a one-dimensional wire cannot directly answer the question. When the conductance of the wire is small, a tunneling Hamiltonian approach is widely used to investigate the effect of interactions in reservoirs on transport Mahan (2013); Larkin and Varlamov (2005); Varlamov and Dorin (1983); Levchenko (2010). However, to discuss the ballistic limit realized in Ref. Krinner et al. (2016), we must go beyond the linear response theory which has been widely used in tunneling experiments with correlated materials Fischer et al. (2007).

In this Letter, motivated by the ETH experiment Krinner et al. (2016) and the theoretical challenge mentioned above, we examine the effects of superfluid fluctuations in reservoirs on transport through a one-dimensional wire. To deal with ballistic transport, we apply the nonlinear response theory Wehrum and Hermeking (1974); Genenko and Ivanchenko (1986); Butcher and Cotter (1991); Rammer (2007); Kamenev (2011) to demonstrate that the breakdown of the quantization of conductance occurs by superfluid fluctuations. We show that transport of preformed pairs induced by superfluid fluctuations is essential to account for the breakdown. We also point out that in addition to the contact resistance, the resistance in the one-dimensional wire plays an important role in spin transport.

The Model.— We consider a system where two macroscopic reservoirs with superfluid fluctuations are connected by a quantum point contact (a one-dimensional channel). In the ETH experiment, the constriction has potential variations that take place over length scales larger than with the Fermi momentum . This implies that when transport near the Fermi energy is concerned, the adiabatic approximation is justified in which the detailed shape in the constriction is irrelevant Imry (2002); sup (). Thus, for the single channel case, we can start with the following Hamiltonian ():

(1) | |||

(2) | |||

(3) |

where is the fermionic annihilation (creation) operator with momentum and spin in reservoir , and the energy is measured from the chemical potential . In addition, is the operator in the real space with its argument () representing a position in the left (right) reservoir. The first and second terms on the right-hand side of Eq. (1) are the single-particle Hamiltonian and the interaction with an attractive coupling (), respectively, and describe the system with a broad-Feshbach resonance used in the ETH experiment Krinner et al. (2016). Below, we focus on the BCS regime ( with the -wave scattering length ) Haussmann et al. (2007); Bloch et al. (2008); Levinsen et al. (2015); sup () above the superfluid transition temperature . To discuss the case of a single conducting channel, we set sup (); Cuevas et al. (1996); Berthod and Giamarchi (2011), where near the Fermi energy the tunneling amplitude can be chosen to be a real constant without loss of generality, and () is the entrance (exit) point in the quantum point contact. In fact, this tunneling Hamiltonian can precisely reproduce the known universal conduction properties in the quantum point contact including ballistic transport with superfluid reservoirs Cuevas et al. (1996); Husmann et al. (2015).

In terms of Eq. (3), the mass and spin current operators are given by , and , where is the number operator with spin in the reservoir . We note that the number operator in each reservoir commutes with the Hamiltonian except for .

In the presence of a chemical-potential difference between the reservoirs, (), the mass (spin) current is induced. Then, the averages of the mass and spin currents at time are given by

(4) |

where , and means the thermal average for the Hamiltonian (1).

In the presence of superfluid fluctuations, we should consider contributions arising from fermionic quasiparticles and fluctuation pairs Larkin and Varlamov (2005). Below, such fluctuations are considered up to the Gaussian level in each propagator, which is reasonable in a regime for the case of three-dimensional reservoirs Larkin and Varlamov (2005) relevant to the ETH experiment.

Fermionic quasiparticle current.—We now examine a steady current induced by fermionic quasiparticles. By the assumption of the steady state, we can put without loss of generality. Then, the mass and spin currents can be expressed as where is the Keldysh Green’s function Kamenev (2011), and we use for the retarded Green’s function Rammer (2007); Kamenev (2011); Berthod and Giamarchi (2011).

As shown in Fig. 1 (b), the single-particle Green’s function is renormalized by the fluctuation pair propagator (Fig. 1 (a)) up to the Gaussian level. As pointed out in Ref. Cuevas et al. (1996), the effect of must be incorporated to all orders in the ballistic limit. By using an analysis similar to the case of noninteracting fermions Genenko and Ivanchenko (1986); Caroli et al. (1971); Cuevas et al. (1996); Berthod and Giamarchi (2011) the fermionic quasiparticle contribution to the conductance per spin is obtained as sup ()

(5) |

where is the Fermi distribution function at temperature , with the retarded Green’s function in the reservoirs, and is the density of state (DOS). We note that the conductance depends neither on nor sup (). To obtain the above result, an expansion up to linear order in is considered, since and no significant deviation from the linear order is found in the ETH experiment. We note that the same expression holds for the mass and spin currents. In the case of small transmittance where in the denominator, Eq. (5) essentially reduces to the Ambegaokar-Baratoff formula Ambegaokar and Baratoff (1963); Larkin and Varlamov (2005). On the other hand, in the absence of the fluctuations, the conductance is reduced to and is equivalent to Landauer’s formula with transmittance , where is the DOS for noninteracting fermions and we use the fact that the change of around the Fermi level is much smaller than that of Genenko and Ivanchenko (1986); Caroli et al. (1971); Cuevas et al. (1996). In this limit, the quantized conductance, , is obtained for .

The superfluid fluctuations renormalize the conductance of fermionic quasiparticles and generate that of preformed pairs. The fermionic quasiparticle conductance, in general, tends to be suppressed due to pseudogap effect Jankó et al. (1997); Tsuchiya et al. (2009); however, in the case of three-dimensional reservoirs, such suppression is found to be negligible in the experimentally relevant regime sup ().

Fluctuation pair current.—We now consider a current carried by the fluctuating (preformed) pairs that makes a dominant contribution to the conductivity in dirty superconductors Larkin and Varlamov (2005). As shown on the left of Fig. 1 (c), the lowest-order diagram already contains a factor . In usual tunneling experiments where Mahan (2013); Larkin and Varlamov (2005); Varlamov and Dorin (1983); Levchenko (2010); Fischer et al. (2007), this contribution is negligible compared with the fermionic quasiparticle current, and has not been considered for realistic situations. However, in the ballistic regime in which , one needs to consider it seriously. To evaluate the fourth-order diagram, we calculate the third-order response function, which is related to an imaginary-time correlation function through analytic continuation Mahan (2013); Wehrum and Hermeking (1974). Up to linear order in , the fluctuation pair current at the fourth order in behaves as . Here, the factor originates from the pair exchange between the reservoirs, and the factor reflects the superfluid fluctuations. We note that as in the case of spin conductivity Enss and Haussmann (2012), the fluctuation pair contribution in the spin current vanishes. This reflects the fact that the pair exchange is not caused by a spin bias. Thus, the enhancement of the current by fluctuation pairs only occurs for mass transport.

We also note that the left-hand side of
Fig. 1 (c)
can be replaced by the right-hand side of Fig. 1 (c) up to linear order in .
Namely, the fluctuation pair contribution can be expressed
in terms of an effective hopping amplitude
^{1}^{1}1Notice that the dimensional of is different from that of , since
the dimension of is different from that of . and
the retarded pair-fluctuation propagator
whose expression in the vicinity of
is given by Larkin and Varlamov (2005)

(6) |

where and with the Fermi velocity and the dimension of the system . Thus, this contribution can be calculated as tunneling of the preformed pairs for a given effective hopping amplitude and bias , and therefore the multiple tunneling processes of the preformed pairs represented by power series in can be systematically evaluated as depicted in Fig. 1 (d). By using the nonlinear response theory, we obtain the conductance contributed from the fluctuation pair per spin up to linear order in as sup ()

We note that the above formula indeed reflects bosonic transport, since the term in the integrand is the derivative of the Bose distribution function with respect to (note that such a term is absent in Eq. (5)). Thus, fluctuation pairs make a positive contribution to the mass conductance. Such a contribution is already visible in the regime sup ().

Conductances in the single-channel regime.—We now compare our theory with the ETH experiment in the ballistic single-channel regime. To this end, one may also consider an interaction effect inside the wire.

For mass transport, the mass current operator commutes with the bulk Hamiltonian containing an interaction in the wire Enss and Haussmann (2012); sup () (except, of course, for the tunneling term), the wire resistance is expected to be negligible, and the conductance calculation obtained above is directly applicable to the ETH experiment. An essential input parameter in the theory is the ratio with the Fermi temperature , which is extracted from the experiment Krinner et al. (2016); sup (). By assigning the ballistic limit , we compare the theory with the experiment and find excellent agreement as shown in Fig. 2 (a). A crucial point here is that the conductance is enhanced due to the bosonic fluctuation-pair contribution. Since our theory is based on an expansion from , some deviation is expected at .

For spin transport, the spin current operator and the Hamiltonian do not commute even in the absence of the tunneling term, giving rise to the wire resistance Sommer et al. (2011a, b); Duine et al. (2010); Bruun (2011); Bruun and Pethick (2011); Enss and Haussmann (2012). In the presence of an attractive interaction, a spin gap, shows up. A typical estimation suggests , where the lower bound is estimated with the Yang-Gaudin model at the density and the upper bound is determined from the binding energy of the confinement-induced resonance , where is the transverse confinement frequency Olshanii (1998); Fuchs et al. (2004). In Fig. 2 (b), we show the spin conductance whose resistance is the sum of the contact resistance and wire resistance. The wire resistance is estimated so as to be compatible with the experiment by assuming Giamarchi (2004). Our result shows that the wire resistance for spin transport is of the order of the contact resistance, implying that even in the ballistic limit a nonnegligible chemical potential drop occurs inside the wire due to the interaction between and spin components.

We also comment on an effect of the spin gap near the contacts. The spin gap in the wire originates from the strong nesting effect allowed in a one-dimension system Giamarchi (2004). On the other hand, near the contact, multiple channels that render the spin gap smeared out through the dimensional crossover are present. Thus, the effect of the spin gap near the contacts is expected to make negligible contributions to the contact resistance.

Effects of the gate potential, trapping, and interaction on particle conductance—Now, we discuss how the particle conductance is affected by the gate potential, trap potential, and interaction. Since the gate and trap potentials shift the energy levels of the conducting channels, these effects can be incorporated by introducing multiple tunneling amplitudes, each of which depends on the gate and trap potentials sup (); Glazman et al. (1988); Büttiker (1990); Krinner et al. (2016). The tunneling amplitudes as the input parameters are determined so as to reproduce the weakest-interaction data in the experiment based on Landauer’s formula for noninteracting fermions, since there, and the superfluid fluctuations are expected to be minuscule sup (). On the other hand, the interaction strength is directly related to how close the system is to , since increasing towards the unitarity leads to an enhancement of Krinner et al. (2016); sup ().

Figure 3 compares the results of our theory with the ETH experiments for different interaction strengths ( in Fig. 3 (a) and in Fig. 3 (b) Krinner et al. (2016)). For the weakest interaction , the grey curves are obtained by fittings of the experimental data by assuming Landauer’s formula. For stronger interaction strengths, where , we use our theory by incorporating the superfluid fluctuations to calculate the particle conductance by assigning tunneling amplitudes estimated from the data at . As shown in Fig. 3, the particle conductance is enhanced and deviates from Landauer’s formula, which is consistent with the experiment Krinner et al. (2016). As in the case of the single-channel regime, enhancement is caused by the preformed pairs. The discrepancy between the theory and the experiment occurring at the larger gate potential in Fig. 3 (b) may be due to an effect of the higher transverse channels that are not treated in the theory.

Summary.—We have shown that superfluid fluctuations cause two competing effects in two-terminal transport through a quantum point contact; suppression of the conductance of fermionic quasiparticles and enhancement due to bosonic preformed pairs. The former is negligible in the ETH experiment, since the depletion of the DOS near the Fermi level is negligible. The latter in the ballistic regime is shown to be significant due to the absence of the Pauli exclusion principle. Thus, the net conductance can exceed the upper bound of the quantized conductance for noninteracting fermions through multiple tunneling processes that are captured with the nonlinear response theory. Such transport is ideally realized in an impurity-free system with perfect transmittance such as cold atoms and high-mobility semiconductors in which diffusive properties in a one-dimensional wire, which tends to suppress the bosonic transport, can be ignored. We have also shown that spin transport is affected by the wire resistance originating from the spin gap whose determination with no ambiguity requires the more precise knowledge of the particle density in the wire.

Acknowledgement.—We thanks Jean-Philippe Brantut, Tilman Esslinger, Sebastian Krinner, Martin Lebrat, and Naoto Tsuji for useful discussions. This work is supported by KAKENHI Grant Nos. JP26287088 and JP15H05855from the Japan Society for the Promotion of Science, the Photon Frontier Network Program from MEXT of Japan, and Grant for Basic Science Research Projects from The Sumitomo Foundation.

Note added.—Recently, we became aware of works by M. Kanász-Nagy et al. Kanász-Nagy et al. (2016) and B. Liu et al. Liu et al. (2016), both of which discuss the anomalous conductances from different perspectives.

## References

- Brantut et al. (2012) J.-P. Brantut, J. Meineke, D. Stadler, S. Krinner, and T. Esslinger, Science 337, 1069 (2012).
- Stadler et al. (2012) D. Stadler, S. Krinner, J. Meineke, J.-P. Brantut, and T. Esslinger, Nature 491, 736 (2012).
- Krinner et al. (2013) S. Krinner, D. Stadler, J. Meineke, J.-P. Brantut, and T. Esslinger, Phys. Rev. Lett. 110, 100601 (2013).
- Brantut et al. (2013) J.-P. Brantut, C. Grenier, J. Meineke, D. Stadler, S. Krinner, C. Kollath, T. Esslinger, and A. Georges, Science 342, 713 (2013).
- Krinner et al. (2015) S. Krinner, D. Stadler, D. Husmann, J.-P. Brantut, and T. Esslinger, Nature 517, 64 (2015).
- Husmann et al. (2015) D. Husmann, S. Uchino, S. Krinner, M. Lebrat, T. Giamarchi, T. Esslinger, and J.-P. Brantut, Science 350, 1498 (2015).
- Krinner et al. (2016) S. Krinner, M. Lebrat, D. Husmann, C. Grenier, J.-P. Brantut, and T. Esslinger, Proc. Natl. Acad. Sci. U.S.A. 113, 201601812 (2016).
- Enss and Haussmann (2012) T. Enss and R. Haussmann, Phys. Rev. Lett. 109, 195303 (2012).
- Imry (2002) Y. Imry, Introduction to mesoscopic physics (Oxford University Press, 2002).
- Tarucha et al. (1995) S. Tarucha, T. Honda, and T. Saku, Solid state communications 94, 413 (1995).
- Maslov and Stone (1995) D. L. Maslov and M. Stone, Physical Review B 52, R5539 (1995).
- Ponomarenko (1995) V. Ponomarenko, Physical Review B 52, R8666 (1995).
- Safi and Schulz (1995) I. Safi and H. Schulz, Physical Review B 52, R17040 (1995).
- Kawabata (1996) A. Kawabata, Journal of the Physical Society of Japan 65, 30 (1996).
- Shimizu (1996) A. Shimizu, Journal of the Physical Society of Japan 65, 1162 (1996).
- Mahan (2013) G. D. Mahan, Many-particle physics (Springer Science & Business Media, 2013).
- Larkin and Varlamov (2005) A. Larkin and A. Varlamov, Theory of fluctuations in superconductors (Oxford University Press, 2005).
- Varlamov and Dorin (1983) A. Varlamov and V. Dorin, Sov. Phys. JETP 57, 1089 (1983).
- Levchenko (2010) A. Levchenko, Phys. Rev. B 81, 012501 (2010).
- Fischer et al. (2007) O. Fischer, M. Kugler, I. Maggio-Aprile, C. Berthod, and C. Renner, Rev. Mod. Phys. 79, 353 (2007).
- Wehrum and Hermeking (1974) R. Wehrum and H. Hermeking, Journal of Physics C: Solid State Physics 7, L107 (1974).
- Genenko and Ivanchenko (1986) Y. A. Genenko and Y. M. Ivanchenko, Theoretical and Mathematical Physics 69, 1056 (1986).
- Butcher and Cotter (1991) P. N. Butcher and D. Cotter, The elements of nonlinear optics, Vol. 9 (Cambridge University press, 1991).
- Rammer (2007) J. Rammer, Quantum field theory of non-equilibrium states (Cambridge University Press, 2007).
- Kamenev (2011) A. Kamenev, Field theory of non-equilibrium systems (Cambridge University Press, 2011).
- (26) See Supplemental Material for the derivation and further discussions .
- Haussmann et al. (2007) R. Haussmann, W. Rantner, S. Cerrito, and W. Zwerger, Phys. Rev. A 75, 023610 (2007).
- Bloch et al. (2008) I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008).
- Levinsen et al. (2015) J. Levinsen, M. Parish, et al., Annual Review of Cold Atoms and Molecules: Volume 3 3, 1 (2015).
- Cuevas et al. (1996) J. C. Cuevas, A. Martín-Rodero, and A. L. Yeyati, Phys. Rev. B 54, 7366 (1996).
- Berthod and Giamarchi (2011) C. Berthod and T. Giamarchi, Phys. Rev. B 84, 155414 (2011).
- Caroli et al. (1971) C. Caroli, R. Combescot, P. Nozieres, and D. Saint-James, Journal of Physics C: Solid State Physics 4, 916 (1971).
- Ambegaokar and Baratoff (1963) V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 486 (1963).
- Jankó et al. (1997) B. Jankó, J. Maly, and K. Levin, Phys. Rev. B 56, R11407 (1997).
- Tsuchiya et al. (2009) S. Tsuchiya, R. Watanabe, and Y. Ohashi, Phys. Rev. A 80, 033613 (2009).
- (36) Notice that the dimensional of is different from that of , since the dimension of is different from that of .
- Sommer et al. (2011a) A. Sommer, M. Ku, G. Roati, and M. W. Zwierlein, Nature 472, 201 (2011a).
- Sommer et al. (2011b) A. Sommer, M. Ku, and M. W. Zwierlein, New Journal of Physics 13, 055009 (2011b).
- Duine et al. (2010) R. A. Duine, M. Polini, H. T. C. Stoof, and G. Vignale, Phys. Rev. Lett. 104, 220403 (2010).
- Bruun (2011) G. M. Bruun, New Journal of Physics 13, 035005 (2011).
- Bruun and Pethick (2011) G. M. Bruun and C. J. Pethick, Phys. Rev. Lett. 107, 255302 (2011).
- Olshanii (1998) M. Olshanii, Phys. Rev. Lett. 81, 938 (1998).
- Fuchs et al. (2004) J. N. Fuchs, A. Recati, and W. Zwerger, Phys. Rev. Lett. 93, 090408 (2004).
- Giamarchi (2004) T. Giamarchi, Quantum physics in one dimension (Oxford university press, 2004).
- Glazman et al. (1988) L. Glazman, G. Lesovik, D. Khmel’Nitskii, and R. Shekhter, Soviet Journal of Experimental and Theoretical Physics Letters 48, 238 (1988).
- Büttiker (1990) M. Büttiker, Phys. Rev. B 41, 7906 (1990).
- Kanász-Nagy et al. (2016) M. Kanász-Nagy, L. Glazman, T. Esslinger, and E. A. Demler, Phys. Rev. Lett. 117, 255302 (2016).
- Liu et al. (2016) B. Liu, H. Zhai, and S. Zhang, Phys. Rev. A 95, 013623 (2016) .