- Letter
- Published:
A cold-atom Fermi–Hubbard antiferromagnet
- Anton Mazurenko1,
- Christie S. Chiu1,
- Geoffrey Ji1,
- Maxwell F. Parsons1,
- Márton Kanász-Nagy1,
- Richard Schmidt1,
- Fabian Grusdt1,
- Eugene Demler1,
- Daniel Greif1 &
- …
- Markus Greiner1
Naturevolume 545, pages462–466 (2017)Cite this article
23kAccesses
87Altmetric
Abstract
Exotic phenomena in systems with strongly correlated electrons emerge from the interplay between spin and motional degrees of freedom. For example, doping an antiferromagnet is expected to give rise to pseudogap states and high-temperature superconductors1. Quantum simulation2,3,4,5,6,7,8 using ultracold fermions in optical lattices could help to answer open questions about the doped Hubbard Hamiltonian9,10,11,12,13,14, and has recently been advanced by quantum gas microscopy15,16,17,18,19,20. Here we report the realization of an antiferromagnet in a repulsively interacting Fermi gas on a two-dimensional square lattice of about 80 sites at a temperature of 0.25 times the tunnelling energy. The antiferromagnetic long-range order manifests through the divergence of the correlation length, which reaches the size of the system, the development of a peak in the spin structure factor and a staggered magnetization that is close to the ground-state value. We hole-dope the system away from half-filling, towards a regime in which complex many-body states are expected, and find that strong magnetic correlations persist at the antiferromagnetic ordering vector up to dopings of about 15 per cent. In this regime, numerical simulations are challenging21 and so experiments provide a valuable benchmark. Our results demonstrate that microscopy of cold atoms in optical lattices can help us to understand the low-temperature Fermi–Hubbard model.
This is a preview of subscription content,access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
9,800 Yen / 30 days
cancel any time
Subscription info for Japanese customers
We have a dedicated website for our Japanese customers. Please go tonatureasia.com to subscribe to this journal.
Prices may be subject to local taxes which are calculated during checkout




