Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Quantum phases of matter on a 256-atom programmable quantum simulator


Motivated by far-reaching applications ranging from quantum simulations of complex processes in physics and chemistry to quantum information processing1, a broad effort is currently underway to build large-scale programmable quantum systems. Such systems provide insights into strongly correlated quantum matter2,3,4,5,6, while at the same time enabling new methods for computation7,8,9,10 and metrology11. Here we demonstrate a programmable quantum simulator based on deterministically prepared two-dimensional arrays of neutral atoms, featuring strong interactions controlled by coherent atomic excitation into Rydberg states12. Using this approach, we realize a quantum spin model with tunable interactions for system sizes ranging from 64 to 256 qubits. We benchmark the system by characterizing high-fidelity antiferromagnetically ordered states and demonstrating quantum critical dynamics consistent with an Ising quantum phase transition in (2 + 1) dimensions13. We then create and study several new quantum phases that arise from the interplay between interactions and coherent laser excitation14, experimentally map the phase diagram and investigate the role of quantum fluctuations. Offering a new lens into the study of complex quantum matter, these observations pave the way for investigations of exotic quantum phases, non-equilibrium entanglement dynamics and hardware-efficient realization of quantum algorithms.

This is a preview of subscription content

Access options

Rent or Buy article

Get time limited or full article access on ReadCube.


All prices are NET prices.

Fig. 1: Programmable two-dimensional arrays of strongly interacting Rydberg atoms.
Fig. 2: Benchmarking of quantum simulator using chequerboard ordering.
Fig. 3: Observation of the (2 + 1)D Ising quantum phase transition on a 16 × 16 array.
Fig. 4: Phase diagram of the two-dimensional square lattice.
Fig. 5: Probing correlations and coherence in the striated phase via quench dynamics.

Data availability

The data that support the findings of this study are available from the corresponding author on reasonable request.


  1. 1.

    Preskill, J. Quantum computing in the NISQ era and beyond. Quantum 2, 79 (2018).

    Google Scholar 

  2. 2.

    Zhang, J. et al. Observation of a many-body dynamical phase transition with a 53-qubit quantum simulator. Nature 551, 601–604 (2017).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Gross, C. & Bloch, I. Quantum simulations with ultracold atoms in optical lattices. Science 357, 995–1001 (2017).

    ADS  CAS  Google Scholar 

  4. 4.

    Choi, J.-y. et al. Exploring the many-body localization transition in two dimensions. Science 352, 1547–1552 (2016).

    ADS  MathSciNet  CAS  PubMed  PubMed Central  MATH  Google Scholar 

  5. 5.

    Mazurenko, A. et al. A cold-atom Fermi–Hubbard antiferromagnet. Nature 545, 462–466 (2017).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Bernien, H. et al. Probing many-body dynamics on a 51-atom quantum simulator. Nature 551, 579–584 (2017).

    ADS  CAS  Google Scholar 

  7. 7.

    Neill, C. et al. Accurately computing electronic properties of a quantum ring. Preprint at (2021).

  8. 8.

    Wang, C. S. et al. Efficient multiphoton sampling of molecular vibronic spectra on a superconducting bosonic processor. Phys. Rev. X 10, 021060 (2020).

    ADS  CAS  Google Scholar 

  9. 9.

    Arute, F. et al. Quantum supremacy using a programmable superconducting processor. Nature 574, 505–510 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Zhong, H.-S. et al. Quantum computational advantage using photons. Science 370, 1460–1463 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Giovannetti, V., Lloyd, S. & Maccone, L. Advances in quantum metrology. Nat. Photon. 5, 222–229 (2011).

    ADS  CAS  Google Scholar 

  12. 12.

    Browaeys, A. & Lahaye, T. Many-body physics with individually controlled Rydberg atoms. Nat. Phys. 16, 132–142 (2020).

    CAS  Google Scholar 

  13. 13.

    Sachdev, S. Quantum Phase Transitions 2nd edn (Cambridge Univ. Press, 2011).

  14. 14.

    Samajdar, R., Ho, W. W., Pichler, H., Lukin, M. D. & Sachdev, S. Complex density wave orders and quantum phase transitions in a model of square-lattice Rydberg atom arrays. Phys. Rev. Lett. 124, 103601 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Pagano, G. et al. Quantum approximate optimization of the long-range Ising model with a trapped-ion quantum simulator. Proc. Natl Acad. Sci. USA 117, 25396–25401 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Friis, N. et al. Observation of entangled states of a fully controlled 20-qubit system. Phys. Rev. X 8, 021012 (2018).

    CAS  Google Scholar 

  17. 17.

    Harris, R. et al. Phase transitions in a programmable quantum spin glass simulator. Science 361, 162–165 (2018).

    ADS  MathSciNet  CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Labuhn, H. et al. Tunable two-dimensional arrays of single Rydberg atoms for realizing quantum Ising models. Nature 534, 667–670 (2016).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Morgado, M. & Whitlock, S. Quantum simulation and computing with Rydberg qubits. AVS Quantum Sci. 3, 023501 (2021).

    ADS  Google Scholar 

  20. 20.

    Lienhard, V. et al. Observing the space- and time-dependent growth of correlations in dynamically tuned synthetic Ising models with antiferromagnetic interactions. Phys. Rev. X 8, 021070 (2018).

    CAS  Google Scholar 

  21. 21.

    Guardado-Sanchez, E. et al. Probing the quench dynamics of antiferromagnetic correlations in a 2D quantum Ising spin system. Phys. Rev. X 8, 021069 (2018).

    CAS  Google Scholar 

  22. 22.

    Kim, H., Park, Y., Kim, K., Sim, H.-S. & Ahn, J. Detailed balance of thermalization dynamics in Rydberg-atom quantum simulators. Phys. Rev. Lett. 120, 180502 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  23. 23.

    de Léséleuc, S. et al. Observation of a symmetry-protected topological phase of interacting bosons with Rydberg atoms. Science 365, 775–780 (2019).

    ADS  MathSciNet  PubMed  PubMed Central  MATH  Google Scholar 

  24. 24.

    Keesling, A. et al. Quantum Kibble–Zurek mechanism and critical dynamics on a programmable Rydberg simulator. Nature 568, 207–211 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Madjarov, I. S. et al. High-fidelity entanglement and detection of alkaline-earth Rydberg atoms. Nat. Phys. 16, 857–861 (2020).

    CAS  Google Scholar 

  26. 26.

    Omran, A. et al. Generation and manipulation of Schrödinger cat states in Rydberg atom arrays. Science 365, 570–574 (2019).

    ADS  MathSciNet  CAS  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Graham, T. M. et al. Rydberg mediated entanglement in a two-dimensional neutral atom qubit array. Phys. Rev. Lett. 123, 230501 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Levine, H. et al. Parallel implementation of high-fidelity multiqubit gates with neutral atoms. Phys. Rev. Lett. 123, 170503 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Madjarov, I. S. et al. An atomic-array optical clock with single-atom readout. Phys. Rev. X 9, 041052 (2019).

    CAS  Google Scholar 

  30. 30.

    Young, A. W. et al. Half-minute-scale atomic coherence and high relative stability in a tweezer clock. Nature 588, 408–413 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Schymik, K.-N. et al. Enhanced atom-by-atom assembly of arbitrary tweezers arrays. Phys. Rev. A 102, 063107 (2020).

    ADS  CAS  Google Scholar 

  32. 32.

    Kumar, A., Wu, T.-Y., Giraldo, F. & Weiss, D. S. Sorting ultracold atoms in a three-dimensional optical lattice in a realization of Maxwell’s demon. Nature 561, 83–87 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Ohl de Mello, D. et al. Defect-free assembly of 2D clusters of more than 100 single-atom quantum systems. Phys. Rev. Lett. 122, 203601 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Barredo, D., Lienhard, V., de Léséleuc, S., Lahaye, T. & Browaeys, A. Synthetic three-dimensional atomic structures assembled atom by atom. Nature 561, 79–82 (2018).

    ADS  CAS  Google Scholar 

  35. 35.

    Barredo, D., de Léséleuc, S., Lienhard, V., Lahaye, T. & Browaeys, A. An atom-by-atom assembler of defect-free arbitrary 2D atomic arrays. Science 354, 1021–1023 (2016).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Jaksch, D. et al. Fast quantum gates for neutral atoms. Phys. Rev. Lett. 85, 2208 (2000).

    ADS  CAS  Google Scholar 

  37. 37.

    Lukin, M. D. et al. Dipole blockade and quantum information processing in mesoscopic atomic ensembles. Phys. Rev. Lett. 87, 037901 (2001).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Gardas, B., Dziarmaga, J., Zurek, W. H. & Zwolak, M. Defects in quantum computers. Sci. Rep. 8, 4539 (2018).

    ADS  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Zurek, W. H., Dorner, U. & Zoller, P. Dynamics of a quantum phase transition. Phys. Rev. Lett. 95, 105701 (2005).

    ADS  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Felser, T., Notarnicola, S. & Montangero, S. Efficient tensor network ansatz for high-dimensional quantum many-body problems. Phys. Rev. Lett. 126, 170603 (2021).

    ADS  MathSciNet  CAS  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Huang, H.-Y., Kueng, R. & Preskill, J. Predicting many properties of a quantum system from very few measurements. Nat. Phys. 16, 1050–1057 (2020).

    CAS  Google Scholar 

  42. 42.

    Verresen, R., Lukin, M. D. & Vishwanath, A. Prediction of toric code topological order from Rydberg blockade. Preprint at (2020).

  43. 43.

    Samajdar, R., Ho, W. W., Pichler, H., Lukin, M. D. & Sachdev, S. Quantum phases of Rydberg atoms on a kagome lattice. Proc. Natl Acad. Sci. USA 118, e2015785118 (2021).

    CAS  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Semeghini, G. et al. Probing topological spin liquids on a programmable quantum simulator. Preprint at (2021).

  45. 45.

    Saskin, S., Wilson, J. T., Grinkemeyer, B. & Thompson, J. D. Narrow-line cooling and imaging of ytterbium atoms in an optical tweezer array. Phys. Rev. Lett. 122, 143002 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Anderegg, L. et al. An optical tweezer array of ultracold molecules. Science 365, 1156–1158 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Liu, L. R. et al. Building one molecule from a reservoir of two atoms. Science 360, 900–903 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Turner, C. J., Michailidis, A. A., Abanin, D. A., Serbyn, M. & Papić, Z. Weak ergodicity breaking from quantum many-body scars. Nat. Phys. 14, 745–749 (2018).

    CAS  Google Scholar 

  49. 49.

    Surace, F. M. et al. Lattice gauge theories and string dynamics in Rydberg atom quantum simulators. Phys. Rev. X 10, 021041 (2020).

    CAS  Google Scholar 

  50. 50.

    Bluvstein, D. et al. Controlling quantum many-body dynamics in driven Rydberg atom arrays. Science 371, 1355–1359 (2021).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Diehl, H. W. The theory of boundary critical phenomena. Int. J. Mod. Phys. B 11, 3503–3523 (1997).

    ADS  Google Scholar 

  52. 52.

    Bañuls, M. C. et al. Simulating lattice gauge theories within quantum technologies. Eur. Phys. J. D 74, 165 (2020).

    ADS  Google Scholar 

  53. 53.

    Notarnicola, S., Collura, M. & Montangero, S. Real-time-dynamics quantum simulation of (1 + 1)-dimensional lattice QED with Rydberg atoms. Phys. Rev. Research 2, 013288 (2020).

    ADS  CAS  Google Scholar 

  54. 54.

    Weimer, H., Müller, M., Lesanovsky, I., Zoller, P. & Büchler, H. P. A Rydberg quantum simulator. Nat. Phys. 6, 382–388 (2010).

    CAS  Google Scholar 

  55. 55.

    Auger, J. M., Bergamini, S. & Browne, D. E. Blueprint for fault-tolerant quantum computation with Rydberg atoms. Phys. Rev. A 96, 052320 (2017).

    Google Scholar 

  56. 56.

    Farhi, E., Goldstone, J. & Gutmann, S. A quantum approximate optimization algorithm. Preprint at (2014).

  57. 57.

    Zhou, L., Wang, S.-T., Choi, S., Pichler, H. & Lukin, M. D. Quantum approximate optimization algorithm: performance, mechanism, and implementation on near-term devices. Phys. Rev. X 10, 021067 (2020).

    CAS  Google Scholar 

  58. 58.

    Wild, D. S., Sels, D., Pichler, H. & Lukin, M. D. Quantum sampling algorithms for near-term devices. Preprint at (2020).

  59. 59.

    Scholl, P. et al. Programmable quantum simulation of 2D antiferromagnets with hundreds of Rydberg atoms. Nature (2021).

  60. 60.

    Kim, D. et al. Large-scale uniform optical focus array generation with a phase spatial light modulator. Opt. Lett. 44, 3178–3181 (2019).

    ADS  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Endres, M. et al. Atom-by-atom assembly of defect-free one-dimensional cold atom arrays. Science 354, 1024–1027 (2016).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Lee, W., Kim, H. & Ahn, J. Defect-free atomic array formation using the Hungarian matching algorithm. Phys. Rev. A 95, 053424 (2017).

    ADS  Google Scholar 

  63. 63.

    Sheng, C. et al. Efficient preparation of 2D defect-free atom arrays with near-fewest sorting-atom moves. Phys. Rev. Research 3, 023008 (2021).

    ADS  Google Scholar 

  64. 64.

    Levine, H. et al. High-fidelity control and entanglement of Rydberg-atom qubits. Phys. Rev. Lett. 121, 123603 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Bowman, D. et al. High-fidelity phase and amplitude control of phase-only computer generated holograms using conjugate gradient minimisation. Opt. Express 25, 11692–11700 (2017).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Zupancic, P. et al. Ultra-precise holographic beam shaping for microscopic quantum control. Opt. Express 24, 13881–13893 (2016).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  67. 67.

    Beterov, I. I., Ryabtsev, I. I., Tretyakov, D. B. & Entin, V. M. Quasiclassical calculations of blackbody-radiation-induced depopulation rates and effective lifetimes of Rydberg nS, nP, and nD alkali-metal atoms with n ≤ 80. Phys. Rev. A 79, 052504 (2009).

    ADS  Google Scholar 

  68. 68.

    Bhattacharjee, S. M. & Seno, F. A measure of data collapse for scaling. J. Phys. Math. Gen. 34, 6375 (2001).

    ADS  CAS  MATH  Google Scholar 

  69. 69.

    Hasenbusch, M. Monte Carlo study of surface critical phenomena: The special point. Phys. Rev. B 84, 134405 (2011).

    ADS  Google Scholar 

  70. 70.

    White, S. R. Density matrix formulation for quantum renormalization groups. Phys. Rev. Lett. 69, 2863 (1992).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  71. 71.

    White, S. R. Density-matrix algorithms for quantum renormalization groups. Phys. Rev. B 48, 10345 (1993).

    ADS  CAS  Google Scholar 

  72. 72.

    Stoudenmire, E. M. & White, S. R. Studying two-dimensional systems with the density matrix renormalization group. Annu. Rev. Condens. Matter Phys. 3, 111–128 (2012).

    Google Scholar 

  73. 73.

    Fishman, M., White, S. R. & Stoudenmire, E. M. The ITensor software library for tensor network calculations. Preprint at (2020).

Download references


We thank many members of the Harvard AMO community, particularly E. Urbach, S. Dakoulas and J. Doyle for enabling safe and productive operation of our laboratories during 2020. We thank H. Bernien, D. Englund, M. Endres, N. Gemelke, D. Kim, P. Stark and A. Zibrov for discussions and experimental help. We acknowledge financial support from the Center for Ultracold Atoms, the US National Science Foundation, the Vannevar Bush Faculty Fellowship, the US Department of Energy (DE-SC0021013 and DOE Quantum Systems Accelerator Center, contract no. 7568717), the Office of Naval Research, the Army Research Office MURI and the DARPA ONISQ programme. T.T.W. acknowledges support from Gordon College. H.L. acknowledges support from the National Defense Science and Engineering Graduate (NDSEG) fellowship. G.S. acknowledges support from a fellowship from the Max Planck/Harvard Research Center for Quantum Optics. D.B. acknowledges support from the NSF Graduate Research Fellowship Program (grant DGE1745303) and the Fannie and John Hertz Foundation. W.W.H. is supported by the Moore Foundation’s EPiQS Initiative grant no. GBMF4306, the NUS development grant AY2019/2020, and the Stanford Institute of Theoretical Physics. S.C. acknowledges support from the Miller Institute for Basic Research in Science. R.S. and S.S. were supported by the US Department of Energy under grant DE-SC0019030. The DMRG calculations were performed using the ITensor Library73. The computations in this paper were run on the FASRC Cannon cluster supported by the FAS Division of Science Research Computing Group at Harvard University.

Author information




S.E., T.T.W., H.L., A.K., G.S, A.O. and D.B. contributed to building the experimental set-up, performed the measurements and analysed the data. Theoretical analysis was performed by R.S., H.P., W.W.H. and S.C. All work was supervised by S.S., M.G., V.V. and M.D.L. All authors discussed the results and contributed to the manuscript.

Corresponding author

Correspondence to Mikhail D. Lukin.

Ethics declarations

Competing interests

M.G., V.V. and M.D.L. are co-founders and shareholders of QuEra Computing. A.K. and A.O. are shareholders of QuEra Computing. All other authors declare no competing interests.

Additional information

Peer review information Nature thanks the 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 Fig. 1 Large arrays of optical tweezers.

The experimental platform produces optical tweezer arrays with up to ~1,000 tweezers and ~50% loading probability per tweezer after 100 ms of magneto-optical trap loading time. a, Camera image of an array of 34 × 30 tweezers (1,020 traps), including aberration correction. b, Sample image of random loading into this tweezer array, with 543 loaded atoms. Atoms are detected on an EMCCD camera with fluorescence imaging.

Extended Data Fig. 2 Correcting for aberrations in the SLM tweezer array.

The aberration correction procedure uses the orthogonality of Zernike polynomials and the fact that correcting aberrations increases tweezer light shifts on the atoms. To independently measure and correct each aberration type, Zernike polynomials are added with variable amplitude to the SLM phase hologram, with values optimized to maximize tweezer light shifts. a, Two common aberration types: horizontal coma (upper) and primary spherical (lower), for which ~50 milliwaves compensation on each reduces aberrations and results in higher-depth traps. b, Correcting for aberrations associated with the 13 lowest-order Zernike polynomials. The sum of all polynomials with their associated coefficients gives the total wavefront correction (RMS ~70 milliwaves) for our optical system, which is summed with the optical tweezer hologram on the SLM. c, Trap depths across a 26 × 13 trap array before and after correction for aberrations. Aberration correction results in tighter focusing (higher trap light shift) and improved homogeneity. Trap depths are measured by probing the light shift of each trap on the \(|5{S}_{1/2},F=2\rangle \to |5{P}_{3/2},F{\prime} =2\rangle \) transition. d, Aberration correction also results in higher and more homogeneous trap frequencies across the array. Trap frequencies are measured by modulating tweezer depths at variable frequencies, resulting in parametric heating and atom loss when the modulation frequency is twice the radial trap frequency. The measurement after correction for aberrations shows a narrower spectrum and higher trap frequencies (averaged over the whole array).

