# Observation of the Efimovian Expansion in Scale Invariant Fermi Gases

###### Abstract

Scale invariance emerges and plays an important role in strongly correlated many-body systems such as critical regimes nearby phase transitions and the unitary Fermi gases. Discrete scaling symmetry also manifests itself in quantum few-body systems such as the Efimov effect. Here we report both theoretical predication and experimental observation of a novel type expansion dynamics for scale invariant quantum gases. When the frequency of the harmonic trap holding the gas decreases continuously as the inverse of time , surprisingly, the expansion of cloud size exhibits a sequence of plateaus. Remarkably, the locations of these plateaus obey a discrete geometric scaling law with a controllable scale factor and the entire expansion dynamics is governed by a log-periodic function. This striking expansion of quantum Fermi gases shares similar scaling laws and same mathematical description as the Efimov effect. Our work demonstrates the first expansion dynamics of a quantum many-body system with the temporal discrete scaling symmetry, which reveals the underlying spatial continuous scaling symmetry of the many-body system.

^{†}

^{†}thanks: They contribute equally to this work

^{†}

^{†}thanks: They contribute equally to this work

Interaction between dilute ultracold atoms is described by the -wave scattering length. For a spin- Fermi gas, when the scattering length diverges at a Feshbach resonance, there is no length scale other than the interparticle spacing in this many-body system, and therefore the system, known as the unitary Fermi gas, becomes scale invariant. The spatial scale invariance leads to universal thermodynamics and transport properties as revealed by many experiments thermodynamics0 (); thermodynamics1 (); thermodynamics2 (); thermodynamics3 (); thermodynamics4 (); thermodynamics5 (); transport0 (); transport1 (); transport2 (); transport3 (); transport4 (); transport5 (); scale (). On the other hand, in a boson system with an infinite scattering length, three-body bound state can form, where an extra length scale of the three-body parameter sets a short-range boundary condition for all three bosons being very close. It turns the continuous scaling symmetry into a discrete scaling symmetry, and gives rise to infinite number of three-body bound states whose energies obey a geometric scaling symmetry. This is well known as the Efimov effect Efimovtheory1 (); Efimovtheory2 (), which has been observed in quite a few cold atom experiments Efimov_exp1 (); Efimov_exp2 (); Efimov_exp3 (); Efimov_exp4 (); Efimov_exp5 (); Efimov_exp6 (); Efimov_exp7 (); Efimov_exp8 (); Efimov_exp9 (), and recent experiments have also confirmed the geometric scaling of the energy spectrum scaling_exp1 (); scaling_exp2 (); scaling_exp3 (); scaling_exp4 (). Both the continuous and the discrete scaling symmetry are interesting emergent phenomena in a strongly interacting system.

For a harmonic trapped gas, the expansion dynamics offers great insight to the property of the gas. Well known example is the anisotropic expansion that proves hydrodynamics behavior due to the Bose condensation BEC1 (); BEC2 () or strong interactions of Fermi gas Thomas (). Other examples are, for instance, slowing down of expansion in a disorder potential provides evidence for localization behaviors localizations1 (); localizations2 () and expansion in the presence of optical lattice reveals correlation effects Bloch (). In this work, we ask a question that, considering a scale invariant quantum gas hold by a harmonic trap, when the trap is gradually opened up with decreasing the trap frequency as ( is a coefficient and is the time), as shown in Fig. 1(a) and (b), how does the gas expand? Naively, by dimensional analysis, one would expect that the cloud size just increases as . Here we show, both theoretically and experimentally, that it is not the case. When is smaller than a critical value, the expansion dynamics displays a discrete scaling symmetry in the time domain. As shown in Fig. 1(c), as a function of displays a sequence of plateaus, which mean that at a set of discrete times the cloud expansion surprisingly stops, despite of the continuous decreasing of the trap frequency. The locations of the plateaus obey a geometric scaling behavior. More interestingly, this striking behavior of a discrete scaling law in the time domain is in fact a consequence of the continuous spatial scaling symmetry.