Similar content being viewed by others
References
Lee, P., Nagaosa, N. & Wen, X.-G. Doping a Mott insulator: physics of high-temperature superconductivity.Rev. Mod. Phys.78, 17–85 (2006)
Friedenauer, A., Schmitz, H., Glueckert, J. T., Porras, D. & Schaetz, T. Simulating a quantum magnet with trapped ions.Nat. Phys.4, 757–761 (2008)
Kim, K. et al. Quantum simulation of frustrated Ising spins with trapped ions.Nature465, 590–593 (2010)
Struck, J. et al. Quantum simulation of frustrated classical magnetism in triangular optical lattices.Science333, 996–999 (2011)
Simon, J. et al. Quantum simulation of antiferromagnetic spin chains in an optical lattice.Nature472, 307–312 (2011)
Yan, B. et al. Observation of dipolar spin-exchange interactions with lattice-confined polar molecules.Nature501, 521–525 (2013)
Drewes, J. H. et al. Antiferromagnetic correlations in two-dimensional fermionic Mott-insulating and metallic phases.Phys. Rev. Lett.118, 170401 (2017)
Murmann, S. et al. Antiferromagnetic Heisenberg spin chain of a few cold atoms in a one-dimensional trap.Phys. Rev. Lett.115, 215301 (2015)
Hofstetter, W., Cirac, J. I., Zoller, P., Demler, E. & Lukin, M. D. High-temperature superfluidity of fermionic atoms in optical lattices.Phys. Rev. Lett.89, 220407 (2002)
Bloch, I., Dalibard, J. & Zwerger, W. Many-body physics with ultracold gases.Rev. Mod. Phys.80, 885–964 (2008)
Jördens, R., Strohmaier, N., Günter, K., Moritz, H. & Esslinger, T. A Mott insulator of fermionic atoms in an optical lattice.Nature455, 204–207 (2008)
Schneider, U. et al. Metallic and insulating phases of repulsively interacting fermions in a 3D optical lattice.Science322, 1520–1525 (2008)
Greif, D., Uehlinger, T., Jotzu, G., Tarruell, L. & Esslinger, T. Short-range quantum magnetism of ultracold fermions in an optical lattice.Science340, 1307–1310 (2013)
Hart, R. A. et al. Observation of antiferromagnetic correlations in the Hubbard model with ultracold atoms.Nature519, 211–214 (2015)
Haller, E. et al. Single-atom imaging of fermions in a quantum-gas microscope.Nat. Phys.11, 738–742 (2015)
Edge, G. J. A. et al. Imaging and addressing of individual fermionic atoms in an optical lattice.Phys. Rev. A92, 063406 (2015)
Parsons, M. F. et al. Site-resolved measurement of the spin-correlation function in the Fermi-Hubbard model.Science353, 1253–1256 (2016)
Boll, M. et al. Spin- and density-resolved microscopy of antiferromagnetic correlations in Fermi-Hubbard chains.Science353, 1257–1260 (2016)
Cheuk, L. W. et al. Observation of spatial charge and spin correlations in the 2D Fermi-Hubbard model.Science353, 1260–1264 (2016)
Brown, P. T. et al. Observation of canted antiferromagnetism with ultracold fermions in an optical lattice. Preprint athttps://arxiv.org/abs/1612.07746 (2016)
Staar, P., Maier, T. & Schulthess, T. C. Dynamical cluster approximation with continuous lattice self energy.Phys. Rev. B88, 115101 (2013)
Manousakis, E. The spin-1/2 Heisenberg antiferromagnet on a square lattice and its application to the cuprous oxides.Rev. Mod. Phys.63, 1–62 (1991)
Sandvik, A. W. Finite-size scaling of the ground-state parameters of the two-dimensional Heisenberg model.Phys. Rev. B56, 11678–11690 (1997)
Liang, J., Kohn, R. N., Becker, M. F. & Heinzen, D. J. High-precision laser beam shaping using a binary-amplitude spatial light modulator.Appl. Opt.49, 1323–1330 (2010)
Ho, T.-L. & Zhou, Q. Universal cooling scheme for quantum simulation. Preprint athttps://arxiv.org/abs/0911.5506 (2009)
Gorelik, E. V. et al. Universal probes for antiferromagnetic correlations and entropy in cold fermions on optical lattices.Phys. Rev. A85, 061602 (2012)
Chakravarty, S., Halperin, B. I. & Nelson, D. R. Two-dimensional quantum Heisenberg antiferromagnet at low temperatures.Phys. Rev. B39, 2344–2371 (1989)
Denteneer, P. J. H. & Van Leeuwen, J. M. J. Spin waves in the half-filled Hubbard Model beyond the random phase approximation.Europhys. Lett.22, 413–418 (1993)
Hofferberth, S. et al. Probing quantum and thermal noise in an interacting many-body system.Nat. Phys.4, 489–495 (2008)
Yamada, K. et al. Doping dependence of the spatially modulated dynamical spin correlations and the superconducting-transition temperature in La2−xSrxCuO4 .Phys. Rev. B57, 6165–6172 (1998)
Zupancic, P. et al. Ultra-precise holographic beam shaping for microscopic quantum control.Opt. Express24, 13881–13893 (2016)
Gaunt, A. L. et al. Robust digital holography for ultracold atom trapping.Sci. Rep.2, 721 (2012)
Hueck, K., Mazurenko, A., Luick, N., Lompe, T. & Moritz, H. Suppression of kHz-frequency switching noise in digital micro-mirror devices.Rev. Sci. Instrum.88, 016103 (2017)
Parsons, M. F. et al. Site-resolved imaging of fermionic6Li in an optical lattice.Phys. Rev. Lett.114, 213002 (2015)
Paiva, T., Scalettar, R., Randeria, M. & Trivedi, N. Fermions in 2D optical lattices: temperature and entropy scales for observing antiferromagnetism and superfluidity.Phys. Rev. Lett.104, 066406 (2010)
Hubbard, J. Electron correlations in narrow energy bands.Proc. R. Soc. Lond. A276, 238–257 (1963)
LeBlanc, J. P. F. & Gull, E. Equation of state of the fermionic two-dimensional Hubbard model.Phys. Rev. B88, 155108 (2013)
LeBlanc, J. P. F. et al. Solutions of the two-dimensional Hubbard model: benchmarks and results from a wide range of numerical algorithms.Phys. Rev. X5, 041041 (2015)
Khatami, E. & Rigol, M. Thermodynamics of strongly interacting fermions in two-dimensional optical lattices.Phys. Rev. A84, 053611 (2011)
Machida, K. Magnetism in La2CuO4 based compounds.Physica C158, 192–196 (1989)
Schulz, H. J. Incommensurate antiferromagnetism in the two-dimensional Hubbard model.Phys. Rev. Lett.64, 1445–1448 (1990)
Schulz, H. J. inThe Hubbard Model: ItsPhysics and Mathematical Physics (eds Baeriswyl, D. et al.) 89–102 (Springer, 1995)
Chubukov, A. V., Sachdev, S. & Ye, J. Theory of two-dimensional quantum Heisenberg antiferromagnets with a nearly critical ground state.Phys. Rev. B49, 11919–11961 (1994)
Caffarel, M. et al. Monte Carlo calculation of the spin stiffness of the two-dimensional Heisenberg model.Europhys. Lett.26, 493–498 (1994)
Chang, C.-C., Gogolenko, S., Perez, J., Bai, Z. & Scalettar, R. T. Recent advances in determinant quantum Monte Carlo.Phil. Mag.95, 1260–1281 (2015)
Gogolenko, S., Bai, Z. & Scalettar, R. Structured orthogonal inversion of block p-cyclic matrices on multicores with GPU accelerators. InEuro-Par 2014 Parallel Processing (eds Silva, F. et al.) 524–535 (Lecture Notes in Computer Science Vol.8632, Springer, 2014)
Jiang, C ., Bai, Z . & Scalettar, R. T. A fast selected inversion algorithm for Green’s function calculations in many-body quantum Monte Carlo simulations. InProc. 30th IEEE International Parallel and Distributed Processing Symposiumhttps://doi.org/10.1109/IPDPS.2016.69 (IEEE, 2016)
Cherng, R. W. & Demler, E. Quantum noise analysis of spin systems realized with cold atoms.New J. Phys.9, 7 (2007)
Braungardt, S., Sen(De), A., Sen, U., Glauber, R. J. & Lewenstein, M. Fermion and spin counting in strongly correlated systems.Phys. Rev. A78, 063613 (2008)
Braungardt, S. et al. Counting of fermions and spins in strongly correlated systems in and out of thermal equilibrium.Phys. Rev. A83, 013601 (2011)
Lamacraft, A. Noise and counting statistics of insulating phases in one-dimensional optical lattices.Phys. Rev. A76, 011603 (2007)
Blanter, Y. & Büttiker, M. Shot noise in mesoscopic conductors.Phys. Rep.336, 1–166 (2000)
de-Picciotto, R. et al. Direct observation of a fractional charge.Nature389, 162–164 (1997)
Saminadayar, L., Glattli, D. C., Jin, Y. & Etienne, B. Observation of the e/3 fractionally charged Laughlin quasiparticle.Phys. Rev. Lett.79, 2526–2529 (1997)
Goldman, V. J. & Su, B. Resonant tunneling in the quantum Hall regime: measurement of fractional charge.Science267, 1010–1012 (1995)
Gritsev, V., Altman, E., Demler, E. & Polkovnikov, A. Full quantum distribution of contrast in interference experiments between interacting one-dimensional Bose liquids.Nat. Phys.2, 705–709 (2006)
Polkovnikov, A., Altman, E. & Demler, E. Interference between independent fluctuating condensates.Proc. Natl Acad. Sci. USA103, 6125–6129 (2006)
Mermin, N. D. & Wagner, H. Absence of ferromagnetism or antiferromagnetism in one- or two-dimensional isotropic Heisenberg models.Phys. Rev. Lett.17, 1133–1136 (1966)
Hohenberg, P. C. Existence of long-range order in one and two dimensions.Phys. Rev.158, 383–386 (1967)
Sandvik, A. W. Stochastic series expansion method with operator-loop update.Phys. Rev. B59, R14157(R) (1999)
Acknowledgements
We thank R. Desbuquois, S. Dickerson, A. Eberlein, A. Kaufman, M. Messer, N. Prokov’ev, S. Sachdev, R. Scalettar, B. Svistunov, W. Zwerger, and M. Zwierlein and his research group for discussions. We thank S. Blatt, D. Cotta, S. Fölling, F. Huber, W. Setiawan and K. Wooley-Brown for early-stage contributions to the experiment. We acknowledge support from AFOSR (MURI), ARO (MURI, NDSEG), the Gordon and Betty Moore foundation EPiQS initiative, HQOC, NSF (CUA, ITAMP, GRFP, SAO) and SNSF.
Author information
Authors and Affiliations
Department of Physics, Harvard University, Cambridge, Massachusetts, USA
Anton Mazurenko, Christie S. Chiu, Geoffrey Ji, Maxwell F. Parsons, Márton Kanász-Nagy, Richard Schmidt, Fabian Grusdt, Eugene Demler, Daniel Greif & Markus Greiner
- Anton Mazurenko
You can also search for this author inPubMed Google Scholar
- Christie S. Chiu
You can also search for this author inPubMed Google Scholar
- Geoffrey Ji
You can also search for this author inPubMed Google Scholar
- Maxwell F. Parsons
You can also search for this author inPubMed Google Scholar
- Márton Kanász-Nagy
You can also search for this author inPubMed Google Scholar
- Richard Schmidt
You can also search for this author inPubMed Google Scholar
- Fabian Grusdt
You can also search for this author inPubMed Google Scholar
- Eugene Demler
You can also search for this author inPubMed Google Scholar
- Daniel Greif
You can also search for this author inPubMed Google Scholar
- Markus Greiner
You can also search for this author inPubMed Google Scholar
Contributions
A.M., C.S.C., G.J., M.F.P. and D.G. performed the experiment and analysed the data. G.J. carried out the determinant quantum Monte Carlo calculations forFig. 2e using the QUEST package. M.K.-N. developed the QMC code for the full-counting statistics and analysed the results together with R.S., F.G. and E.D. M.G. supervised the work. All authors contributed extensively to the writing of the manuscript and to discussions.
Corresponding author
Correspondence toMarkus Greiner.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Additional information
Reviewer InformationNature thanks T. Giamarchi and the other anonymous reviewer(s) for their contribution to the peer review of this work.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Figure 1 Amplitude of light fields applied to atoms.
a, The computed light field generated by the DMD, applied to the atoms for half-filled samples. A gradient compensates residual gradients in the lattice. The rim of the doughnut provides sharp walls for the inner subsystem. A small peak in the centre flattens the potential when combined with the optical lattice. The plot on the right shows a schematic of a radial cut of the potential, including the contribution of the lattice.b, The amplitude of the light field with an offset in the centre of the trap, used to dope the system with a finite population of holes.
Extended Data Figure 2 Average density profile in the system.
a, The average single-particle density map for a sample at half-filling shows a central region of uniform density, surrounded by a doughnut-shaped ring of low density. The dotted white circle indicates our system size, excluding edge effects.b, The azimuthal average of the single-particle densityns shown ina, for the system and for the inner edge of the doughnut where the density drops off to the reservoir density. The vertical dotted lines denote the boundary of the system.c, Azimuthal average of the single-particle densityns for three values of the hole dopingδ used in the experiment, indicating uniformity of atom number across our system to within 4%. The horizontal lines are at the system-wide average densities. Error bars inc are one standard deviation of the sample mean. The figure is based on 2,105 experimental realizations.
Extended Data Figure 3 Comparison of staggered magnetizations obtained directly through single-spin images and from spin correlations.
We calculate the corrected staggered magnetization from images with one spin state removed (main text). It can also be calculated from the spin correlator (Methods), with the two methods being identical in the limit of no noise and exactly one particle per site. Plotting these two quantities against each other, we find very good agreement with the liney = x (dotted line), indicating that any error due to deviation from one particle per site is small. The comparison is performed for the datasets used inFig. 2 (labelled temperature) andFig. 4 (labelled density). Error bars are computed as described in Methods.
Extended Data Figure 4 Alternative basis measurement.
We optionally apply a π/2 or π microwave pulse before the spin removal pulse and correlation measurement. The sign-corrected spin correlation functions (−1)iCd are insensitive to the presence and duration of this microwave pulse, consistent with an SU(2) symmetry of the state. The error bars are computed as described in Methods. This figure is based on 667 experimental realizations.
Extended Data Figure 5 Staggered magnetization obtained from spin correlations, with and without the nearest-neighbour contribution included.
To investigate the contributions to the corrected staggered magnetization at high dopingsδ, we consider the value calculated from the spin correlator (blue circles). We then omit the longest-range correlations, which have the greatest level of noise owing to the low number of pairs of sites extending across the cloud, as well as the nearest-neighbour correlations, which are essentially the only non-zero correlator outside of the antiferromagnetic phase (red circles). In the high-doping regime, we see that the greatest contribution to the staggered magnetization is the nearest-neighbour correlation, followed by the noisy longest-range correlations. Error bars are one standard deviation of the sample mean.
Rights and permissions
About this article
Cite this article
Mazurenko, A., Chiu, C., Ji, G.et al. A cold-atom Fermi–Hubbard antiferromagnet.Nature545, 462–466 (2017). https://doi.org/10.1038/nature22362
Received:
Accepted:
Published:
Issue Date: