Bosonic thermoelectric transport and breakdown of universality

In this paper we compare Bose transport in normal phase atomic gases with its counterpart in Fermi gases, illustrating the non-universality of two dimensional bosonic transport associated with different dissipation mechanisms. Near the superfluid transition temperature $T_c$, a striking similarity between the fermionic and bosonic transport emerges because super-conducting(fluid) fluctuation transport for Fermi gases is dominated by the bosonic, Cooper pair component. As in fluctuation theory, one finds that the Seebeck coefficient changes sign at $T_c$ and the Lorenz number approaches zero at $T_c$. Our findings appear semi-quantitatively consistent with recent Bose gas experiments.

We discuss the general principles of transport in normal phase atomic gases, comparing Bose and Fermi systems. Our study shows that two dimensional bosonic transport is non-universal with respect to different dissipation mechanisms. Near the superfluid transition temperature Tc, a striking similarity between the fermionic and bosonic transport emerges because super-conducting (fluid) fluctuation transport for Fermi gases is dominated by the bosonic, Cooper pair component. As in fluctuation theory, one finds that the Seebeck coefficient changes sign at Tc and the Lorenz number approaches zero at Tc. Our findings appear quantitatively consistent with recent Bose gas experiments. Cold atom samples differ from electronic systems in significant ways that offer new opportunities to investigate transport phenomena [1][2][3][4][5]. Cold atom systems are highly versatile: many of them have tunable interactions, can be confined in a single potential well or in lattices, as well as in different dimensions. Furthermore, at nanoKelvin temperatures, the dynamics of atoms is very slow (milliseconds), in contrast to the fast dynamics (picoseconds) of electrons in materials. This allows a detailed scrutiny of atomic motion. One powerful tool in recent cold atom experiments is in situ imaging of atoms, which reveals high space-time resolution images of atomic distributions in snapshots [6][7][8]. In the condensed matter analogy, this technique is equivalent to following every electron with a femtosecond temporal resolution.
In further pursuit of understanding both these dynamics and the analogy between electrons and cold atoms, in this paper we apply the theory of dissipative transport (developed for electrons) to ultracold trapped atoms. In the last section of the paper we address recent experimental data [1].
Dissipative transport is to be contrasted with ballistic transport where there are no relevant scattering processes to limit the transport lifetime. Thus, ballistic transport is observable when the mean free path is long compared to the relevant dimensions of the system. We stress that the analogy between electrons and cold atoms requires careful consideration. Universal transport laws developed for electrons need to be revisited when applied to bosonic atoms. In addition, cold atom samples are trapped in a conservative potential and isolated in vacuum. The lack of thermal and particle reservoirs means total particle and total energy are usually conserved quantities. A statistical description based on a grand-canonical ensemble may apply only to local observables [9].
Furthermore, atoms are neutral and the analogue of electrical current will be the particle flow or mass flow driven, not by electric field, but by a chemical potential gradient ∇µ or temperature gradient ∇T . Interestingly, unlike neutral liquid Helium, cold gas superfluids allow the imposition of a non-zero ∇µ and thus are rather uniquely amenable to the transport studies we present. Finally, to employ dissipative, as distinguished from ballistic transport theory, the samples should be in the hydrodynamic regime with coherence length much shorter than the sample and relaxation time less than the measurement time resolution. When these criteria are satisfied the atoms can reach a local equilibrium associated with coarse graining the system over a proper length and time scale. Thermodynamic quantities such as temperature and chemical potential can therefore be defined locally, and their gradients constitute the thermodynamic forces. This approximation, which we will call the local equilibrium approximation (LEA), limits the resolution of both temporal and spatial measurements. It should not be confused with the local density approximation that describes inhomogeneously trapped atoms in equilibrium.
Since cold atoms are typically confined in conservative potentials, they are usually free from impurities and background ionic lattices, including their phononic excitations. The fact that atoms are in a clean environment, on the one hand, may simplify many-body calculations; this removes some complexity encountered in electron transport in materials. On the other hand, there must necessarily be a relaxation mechanism in order to achieve a steady state and apply the fundamentals of transport theory.
In this paper, we discuss the transport behavior of ultracold gases in the normal phase under the assumption of the LEA, the applicability of which is described in more detail below. We characterize transport in two limits where theory is relatively straightforward and well established. Thus we consider non-interacting particles as well as (generally interacting) particles but near su-perfluid condensation, emphasizing the similarities and differences between Fermi and Bose systems. When the interactions are essentially negligible, a normal degenerate Fermi gas in the dissipative transport regime has universal thermoelectric properties. By contrast universality is absent for the bosonic counterpart. As the system approaches the superfluid critical temperature, the transport tends to be dominated by condensate (superconducting and superfluid) fluctuations in which regime, the Fermi and Bose gases are more similar. It is in this latter regime where there are now (Bose gas) experiments [1]. Importantly these appear to be reasonably consistent with the theory presented here.

A. Transport Equations
We define the transport coefficients in terms of the linear response of the particle current J p and the heat current J Q to temperature gradient ∇T and and chemical potential gradient ∇µ as where ∇µ is analogous to the electric field for a charged system. It is convenient to introduce dimensionless ratios to write this equation in the form (We work with units such that = k B = 1.) where the conductivity σ ≡ L 11 , the Seebeck coefficient (or thermopower) S = P/T ≡ L 12 /L 11 and the Lorenz number (which is given by ratio of thermal to mass conductivity, divided by temperature) L ≡ (L 22 L 11 − L 12 L 21 )/T L 2 11 . Note that L 12 and L 21 satisfy the Onsager relation L 12 = L 21 /T .

B. Kubo versus Boltzmann Approaches
There are two rather complementary approaches to addressing dissipative transport in quantum systems. One can adopt the approach of linear response theory and apply the Kubo formula, or alternatively one can employ the Boltzmann transport equation which is based on kinetic theory. Each has its strengths in accomodating different physical aspects of the system at hand and the choice for which to apply is governed by the underlying goals. The Kubo formalism is more suitable for imposing conservation principles and sum rules, but it does not generally build in the specifics of the dissipation mechanism. By contrast, the Boltzmann transport equation focuses on the details of the collision processes leading to dissipation, but is not as well suited to impose or to verify conservation principles. For non-interacting gases treated in the Kubo formalism, the results are essentially identical to those based on Boltzmann transport equation at the level of simple relaxation time approximations.
In this paper, we show that Kubo based approaches are well suited to addressing transport of both noninteracting as well as interacting atomic gases in the near condensation regime. Since it is likely that dissipation in the ultracold gases is linked to the details of the experimental set up, we will introduce dissipation via a phenomenological parameterization. The Kubo approach in this form is also suitable for delivering one of the messages of this paper: that in bosonic systems (despite strong evidence for universality in many contexts), transport is highly non-universal and depends on the details of the dissipation. The philosophy behind our phenomenological approach to dissipation is rather similar to that articulated by Kadanoff and Martin who in a series of papers emphasized the importance of the Kubo-based correlation functions and their symmetries [10]. In related work on superfluids [11], they argued for the suitability of introducing a phenomenological parameterization of the lifetimes associated with transport.
In this paper we parameterize the relaxation time. We introduce the quantity τ (ǫ) where ǫ effectively represents the energy; presumably at low temperatures ǫ → 0. This is treated phenomenologically as τ (ǫ) = τ 0 ǫ η/2 . While η appears to be arbitrary, in this paper we consider the two different values chosen because they have been addressed in the solid state literature: η = −2 and η = 0. The latter is associated with impurity scattering models [12], while the former is more naturally associated with strictly bosonic transport as, for example in superconducting fluctuation theories of transport [13,14]. Then, reflecting the odd energy or frequency dependene of the boson self energy, to leading order (in small ǫ) this constrains τ (ǫ) and thus η = −2.

C. Role of Dissipation
Equilibration is necessary but not sufficient for establishing steady state transport. A relaxation mechanism is needed for equilibration. Such a mechanism is also central to many non-equilibrium studies, involving for example, interaction quenches in the ultracold gases. This was demonstrated [15] by our group and in recent experimental studies [16,17]. It is, however, important to stress that mass transport, for example, as in the particle conductivity, cannot exclusively depend on inter-particle collisions as total momentum is conserved in the presence of Galilean invariance. A well behaved transport requires that there be a source of momentum relaxation.
In cold atom experiments dissipation depends on the experimental design, and will arise in a variety of ways. It will be present in thin channel transport and in experiments with optical lattices or speckle potentials. In addition, there are systematic [18] and intrinsic [19] fluctuations of the optical potentials that contribute to re-laxation. The presence of a harmonic trap alone, which breaks Galilean invariance, is insufficient to avoid an unphysical infinite mass conductivity if there were otherwise no source of relaxation [20].
A recent experiment by one of us [1] has measured the thermoelectric transport coefficients using three-body recombination loss and heating. These processes lead to locally measurable mass and energy flow or currents. This is in contrast to condensed matter experiments where the temperature and chemical potential (or scalar EM potential) are externally applied. With simple models one can quantify the associated ∇µ and ∇T and thereby deduce transport coefficient ratios. This paper aims to address transport more generally applying to a variety of different experiments, not just to the three-body recombination loss mechanism. However, it should be stressed that, for this particular transport scheme, the LEA will be demonstrated to be appropriate. Most importantly, along with particle and energy loss, momentum relaxation is present, through the loss of particles.

D. Local Equilibrium Approximation
Indeed, a central assumption of this paper is the applicability of the local equilibrium approximation. In dissipative transport theory developed for electrons in solids it is assumed that in the presence of external fields or perturbations, a steady state can be achieved. This is possible when the atoms are in the hydrodynamic regime.
Here the coherence length l c is significantly less than the system size l s and thermal relaxation time τ is much shorter than the system lifetime τ s . One can consider the local chemical potential gradients ∇µ(x; t) and temperature gradients ∇T (x; t) which are essential quantities in transport theory. That the dynamics of an atomic sample can be completely described by local thermodynamic variables relies on the rapid and short-ranged scattering by an external potential, which leads to relaxation and serves to establish a steady state. Note that ∇µ(x; t) and ∇T (x; t) can only be defined by coarse graining local variables over a proper length scale l m and time scale τ m . This leads to the constraint l c < l m ≪ l s and τ < τ m ≪ τ s where the subscript m represents the experimental or measurement variables. When these inequalities hold one is in the regime of validity of the Local Equilibrium Approximation. This LEA is a stronger condition than a hydrodynamic approximation since it assumes a history-independent local and short-time thermodynamic equilibrium. The LEA, which we assume to be valid throughout this study, will break down for atoms in the collisionless regime or in a smooth trapping potential, where we have l c ≫ l s and τ m ≫ τ s , respectively.
In support of our use of the LEA is the observation that the density profiles are well described by an equilibrium equation of state [6]. Also important is the fact that measurements of the coherence times and correlation lengths (which are discussed in Section V) are consistent with the constraints, as outlined above.

II. TRANSPORT FOR NON-INTERACTING BOSONS AND FERMIONS
The transport properties of weakly interacting, normal Bose and Fermi gases are strikingly different. In a degenerate Fermi gas, it is well known that the phase space contributions to transport are confined to a narrow energy range around the Fermi energy. By contrast, for Bose gases, there is no such constraint. As a consequence, the magnitude of transport coefficients L ij for bosons tends to be much larger and much more sensitive to the detailed assumptions about the nature of the dissipation.
For a Hamiltonian with only one-body terms, we are able to compare the transport properties of bosons and fermions using an exact expression based on the Kubo formula first derived for fermions in Ref. 21 but readily generalized to bosons, ± e ω/T ± 1) −1 the Fermi/Bose distribution, and the fugacity z ± is defined as z + = e µ/T for fermions, and z − = e (µ−µc)/T for bosons below critical chemical potential µ < µ c . We introduce ρ(ǫ) = (m/2π) d/2 ǫ d/2−1 /Γ(d/2) the density of states for free particles of mass m in d dimensions (we use units such that = k B = 1). In order to describe a bosonic system at positive chemical potential, one has to include the (non-dissipative) inter-boson interactions, so that the chemical potential term is given by µ − µ c . See for instance [22] for an estimate of µ c in weak coupling.
The relaxation time introduced earlier appears as τ (ǫ) in Eq. (4). This parameter plays the role of the relaxation time in Boltzmann theories of transport [40].
The transport coefficients for non-interacting fermions/bosons can be derived from Eq. (4) as where ζ = (d + η)/2, λ dB = 2π mT is the de Brogile length, and are the polylogarithm functions. For degenerate fermions z + ≫ 1, the factor b (1) + is strongly peaked at the Fermi energy. Thus the energy dependence of τ is not relevant and the Lorenz number in this regime is a universal number L f = π 2 /3 (independent of τ and η), known as the Wiedermann-Franz law. By contrast, for bosons all energies ǫ contribute to the integrals, so that L depends on the details of the dissipation and has a non-universal behavior. Figure 1 presents a comparison of the behavior of fermionic and bosonic transport for one particular transport coefficient ratio (the Lorenz number) L, for d = 2 as a function of the log of the fugacity. The figure illustrates how different values of the dissipation exponents (η = −2 and η = 0) influence the transport. One sees that for very negative ln(z), the transport is independent of statistics and also sensitive to the details of the dissipation. Each component of L is larger for the bosonic case, but this effect is not apparent when plotted as a ratio. One sees that as the fermions cross over to the large positive fugacities, the Lorenz number approaches the universal Wiedermann-Franz law.
While the universality in the fermionic case is evident, this is clearly not the case for bosons. For Bose systems, we observe that even in the quantum regime the behavior depends on the dissipation mechanism, implying a nonuniversality of the transport.

III. SCALE INVARIANCE OF TWO-DIMENSIONAL BOSONIC TRANSPORT
Although, experimentally it has not been tested in Ref. 1, here we argue that an scale invariance should be observed in transport in a two-dimensional Bose gas [41]. This scaling arises from physics similar to that found in the thermodynamics of two-dimensional dilute Bose gases [16,23]. This scale invariance reflects the µ = T = 0 quantum phase transition (QPT) between the vacuum and superfluid phase [24]. This zero-temperature phase transition must not be confused with the finitetemperature phase transition between the normal and superfluid phase discussed in the rest of the paper. The presence of this QCP implies that a thermodynamic property such as the pressure is a universal function of the form Here P(x, y) is a universal function (independent of the microscopic details of the system) andg(x) is a renormalized interaction depending logarithmically on x, which for weakly-interacting bosons can be taken as constant, which we assume in the following. (Ifg(T ) varies significantly, the scale invariance is then lost.) Under these circumstances there is scale invariance, so that all thermodynamic functions in the dilute two-dimensional Bose gas depend only on µ/T once the interaction strength has been fixed. The scaling of the pressure is valid in the critical regime close to the QCP, that is µma 2 , T ma 2 ≪ 1 [25] (a is the s-wave scattering length). For dilute gases, this corresponds to the whole T -µ plane in the range of µ and T relevant for the experiments. From these results at equilibrium, we can infer that due to this same QCP, the transport coefficients also obey scaling relations. Following Ref. [26], one shows that close to a critical point, the conductivity scale as σ(µ, T, g) = s 2−d σ µs 1/ν , T s z ,g(s) , where S d (x, y) = σ(x, 1, y) is a universal function. This reasoning is easily generalized to the other transport coefficients and explicitly in dimension two with A d , K d and L ij universal functions. Note that the scaling arguments are in principle valid only for the singular part of the transport coefficient, and one could expect the presence of a regular part, which would not scale accordingly. However, since all four transport coefficients vanish for all µ ≤ 0 at T = 0 (the system is empty), we can infer that the regular parts are identically zero in the vicinity of the QCP (i.e. for µma 2 , T ma 2 ≪ 1). The normal-superfluid phase transition at finite T (as distinguished from the above vacuum-superfluid zerotemperature phase transition) is characterized by a singularity of the pressure, which translates into a singularity of P(x, y) at a given x c (y), where P(x, y) is universal. One finds that the critical chemical potential (at a given temperature) is given by µ c (T ) = T M g , where M(x) is also universal. Building on the universality of µ c (T )/T , we arrive at Here, we have shown explicitly the dependence of the transport coefficients on the dissipation mechanism through the dissipative coefficient η introduced earlier.
In particular, this scaling holds far from the normalsuperfluid critical regime [42]. The dependence of the scaling functions on η, in turn, implies that even though for a given kind of dissipation (defined by η) the transport coefficients of a dilute Bose gas are universal (i.e. described by a function L ij ) in the whole µ-T plane, it will be different for different kinds of dissipation mechanisms (i.e. different power laws), thus defining different universality classes [43]. It follows directly from Eq. (12) that the Seebeck coefficient (or thermopower) S and the Lorenz ratio L are also scale invariant.

A. Fluctuation theory
In Section II we discussed transport in non-interacting normal gases. In this context we have seen that there is a dissipative constant parameterized by η which is generally unconstrained. Moreover η played a central role in determining the transport coefficients of Bose gases, leading to highly non-universal behavior. We turn now to Bose and Fermi systems which are normal but near condensation (Bose-Einstein condensation for bosons and the superconducting transition for fermions). This focus stems from recent Bose gas experiments [1]. For the Fermi case, interactions are essential and one cannot apply the theory of Section II. For the Bose case, the near-condensation regime provides constraints on the character of dissipation. In this regime we can appeal to theories of fluctuation dynamics (such as time dependent Landau Ginsburg theory (TDGL)) to constrain the time dependences and hence frequency dependences of the dissipation. This, in turn, serves to constrain η, as these dynamical theories contain linear time derivatives leading to lifetimes which vary linearly with frequency or energy. Relating to the discussion in Section II, interestingly one arrives at the same limiting behavior of the Bose transport coefficients using either the fluctuation scheme or Eq. (4) with the equivalent η = −2.
Our focus on normal state transport in the ultracold gases near their condensation temperatures stems both from recent Bose gas experiments [1] and also from the observation that Bose and Fermi superfluids have rather similar properties in this regime. Both are dominated by a bosonic condensate fluctuation contribution. In this context bosonic transport contributions in superconductors (in the narrow critical regime) have been of great interest to the condensed matter community over many decades [13,14]. Here one attributes the transport coefficients to fermionic pairs ("composite" bosons) in low but non-zero momentum states. These pre-formed pairs yield a greatly enhanced transport in the normal state, as compared with the behavior deduced from weakly interacting fermion theories such as those in the previous section. More recently there has been a focus on the normal state of high temperature superconductors where it appears that condensate fluctuation transport may set in at much higher temperatures T * well above the transition at T c . This reflects the so-called "pseudogap" physics [27] and this thought to be related to pre-formed pairs. We [27,28] and others have argued that this pseudogap should be relevant to Fermi gases at unitarity, where pairing is strong, but the system is still fermionic. Thus, if indeed, there is a pseudogap at unitarity, the behavior of transport in the Fermi gases may serve to reveal it in future experiments.
Transport properties associated with condensate fluctuations are based on two assumptions. (i) The fluctuations (often called "Gaussian fluctuations") represent relatively independent bosons. At most one introduces Hartree correlations, and (ii) The dynamics of the fluctuations incorporates linear time derivatives. Higher order derivatives are less important due to critical slowing down effects. This, importantly, has implications for the frequency dependence of the dissipation which ultimately affects the behavior of transport.
With these assumptions, the transport calculations are straightforward and one can evaluate [15,29] the various transport coefficients following Gaussian fluctuation models [13,30,31]. Entering are two factors of the bosonic spectral function, called A ≡ A(k, ω): Here, b (1) ≡ − ∂ ∂ω b(ω), with b(ω) = (e ω/T − 1) −1 the Bose distribution. We characterize the spectral function in terms of the bosonic propagator so that where the critical chemical potential µ c is a phenomenological parameter, and we take Σ 1 (0, 0) = 0. These constraints are convenient, rather than reflective of any deep physics. We view the effective Hamiltonian (including self energy and dissipation) as a Hartree-like theory [22] where by definition Σ 1 (0, 0) = 0 and in this way the superfluid transition occurs at a finite (positive) chemical potential.
A crucial feature of bosons is that Σ 2 changes sign at ω = 0, so that the spectral function A has the sign of ω. This is satisfied in the usual TDGL dynamics, where at low ω one has Σ 2 ∝ ω n with n = 1, corresponding to Ohmic dissipation. It should be pointed out, the apparent absence of vertex correction terms in Eq. (13) has been validated in the TDGL literature (where the bosons only experience Hartree interactions) and more directly in Ref. 29 it can be seen to be consistent with the current conservation constraint which provides the framework for including vertex terms in a Kubo formalism. For the purposes of completeness we briefly revisit the transport nomenclature found in superconducting fluctuation theories, [13,14,32], where the fluctuation propagator is found to be proportional to (17) Note that Eq. (17) can be seen to be equivalent to Eq. (14), when the dissipation is Ohmic. Here D is the so-called diffusion coefficient and λ = [2πT /gN (0)][N ′ (0)/N (0)] depends on the density of states at the Fermi energy N (0) and its derivative. By way of further comparison, we note that in solid state physics one usually describes the approach of the transition by T − T c , whereas µ − µ c is the variable associated with cold atom experiments, although these two parametrizations are equivalent.

B. Divergences in Transport
As δμ ≡ (µ − µ c )/T goes to zero (but away from the strictly critical regime) we can deduce the transport coefficients for the two-dimensional case from the integrals in Eqs. (13) and their counterparts. The literature is based on the Ohmic case n = 1 which defines the nature of the divergences in the L ij . These can be summarized [13,14,32] in terms of proportionality relations From Eqs. (18) we deduce that the Seebeck coefficient (thermopower) and Lorenz number, which involve ratios of these coefficients, behave as S ∝ −|δμ| ln |δμ|, L ∝ |δμ|.
Importantly, the thermopower changes sign at condensation and the Lorenz number tends to zero linearly. Recall that L must be greater or equal to zero for thermodynamic stability [44].
In Figure 2 we plot the transport coefficients for d = 2 as obtained from the fluctuation theory (Eq. (13)) as a function of scaled chemical potential, for a typical Γ = 0.1 [15]. We have checked that the value of Γ does not change qualitatively the transport coefficients or power laws for |δμ| 1. This figure shows that the thermal conductivity (related to L 22 ) and the mass conductivity (related to L 11 ) as well as the off-diagonal conductivity satisfying the Onsager relation (L 12 = L 21 /T ) all diverge at the transition. It should be clear, however, that this divergence is strongest for the mass conductivity. Important here is that, because of the divergence of σ = L 11 , as µ goes to µ c , L (and S) vanish at T c .

V. RESULTS: COMPARISON WITH EXPERIMENT NEAR CONDENSATION
We turn now to a comparison of theory and the experiments of Ref. 1. In Ref. 1, a Bose gas was confined in an optical trap with peak density n = 5 × 10 13 cm −3 , a large scattering length a = 22 nm, and the temperature T = 35 nK. The trap was highly oblate with very strong confinement in the z−direction and weak confinement in the two radial directions. The sample is thus in the quasi-2D regime. For these parameters we are able to address the validity of the LEA. In the normal gas regime, the scattering rate determines the thermalization time scale τ c = α/(4πna 2 v) = 2 ms where α = 2.7 is a constant [33], and the coherence length scale is given by the mean-free path l c = 3.5µm. An optical potential corrugation of up to 3 nK leads to a measured relaxation time of 20 ms of collective excitation. Both these length and time scales are small compared to the the sample size l s = 60µm≫ l c and sample lifetime of τ s = 1s≫ τ c , which suggests that the LEA is valid over a coarse grained spatial scale of 3.5µm and time scale 20 ms.
Our comparison between theory and experiment is plotted in Figs. 3. Here we show the thermopower (upper panel) and the Lorenz number (lower panel) for a twodimensional Bose gas, obtained from Eq. (13). In the simplest terms the data suggests that the thermopower changes sign at or near the transition. Experimentally, the Lorenz number approaches zero at or near T c . The error bars are sufficiently large and one should not infer that L becomes negative on the other side of the transition, as this is unphysical. We stress that our theory is applicable only to the normal phase. Indeed, one does not have many reliable theoretical tools for treating bosons (in any capacity, not just transport) near, but below the transition. In the normal phase, the figure shows that the theory and experiment of Ref. 1 coincide rather well for the thermopower S, particularly near the condensation regime. The Lorenz number is off by about a factor of 2. However, given the simplicity of the model, a disparity of order unity is not unreasonable. It should be stressed that the behavior of S and L are very strongly reflective of the divergence of the mass conductivity at the phase transition. The fact that both are found in the data to vanish at or near the transition in the experiments serves to help validate this experimental methodology [1].

VI. CONCLUSIONS
Many communities from nuclear physics, astrophysics and condensed matter are interested in transport in the cold gases. The shear viscosity has been a rather recent focus [34][35][36][37][38] and, as in this paper, there is an emerging interest in thermoelectric properties [1,3]. There will undoubtedly be many more such experimental transport studies with time. The theory in this paper is directed towards such experiments. We have emphasized that the mechanisms for momentum dissipation are very sensitive to the experimental cold atom set up. This has led us to adopt a Kubo formula approach which treats dissipation at a non-specific phenomenological level. An important consequence of our studies is the observation of quantum critical scaling and of the non-universality of (bosonic) transport, despite universality in thermodynamics. This is in contrast to Fermi liquid or gas systems where the processes contributing to transport are localized around the Fermi energy, and one finds, for example a Wiedermann-Franz law for the Lorenz ratio.
Another major finding of this work is the observation that bosonic transport in the normal state near condensation appears to be semi-quantitatively consistent with observations and theory in condensed matter studies of (bosonic) superconducting fluctuations [13,14]. In this way our paper establishes a connection between transport properties in condensed matter and in normal-phase cold atom systems. This connection has other implications for future experiments. Because ultracold gases are clean and well controlled, they may help elucidate transport in (fermionic) high T c systems [27] where bosonic degrees of freedom (through pre-formed pairs, or vortices or bipolarons, etc.) are believed to dominate transport. In particular, high temperature superconductors in condensed matter appear to have an anomalously large critical regime where these fluctuation effects are observed which is associated with the famous "pseudogap". In the cold Fermi gases near unitarity there is some debate about whether such a pseudogap exists or not. This would suggest that future transport studies of fermionic superfluids may help elucidate this issue.
In contrast to the behavior of degenerate Fermi gases, we reiterate the strong sensitivity of Bose gases to the dissipation mechanism. For bosons, different dissipation mechanisms imply different scaling laws, possibly defining different universality classes. For the future, the universality found in the thermodynamics of dilute Bose gases, either loaded or not in a lattice, [16,23,25,39], will also need to be addressed in the context of transport. This would allow a direct verification of our prediction of the breakdown of universality in transport. This work is supported by NSF-MRSEC Grant 0820054; CC acknowledges support from NSF PHY-1206095 and ARO-MURI 63834-PH-MUR. We thank Chih-Chun Chien and Yan He for helpful conversations. CC thanks Grenier Charles and Antoine Georges for useful discussion.