Extended Data Fig. 3 Rearrangement protocol.

a, Sample sequence of individual rearrangement steps. There are two pre-sorting moves (1, 2). Move 3 is the single ejection move. Moves 4–6 consist of parallel vertical sorting within each column, including both upward and downward moves. The upper panel illustrates the frequency spectrum of the waveform in the vertical and horizontal AODs during these moves, with the underlying grid corresponding to the calibrated frequencies that map to SLM array rows and columns. b, Spectrograms representing the horizontal and vertical AOD waveforms over the duration of a single vertical frequency scan during a realistic rearrangement procedure for a 26 × 13 array. The heat-maps show frequency spectra of the AOD waveforms over small time intervals during the scan.

Extended Data Fig. 4 Generating homogeneous Rydberg beams.

a, Measured Gaussian-beam illumination on the SLM for shaping the 420-nm Rydberg beam. A Gaussian fit to these data is used as an input for the hologram optimization algorithm. b, Measured wavefront error through our optical system (after correction), showing a reduction of aberrations to λ/100. c, Computer-generated hologram for creating the 420-nm top-hat beam. d, Measured light intensity of the 420-nm top-hat beam (top), and the cross-section along where atoms will be positioned (bottom). Vertical lines denote the 105-μm region where the beam should be flat. e, Using the measured top-hat intensity, a phase correction is calculated for adding to the initial hologram. f, Resulting top-hat beam after feedback shows considerably improved homogeneity. pk–pk, peak to peak.

Extended Data Fig. 5 Characterizing microwave-enhanced Rydberg detection fidelity.