To explain this intriguing dynamics, we shall first point out the insight that why is so special. For simplicity, we first consider a three-dimensional isotropic trap . Let us consider that a many-body system, if in the absence of a trapping potential, is invariant under a scale transformation , while in the presence of a static harmonic trap, the fixed harmonic length introduces an additional length scale that could break this spatial scale invariance. Nevertheless, if changes as , the time-dependent Schrödinger equation exhibits a space-time scaling symmetry under the transformation and . With the scaling symmetry, it is straightforward to derive that the equation-of-motion for the cloud size is given by (see supplementary material for the detailed derivation):

(1) |

Obviously, the differential equation is invariant under a continuous scaling of time . However, in practices, one should always start with a finite initial trap frequency before turning it down, which corresponds to an initial time with , as shown in Fig. 1(b). The system is at equilibrium for , and at , and for . This sets a boundary condition for Eq. 1 which can turn the continuous scaling symmetry in the time domain into a discrete one. The solution of this differential equation depends on the value of . When , the solution is log-periodic function as

(2) |

where and . Eq. 2 clearly reveals the discrete scaling symmetry, i.e. when , and for all the -th order derivatives. Therefore, at time , the first- and the second-order time derivatives for become zero and the cloud expansion is strongly suppressed, that is to say, the expansion dynamics shows a series of plateaus around each . While when , simply follows a power law as for , where .

Here we shall emphasize that this intriguing expansion dynamics is a universal phenomenon for all scale invariant quantum gases, which include, for examples, non-interacting gas, the unitary Fermi gas in three-dimension, weakly-interacting gases in two-dimension (when anomaly can be ignored), and the Tonks gas in one-dimension. It is also independent of temperature.

Before proceeding to experimental observation of such dynamical expansion in the ultracold Fermi gases, we would also like to bring out the analogy to the Efimov effect. First, when solving the three-body problem of bosons in the hyper-spherical coordinate, one finally reaches an effective potential as ( is the so-called hyper-radius) that scales the same way as the kinetic energy, and the Schrödinger equation then reduces to a one-dimensional scale invariant differential equation Efimovtheory1 (); Efimovtheory2 (). Similarly, our equation-of-motion for is also a scale invariant differential equation. Secondly, the short-range boundary condition (i.e. the three-body parameter) plays the similar role as the initial trap frequency here, which sets a boundary condition for the differential equation and turns the symmetry into a discrete scaling symmetry. Thirdly, the solution for the three-body wave function is also a log-periodic function, which is the same function as our solution Eq. 2. Finally, in our case plays the similar role as the mass ratio in the Efimov problem that controls the scaling factor, as well as whether the effect will occur. Hence, this particular dynamical expansion of a scale invariant gas shares the same symmetry property and similar mathematical description as the three-boson problem of bosons with infinite scattering length. It can be viewed as the counterpart of the Efimov effect in the time domain, and thus is named as “the Efimovian expansion” here.

In our experiment, we use a balanced mixture of Li fermions in the lowest two hyperfine states and . Fermionic atoms are loaded into a cross-dipole trap to perform evaporative cooling. The resulting potential has a cylindrical symmetry around the propagation axis of the laser and the trap anisotropic frequency ratio is about . The trap anisotropy causes an additional complication compared with the isotropic case discussed above. Nevertheless, as shown in the supplementary material, similar results can be obtained for the anisotropic case. Starting at the initial time , the trap potential is lowered as

(3) |

Since and the effect is more pronounced along the axial direction than in the transverse direction. (See supplementary material for detailed reason.) Therefore, hereafter we focus on the cloud expansion along the axial direction. Theory shows that the axial cloud square size obeys the same form as Eq. 2 except

(4) |

where is a factor related to the breathing mode frequency, and for the non-interacting gas and for the unitary Fermi gas along the axial direction. Feshbach resonance is used to tune the interaction of the atoms either to the non-interacting regime with the magnetic field G or to the unitary regime with G. The trap frequency is lowered by decreasing the laser intensity, and is controlled by the speed of how fast the laser intensity decreases, with the initial axial trap depth always fixed at where is the full trap potential. Thus, different corresponds to different , where is the initial axial trap frequency. Finally, after certain expansion time with the trap, the trap is completely turned off and the cloud is probed by standard resonant absorption imaging techniques after a time-of-flight expansion time . Each data point is an average of shots of the measurements at identical parameters.

The time-of-flight density profile along the axial direction is fitted by a Gaussian function as , from which we obtain . is related to the in-situ cloud size by a scale factor via . can be obtained by either hydrodynamic or ballistic expansion equation with the time-of-flight time (see the supplementary material for detail), with which we can deduce from . Since the trap is quite anisotropic, the cloud size expands slowly along the axial direction during a short time-of-flight, and the expansion factor only gives a quantitive correction to the results. Fig. 2 shows the typical measurements of with different for both the non-interacting and the unitary Fermi gases. For instance, for , we decrease the trap frequency from Hz to Hz within ms. Dots are the measured data, and the solid and the dashed lines are both theoretical curves based on Eq. 2, taking as fitting parameter or using given by Eq. 4, respectively. Since is obtained by a Gaussian fit to the density profile, and thus the theoretical expression for is simply a square root of Eq. 2. Fig. 2 clearly shows the plateaus for the expansion dynamics and the excellent agreement between theory and experiment. In the inset of Fig. 2(b) we also show density profiles for three successive measurement times inside a plateau. The entire density profiles almost perfectly overlap with each other, which clearly confirms that the expansion stops at the plateau.

For smaller , the trap frequency decreases slower, we find that the plateaus becomes denser and the difference in height between two adjacent plateaus becomes smaller. Eventually in the limit , one reaches the adiabatic limit and Eq. 2 reduces to . (See the supplementary material for discussions) For larger , the trap frequency drops faster and the plateau appears in a later time. In practice the experiment has a limit, that is, when the trap frequency is too below (say, below Hz), it is hard to get high quality data due to the low resolution of the images after the time-of-flight expansion. It prohibits us from observing plateaus for large and investigating the critical regime where the Efimovian expansion behavior disappears.

Now we further demonstrate that this dynamics is universal. First, we should verify that relates to via Eq. 4. In the experiment, is determined by the trap frequencies measured by the parametric resonance and is extracted from the best fit of the expansion data in Fig. 2. The universal relation between and is plotted in Fig. 3(a). can fit very well with a linear function which gives the slope for the non-interacting case and for the unitary Fermi gas. These are in good agreement with for the non-interacting case and for the unitary case. The Efimovian expansion is also robust and insensitive to the temperature and atoms number of the Fermi gas. Fig. 3 (b) shows the dependence of of the unitary Fermi gas on these parameters, where is fixed and both the temperature and the atom number are varied by controlling the evaporative cooling process. It shows that is almost a constant holding for the changing of the atoms number of each spin component from to and the dimensionless temperature varying from to , where is the Fermi temperature.

Second, we notice that the non-interacting and the unitary cases only differ in the relation between and , and once is given to be the same, the dynamics is the same for these two different systems. In another word, as a function of (or ) and is a universal function for all scale invariant state-of-matters. To better show this, we choose an expansion curve for the unitary Fermi gas and another expansion curve for the non-interacting gas with the same value of , and we plot these two expansion curves together, as shown in Fig. 4. For these two cases we have examined (two different ), the entire expansion curves perfectly coincide with each other. This proves experimentally that the Efimovian dynamics is universal for different scale invariant quantum gases.

At last, we study a time-reversal compression process. Consider an expansion process from to , the trap frequency decreases from to . Now we consider an inverted process of increasing the trap frequency as , the trap frequency increases from to when changes from to . In order for the compression dynamics to really invert the expansion dynamics, has to be carefully chosen to satisfy . We perform such an experiment, and both the expansion and the compression dynamics are presented in Fig. 5 for comparison, which shows that the dynamical process with a carefully chosen boundary is time-reversal symmetric. The little asymmetry is because that the lowering trap (black dots) also plays a role as the evaporative cooling, which causes the atoms cooler and the cloud sizes smaller correspondingly.

In summary, we have discovered a novel type expansion dynamics for scale invariant quantum gases. It exhibits geometric scaling behavior, which is mathematically described by a log-periodic function, and displays plateaus features, which counter-intuitively shows that the cloud nearly stops expanding despite of the continuously decreasing of the confinement potential. Moreover, it is universal for all scale invariant quantum gases. Therefore, it can serve as a tool to detect the emergent scale invariance, similar as the anisotropic expansion now serving as a tool to probe the emergent hydrodynamics behavior. It can be used in future studies of dynamics in the quantum critical regime critical1 (); critical2 (); critical3 (); critical4 (); critical5 () and to calibrate scale symmetry anomaly in a two-dimensional quantum gases anomaly1 (); anomaly2 (); anomaly3 (); anomaly4 ().

## References

- (1) J. Kinast, A. Turlapov, J. E. Thomas, Q. Chen, J. Stajic, and K. Levin, Science 307, 1296 (2005).
- (2) S. Nascimbne, N. Navon, K. J. Jiang, F. Chevy, and C. Salomon, Nature 463, 1057 (2010).
- (3) M. Horikoshi, S. Nakajima, M. Ueda, and T. Mukaiyama, Science 327, 442 (2010).
- (4) N. Navon, S. Nascimbne, F. Chevy, and C. Salomon, Science 328, 729 (2010).
- (5) M. J. H. Ku, A. T. Sommer, L. W. Cheuk, and M. W. Zwierlein, Science 335, 563 (2012).
- (6) L. A. Sidorenkov, M. K. Tey, R. Grimm, Y.-H. Hou, L. Pitaevskii, and S. Stringari, Nature 498, 78 (2013).
- (7) C. Cao, E. Elliott, J. Joseph, H. Wu, J. Petricka, T. Schaefer, and J. E. Thomas, Science 58, 33101 (2011).
- (8) A. Sommer, M. Ku, G. Roati, and M. W. Zwierlein, Nature 472, 201 (2011).
- (9) D. Stadler, S. Krinner1, J. Meineke, J.-P. Brantut, and T. Esslinger, Nature 491, 736 (2012).
- (10) J.-P. Brantut, C. Grenier, J. Meineke, D. Stadler, S. Krinner, C. Kollath, T. Esslinger, and A. Georges, Science 342, 713 (2013).
- (11) M. Koschorreck, D. Pertot, E. Vogt, and M. Khl, Nat. Phys. 9, 405 (2013).
- (12) A. B. Bardon, S. Beattie, C. Luciuk, W. Cairncross, D. Fine, N. S. Cheng, G. J. A. Edge, E. Taylor, S. Zhang, S. Trotzky, and J. H. Thywissen, Science 344, 722 (2014).
- (13) E. Elliott, J. Joseph, and J. E. Thomas, Phys. Rev. Lett. 112, 040405 (2014).
- (14) V. Efimov, Phys. Lett. 33B, 563 (1970)
- (15) E. Braaten and H.-W. Hammer, Phys. Rep. 428, 259 (2006).
- (16) T. Kraemer, M. Mark, P. Waldburger, J. Danzl, C. Chin, B. Engeser, A. Lange, K. Pilch, A. Jaakkola, H. Ngerl, and R. Grimm, Nature 440, 315 (2006).
- (17) S. Knoop, F. Ferlaino, M. Mark, M. Berninger, H. Schobel, H. Nagerl, and R. Grimm, Nat. Phys. 5, 227 (2009).
- (18) G. Barontini, C. Weber, F. Rabatti, J. Catani, G. Thalhammer, M. Inguscio, and F. Minardi, Phys. Rev. Lett. 103, 043201 (2009).
- (19) J. R. Williams, E. L. Hazlett, J. H. Huckans, R. W. Stites, Y. Zhang, and K. M. O¡¯Hara, Phys. Rev. Lett. 103, 130404 (2009).
- (20) S. Pollack, D. Dries, and R. Hulet, Science 326, 1683 (2009).
- (21) M. Berninger, A. Zenesini, B. Huang, W. Harm, H.-C. Ngerl, F. Ferlaino, R. Grimm, P. S. Julienne, and J. M. Hutson, Phys. Rev. Lett. 107, 120401 (2011).
- (22) S. Roy, M. Landini, A. Trenkwalder, G. Semeghini, G. Spagnolli, A. Simoni, M. Fattori, M. Inguscio, and G. Modugno, Phys. Rev. Lett. 111, 053202 (2013).
- (23) R. S. Bloom, M.-G. Hu, T. D. Cumby, and D. S. Jin, Phys. Rev. Lett. 111, 105301 (2013).
- (24) M. Kunitski, S. Zeller, J. Voigtsberger, A. Kalinin, L. Ph. H. Schmidt, M. Schffler, A. Czasch, W. Schllkopf, R. E. Grisenti, T. Jahnke, D. Blume, and R. Drner, Science 348, 551 (2015).
- (25) M. Zaccanti, B. Deissler, C. D’Errico, M. Fattori, M. Jona-Lasinio, S. Mller, G. Roati, M. Inguscio, and G. Modugno, Nat. Phys. 5, 586 (2009).
- (26) B. Huang, L. A. Sidorenkov, R. Grimm, and J. M. Hutson, Phys. Rev. Lett. 112, 190401 (2014).
- (27) S.-K. Tung, K. Jimnez-Garca, J. Johansen, C. V. Parker, and C. Chin, Phys. Rev. Lett. 113, 240402 (2014).
- (28) R. Pires, J. Ulmanis, S. Häfner, M. Repp, A. Arias, E. D. Kuhnle, and M. Weidemüller, Phys. Rev. Lett. 112, 250404 (2014)
- (29) M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Science 269, 198 (1995).
- (30) K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995).
- (31) K. M. O’Hara, S. L. Hemmer, M. E. Gehm, S. R. Granade, and J. E. Thomas, Science 298, 2179 (2002).
- (32) C. Fort, L. Fallani, V. Guarrera, J. E. Lye, M. Modugno, D. S. Wiersma, and M. Inguscio, Phys. Rev. Lett. 95, 170410 (2005).
- (33) L. Sanchez-Palencia, D. Clment, P. Lugan, P. Bouyer, G. V. Shlyapnikov, and A. Aspect, Phys. Rev. Lett. 98, 210401 (2007).
- (34) M. Greiner, O. Mandel, T. Esslinger, T. W. Hnsch, and I. Bloch, Nature 415, 39 (2002).
- (35) P. Coleman and A. J. Schofield, Nature 433, 226 (2005).
- (36) T. Donner , S. Ritter, T. Bourdel, A. ttl, M. Khl, and T. Esslinger, Science 315, 1556 (2007).
- (37) C.-L. Hung, X. Zhang, N. Gemelke, and C. Chin, Nature 470, 236 (2011).
- (38) J. Bork, Y. Zhang, L. Diekhner, L. Borda, P. Simon, J. Kroha, P. Wahl, and K. Kern, Nat. Phys. 7, 901 (2011).
- (39) X. Zhang, C.-L. Hung, S.-K. Tung, and C. Chin, Science 335, 6072 (2012).
- (40) E. Vogt, M. Feld, B. Frhlich, D. Pertot, M. Koschorreck, and M. Khl, Phys. Rev. Lett. 108, 070404 (2012).
- (41) J. Hofmann, Phys. Rev. Lett. 108, 185303 (2012).
- (42) E. Taylor and M. Randeria, Phys. Rev. Lett. 109, 135301 (2012).
- (43) C. Gao and Z. Yu, Phys. Rev. A 86, 043609 (2012).

## Appendix A Supplemental Material

### a.1 The dynamical scaling solution

In the following, we consider a unitary Fermi gas trapped in a harmonic potential whose frequency is time dependent. The Hamiltonian is given by

(5) |

where the non-interacting part

(6) |

represents the kinetic energy plus a time dependent harmonic trap, is total particle number and represents the short range interaction between spin up and down fermions. Due to the divergence of the -wave scattering length , is scale invariant in the zero interaction range limit such that . Thanks to this property of and the specific choice of time dependence in , the total Hamiltonian has continuous scaling symmetry in both configurational and temporal spaces.

Isotropic trapping. To illustrate the basic idea and get a better physical intuition, we first consider the simpler isotropic trap with . In this case, we can calculate the cloud size directly from its equation of motion. The first derivative of is given as

(7) |

where is the generator of a spacial scaling transformation. On the other hand, the equation of motion of is

(8) |

where we have used the fact that is scale invariant such that . Combining these two equations we obtain

(9) |

To make the equation closed, one still needs to calculate which can be obtain by the Hellmann-Feynman theorem:

(10) |

Combing (9) and (10), we finally obtain the following equation of motion for

(11) |

Since this is a third order differential equation, one needs three initial conditions, which are the values of and at the starting point , to fix the solution. On one hand, since the system remains at static before the expansion, arbitrary order of time derivative of remains 0 for all . On the other hand, since is continuous while is discontinuous at , it is easy to check that only the first and second derivative of is continuous across while higher derivatives are discontinuous. Thus we will apply the following initial conditions

(12) | |||||

(13) |

The solution of (11) has very different behavior for small and large value of . For , we have

(14) |

where . One can see that in the limit , the cloud size follows a simple power law , where can be seen as an anomalous dimension in the time domain. As a result, the continuous scaling symmetry is still preserved for and there is no dynamic Efimov effect in this case.

The situation is much more interesting when . In this case, we have

(15) |

where and . Instead of a simple power law form, now contains a logarithmic periodic part which breaks the continuous scaling symmetry down to a discrete one in the time domain and satisfy

(16) |

for all integers of .

Adiabatic limit. Now let us consider the adiabatic limit() of the expansion. Physically, it represents the situation that the trap expands extremely slow. In this limit, Eq. 15 can be expressed as

(17) |

We find that the cloud size follows the size of the trap and grows proportional to . This behavior is consistent with the adiabatic theorem which claims that the system remains in the instantaneous ground state of .

Anisotropic trapping. Above approach, although elegant and rigorous, is only valid for an isotropic harmonic trap. For more general systems with anisotropic harmonic trap, there is no known exact solutions for the dynamic expansion. In this case, inspired by the logarithmic time dependence, we first perform the following scaling transformation

(18) | |||||

(19) | |||||

(20) |

where is the new position vector, and is the total number of fermions. This new set of space and time coordinates are now all dimensionless variables. After this change of variables, we obtain a new Hamiltonian for the transformed wave function

(21) |

Due to the scaling symmetry of the system, the original dynamic expansion problem is now transformed into a quantum quench problem under the new static Hamiltonian given by Eq. 21. This Hamiltonian corresponds to a unitary Fermi gas in a static harmonic trap with trapping frequency where and with a single particle gauge field . The initial state of this quench process is an eigen-state of another Hamiltonian

(22) |

Since the gauge field satisfies , it can be eliminated by a local gauge transformation which will introduce an initial current to through the phase factor . When , the quench will excite one or multiple breathing modes with small oscillation amplitude. In the axial symmetric case when and , this problem can be solved with a hydrodynamic approach proposed in Hou et al. (2013). In this method, the system is described by a set of density field such as the number density , entropy density . An exact solution to the two-fluid hydrodynamic equations was found in Ref. Hou et al. (2013) with a scaling ansatz

(23) |

where and the similar form of solution is also applied for . The axial size of the cloud for transformed wave function is thus given as

(24) |

where is a constant. After transforming back to the original coordinates and when , one finally gets

(25) |

As a result, the axial expansion of cloud is fully determined by the function. In the small perturbation regime, the functions satisfy the following linearized coupled differential equations Hou et al. (2013)

(26) | |||||

(27) |

It is obvious that the solution of takes a simple form , with and to be fixed by initial conditions. Here is the eigen-frequency of (26) and (27) whose explicit expression was given in Hou et al. (2013). In the limit of , becomes and . Moreover, the amplitude vanishes in this limit, thus we can subsitute this solution into (25) and implement the same initial conditions given in (12) and (13). We find that are given by exactly the same expression as in (15) but with a different .

### a.2 Experimental methods

Our experimental setup was presented in Deng et al. (2015). A large atom number magneto-optical trap (MOT) is realized by a laser system of 2.5-watts laser output with Raman fiber amplifier and intracavity-frequency-doubler, as shown in Fig. 1(a). After the loading stage, cold atoms are obtained. Then MOT is compressed to obtain atoms of near Doppler-limited temperature of 280 and a density of cm. Following this compressed stage, the MOT gradient magnets are extinguished and repumping beams are switched off faster than the cooling beams. By optical pumping, a balanced mixture of atoms in the two lowest hyperfine states and is prepared.

The optical crossed dipole trap is consisted by a 200 W Ytterbium fiber laser at 1070nm. The resulting potential has a cylindrical symmetry around the propagation axis of the laser. The evaporative cooling is performed at the Feshbach resonance of the magnetic field Gauss. The trap is turned on 100 ms before the MOT compressed stage. After the atoms are loaded into the dipole trap we hold the atoms 200 ms on the trap and then forced evaporative cooling is followed by lowering the trap laser power, as shown in Fig. 1(b). First a simple exponential ramp [] is used as a lowering curve, where is the full trap depth. The time constant is selected to control the trapping potential down to a variable value, which allows us to vary the final temperature and atom number of the cloud. After the forced evaporative cooling the dipole trap is recompressed to a value . Then the trap depth is held 0.5 s for the equilibrium and is lowered with the relation . After the time , the trap is turned off and a standard absorption image is followed to extract the cloud profiles. In the analysis of the images we use a Gaussian function to fit the column density and get the mean cloud size.

### a.3 The expansion factor

The dynamical evolution is investigated by measuring the cloud radii after the cloud is released from the trap. So the size of the observed expanded cloud is related to that of the trapped gas by a scale factor , such that . The methods to analyze of the non-interacting and the unitary Fermi gas are presented as follows:

Non-interacting Fermi gas. For the non-interacting Fermi gas, we can solve the dynamics of the time-of-flight exactly. Applying the Heisenberg equation to , we obtain

(28) |

Here is the kinetic energy along the -direction, which is a constant during the time-of-flight. We assume that the system is stationary before we release the trap. Then according to the Virial theorem,

(29) |

where and is the trapping frequency and the cloud size just before the time-of-flight. Eq. 28 can be solved exactly and we obtain the scale factor

(30) |

where is the flight time.

Unitary Fermi gas. To calculate the scale factor of the unitary Fermi gas, we follow the hydrodynamic approach used in Ref. Elliott et al. (2014). First the density is given by

(31) |

where , is a time dependent scale factor. is the density profile of the trapped cloud before we start the Efimovian expansion. Here, is the volume scale factor, which is independent of the spatial coordinates in the scaling approximation. With Eq. 31 and a velocity field that is linear in the spatial coordinates, . Here , and , where is the mean-square cloud radius of the trapped cloud in the -direction, just before releasing. Then, with these scaling assumptions and for the unitary Fermi gas, the fluids dynamics equations yield

(32) |

where the heating from the shear viscosity is neglected. To add this term into the equation is straightforward, as shown in Ref. Cao et al. (2011).

Before the Efimovian expansion, the balance of the forces arising from the pressure and the trapping potential yields

(33) |

Substituting Eq. 33 into Eq. 32, we get the evolution equations of the expansion factors,

(34) |

where

(35) |

Here, is the total trap potential including the optical trapping potential and the residual magnetic potential from the Feshbach coils. For the harmonic trap potentials, Eq. 35 becomes

(36) |

where is measured from the cloud profile and is

(37) |

Here is the time when we start the Efimovian expansion, is the time we turn off the trap and is the time-of-flight. We numerically solve Eq. 36 to determine with the initial conditions and . The scale factor of the time-of-flight is then . In Fig. 7, we show and of the unitary Fermi gas for flight time .

## References

- Hou et al. (2013) Y. Hou, L. Pitaevskii, and S. Stringari, Phys. Rev. A 87, 033620 (2013).
- Deng et al. (2015) S. Deng, P. Diao, Q. Yu, and H. Wu, Chin. Phys. Lett. 32, 053401 (2015).
- Elliott et al. (2014) E. Elliott, J. Joseph, and J. Thomas, Phys. Rev. Lett. 112, 040405 (2014).
- Cao et al. (2011) C. Cao, E. Elliott, J. Joseph, H. Wu, J. Petricka, T. Schaefer, and J. Thomas, Science. 331, 58 (2011).