# Detecting Friedel oscillations in ultracold Fermi gases

###### Abstract

Investigating Friedel oscillations in ultracold gases would complement the studies performed on solid state samples with scanning-tunneling microscopes. In atomic quantum gases interactions and external potentials can be tuned freely and the inherently slower dynamics allow to access non-equilibrium dynamics following a potential or interaction quench. Here, we examine how Friedel oscillations can be observed in current ultracold gas experiments under realistic conditions. To this aim we numerically calculate the amplitude of the Friedel oscillations which a potential barrier provokes in a 1D Fermi gas and compare it to the expected atomic and photonic shot noise in a density measurement. We find that to detect Friedel oscillations the signal from several thousand one-dimensional systems has to be averaged. However, as up to 100 parallel one-dimensional systems can be prepared in a single run with present experiments, averaging over about 100 images is sufficient.

## 1 Introduction

Disturbing a homogeneous Fermi gas with an impurity gives rise to Friedel oscillations Friedel1958 (); Villain2016 (). The density distribution close to the impurity shows a spatially oscillating structure which decays with increasing distance and whose periodicity is given by half the Fermi wavelength. Friedel oscillations occur e.g. in metals when the free electron gas is disturbed by the potential associated with impurity atoms. They mediate long range interactions between individual impurities, which can give rise to the formation of ordered superstructures in adsorbates Tsong1973 (); Lau1978 () and are relevant for the interactions between magnetic impurities Briner1998 (). Using scanning tunneling microscopy (STM) Friedel oscillations have been observed in two-dimensional and one-dimensional electron gases at surfaces of solids and have served as a tool for the measurement of bandstructures and Fermi surfaces Crommie1993 (); Hasegawa1993 (); Sprunger1997 (); Hofmann1997 (); Yokoyama1998 ().

While Friedel oscillations in non-interacting systems are fully understood, the precise impact of interactions remains an open issue. Theoretical and experimental results suggest an enhancement of the oscillation amplitude for repulsive interactions Sprunger1997 (); Egger1995 (); Simion2005 (), but systematic experimental studies remain to be done. Furthermore, all observations so far have reported on static Friedel oscillations since STM cannot resolve the dynamics of electronic systems. Ultracold Fermi gases Stringari2008 (); Ketterle2008 (); Zwerger2008 () have the potential to contribute to both aspects of the topic: Interactions can easily be tuned via Feshbach resonances and the dynamics of these systems can be resolved due to their much longer intrinsic timescales. Yet so far, Friedel oscillations in ultracold gases Zwerger2005 (); Wonneberger2004 (); Eggert2009 () have not been observed. Motivated by recent advances in the generation and study of ultracold Fermi gases we investigate the feasibility of observing Friedel oscillations in a one-dimensional gas of ultracold non-interacting fermions Moritz2005 (); Zimmermann2011 (); Zwierlein2015 (); Kuhr2015 (); Greiner2015 (); Thywissen2015 (); Bloch2016 (); Hueck2017b ().

## 2 Density distribution around an impurity potential

In homogeneous non-interacting fermionic systems the one-particle eigenstates of the Hamiltonian are given by plane waves characterized by a well defined momentum . When inserting a localized impurity potential the plane waves are scattered, giving rise to standing wave patterns in each single-particle wavefunction. Close to an abrupt potential change the standing waves corresponding to different occupied momenta are in phase and add up to an oscillatory modification

(1) |

of the many-body density with respect to the unperturbed density. Here, denotes the Fermi vector and the dimensionality of the system. The phase shift depends on the precise shape of the impurity potential and the dimensionality. Since the decay of Friedel oscillations is weakest in one dimension we restrict our investigations to this case.

As a first step it is instructive to see how Friedel oscillations emerge for the simplest case, i.e. in a box potential of length bound by infinitely high walls for . Filling fermions of identical spin into the lowest eigenstates yields:

(2) |

where is the mean density far from the impurity and . This formula allows to determine the maximum amplitude of Friedel oscillations in a non-interacting system. The peak-valley amplitudes of the first and second Friedel oscillation as defined in Fig. 1 are and .

However, due to the finite slope of the impurity potential in an experimental realization the plane waves with different momenta are reflected with different phase shifts. The standing waves are hence not in phase even close to the impurity and the amplitude of the Friedel oscillations is decreased. In order to quantify this effect we perform a numerical study for a Gaussian impurity potential having a height and a -radius :

(3) |

This barrier is placed at the center of a finite size system which is limited by narrow and high potential walls outside the region of interest. We obtain the one-particle orbitals by numerically solving the discretized Schrödinger equation. The expectation value for the particle density operator of a Fermi gas at temperature and chemical potential is given by Riechers2017 ()

(4) |

Here, is the Fermi distribution function and the energy of the orbital with wavevector .

According to equation (2) the wavelength of the Friedel oscillations is given by the inverse density and therefore equals the average particle distance. Because Friedel oscillations on scales below the resolution of the imaging system cannot be observed the maximal density in a possible experiment is constrained. Accordingly we ensure that the Friedel wavelength is 4 times larger than typical resolutions of by choosing the chemical potential such that the average density is .

Figure 1 shows the density distribution around a gaussian barrier at zero temperature calculated with the approach outlined above. The peak-valley amplitudes of the first and second oscillation with respect to the average particle density are and , respectively.

Since an observation of at least two density maxima is crucial for an experimental determination of we study the impact of the impurity potential height and radius on . As shown in Fig. 2 the results range from the absence of significant Friedel oscillations for to for and , where is defined as the Fermi energy of the unperturbed system with a density . The amplitude of the oscillation is larger the more abrupt and pronounced the change in the potential is. Very narrow barriers allow for tunneling and therefore do not provoke strong Friedel oscillations.

In the second step of our analysis we quantitatively study the influence of finite temperature. An increase in temperature is accompanied by a loss of coherence and therefore a decrease in the amplitude of the Friedel oscillations is expected. In Fig. 3 results on the temperature dependence are shown for the parameters and used also for Fig. 1. As expected decreases monotonously with increasing temperature. For a temperature of that can reliably achieved in quantum gas experiments has decreased to . Hence, the expected peak-valley amplitude of the Friedel oscillations in current experiments will most likely be limited by the finite temperature rather than the finite barrier width and height.

## 3 Experimental scenario

The experimental observation of Friedel oscillations in ultracold Fermi gases requires the ability to create impurity potentials as well as to image the atomic density with a resolution on the order of half the Fermi wavelength. This has become possible in recent years with the introduction of high resolution imaging using microscope objectives Zimmermann2011 (); Zwierlein2015 (); Kuhr2015 (); Greiner2015 (); Thywissen2015 (). Quantum gas microscopy has already enabled the study of 1D fermionic lattice systems with single atom sensitivity Bloch2016 ().

In the following we describe a promising experimental scenario in which Friedel oscillations should be observable. A 2D Fermi gas of e.g. several hundred fermionic atoms (e.g.Li, K or Yb) is prepared in a single 2D layer, which is sliced into about one hundred 1D tubes by imposing an optical lattice. The distance between the 1D tubes is given by the lattice spacing. A repulsive barrier can then be projected onto the atoms using a repulsive optical potential shaped by means of spatial light modulators (see e.g. Boyer2006 (); Greiner2016 (); Hueck2017 ()) providing diffraction limited feature sizes below . It will remain a technical challenge to keep density deviations caused by imperfections in the potential landscape significantly smaller than the amplitude of the Friedel oscillations.

## 4 Expected signal to noise ratio

As shown in Sect. 2 the observation of Friedel oscillations in ultracold Fermi gases requires the detection of signal amplitudes as small as of the 1D density. This is only possible if the signal to noise ratio (SNR) of the density measurements exceeds . The two most important sources of noise in density images are the atomic shot noise and the detection noise. In the following we first focus on the atomic shot noise, and consider only a single one-dimensional system in order to discuss the signal to noise ratio in single atom sensitive fluorescence imaging. Finally, we calculate the signal to noise ratios achievable in absorption imaging.

In order to be able to resolve the Friedel oscillations spatially, the number of detection bins per Friedel wavelength should be larger than 4, otherwise the wavelength cannot be determined. The linear size of the bins is given by

(5) |

As the average interatomic distance equals the wavelength of the Friedel oscillation, the average number of atoms located within the area of a single bin is . The atomic shot noise is approximately NoiseSuppression (), yielding a relative atomic shot noise per bin of which is independent of the density. Even for a minimal number of bins of and no further noise sources the relative noise would be . This shows that suppressing the atomic shot noise to a relative level of requires an average over 1600 measurements from individual 1D systems. For state of the art fluorescence imaging no further significant detection noise is added. Here, a very deep optical lattice is used to pin the atoms to one site during detection and single atom, single site sensitive detection is achieved Zwierlein2015 (); Kuhr2015 (); Greiner2015 (); Thywissen2015 (); Bloch2016 (). We note that for this detection method the mean interparticle distance along the tubes must be at least times larger than the optical pinning lattice spacing of typically in order to be able to spatially resolve the Friedel oscillations. Since in the proposed experimental setup up to 100 parallel tubes can be prepared and imaged in each realization, only data from 10 to 100 separate runs would have to be averaged. Density measurements with such sensitivity have already been performed by averaging over 1000 tubes in a bosonic quantum gas microscope setup Bloch2012 ().

Most quantum gas experiments rely on measuring two-dimensional column densities via absorption imaging. Here, the detection noise becomes relevant and is mainly caused by photon shot noise. In the following we perform a calculation of the signal to noise ratio including photon shot noise. We find that the major limitation is still given by the atomic shot noise and that photon shot noise reduces the signal to noise ratio by at most a factor of 1.25.

In absorption imaging the 2D density is measured rather than the 1D density. For the calculation we choose the detection bins to be two-dimensional pixels with pixel lengths along the tube direction and perpendicular to it, where is the distance between tubes. This ensures that effectively only one 1D system is measured per row of pixels, despite the fact that the tube structure proposed in Sect. 3 cannot be resolved for typical lattice spacings between the tubes. Enlarging the pixel size perpendicular to the tube direction would be analogous to averaging over several parallel 1D systems. The average number of atoms per pixel is and the atomic shot noise approximately yielding per pixel.

In absorption imaging the 2D density is measured indirectly by determining the number of photons scattered by the atoms

(6) |

Here denotes the number of photons transmitted by the atoms when illuminated by photons and is the number transmitted in the absence of atoms. and originate from identically prepared laser pulses and have the same mean value, but are stochastically independent with . The transmission coefficient can be a approximated as in the limit of low optical densities which is relevant here. denotes the scattering cross section of the corresponding atomic transition.

Eq. 6 shows that it is convenient to regard the number of scattered photons as the relevant signal and to compare it with the corresponding standard deviation to determine the signal to noise ratio of the density measurement . Gaussian error propagation yields

(7) |

for the variance of the scattered photons. The signal to noise ratio then reads

(8) |

For very high numbers of incoming photons - i.e. in the limit of vanishing relative photon shot noise - the SNR is ultimately limited by

(9) |

and hence limited by atomic shot noise as in the case of single atom sensitive fluorescence imaging.

It is therefore desirable to work with high intensities and long illumination times . However, the number of photons that can be scattered by an individual atom is limited by motional blurring. When scattering photons for an extended time at an intensity an atom performs a random walk in momentum space. This leads to a motional blurring of its position with respect to its original position. To ensure that the density distribution is not altered significantly during the imaging process the condition should be fulfilled. The full calcuation Horikoshi2016 () yields an upper bound for the illumination time

(10) |

Here the saturation parameter is used where refers to the saturation intensity and to the atomic mass. , and are the linewidth, wavelength and frequency of the atomic transition. The optimal signal to noise ratio is achieved for the maximal illumination time and its dependence on and can be calculated by using Eq. (8) and

(11) |

We evaluate the optimal signal to noise ratio for an experimentally accessible configuration, i.e. for atoms, a lattice spacing of and . The results for varying density and saturation are presented in Fig. 4. For any saturation the signal to noise ratio improves with increasing density but care must be taken that the 1D density is such that the Friedel wavelength remains times larger than the optical resolution. Within the parameter range considered here a maximal SNR per pixel of can be achieved for a light atom such as Li. For heavier atoms such as K the signal to noise ratio approaches the atomic shot noise limit of .

## 5 Conclusion

In this article we study the feasibility of observing Friedel oscillations in ultracold one-dimensional Fermi gases. We numerically calculate the amplitude of the density oscillations for a suitable experimental setup and find that for currently achievable temperatures of it is on the order of of the total density. We then calculate the expected noise for a density measurement on a single 1D system, which is limited by atomic shot noise and exceeds the amplitude of the Friedel oscillations by a factor of 20. Nevertheless, since many 1D systems can be observed in a single run the noise amplitude can be sufficiently reduced by averaging over 100 images. Therefore we conclude that an observation of Friedel oscillations is experimentally feasible. This would open up the possibility to investigate their non-equilibrium dynamics and to use them to probe Fermi liquids with attractive and repulsive interactions.

## 6 Authors contributions

All authors were involved in the discussion of the physical setup and the interpretation of the numerical results. The calculations were mainly performed by K.R.. All the authors were involved in the preparation of the manuscript.

###### Acknowledgements.

The research leading to these results has received funding from the European Unionâs Seventh Framework Programme (FP7/ 2007-2013) under grant agreement No. 335431 and by the DFG in the framework of SFB 925 and the excellence cluster the Hamburg Centre for Ultrafast Imaging CUI. We thank T. Giamarchi, L. Mathey and F. Werner for stimulating discussions.## References

- (1) J. Friedel, Del Nuovo Cimento 2, 287 (1958)
- (2) J. Villain, M. Lavagna, P. Bruno, Comptes Rendus Physique 17, 302 (2016)
- (3) A. Recati, J. N. Fuchs, C. S. Peca, W. Zwerger. Phys. Rev. A 72, 023616 (2005)
- (4) T. T. Tsong, Phys. Rev. Lett. 31, 1207 (1973)
- (5) K. H. Lau, W. Kohn, Surface Science 75, 69 (1978)
- (6) B. Briner, P. Hofmann, M. Doering, H. P. Rust, E. Plummer, A. Bradshaw, Phys. Rev. B 58, 13931 (1998)
- (7) M. F. Crommie, C. P. Lutz, D. M. Eigler, Nature 363, 524 (1993)
- (8) Y. Hasegawa, P. Avouris, Phys. Rev. Lett. 71, 1071 (1993)
- (9) P. T. Sprunger, L. Petersen, E. W. Plummer, E. Laegsgaard, F. Besenbacher, Science 275, 1764 (1997)
- (10) P. Hofmann, B. Briner, M. Doering, H. P. Rust, E. Plummer, A. Bradshaw, Phys. Rev. Lett. 79, 265 (1997)
- (11) T. Yokoyama, M. Okamoto, K. Takayanagi, Phys. Rev. Lett. 80, 3423 (1998)
- (12) R. Egger, H. Grabert, Phys. Rev. Lett. 75, 3505 (1995)
- (13) G. E. Simion, G. F. Giuliani, Phys. Rev. B 72, 045127 (2005)
- (14) S. Giorgini, L. P. Pitaevski, S. Stringari, Rev. Mod. Phys. 80, 1215 (2008)
- (15) W. Ketterle, M. W. Zwierlein, in Proceedings of the International School of Physics: Enrico Fermi, Course CLXIV, Varenna 2006, edited by M. Inguscio, W. Ketterle, C. Salomon (IOS Press, Amsterdam, 2008), p. 95
- (16) I. Bloch, J. Dalibard, W. Zwerger, Rev. Mod. Phys. 80, 885 (2008)
- (17) S. N. Artemenko, G. Xianlong, W. Wonneberger, J. Phys. B 37, S49 (2004)
- (18) S. A. Söffing, M. Bortz, I. Schneider, A. Struck, M. Fleischhauer, S. Eggert, Phys. Rev. B 79, 195114 (2009)
- (19) H. Moritz, T. Stöferle, K. Günter, M. Köhl, T. Esslinger, Phys. Rev. Lett. 94, 210401 (2005)
- (20) B. Zimmermann, T. Mueller, J. Meineke, T. Esslinger, H. Moritz, New J. Phys. 13, 043007 (2011)
- (21) L. W. Cheuk, M. A. Nichols, M. Okan, T. Gersdorf, V. V. Ramasesh, W. S. Bakr, T. Lompe, M. W. Zwierlein, Phys. Rev. Lett. 114, 193001 (2015)
- (22) E. Haller, J. Hudson, A. Kelly, D. A. Cotta, B. Peaudecerf, G. D. Bruce, S. Kuhr, Nature Phys. 11, 738 (2015)
- (23) M. F. Parsons, F. Huber, A. Mazurenko, C. S. Chiu, W. Setiawan, K. Wooley-Brown, S. Blatt, M. Greiner, Phys. Rev. Lett. 114, 213002 (2015)
- (24) G. J. A. Edge, R. Anderson, D. Jervis, D. C. McKay, R. Day, S. Trotzky, J. H. Thywissen, Phys. Rev. A 92, 063406 (2015)
- (25) M. Boll, T. A. Hilker, G. Salomon, A. Omran, J. Nespolo, L. Pollet, I. Bloch, C. Gross, Science 353, 1257 (2016)
- (26) K. Hueck, N. Luick, L. Sobirey, J. Siegl, T. Lompe, H. Moritz, manuscript in preparation (2017)
- (27) K. Riechers, M. Sc. thesis, University of Hamburg (2017)
- (28) V. Boyer, R. M. Godun, G. Smirne, D. Cassettari, C. M. Chandrashekar, A. B. Deb, Z. J. Laczik, C. J. Foot, Phys. Rev. A 73, 031402 (2006)
- (29) P. Zupancic, P. M. Preiss, R. Ma, A. Lukin, M. E. Tai, M. Rispoli, R. Islam, M. Greiner, Opt. Exp. 24, 13881 (2016)
- (30) K. Hueck, A. Mazurenko, N. Luick, T. Lompe, H. Moritz, Rev. Sci. Instrum. 88, 016103 (2017)
- (31) The suppression of atomic shot noise for low temperature Fermi gases due to antibunching Mueller2010 () is only significant for detection volumes larger than the correlation length, which is not the case here.
- (32) T. Müller, B. Zimmermann, J. Meineke, J.-P. Brantut, T. Esslinger, H. Moritz, Phys. Rev. Lett. 105, 040401 (2010)
- (33) M. Cheneau, P. Barmettler, D. Poletti, M. Endres, P. Schauß, T. Fukuhara, C. Gross, I. Bloch, C. Kollath, S. Kuhr, Nature 481, 484 (2012)
- (34) M. Horikoshi, A. Ito, T. Ikemachi, Y. Aratake, M. Kuwata-Gonokami, M. Koashi, arXiv:1608.07152 (2016)