The effect of strong microwave (MW) pulses on Rydberg atoms is measured by preparing atoms in \(|g\rangle \), exciting to \(|r\rangle \) with a Rydberg π-pulse, and then applying the microwave pulse before de-exciting residual Rydberg atoms with a final Rydberg π-pulse. (The entire sequence occurs while tweezers are briefly turned off.) a, Broad resonances are observed with varying microwave frequency, corresponding to transitions from \(|r\rangle \) = \(|70S\rangle \) to other Rydberg states. Note that the transitions to \(|69P\rangle \) and \(|70P\rangle \) are in the range of 10–12 GHz, and over this entire range there is strong transfer out of  \(|r\rangle \). Other resonances might be due to multiphoton effects. b, With fixed 6.9-GHz microwave frequency and varying pulse time, there is a rapid transfer out of the Rydberg state on the timescale of several nanoseconds. Over short timescales, there may be coherent oscillations that return population back to \(|r\rangle \), so a 100-ns pulse is used for enhancement of loss signal of \(|r\rangle \) in the experiment.

Extended Data Fig. 6 Coarse-grained local staggered magnetization.

a, Examples of Rydberg populations ni after a faster (top) and slower (bottom) linear sweep. b, Corresponding coarse-grained local staggered magnetizations mi clearly show larger extents of antiferromagnetically ordered domains (dark blue or dark red) for the slower sweep (bottom) than for the faster sweep (top), as expected from the Kibble–Zurek mechanism. c, Isotropic correlation functions \({G}_{m}^{(2)}\) for the corresponding coarse-grained local staggered magnetizations after a faster (top) or a slower (bottom) sweep. d, As a function of radial distance, correlations \({G}_{m}^{(2)}\) decay exponentially with a length scale corresponding to the correlation length ξ. The two decay curves correspond to faster (orange) and slower (blue) sweeps.

Extended Data Fig. 7 Extracting the quantum critical point.

a, The mean Rydberg excitation density \(\langle n\rangle \) versus detuning Δ/Ω on a 16 × 16 array. The data are fitted within a window (dashed lines) to a cubic polynomial (red curve) as a means of smoothing the data. b, The peak in the numerical derivative of the fitted data (red curve) corresponds to the critical point Δc/Ω = 1.12(4) (red shaded regions show uncertainty ranges, obtained from varying fit windows). In contrast, the point-by-point slope of the data (grey) is too noisy to be useful. c, Order parameter \(\mathop{ {\mathcal F} }\limits^{ \sim }({\rm{\pi }},{\rm{\pi }})\) for the chequerboard phase versus Δ/Ω measured on a 16 × 16 array with the value of the critical point from b superimposed (red line), showing the clear growth of the order parameter after the critical point. d, DMRG simulations of  \(\langle n\rangle \) versus Δ/Ω on a 10 × 10 array. For comparison against the experimental fitting procedure, the data from numerics are also fitted to a cubic polynomial within the indicated window (dashed lines). e, The point-by-point slope of the numerical data (blue curve) has a peak at Δc/Ω = 1.18 (blue dashed line), in good agreement with the results (red dashed line) from both the numerical derivative of the cubic fit on the same data (red curve) and the result of the experiment. f, DMRG simulation of \(\mathop{ {\mathcal F} }\limits^{ \sim }({\rm{\pi }},{\rm{\pi }})\) versus Δ/Ω, with the exact quantum critical point from numerics shown (red line).

Extended Data Fig. 8 Optimization of data collapse.

a, Distance D between rescaled correlation length \(\mathop{\xi }\limits^{ \sim }\) versus \(\mathop{\varDelta }\limits^{ \sim }\) curves depends on both the location of the quantum critical point Δc/Ω and on the correlation length critical exponent ν. The independently determined Δc/Ω (blue line, with uncertainty range in grey) and the experimentally extracted value of ν (dashed red line, with uncertainty range corresponding to the red shaded region) are marked on the plot. b, Our determination of ν (red) from data collapse around the independently determined Δc/Ω (blue) is consistent across arrays of different sizes. ce, Data collapse is clearly better at the experimentally determined value (ν = 0.62) as compared with the mean-field (ν = 0.5) or the (1 + 1)D (ν = 1) values. The horizontal extent of the data corresponds to the region of overlap of all rescaled datasets.

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ebadi, S., Wang, T.T., Levine, H. et al. Quantum phases of matter on a 256-atom programmable quantum simulator. Nature 595, 227–232 (2021).

Download citation

Further reading


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links