Showing 1–48 of 48 results for author: Carrasquilla, J

Observation of a non-Hermitian supersonic mode

Authors: Yuxuan Zhang, Juan Carrasquilla, Yong Baek Kim

Abstract: Quantum computers have long been anticipated to excel in simulating quantum many-body physics. While most previous work has focused on Hermitian physics, we demonstrate the power of variational quantum circuits for resource-efficient simulations of dynamical and equilibrium physics in non-Hermitian systems, revealing new phenomena beyond standard Hermitian quantum machines. Using a variational qua… ▽ More Quantum computers have long been anticipated to excel in simulating quantum many-body physics. While most previous work has focused on Hermitian physics, we demonstrate the power of variational quantum circuits for resource-efficient simulations of dynamical and equilibrium physics in non-Hermitian systems, revealing new phenomena beyond standard Hermitian quantum machines. Using a variational quantum compilation scheme for fermionic systems, we reduce gate count, save qubits, and eliminate the need for postselection, a major challenge in simulating non-Hermitian dynamics via standard Trotterization. Experimentally, we observed a supersonic mode in the connected density-density correlation function on an $ n = 18 $ fermionic chain after a non-Hermitian, locally interacting quench, which would otherwise be forbidden by the Lieb-Robinson bound in a Hermitian system. Additionally, we investigate sequential quantum circuits generated by tensor networks for ground state preparation, here defined as the eigenstate with the lowest real part eigenvalue, using a variance minimization scheme. Through a trapped-ion implementation on the Quantinuum H1 quantum processor, we accurately capture correlation functions and energies across an exceptional point on a dissipative spin chain up to length $ n = 20 $ using only 3 qubits. Motivated by these advancements, we provide an analytical example demonstrating that simulating single-qubit non-Hermitian dynamics for $Θ(\log(n))$ time from certain initial states is exponentially hard on a quantum computer, offering insights into the opportunities and limitations of using quantum computation for simulating non-Hermitian physics. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Recurrent neural network wave functions for Rydberg atom arrays on kagome lattice

Authors: Mohamed Hibat-Allah, Ejaaz Merali, Giacomo Torlai, Roger G Melko, Juan Carrasquilla

Abstract: Rydberg atom array experiments have demonstrated the ability to act as powerful quantum simulators, preparing strongly-correlated phases of matter which are challenging to study for conventional computer simulations. A key direction has been the implementation of interactions on frustrated geometries, in an effort to prepare exotic many-body states such as spin liquids and glasses. In this paper,… ▽ More Rydberg atom array experiments have demonstrated the ability to act as powerful quantum simulators, preparing strongly-correlated phases of matter which are challenging to study for conventional computer simulations. A key direction has been the implementation of interactions on frustrated geometries, in an effort to prepare exotic many-body states such as spin liquids and glasses. In this paper, we apply two-dimensional recurrent neural network (RNN) wave functions to study the ground states of Rydberg atom arrays on the kagome lattice. We implement an annealing scheme to find the RNN variational parameters in regions of the phase diagram where exotic phases may occur, corresponding to rough optimization landscapes. For Rydberg atom array Hamiltonians studied previously on the kagome lattice, our RNN ground states show no evidence of exotic spin liquid or emergent glassy behavior. In the latter case, we argue that the presence of a non-zero Edwards-Anderson order parameter is an artifact of the long autocorrelations times experienced with quantum Monte Carlo simulations. This result emphasizes the utility of autoregressive models, such as RNNs, to explore Rydberg atom array physics on frustrated lattices and beyond. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 13 pages, 5 figures, 3 tables. Link to GitHub repository: https://github.com/mhibatallah/RNNWavefunctions

The statistical mechanics and machine learning of the $α$-Rényi ensemble

Authors: Andrew Jreissaty, Juan Carrasquilla

Abstract: We study the statistical physics of the classical Ising model in the so-called $α$-Rényi ensemble, a finite-temperature thermal state approximation that minimizes a modified free energy based on the $α$-Rényi entropy. We begin by characterizing its critical behavior in mean-field theory in different regimes of the Rényi index $α$. Next, we re-introduce correlations and consider the model in one an… ▽ More We study the statistical physics of the classical Ising model in the so-called $α$-Rényi ensemble, a finite-temperature thermal state approximation that minimizes a modified free energy based on the $α$-Rényi entropy. We begin by characterizing its critical behavior in mean-field theory in different regimes of the Rényi index $α$. Next, we re-introduce correlations and consider the model in one and two dimensions, presenting an exact analysis of the former and devising an unconventional Monte Carlo approach to the study of the latter. Remarkably, we find that while mean-field predicts a continuous phase transition below a threshold index value of $α\sim 1.303$ and a first-order transition above it, the Monte Carlo results in two dimensions point to a continuous transition at all $α$. We conclude by performing a variational minimization of the $α$-Rényi free energy using a recurrent neural network (RNN) ansatz where we find that the RNN performs well in two dimensions when compared to the Monte Carlo simulations. Our work highlights the potential opportunities and limitations associated with the use of the $α$-Rényi ensemble formalism in probing the thermodynamic equilibrium properties of classical and quantum systems. △ Less

Submitted 18 April, 2024; v1 submitted 5 April, 2024; originally announced April 2024.

Computational supremacy in quantum simulation

Authors: Andrew D. King, Alberto Nocera, Marek M. Rams, Jacek Dziarmaga, Roeland Wiersema, William Bernoudy, Jack Raymond, Nitin Kaushal, Niclas Heinsdorf, Richard Harris, Kelly Boothby, Fabio Altomare, Andrew J. Berkley, Martin Boschnak, Kevin Chern, Holly Christiani, Samantha Cibere, Jake Connor, Martin H. Dehn, Rahul Deshpande, Sara Ejtemaee, Pau Farré, Kelsey Hamer, Emile Hoskinson, Shuiyuan Huang , et al. (37 additional authors not shown)

Abstract: Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. Establishing this capability, especially for impactful and meaningful problems, remains a central challenge. One such problem is the simulation of nonequilibrium dynamics of a magnetic spin system quenched through a quantum phase transition. State-of-the-art classical simulations dem… ▽ More Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. Establishing this capability, especially for impactful and meaningful problems, remains a central challenge. One such problem is the simulation of nonequilibrium dynamics of a magnetic spin system quenched through a quantum phase transition. State-of-the-art classical simulations demand resources that grow exponentially with system size. Here we show that superconducting quantum annealing processors can rapidly generate samples in close agreement with solutions of the Schrödinger equation. We demonstrate area-law scaling of entanglement in the model quench in two-, three- and infinite-dimensional spin glasses, supporting the observed stretched-exponential scaling of effort for classical approaches. We assess approximate methods based on tensor networks and neural networks and conclude that no known approach can achieve the same accuracy as the quantum annealer within a reasonable timeframe. Thus quantum annealers can answer questions of practical importance that classical computers cannot. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Neural network approach to quasiparticle dispersions in doped antiferromagnets

Authors: Hannah Lange, Fabian Döschl, Juan Carrasquilla, Annabelle Bohrdt

Abstract: Numerically simulating spinful, fermionic systems is of great interest from the perspective of condensed matter physics. However, the exponential growth of the Hilbert space dimension with system size renders an exact parameterization of large quantum systems prohibitively demanding. This is a perfect playground for neural networks, owing to their immense representative power that often allows to… ▽ More Numerically simulating spinful, fermionic systems is of great interest from the perspective of condensed matter physics. However, the exponential growth of the Hilbert space dimension with system size renders an exact parameterization of large quantum systems prohibitively demanding. This is a perfect playground for neural networks, owing to their immense representative power that often allows to use only a fraction of the parameters that are needed in exact methods. Here, we investigate the ability of neural quantum states (NQS) to represent the bosonic and fermionic $t-J$ model - the high interaction limit of the Fermi-Hubbard model - on different 1D and 2D lattices. Using autoregressive recurrent neural networks (RNNs) with 2D tensorized gated recurrent units, we study the ground state representations upon doping the half-filled system with holes. Moreover, we present a method to calculate dispersion relations from the neural network state representation, applicable to any neural network architecture and any lattice geometry, that allows to infer the low-energy physics from the NQS. To demonstrate our approach, we calculate the dispersion of a single hole in the $t-J$ model on different 1D and 2D square and triangular lattices. Furthermore, we analyze the strengths and weaknesses of the RNN approach for fermionic systems, pointing the way for an accurate and efficient parameterization of fermionic quantum systems using neural quantum states. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 14 pages + SM, 12 figures

A Framework for Demonstrating Practical Quantum Advantage: Racing Quantum against Classical Generative Models

Authors: Mohamed Hibat-Allah, Marta Mauri, Juan Carrasquilla, Alejandro Perdomo-Ortiz

Abstract: Generative modeling has seen a rising interest in both classical and quantum machine learning, and it represents a promising candidate to obtain a practical quantum advantage in the near term. In this study, we build over a proposed framework for evaluating the generalization performance of generative models, and we establish the first quantitative comparative race towards practical quantum advant… ▽ More Generative modeling has seen a rising interest in both classical and quantum machine learning, and it represents a promising candidate to obtain a practical quantum advantage in the near term. In this study, we build over a proposed framework for evaluating the generalization performance of generative models, and we establish the first quantitative comparative race towards practical quantum advantage (PQA) between classical and quantum generative models, namely Quantum Circuit Born Machines (QCBMs), Transformers (TFs), Recurrent Neural Networks (RNNs), Variational Autoencoders (VAEs), and Wasserstein Generative Adversarial Networks (WGANs). After defining four types of PQAs scenarios, we focus on what we refer to as potential PQA, aiming to compare quantum models with the best-known classical algorithms for the task at hand. We let the models race on a well-defined and application-relevant competition setting, where we illustrate and demonstrate our framework on 20 variables (qubits) generative modeling task. Our results suggest that QCBMs are more efficient in the data-limited regime than the other state-of-the-art classical generative models. Such a feature is highly desirable in a wide range of real-world applications where the available data is scarce. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 17 pages, 5 figures, 3 tables

Investigating Topological Order using Recurrent Neural Networks

Authors: Mohamed Hibat-Allah, Roger G. Melko, Juan Carrasquilla

Abstract: Recurrent neural networks (RNNs), originally developed for natural language processing, hold great promise for accurately describing strongly correlated quantum many-body systems. Here, we employ 2D RNNs to investigate two prototypical quantum many-body Hamiltonians exhibiting topological order. Specifically, we demonstrate that RNN wave functions can effectively capture the topological order of t… ▽ More Recurrent neural networks (RNNs), originally developed for natural language processing, hold great promise for accurately describing strongly correlated quantum many-body systems. Here, we employ 2D RNNs to investigate two prototypical quantum many-body Hamiltonians exhibiting topological order. Specifically, we demonstrate that RNN wave functions can effectively capture the topological order of the toric code and a Bose-Hubbard spin liquid on the kagome lattice by estimating their topological entanglement entropies. We also find that RNNs favor coherent superpositions of minimally-entangled states over minimally-entangled states themselves. Overall, our findings demonstrate that RNN wave functions constitute a powerful tool to study phases of matter beyond Landau's symmetry-breaking paradigm. △ Less

Submitted 25 October, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: 15 pages, 6 figures, 2 tables. Published version in Physical Review B

Journal ref: Phys. Rev. B 108, 075152, August 2023

Variational Benchmarks for Quantum Many-Body Problems

Authors: Dian Wu, Riccardo Rossi, Filippo Vicentini, Nikita Astrakhantsev, Federico Becca, Xiaodong Cao, Juan Carrasquilla, Francesco Ferrari, Antoine Georges, Mohamed Hibat-Allah, Masatoshi Imada, Andreas M. Läuchli, Guglielmo Mazzola, Antonio Mezzacapo, Andrew Millis, Javier Robledo Moreno, Titus Neupert, Yusuke Nomura, Jannes Nys, Olivier Parcollet, Rico Pohle, Imelda Romero, Michael Schmid, J. Maxwell Silvester, Sandro Sorella , et al. (8 additional authors not shown)

Abstract: The continued development of novel many-body approaches to ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. Here we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide the most extensive curated dataset of variational calculations of many-body quantum systems to date, ide… ▽ More The continued development of novel many-body approaches to ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. Here we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide the most extensive curated dataset of variational calculations of many-body quantum systems to date, identifying cases where state-of-the-art numerical approaches show limited accuracy, and novel algorithms or computational platforms, such as quantum computing, could provide improved accuracy. The V-score can be used as a metric to assess the progress of quantum variational methods towards quantum advantage for ground-state problems, especially in regimes where classical verifiability is impossible. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 25 pages, 5 figures

Kagome qubit ice

Authors: Alejandro Lopez-Bezanilla, Jack Raymond, Kelly Boothby, Juan Carrasquilla, Cristiano Nisoli, Andrew D. King

Abstract: Topological phases of spin liquids with constrained disorder can host a kinetics of fractionalized excitations. However, spin-liquid phases with distinct kinetic regimes have proven difficult to observe experimentally. Here we present a realization of kagome spin ice in the superconducting qubits of a quantum annealer, and use it to demonstrate a field-induced kinetic crossover between spin-liquid… ▽ More Topological phases of spin liquids with constrained disorder can host a kinetics of fractionalized excitations. However, spin-liquid phases with distinct kinetic regimes have proven difficult to observe experimentally. Here we present a realization of kagome spin ice in the superconducting qubits of a quantum annealer, and use it to demonstrate a field-induced kinetic crossover between spin-liquid phases. Employing fine control over local magnetic fields, we show evidence of both the Ice-I phase and an unconventional field-induced Ice-II phase. In the latter, a charge-ordered yet spin-disordered topological phase, the kinetics proceeds via pair creation and annihilation of strongly correlated, charge conserving, fractionalized excitations. As these kinetic regimes have resisted characterization in other artificial spin ice realizations, our results demonstrate the utility of quantum-driven kinetics in advancing the study of topological phases of spin liquids. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Supplementing Recurrent Neural Network Wave Functions with Symmetry and Annealing to Improve Accuracy

Authors: Mohamed Hibat-Allah, Roger G. Melko, Juan Carrasquilla

Abstract: Recurrent neural networks (RNNs) are a class of neural networks that have emerged from the paradigm of artificial intelligence and has enabled lots of interesting advances in the field of natural language processing. Interestingly, these architectures were shown to be powerful ansatze to approximate the ground state of quantum systems. Here, we build over the results of [Phys. Rev. Research 2, 023… ▽ More Recurrent neural networks (RNNs) are a class of neural networks that have emerged from the paradigm of artificial intelligence and has enabled lots of interesting advances in the field of natural language processing. Interestingly, these architectures were shown to be powerful ansatze to approximate the ground state of quantum systems. Here, we build over the results of [Phys. Rev. Research 2, 023358 (2020)] and construct a more powerful RNN wave function ansatz in two dimensions. We use symmetry and annealing to obtain accurate estimates of ground state energies of the two-dimensional (2D) Heisenberg model, on the square lattice and on the triangular lattice. We show that our method is superior to Density Matrix Renormalisation Group (DMRG) for system sizes larger than or equal to $14 \times 14$ on the triangular lattice. △ Less

Submitted 12 January, 2024; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: 11 pages, 4 figures, 1 table. Corrected typos. Originally published in Machine Learning and the Physical Sciences Workshop (NeurIPS 2021), see: https://ml4physicalsciences.github.io/2021/files/NeurIPS_ML4PS_2021_92.pdf. Our reproducibility code can be found at https://github.com/mhibatallah/RNNWavefunctions

Journal ref: Machine Learning and the Physical Sciences, NeurIPS 2021

Modern applications of machine learning in quantum sciences

Authors: Anna Dawid, Julian Arnold, Borja Requena, Alexander Gresch, Marcin Płodzień, Kaelan Donatella, Kim A. Nicoli, Paolo Stornati, Rouven Koch, Miriam Büttner, Robert Okuła, Gorka Muñoz-Gil, Rodrigo A. Vargas-Hernández, Alba Cervera-Lierta, Juan Carrasquilla, Vedran Dunjko, Marylou Gabrié, Patrick Huembeli, Evert van Nieuwenburg, Filippo Vicentini, Lei Wang, Sebastian J. Wetzel, Giuseppe Carleo, Eliška Greplová, Roman Krems , et al. (4 additional authors not shown)

Abstract: In this book, we provide a comprehensive introduction to the most recent advances in the application of machine learning methods in quantum sciences. We cover the use of deep learning and kernel methods in supervised, unsupervised, and reinforcement learning algorithms for phase classification, representation of many-body quantum states, quantum feedback control, and quantum circuits optimization.… ▽ More In this book, we provide a comprehensive introduction to the most recent advances in the application of machine learning methods in quantum sciences. We cover the use of deep learning and kernel methods in supervised, unsupervised, and reinforcement learning algorithms for phase classification, representation of many-body quantum states, quantum feedback control, and quantum circuits optimization. Moreover, we introduce and discuss more specialized topics such as differentiable programming, generative models, statistical approach to machine learning, and quantum machine learning. △ Less

Submitted 15 November, 2023; v1 submitted 8 April, 2022; originally announced April 2022.

Comments: 288 pages, 92 figures. We have a publishing contract with Cambridge University Press. Figures and tex files are available at https://github.com/Shmoo137/Lecture-Notes

Measurement-induced entanglement phase transitions in variational quantum circuits

Authors: Roeland Wiersema, Cunlu Zhou, Juan Felipe Carrasquilla, Yong Baek Kim

Abstract: Variational quantum algorithms (VQAs), which classically optimize a parametrized quantum circuit to solve a computational task, promise to advance our understanding of quantum many-body systems and improve machine learning algorithms using near-term quantum computers. Prominent challenges associated with this family of quantum-classical hybrid algorithms are the control of quantum entanglement and… ▽ More Variational quantum algorithms (VQAs), which classically optimize a parametrized quantum circuit to solve a computational task, promise to advance our understanding of quantum many-body systems and improve machine learning algorithms using near-term quantum computers. Prominent challenges associated with this family of quantum-classical hybrid algorithms are the control of quantum entanglement and quantum gradients linked to their classical optimization. Known as the barren plateau phenomenon, these quantum gradients may rapidly vanish in the presence of volume-law entanglement growth, which poses a serious obstacle to the practical utility of VQAs. Inspired by recent studies of measurement-induced entanglement transition in random circuits, we investigate the entanglement transition in variational quantum circuits endowed with intermediate projective measurements. Considering the Hamiltonian Variational Ansatz (HVA) for the XXZ model and the Hardware Efficient Ansatz (HEA), we observe a measurement-induced entanglement transition from volume-law to area-law with increasing measurement rate. Moreover, we provide evidence that the transition belongs to the same universality class of random unitary circuits. Importantly, the transition coincides with a "landscape transition" from severe to mild/no barren plateaus in the classical optimization. Our work paves an avenue for greatly improving the trainability of quantum circuits by incorporating intermediate measurement protocols in currently available quantum hardware. △ Less

Submitted 15 November, 2021; originally announced November 2021.

Optimal control of quantum thermal machines using machine learning

Authors: Ilia Khait, Juan Carrasquilla, Dvira Segal

Abstract: Identifying optimal thermodynamical processes has been the essence of thermodynamics since its inception. Here, we show that differentiable programming (DP), a machine learning (ML) tool, can be employed to optimize finite-time thermodynamical processes in a quantum thermal machine. We consider the paradigmatic quantum Otto engine with a time-dependent harmonic oscillator as its working fluid, and… ▽ More Identifying optimal thermodynamical processes has been the essence of thermodynamics since its inception. Here, we show that differentiable programming (DP), a machine learning (ML) tool, can be employed to optimize finite-time thermodynamical processes in a quantum thermal machine. We consider the paradigmatic quantum Otto engine with a time-dependent harmonic oscillator as its working fluid, and build upon shortcut-to-adiabaticity (STA) protocols. We formulate the STA driving protocol as a constrained optimization task and apply DP to find optimal driving profiles for an appropriate figure of merit. Our ML scheme discovers profiles for the compression and expansion strokes that are superior to previously-suggested protocols. Moreover, using our ML algorithm we show that a previously-employed, intuitive energetic cost of the STA driving suffers from a fundamental flaw, which we resolve with an alternative construction for the cost function. Our method and results demonstrate that ML is beneficial both for solving hard-constrained quantum control problems and for devising and assessing their theoretical groundwork. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Comments: 5 pages + 5 figures, SM: 3 page + 1 figure. Comments are welcomed!

Journal ref: Phys. Rev. Research 4, L012029 2022

Neural networks in quantum many-body physics: a hands-on tutorial

Authors: Juan Carrasquilla, Giacomo Torlai

Abstract: Over the past years, machine learning has emerged as a powerful computational tool to tackle complex problems over a broad range of scientific disciplines. In particular, artificial neural networks have been successfully deployed to mitigate the exponential complexity often encountered in quantum many-body physics, the study of properties of quantum systems built out of a large number of interacti… ▽ More Over the past years, machine learning has emerged as a powerful computational tool to tackle complex problems over a broad range of scientific disciplines. In particular, artificial neural networks have been successfully deployed to mitigate the exponential complexity often encountered in quantum many-body physics, the study of properties of quantum systems built out of a large number of interacting particles. In this Article, we overview some applications of machine learning in condensed matter physics and quantum information, with particular emphasis on hands-on tutorials serving as a quick-start for a newcomer to the field. We present supervised machine learning with convolutional neural networks to learn a phase transition, unsupervised learning with restricted Boltzmann machines to perform quantum tomography, and variational Monte Carlo with recurrent neural-networks for approximating the ground state of a many-body Hamiltonian. We briefly review the key ingredients of each algorithm and their corresponding neural-network implementation, and show numerical experiments for a system of interacting Rydberg atoms in two dimensions. △ Less

Submitted 26 January, 2021; originally announced January 2021.

Comments: 21 pages, 7 figures

Variational Neural Annealing

Authors: Mohamed Hibat-Allah, Estelle M. Inack, Roeland Wiersema, Roger G. Melko, Juan Carrasquilla

Abstract: Many important challenges in science and technology can be cast as optimization problems. When viewed in a statistical physics framework, these can be tackled by simulated annealing, where a gradual cooling procedure helps search for groundstate solutions of a target Hamiltonian. While powerful, simulated annealing is known to have prohibitively slow sampling dynamics when the optimization landsca… ▽ More Many important challenges in science and technology can be cast as optimization problems. When viewed in a statistical physics framework, these can be tackled by simulated annealing, where a gradual cooling procedure helps search for groundstate solutions of a target Hamiltonian. While powerful, simulated annealing is known to have prohibitively slow sampling dynamics when the optimization landscape is rough or glassy. Here we show that by generalizing the target distribution with a parameterized model, an analogous annealing framework based on the variational principle can be used to search for groundstate solutions. Modern autoregressive models such as recurrent neural networks provide ideal parameterizations since they can be exactly sampled without slow dynamics even when the model encodes a rough landscape. We implement this procedure in the classical and quantum settings on several prototypical spin glass Hamiltonians, and find that it significantly outperforms traditional simulated annealing in the asymptotic limit, illustrating the potential power of this yet unexplored route to optimization. △ Less

Submitted 25 January, 2021; originally announced January 2021.

Comments: 19 pages, 9 figures, 1 table

Generative models for sampling of lattice field theories

Authors: Matija Medvidovic, Juan Carrasquilla, Lauren E. Hayward, Bohdan Kulchytskyy

Abstract: We explore a self-learning Markov chain Monte Carlo method based on the Adversarial Non-linear Independent Components Estimation Monte Carlo, which utilizes generative models and artificial neural networks. We apply this method to the scalar $\varphi^4$ lattice field theory in the weak-coupling regime and, in doing so, greatly increase the system sizes explored to date with this self-learning tech… ▽ More We explore a self-learning Markov chain Monte Carlo method based on the Adversarial Non-linear Independent Components Estimation Monte Carlo, which utilizes generative models and artificial neural networks. We apply this method to the scalar $\varphi^4$ lattice field theory in the weak-coupling regime and, in doing so, greatly increase the system sizes explored to date with this self-learning technique. Our approach does not rely on a pre-existing training set of samples, as the agent systematically improves its performance by bootstrapping samples collected by the model itself. We evaluate the performance of the trained model by examining its mixing time and study the ergodicity of generated samples. When compared to methods such as Hamiltonian Monte Carlo, this approach provides unique advantages such as the speed of inference and a compressed representation of Monte Carlo proposals for potential use in downstream tasks. △ Less

Submitted 5 January, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

Comments: (7 pages, 3 figures), Updated references. Appeared at the 3rd NeurIPS Workshop on Machine Learning and the Physical Sciences (NeurIPS 2020), Vancouver, Canada

U(1) symmetric recurrent neural networks for quantum state reconstruction

Authors: Stewart Morawetz, Isaac J. S. De Vlugt, Juan Carrasquilla, Roger G. Melko

Abstract: Generative models are a promising technology for the enhancement of quantum simulators. These machine learning methods are capable of reconstructing a quantum state from experimental measurements, and can aid in the calculation of physical observables. In this paper, we employ a recurrent neural network (RNN) to reconstruct the ground state of the spin-1/2 XY model, a prototypical Hamiltonian expl… ▽ More Generative models are a promising technology for the enhancement of quantum simulators. These machine learning methods are capable of reconstructing a quantum state from experimental measurements, and can aid in the calculation of physical observables. In this paper, we employ a recurrent neural network (RNN) to reconstruct the ground state of the spin-1/2 XY model, a prototypical Hamiltonian explored in trapped ion simulators. We explore its performance after enforcing a U(1) symmetry, which was recently shown by Hibat-Allah et al. [Phys. Rev. Research 2, 023358 (2020)] to preserve the autoregressive nature of the RNN. By studying the reconstruction of the XY model ground state from projective measurement data, we show that imposing U(1) symmetry on the RNN significantly increases the efficiency of learning, particularly in the early epoch regime. We argue that this performance increase may result from the tendency of the enforced symmetry to alleviate vanishing and exploding gradients, which helps stabilize the training process. Thus, symmetry-enforced RNNs may be particularly useful for applications of quantum simulators where a rapid feedback between optimization and circuit preparation is necessary, such as in hybrid classical-quantum algorithms. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Comments: 10 pages, 8 figures

Journal ref: Phys. Rev. A 104, 012401 (2021)

Autoregressive Transformer Neural Network for Simulating Open Quantum Systems via a Probabilistic Formulation

Authors: Di Luo, Zhuo Chen, Juan Carrasquilla, Bryan K. Clark

Abstract: The theory of open quantum systems lays the foundations for a substantial part of modern research in quantum science and engineering. Rooted in the dimensionality of their extended Hilbert spaces, the high computational complexity of simulating open quantum systems calls for the development of strategies to approximate their dynamics. In this paper, we present an approach for tackling open quantum… ▽ More The theory of open quantum systems lays the foundations for a substantial part of modern research in quantum science and engineering. Rooted in the dimensionality of their extended Hilbert spaces, the high computational complexity of simulating open quantum systems calls for the development of strategies to approximate their dynamics. In this paper, we present an approach for tackling open quantum system dynamics. Using an exact probabilistic formulation of quantum physics based on positive operator-valued measure (POVM), we compactly represent quantum states with autoregressive transformer neural networks; such networks bring significant algorithmic flexibility due to efficient exact sampling and tractable density. We further introduce the concept of String States to partially restore the symmetry of the autoregressive transformer neural network and improve the description of local correlations. Efficient algorithms have been developed to simulate the dynamics of the Liouvillian superoperator using a forward-backward trapezoid method and find the steady state via a variational formulation. Our approach is benchmarked on prototypical one and two-dimensional systems, finding results which closely track the exact solution and achieve higher accuracy than alternative approaches based on using Markov chain Monte Carlo to sample restricted Boltzmann machines. Our work provides general methods for understanding quantum dynamics in various contexts, as well as techniques for solving high-dimensional probabilistic differential equations in classical setups. △ Less

Submitted 7 June, 2024; v1 submitted 11 September, 2020; originally announced September 2020.

Protocol Discovery for the Quantum Control of Majoranas by Differentiable Programming and Natural Evolution Strategies

Authors: Luuk Coopmans, Di Luo, Graham Kells, Bryan K. Clark, Juan Carrasquilla

Abstract: Quantum control, which refers to the active manipulation of physical systems described by the laws of quantum mechanics, constitutes an essential ingredient for the development of quantum technology. Here we apply Differentiable Programming (DP) and Natural Evolution Strategies (NES) to the optimal transport of Majorana zero modes in superconducting nanowires, a key element to the success of Major… ▽ More Quantum control, which refers to the active manipulation of physical systems described by the laws of quantum mechanics, constitutes an essential ingredient for the development of quantum technology. Here we apply Differentiable Programming (DP) and Natural Evolution Strategies (NES) to the optimal transport of Majorana zero modes in superconducting nanowires, a key element to the success of Majorana-based topological quantum computation. We formulate the motion control of Majorana zero modes as an optimization problem for which we propose a new categorization of four different regimes with respect to the critical velocity of the system and the total transport time. In addition to correctly recovering the anticipated smooth protocols in the adiabatic regime, our algorithms uncover efficient but strikingly counter-intuitive motion strategies in the non-adiabatic regime. The emergent picture reveals a simple but high fidelity strategy that makes use of pulse-like jumps at the beginning and the end of the protocol with a period of constant velocity in between the jumps, which we dub the jump-move-jump protocol. We provide a transparent semi-analytical picture, which uses the sudden approximation and a reformulation of the Majorana motion in a moving frame, to illuminate the key characteristics of the jump-move-jump control strategy. We verify that the jump-move-jump protocol remains robust against the presence of interactions or disorder, and corroborate its high efficacy on a realistic proximity coupled nanowire model. Our results demonstrate that machine learning for quantum control can be applied efficiently to quantum many-body dynamical systems with performance levels that make it relevant to the realization of large-scale quantum technology. △ Less

Submitted 9 April, 2021; v1 submitted 20 August, 2020; originally announced August 2020.

Comments: Revised manuscript: Addition of results on proximity coupled semiconducting nanowire, disorder and interactions. Extra discussions and new conclusion. 18 pages, 8 figures + appendix (12 pages, 8 figures)

Journal ref: PRX Quantum 2, 020332 (2021)

Exploring entanglement and optimization within the Hamiltonian Variational Ansatz

Authors: Roeland Wiersema, Cunlu Zhou, Yvette de Sereville, Juan Felipe Carrasquilla, Yong Baek Kim, Henry Yuen

Abstract: Quantum variational algorithms are one of the most promising applications of near-term quantum computers; however, recent studies have demonstrated that unless the variational quantum circuits are configured in a problem-specific manner, optimization of such circuits will most likely fail. In this paper, we focus on a special family of quantum circuits called the Hamiltonian Variational Ansatz (HV… ▽ More Quantum variational algorithms are one of the most promising applications of near-term quantum computers; however, recent studies have demonstrated that unless the variational quantum circuits are configured in a problem-specific manner, optimization of such circuits will most likely fail. In this paper, we focus on a special family of quantum circuits called the Hamiltonian Variational Ansatz (HVA), which takes inspiration from the quantum approximation optimization algorithm and adiabatic quantum computation. Through the study of its entanglement spectrum and energy gradient statistics, we find that HVA exhibits favorable structural properties such as mild or entirely absent barren plateaus and a restricted state space that eases their optimization in comparison to the well-studied "hardware-efficient ansatz." We also numerically observe that the optimization landscape of HVA becomes almost trap free when the ansatz is over-parameterized. We observe a size-dependent "computational phase transition" as the number of layers in the HVA circuit is increased where the optimization crosses over from a hard to an easy region in terms of the quality of the approximations and speed of convergence to a good solution. In contrast with the analogous transitions observed in the learning of random unitaries which occur at a number of layers that grows exponentially with the number of qubits, our Variational Quantum Eigensolver experiments suggest that the threshold to achieve the over-parameterization phenomenon scales at most polynomially in the number of qubits for the transverse field Ising and XXZ models. Lastly, as a demonstration of its entangling power and effectiveness, we show that HVA can find accurate approximations to the ground states of a modified Haldane-Shastry Hamiltonian on a ring, which has long-range interactions and has a power-law entanglement scaling. △ Less

Submitted 16 November, 2020; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: Updated figure 6, 7 and 11. Other minor changes

Journal ref: PRX Quantum 1, 020319 (2020)

Machine Learning for Quantum Matter

Authors: Juan Carrasquilla

Abstract: Quantum matter, the research field studying phases of matter whose properties are intrinsically quantum mechanical, draws from areas as diverse as hard condensed matter physics, materials science, statistical mechanics, quantum information, quantum gravity, and large-scale numerical simulations. Recently, researchers interested quantum matter and strongly correlated quantum systems have turned the… ▽ More Quantum matter, the research field studying phases of matter whose properties are intrinsically quantum mechanical, draws from areas as diverse as hard condensed matter physics, materials science, statistical mechanics, quantum information, quantum gravity, and large-scale numerical simulations. Recently, researchers interested quantum matter and strongly correlated quantum systems have turned their attention to the algorithms underlying modern machine learning with an eye on making progress in their fields. Here we provide a short review on the recent development and adaptation of machine learning ideas for the purpose advancing research in quantum matter, including ideas ranging from algorithms that recognize conventional and topological states of matter in synthetic an experimental data, to representations of quantum states in terms of neural networks and their applications to the simulation and control of quantum systems. We discuss the outlook for future developments in areas at the intersection between machine learning and quantum many-body physics. △ Less

Submitted 19 August, 2020; v1 submitted 24 March, 2020; originally announced March 2020.

Comments: 34 pages, 4 figures, 290 references. Review paper

Journal ref: Advances in Physics: X, 5:1 (2020)

Visualizing Strange Metallic Correlations in the 2D Fermi-Hubbard Model with AI

Authors: Ehsan Khatami, Elmer Guardado-Sanchez, Benjamin M. Spar, Juan Felipe Carrasquilla, Waseem S. Bakr, Richard T. Scalettar

Abstract: Strongly correlated phases of matter are often described in terms of straightforward electronic patterns. This has so far been the basis for studying the Fermi-Hubbard model realized with ultracold atoms. Here, we show that artificial intelligence (AI) can provide an unbiased alternative to this paradigm for phases with subtle, or even unknown, patterns. Long- and short-range spin correlations spo… ▽ More Strongly correlated phases of matter are often described in terms of straightforward electronic patterns. This has so far been the basis for studying the Fermi-Hubbard model realized with ultracold atoms. Here, we show that artificial intelligence (AI) can provide an unbiased alternative to this paradigm for phases with subtle, or even unknown, patterns. Long- and short-range spin correlations spontaneously emerge in filters of a convolutional neural network trained on snapshots of single atomic species. In the less well-understood strange metallic phase of the model, we find that a more complex network trained on snapshots of local moments produces an effective order parameter for the non-Fermi-liquid behavior. Our technique can be employed to characterize correlations unique to other phases with no obvious order parameters or signatures in projective measurements, and has implications for science discovery through AI beyond strongly correlated systems. △ Less

Submitted 23 December, 2020; v1 submitted 27 February, 2020; originally announced February 2020.

Comments: 12 pages, 9 figures; updated in accord with the published version

Journal ref: Phys. Rev. A 102, 033326 (2020)

Recurrent Neural Network Wave Functions

Authors: Mohamed Hibat-Allah, Martin Ganahl, Lauren E. Hayward, Roger G. Melko, Juan Carrasquilla

Abstract: A core technology that has emerged from the artificial intelligence revolution is the recurrent neural network (RNN). Its unique sequence-based architecture provides a tractable likelihood estimate with stable training paradigms, a combination that has precipitated many spectacular advances in natural language processing and neural machine translation. This architecture also makes a good candidate… ▽ More A core technology that has emerged from the artificial intelligence revolution is the recurrent neural network (RNN). Its unique sequence-based architecture provides a tractable likelihood estimate with stable training paradigms, a combination that has precipitated many spectacular advances in natural language processing and neural machine translation. This architecture also makes a good candidate for a variational wave function, where the RNN parameters are tuned to learn the approximate ground state of a quantum Hamiltonian. In this paper, we demonstrate the ability of RNNs to represent several many-body wave functions, optimizing the variational parameters using a stochastic approach. Among other attractive features of these variational wave functions, their autoregressive nature allows for the efficient calculation of physical estimators by providing independent samples. We demonstrate the effectiveness of RNN wave functions by calculating ground state energies, correlation functions, and entanglement entropies for several quantum spin models of interest to condensed matter physicists in one and two spatial dimensions. △ Less

Submitted 20 June, 2020; v1 submitted 7 February, 2020; originally announced February 2020.

Comments: The GitHub link to the open-source code is fixed

Journal ref: Phys. Rev. Research 2, 023358 (2020)

Probabilistic Simulation of Quantum Circuits with the Transformer

Authors: Juan Carrasquilla, Di Luo, Felipe Pérez, Ashley Milsted, Bryan K. Clark, Maksims Volkovs, Leandro Aolita

Abstract: The fundamental question of how to best simulate quantum systems using conventional computational resources lies at the forefront of condensed matter and quantum computation. It impacts both our understanding of quantum materials and our ability to emulate quantum circuits. Here we present an exact formulation of quantum dynamics via factorized generalized measurements which maps quantum states to… ▽ More The fundamental question of how to best simulate quantum systems using conventional computational resources lies at the forefront of condensed matter and quantum computation. It impacts both our understanding of quantum materials and our ability to emulate quantum circuits. Here we present an exact formulation of quantum dynamics via factorized generalized measurements which maps quantum states to probability distributions with the advantage that local unitary dynamics and quantum channels map to local quasi-stochastic matrices. This representation provides a general framework for using state-of-the-art probabilistic models in machine learning for the simulation of quantum many-body dynamics. Using this framework, we have developed a practical algorithm to simulate quantum circuits with the Transformer, a powerful ansatz responsible for the most recent breakthroughs in natural language processing. We demonstrate our approach by simulating circuits which build GHZ and linear graph states of up to 60 qubits, as well as a variational quantum eigensolver circuit for preparing the ground state of the transverse field Ising model on six qubits. Our methodology constitutes a modern machine learning approach to the simulation of quantum physics with applicability both to quantum circuits as well as other quantum many-body systems. △ Less

Submitted 23 December, 2019; originally announced December 2019.

Comments: 12 pages, 11 figures including supplementary materials

Journal ref: Phys. Rev. A 104, 032610 (2021)

Wavefunction positivization via automatic differentiation

Authors: Giacomo Torlai, Juan Carrasquilla, Matthew T. Fishman, Roger G. Melko, Matthew P. A. Fisher

Abstract: We introduce a procedure to systematically search for a local unitary transformation that maps a wavefunction with a non-trivial sign structure into a positive-real form. The transformation is parametrized as a quantum circuit compiled into a set of one and two qubit gates. We design a cost function that maximizes the average sign of the output state and removes its complex phases. The optimizatio… ▽ More We introduce a procedure to systematically search for a local unitary transformation that maps a wavefunction with a non-trivial sign structure into a positive-real form. The transformation is parametrized as a quantum circuit compiled into a set of one and two qubit gates. We design a cost function that maximizes the average sign of the output state and removes its complex phases. The optimization of the gates is performed through automatic differentiation algorithms, widely used in the machine learning community. We provide numerical evidence for significant improvements in the average sign for a two-leg triangular Heisenberg ladder with next-to-nearest neighbour and ring-exchange interactions. This model exhibits phases where the sign structure can be removed by simple local one-qubit unitaries, but also an exotic Bose-metal phase whose sign structure induces "Bose surfaces" with a fermionic character and a higher entanglement that requires deeper circuits. △ Less

Submitted 10 February, 2020; v1 submitted 11 June, 2019; originally announced June 2019.

Comments: 9 pages, 5 figures

Journal ref: Phys. Rev. Research 2, 032060 (2020)

Neural Gutzwiller-projected variational wave functions

Authors: Francesco Ferrari, Federico Becca, Juan Carrasquilla

Abstract: Variational wave functions have enabled exceptional scientific breakthroughs related to the understanding of novel phases of matter. Examples include the Bardeen-Cooper-Schrieffer theory of superconductivity, the description of the fractional quantum Hall effect through the Laughlin state, and Feynman's variational understanding of large-scale quantum effects in liquid Helium. More recently, Gutzw… ▽ More Variational wave functions have enabled exceptional scientific breakthroughs related to the understanding of novel phases of matter. Examples include the Bardeen-Cooper-Schrieffer theory of superconductivity, the description of the fractional quantum Hall effect through the Laughlin state, and Feynman's variational understanding of large-scale quantum effects in liquid Helium. More recently, Gutzwiller-projected wave functions, typically constructed from fermionic degrees of freedom, have been employed to examine quantum spin models in the presence of competing interactions, where exotic phases with no spontaneous symmetry breaking and fractional excitations may exist. In this work, we investigate the aforementioned fermionic wave functions supplemented with neural networks, specifically with the so-called restricted Boltzmann machine (RBM), to boost their accuracy and obtain reliable approximations to the ground state of generic spin models. In particular, we apply our neural augmented fermionic construction to the description of both magnetically ordered and disordered phases of increasing complexity, including cases where the ground state displays a non-trivial sign structure. Even though the RBM state is by far more effective for Néel states endowed with a particularly simple sign structure, it provides a significant improvement over the original fermionic state in highly frustrated regimes where a complex sign structure is anticipated, thus marking the path to an understanding of strongly-correlated spin models on the lattice via neural Gutzwiller-projected variational wave functions. △ Less

Submitted 13 September, 2019; v1 submitted 2 June, 2019; originally announced June 2019.

Comments: 13 pages, 10 figures

Journal ref: Phys. Rev. B 100, 125131 (2019)

Reconstructing quantum states with generative models

Authors: Juan Carrasquilla, Giacomo Torlai, Roger G. Melko, Leandro Aolita

Abstract: A major bottleneck in the quest for scalable many-body quantum technologies is the difficulty in benchmarking their preparations, which suffer from an exponential `curse of dimensionality' inherent to their quantum states. We present an experimentally friendly method for density matrix reconstruction based on deep neural-network generative models. The learning procedure comes with a built-in appro… ▽ More A major bottleneck in the quest for scalable many-body quantum technologies is the difficulty in benchmarking their preparations, which suffer from an exponential `curse of dimensionality' inherent to their quantum states. We present an experimentally friendly method for density matrix reconstruction based on deep neural-network generative models. The learning procedure comes with a built-in approximate certificate of the reconstruction and makes no assumptions on the state under scrutiny, making it both reliable and unconditional. It can efficiently handle a broad class of complex systems including prototypical states in quantum information, as well as ground states of local spin models common to condensed matter physics. The key insight is to reduce the state tomography task to an unsupervised learning problem of the statistics of an informationally complete set of quantum measurements. This constitutes a modern machine learning approach to the validation of large quantum devices, which may prove relevant as a neural-network ansatz over mixed states suitable for variational optimization. △ Less

Submitted 29 July, 2019; v1 submitted 24 October, 2018; originally announced October 2018.

Comments: 29 pages, 9 figures. Includes appendix. Some typos in the description of the POVMs have been fixed. The numerically generated measurements, the implementation of the generative models, and the code to numerically generate the data sets used in this manuscript are available at https://github.com/carrasqu/POVM_GENMODEL

Journal ref: Nature Machine Intelligence, vol. 1, 155-161 (2019)

Observation of topological phenomena in a programmable lattice of 1,800 qubits

Authors: Andrew D. King, Juan Carrasquilla, Isil Ozfidan, Jack Raymond, Evgeny Andriyash, Andrew Berkley, Mauricio Reis, Trevor M. Lanting, Richard Harris, Gabriel Poulin-Lamarre, Anatoly Yu. Smirnov, Christopher Rich, Fabio Altomare, Paul Bunyk, Jed Whittaker, Loren Swenson, Emile Hoskinson, Yuki Sato, Mark Volkmann, Eric Ladizinsky, Mark Johnson, Jeremy Hilton, Mohammad H. Amin

Abstract: The celebrated work of Berezinskii, Kosterlitz and Thouless in the 1970s revealed exotic phases of matter governed by topological properties of low-dimensional materials such as thin films of superfluids and superconductors. Key to this phenomenon is the appearance and interaction of vortices and antivortices in an angular degree of freedom---typified by the classical XY model---due to thermal flu… ▽ More The celebrated work of Berezinskii, Kosterlitz and Thouless in the 1970s revealed exotic phases of matter governed by topological properties of low-dimensional materials such as thin films of superfluids and superconductors. Key to this phenomenon is the appearance and interaction of vortices and antivortices in an angular degree of freedom---typified by the classical XY model---due to thermal fluctuations. In the 2D Ising model this angular degree of freedom is absent in the classical case, but with the addition of a transverse field it can emerge from the interplay between frustration and quantum fluctuations. Consequently a Kosterlitz-Thouless (KT) phase transition has been predicted in the quantum system by theory and simulation. Here we demonstrate a large-scale quantum simulation of this phenomenon in a network of 1,800 in situ programmable superconducting flux qubits arranged in a fully-frustrated square-octagonal lattice. Essential to the critical behavior, we observe the emergence of a complex order parameter with continuous rotational symmetry, and the onset of quasi-long-range order as the system approaches a critical temperature. We use a simple but previously undemonstrated approach to statistical estimation with an annealing-based quantum processor, performing Monte Carlo sampling in a chain of reverse quantum annealing protocols. Observations are consistent with classical simulations across a range of Hamiltonian parameters. We anticipate that our approach of using a quantum processor as a programmable magnetic lattice will find widespread use in the simulation and development of exotic materials. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: 17 pages, 15 figures

Journal ref: Nature, volume 560, pages 456-460 (2018)

Tripartite entangled plaquette state in a cluster magnet

Authors: Juan Carrasquilla, Gang Chen, Roger G. Melko

Abstract: Using large-scale quantum Monte Carlo simulations we show that a spin-$1/2$ XXZ model on a two-dimensional anisotropic Kagome lattice exhibits a tripartite entangled plaquette state that preserves all of the Hamiltonian symmetries. It is connected via phase boundaries to a ferromagnet and a valence-bond solid that break U(1) and lattice translation symmetries, respectively. We study the phase diag… ▽ More Using large-scale quantum Monte Carlo simulations we show that a spin-$1/2$ XXZ model on a two-dimensional anisotropic Kagome lattice exhibits a tripartite entangled plaquette state that preserves all of the Hamiltonian symmetries. It is connected via phase boundaries to a ferromagnet and a valence-bond solid that break U(1) and lattice translation symmetries, respectively. We study the phase diagram of the model in detail, in particular the transitions to the tripartite entangled plaquette state, which are consistent with conventional order-disorder transitions. Our results can be interpreted as a description of the charge sector dynamics of a Hubbard model applied to the description of the spin liquid candidate ${\mathrm{LiZn}}_{2}{\mathrm{Mo}}_{3}{\mathrm{O}}_{8}$, as well as a model of strongly correlated bosonic atoms loaded onto highly tunable {\it trimerized} optical Kagome lattices. △ Less

Submitted 11 April, 2017; originally announced April 2017.

Comments: 10 pages, 10 figures

Journal ref: Phys. Rev. B 96, 054405 (2017)

Many-body quantum state tomography with neural networks

Authors: Giacomo Torlai, Guglielmo Mazzola, Juan Carrasquilla, Matthias Troyer, Roger Melko, Giuseppe Carleo

Abstract: The experimental realization of increasingly complex synthetic quantum systems calls for the development of general theoretical methods, to validate and fully exploit quantum resources. Quantum-state tomography (QST) aims at reconstructing the full quantum state from simple measurements, and therefore provides a key tool to obtain reliable analytics. Brute-force approaches to QST, however, demand… ▽ More The experimental realization of increasingly complex synthetic quantum systems calls for the development of general theoretical methods, to validate and fully exploit quantum resources. Quantum-state tomography (QST) aims at reconstructing the full quantum state from simple measurements, and therefore provides a key tool to obtain reliable analytics. Brute-force approaches to QST, however, demand resources growing exponentially with the number of constituents, making it unfeasible except for small systems. Here we show that machine learning techniques can be efficiently used for QST of highly-entangled states, in both one and two dimensions. Remarkably, the resulting approach allows one to reconstruct traditionally challenging many-body quantities - such as the entanglement entropy - from simple, experimentally accessible measurements. This approach can benefit existing and future generations of devices ranging from quantum computers to ultra-cold atom quantum simulators. △ Less

Submitted 23 October, 2017; v1 submitted 15 March, 2017; originally announced March 2017.

Comments: Update version and method, now discussing how to reconstruct the complex amplitudes of the wave function

Journal ref: Nature Physics 14, 447-450 (2018)

Machine Learning Phases of Strongly Correlated Fermions

Authors: Kelvin Ch'ng, Juan Carrasquilla, Roger G. Melko, Ehsan Khatami

Abstract: Machine learning offers an unprecedented perspective for the problem of classifying phases in condensed matter physics. We employ neural-network machine learning techniques to distinguish finite-temperature phases of the strongly correlated fermions on cubic lattices. We show that a three dimensional convolutional network trained on auxiliary field configurations produced by quantum Monte Carlo si… ▽ More Machine learning offers an unprecedented perspective for the problem of classifying phases in condensed matter physics. We employ neural-network machine learning techniques to distinguish finite-temperature phases of the strongly correlated fermions on cubic lattices. We show that a three dimensional convolutional network trained on auxiliary field configurations produced by quantum Monte Carlo simulations of the Hubbard model can correctly predict the magnetic phase diagram of the model at the average density of one (half filling). We then use the network, trained at half filling, to explore the trend in the transition temperature as the system is doped away from half filling. This transfer learning approach predicts that the instability to the magnetic phase extends to at least 5% doping in this region. Our results pave the way for other machine learning applications in correlated quantum many-body systems. △ Less

Submitted 10 September, 2017; v1 submitted 8 September, 2016; originally announced September 2016.

Comments: 9 pages, 8 figures

Journal ref: Phys. Rev. X 7, 031038 (2017)

Machine learning quantum phases of matter beyond the fermion sign problem

Authors: Peter Broecker, Juan Carrasquilla, Roger G. Melko, Simon Trebst

Abstract: State-of-the-art machine learning techniques promise to become a powerful tool in statistical mechanics via their capacity to distinguish different phases of matter in an automated way. Here we demonstrate that convolutional neural networks (CNN) can be optimized for quantum many-fermion systems such that they correctly identify and locate quantum phase transitions in such systems. Using auxiliary… ▽ More State-of-the-art machine learning techniques promise to become a powerful tool in statistical mechanics via their capacity to distinguish different phases of matter in an automated way. Here we demonstrate that convolutional neural networks (CNN) can be optimized for quantum many-fermion systems such that they correctly identify and locate quantum phase transitions in such systems. Using auxiliary-field quantum Monte Carlo (QMC) simulations to sample the many-fermion system, we show that the Green's function (but not the auxiliary field) holds sufficient information to allow for the distinction of different fermionic phases via a CNN. We demonstrate that this QMC + machine learning approach works even for systems exhibiting a severe fermion sign problem where conventional approaches to extract information from the Green's function, e.g.~in the form of equal-time correlation functions, fail. We expect that this capacity of hierarchical machine learning techniques to circumvent the fermion sign problem will drive novel insights into some of the most fundamental problems in statistical physics. △ Less

Submitted 28 August, 2016; originally announced August 2016.

Comments: Comments: 8 pages, 6 figures

Journal ref: Scientific Reports 7, 8823 (2017)

Machine learning phases of matter

Authors: Juan Carrasquilla, Roger G. Melko

Abstract: Neural networks can be used to identify phases and phase transitions in condensed matter systems via supervised machine learning. Readily programmable through modern software libraries, we show that a standard feed-forward neural network can be trained to detect multiple types of order parameter directly from raw state configurations sampled with Monte Carlo. In addition, they can detect highly no… ▽ More Neural networks can be used to identify phases and phase transitions in condensed matter systems via supervised machine learning. Readily programmable through modern software libraries, we show that a standard feed-forward neural network can be trained to detect multiple types of order parameter directly from raw state configurations sampled with Monte Carlo. In addition, they can detect highly non-trivial states such as Coulomb phases, and if modified to a convolutional neural network, topological phases with no conventional order parameter. We show that this classification occurs within the neural network without knowledge of the Hamiltonian or even the general locality of interactions. These results demonstrate the power of machine learning as a basic research tool in the field of condensed matter and statistical physics. △ Less

Submitted 5 May, 2016; originally announced May 2016.

Comments: 18 pages, 8 figures, 1 table

Spinon walk in quantum spin ice

Authors: Yuan Wan, Juan Carrasquilla, Roger G. Melko

Abstract: We study a minimal model for the dynamics of spinons in quantum spin ice. The model captures the essential strong coupling between the spinon and the disordered background spins. We demonstrate that the spinon motion can be mapped to a random walk with an entropy-induced memory in imaginary time. Our numerical simulation of the spinon walk indicates that the spinon propagates as a massive quasipar… ▽ More We study a minimal model for the dynamics of spinons in quantum spin ice. The model captures the essential strong coupling between the spinon and the disordered background spins. We demonstrate that the spinon motion can be mapped to a random walk with an entropy-induced memory in imaginary time. Our numerical simulation of the spinon walk indicates that the spinon propagates as a massive quasiparticle at low energy despite its strong coupling to the spin background at the microscopic energy scale. We discuss the experimental implications of our findings. △ Less

Submitted 22 April, 2016; v1 submitted 4 October, 2015; originally announced October 2015.

Comments: 12 pages, 10 figures

Journal ref: Phys. Rev. Lett. 116, 167202 (2016)

Gap scaling at Berezinskii-Kosterlitz-Thouless quantum critical points in one-dimensional Hubbard and Heisenberg models

Authors: M. Dalmonte, J. Carrasquilla, L. Taddia, E. Ercolessi, M. Rigol

Abstract: We discuss how to locate critical points in the Berezinskii-Kosterlitz-Thouless (BKT) universality class by means of gap-scaling analyses. While accurately determining such points using gap extrapolation procedures is usually challenging and inaccurate due to the exponentially small value of the gap in the vicinity of the critical point, we show that a generic gap-scaling analysis, including the e… ▽ More We discuss how to locate critical points in the Berezinskii-Kosterlitz-Thouless (BKT) universality class by means of gap-scaling analyses. While accurately determining such points using gap extrapolation procedures is usually challenging and inaccurate due to the exponentially small value of the gap in the vicinity of the critical point, we show that a generic gap-scaling analysis, including the effects of logarithmic corrections, provides very accurate estimates of BKT transition points in a variety of spin and fermionic models. As a first example, we show how the scaling procedure, combined with density-matrix-renormalization-group simulations, performs extremely well in a non-integrable spin-$3/2$ XXZ model, which is known to exhibit strong finite-size effects. We then analyze the extended Hubbard model, whose BKT transition has been debated, finding results that are consistent with previous studies based on the scaling of the Luttinger-liquid parameter. Finally, we investigate an anisotropic extended Hubbard model, for which we present the first estimates of the BKT transition line based on large-scale density-matrix-renormalization-group simulations. Our work demonstrates how gap-scaling analyses can help to locate accurately and efficiently BKT critical points, without relying on model-dependent scaling assumptions. △ Less

Submitted 30 April, 2015; v1 submitted 17 December, 2014; originally announced December 2014.

Comments: 8 pages, 7 figures

Journal ref: Phys. Rev. B 91, 165136 (2015)

Spectral tensor networks for many-body localization

Authors: A. Chandran, J. Carrasquilla, I. H. Kim, D. A. Abanin, G. Vidal

Abstract: Subsystems of strongly disordered, interacting quantum systems can fail to thermalize because of the phenomenon of many-body localization (MBL). In this article, we explore a tensor network description of the eigenspectra of such systems. Specifically, we will argue that the presence of a complete set of local integrals of motion in MBL implies an efficient representation of the entire spectrum of… ▽ More Subsystems of strongly disordered, interacting quantum systems can fail to thermalize because of the phenomenon of many-body localization (MBL). In this article, we explore a tensor network description of the eigenspectra of such systems. Specifically, we will argue that the presence of a complete set of local integrals of motion in MBL implies an efficient representation of the entire spectrum of energy eigenstates with a single tensor network, a \emph{spectral} tensor network. Our results are rigorous for a class of idealized systems related to MBL with integrals of motion of finite support. In one spatial dimension, the spectral tensor network allows for the efficient computation of expectation values of a large class of operators (including local operators and string operators) in individual energy eigenstates and in ensembles. △ Less

Submitted 2 October, 2014; originally announced October 2014.

Comments: 7 pages, 6 figures

Journal ref: Phys. Rev. B 92, 024201 (2015)

Quantum distillation and confinement of vacancies in a doublon sea

Authors: Lin Xia, Laura A. Zundel, Juan Carrasquilla, Aaron Reinhard, Josh M. Wilson, Marcos Rigol, David S. Weiss

Abstract: Ultracold atomic gases have revolutionized the study of non-equilibrium dynamics in quantum many-body systems. Many counterintuitive non-equilibrium effects have been observed, such as suppressed thermalization in a one-dimensional (1D) gas, the formation of repulsive self-bound dimers, and identical behaviors for attractive and repulsive interactions. Here, we observe the expansion of a bundle of… ▽ More Ultracold atomic gases have revolutionized the study of non-equilibrium dynamics in quantum many-body systems. Many counterintuitive non-equilibrium effects have been observed, such as suppressed thermalization in a one-dimensional (1D) gas, the formation of repulsive self-bound dimers, and identical behaviors for attractive and repulsive interactions. Here, we observe the expansion of a bundle of ultracold 1D Bose gases in a flat-bottomed optical lattice potential. By combining in situ measurements with photoassociation, we follow the spatial dynamics of singly, doubly, and triply occupied lattice sites. The system sheds interaction energy by dissolving some doublons and triplons. Some singlons quantum distill out of the doublon center, while others remain confined. Our Gutzwiller mean-field model captures these experimental features in a physically clear way. These experiments might be used to study thermalization in systems with particle losses or the evolution of quantum entanglement, or if applied to fermions, to prepare very low entropy states. △ Less

Submitted 9 September, 2015; v1 submitted 9 September, 2014; originally announced September 2014.

Comments: 22 pages, 8 figures

Journal ref: Nature Physics 11, 316-320 (2015)

Quantum Kagome Ice

Authors: Juan Carrasquilla, Zhihao Hao, Roger G. Melko

Abstract: Actively shought since the turn of the century, two-dimensional quantum spin liquids (QSLs) are exotic phases of matter where magnetic moments remain disordered even at extremely low temperatures. Despite ongoing searches, QSLs remain elusive, due to a lack of concrete knowledge of the microscopic mechanisms that inhibit magnetic order in real materials. Here, we study a theoretical model for a br… ▽ More Actively shought since the turn of the century, two-dimensional quantum spin liquids (QSLs) are exotic phases of matter where magnetic moments remain disordered even at extremely low temperatures. Despite ongoing searches, QSLs remain elusive, due to a lack of concrete knowledge of the microscopic mechanisms that inhibit magnetic order in real materials. Here, we study a theoretical model for a broad class of frustrated magnetic rare-earth pyrochlore materials called "quantum spin ices". When subject to an external magnetic field along the [111] crystallographic direction, the resulting spin interactions contain a mix of geometric frustration and quantum fluctuations in decoupled two-dimensional kagome planes. Using large-scale quantum Monte Carlo simulations, we identify a simple set of interactions sufficient to promote a groundstate with no magnetic long-range order, and a gap to excitations, consistent with a $Z_2$ spin liquid phase. This suggests a systematic experimental procedure to search for two-dimensional QSLs within the broader class of three-dimensional pyrochlore quantum spin ice materials. △ Less

Submitted 30 June, 2014; originally announced July 2014.

Comments: 11 pages, 13 figures. Includes supplemental information

Journal ref: Nature Communications 6, Article number: 7421 2015

Spiral antiferromagnets beyond the spin-wave approximation: frustrated $XY$ and Heisenberg models in the honeycomb lattice

Authors: Andrea Di Ciolo, Juan Carrasquilla, Federico Becca, Marcos Rigol, Victor Galitski

Abstract: We examine the stability of classical states with a generic incommensurate spiral order against quantum fluctuations. Specifically, we focus on the frustrated spin-1/2 $XY$ and Heisenberg models on the honeycomb lattice with nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ antiferromagnetic couplings. Our variational approach is based on the Jastrow wave functions, which include quantum corr… ▽ More We examine the stability of classical states with a generic incommensurate spiral order against quantum fluctuations. Specifically, we focus on the frustrated spin-1/2 $XY$ and Heisenberg models on the honeycomb lattice with nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ antiferromagnetic couplings. Our variational approach is based on the Jastrow wave functions, which include quantum correlations on top of classical spin waves. We perform a systematic optimization of wave vectors and Jastrow pseudo-potentials within this class of variational states and find that quantum fluctuations favor collinear states over generic coplanar spirals. The Néel state with ${\bf Q}=(0,0)$ extends its stability well beyond the classical value $J_2/J_1=1/6$. Most importantly, the collinear states with ${\bf Q}=(0,2π/\sqrt{3})$ (and the two symmetry-related states) are found to be stable in a large regime with intermediate frustration, while at the classical level they are limited to the point $J_2/J_1=0.5$. For large frustration, the $120^\circ$ state is stabilized for finite values of $J_2/J_1$ in both models. △ Less

Submitted 28 March, 2014; v1 submitted 16 December, 2013; originally announced December 2013.

Comments: 10 pages, 10 figures, 4 tables

Journal ref: Phys. Rev. B 89, 094413 (2014)

Nature of the phases in the frustrated XY model on the honeycomb lattice

Authors: Juan Carrasquilla, Andrea Di Ciolo, Federico Becca, Victor Galitski, Marcos Rigol

Abstract: We study the phase diagram of the frustrated XY model on the honeycomb lattice by using accurate correlated wave functions and variational Monte Carlo simulations. Our results suggest that a spin-liquid state is energetically favorable in the region of intermediate frustration, intervening between two magnetically ordered phases. The latter ones are represented by classically ordered states supple… ▽ More We study the phase diagram of the frustrated XY model on the honeycomb lattice by using accurate correlated wave functions and variational Monte Carlo simulations. Our results suggest that a spin-liquid state is energetically favorable in the region of intermediate frustration, intervening between two magnetically ordered phases. The latter ones are represented by classically ordered states supplemented with a long-range Jastrow factor, which includes relevant correlations and dramatically improves the description provided by the purely classical solution of the model. The construction of the spin-liquid state is based on a decomposition of the underlying bosonic particles in terms of spin-1/2 fermions (partons), with a Gutzwiller projection enforcing no single occupancy, as well as a long-range Jastrow factor. △ Less

Submitted 5 January, 2014; v1 submitted 8 July, 2013; originally announced July 2013.

Comments: 9 pages, 6 figures, 3 tables

Journal ref: Phys. Rev. B 88, 241109(R) (2013)

Self trapping in the two-dimensional Bose-Hubbard model

Authors: Andrew Jreissaty, Juan Carrasquilla, Marcos Rigol

Abstract: We study the expansion of harmonically trapped bosons in a two-dimensional lattice after suddenly turning off the confining potential. We show that, in the presence of multiple occupancies per lattice site and strong interactions, the system exhibits a clear dynamical separation into slowly and rapidly expanding clouds. We discuss how this effect can be understood within a simple picture by invoki… ▽ More We study the expansion of harmonically trapped bosons in a two-dimensional lattice after suddenly turning off the confining potential. We show that, in the presence of multiple occupancies per lattice site and strong interactions, the system exhibits a clear dynamical separation into slowly and rapidly expanding clouds. We discuss how this effect can be understood within a simple picture by invoking doublons and Bose enhancement. This picture is corroborated by an analysis of the momentum distribution function in the regions with slowly and rapidly expanding bosons. △ Less

Submitted 4 June, 2014; v1 submitted 2 July, 2013; originally announced July 2013.

Comments: 7 pages, 8 figures, as published

Journal ref: Phys. Rev. A 88, 031606(R) (2013)

Scaling of the gap, fidelity susceptibility, and Bloch oscillations across the superfluid to Mott insulator transition in the one-dimensional Bose-Hubbard model

Authors: Juan Carrasquilla, Salvatore R. Manmana, Marcos Rigol

Abstract: We investigate the interaction-induced superfluid-to-Mott insulator transition in the one-dimensional Bose-Hubbard model (BHM) for fillings $n=1$, $n=2$, and $n=3$ by studying the single-particle gap, the fidelity susceptibility, and the amplitude of Bloch oscillations via density-matrix renormalization-group methods. We apply a generic scaling procedure for the gap, which allows us to determine t… ▽ More We investigate the interaction-induced superfluid-to-Mott insulator transition in the one-dimensional Bose-Hubbard model (BHM) for fillings $n=1$, $n=2$, and $n=3$ by studying the single-particle gap, the fidelity susceptibility, and the amplitude of Bloch oscillations via density-matrix renormalization-group methods. We apply a generic scaling procedure for the gap, which allows us to determine the critical points with very high accuracy. We also study how the fidelity susceptibility behaves across the phase transition. Furthermore, we show that in the BHM, and in a system of spinless fermions, the amplitude of Bloch oscillations after a tilt of the lattice vanishes at the critical points. This indicates that Bloch oscillations can serve as a tool to detect the transition point in ongoing experiments with ultracold gases. △ Less

Submitted 19 April, 2013; v1 submitted 10 December, 2012; originally announced December 2012.

Comments: 6 pages and 4 figures, URL: http://link.aps.org/doi/10.1103/PhysRevA.87.043606

Journal ref: Phys. Rev. A 87, 043606 (2013)

Superfluid to normal phase transition in strongly correlated bosons in two and three dimensions

Authors: Juan Carrasquilla, Marcos Rigol

Abstract: Using quantum Monte Carlo simulations, we investigate the finite-temperature phase diagram of hard-core bosons (XY model) in two- and three-dimensional lattices. To determine the phase boundaries, we perform a finite-size-scaling analysis of the condensate fraction and/or the superfluid stiffness. We then discuss how these phase diagrams can be measured in experiments with trapped ultracold gases,… ▽ More Using quantum Monte Carlo simulations, we investigate the finite-temperature phase diagram of hard-core bosons (XY model) in two- and three-dimensional lattices. To determine the phase boundaries, we perform a finite-size-scaling analysis of the condensate fraction and/or the superfluid stiffness. We then discuss how these phase diagrams can be measured in experiments with trapped ultracold gases, where the systems are inhomogeneous. For that, we introduce a method based on the measurement of the zero-momentum occupation, which is adequate for experiments dealing with both homogeneous and trapped systems, and compare it with previously proposed approaches. △ Less

Submitted 26 October, 2012; v1 submitted 29 May, 2012; originally announced May 2012.

Comments: 13 pages, 11 figures. http://link.aps.org/doi/10.1103/PhysRevA.86.043629

Journal ref: Phys. Rev. A 86, 043629 (2012)

Phase diagram of the half-filled one-dimensional t-V-V' model

Authors: Tapan Mishra, Juan Carrasquilla, Marcos Rigol

Abstract: We study the phase diagram of spinless fermions with nearest and next-nearest-neighbor interactions in one dimension utilizing the (finite-size) density-matrix renormalization group (DMRG) method. The competition between nearest and next-nearest-neighbor interactions and nearest-neighbor hopping generates four phases in this model: two charge-density-wave insulators, a Luttinger liquid phase, and… ▽ More We study the phase diagram of spinless fermions with nearest and next-nearest-neighbor interactions in one dimension utilizing the (finite-size) density-matrix renormalization group (DMRG) method. The competition between nearest and next-nearest-neighbor interactions and nearest-neighbor hopping generates four phases in this model: two charge-density-wave insulators, a Luttinger liquid phase, and a bond-order phase. We use finite-size scaling of the gap and various structure factors to determine the phase diagram. △ Less

Submitted 26 September, 2011; v1 submitted 10 August, 2011; originally announced August 2011.

Comments: 7 pages, 9 figures

Journal ref: Phys. Rev. B 84, 115135 (2011)

Expansion of Bose-Hubbard Mott insulators in optical lattices

Authors: Mark Jreissaty, Juan Carrasquilla, F. Alexander Wolf, Marcos Rigol

Abstract: We study the expansion of bosonic Mott insulators in the presence of an optical lattice after switching off a confining potential. We use the Gutzwiller mean-field approximation and consider two different setups. In the first one, the expansion is restricted to one direction. We show that this leads to the emergence of two condensates with well defined momenta, and argue that such a construct can… ▽ More We study the expansion of bosonic Mott insulators in the presence of an optical lattice after switching off a confining potential. We use the Gutzwiller mean-field approximation and consider two different setups. In the first one, the expansion is restricted to one direction. We show that this leads to the emergence of two condensates with well defined momenta, and argue that such a construct can be used to create atom lasers in optical lattices. In the second setup, we study Mott insulators that are allowed to expand in all directions in the lattice. In this case, a simple condensate is seen to develop within the mean-field approximation. However, its constituent bosons are found to populate many nonzero momentum modes. An analytic understanding of both phenomena in terms of the exact dispersion relation in the hard-core limit is presented. △ Less

Submitted 9 July, 2012; v1 submitted 4 August, 2011; originally announced August 2011.

Comments: 11 pages, 9 figures. Figures 2,3,4 corrected

Journal ref: Phys. Rev. A 84, 043610 (2011)

Bose-glass, superfluid, and rung-Mott phases of hard-core bosons in disordered two-leg ladders

Authors: Juan Carrasquilla, Federico Becca, Michele Fabrizio

Abstract: By means of Monte Carlo techniques, we study the role of disorder on a system of hard-core bosons in a two-leg ladder with both intra-chain ($t$) and inter-chain ($t^\prime$) hoppings. We find that the phase diagram as a function of the boson density, disorder strength, and $t^\prime/t$ is far from being trivial. This contrasts the case of spin-less fermions where standard localization arguments a… ▽ More By means of Monte Carlo techniques, we study the role of disorder on a system of hard-core bosons in a two-leg ladder with both intra-chain ($t$) and inter-chain ($t^\prime$) hoppings. We find that the phase diagram as a function of the boson density, disorder strength, and $t^\prime/t$ is far from being trivial. This contrasts the case of spin-less fermions where standard localization arguments apply and an Anderson-localized phase pervades the whole phase diagram. A compressible Bose-glass phase always intrudes between the Mott insulator with zero (or one) bosons per site and the superfluid that is stabilized for weak disorder. At half filling, there is a direct transition between a (gapped) rung-Mott insulator and a Bose glass, which is driven by exponentially rare regions where disorder is suppressed. Finally, by doping the rung-Mott insulator, a direct transition to the superfluid is possible only in the clean system, whereas the Mott phase is always surrounded by the a Bose glass when disorder is present. The phase diagram based on our numerical evidence is finally reported. △ Less

Submitted 16 February, 2011; originally announced February 2011.

Comments: 8 pages, 7 figures

Journal ref: Phys. Rev. B 83, 245101 (2011)

Extracting the Mott gap from energy measurements in trapped atomic gases

Authors: Juan Carrasquilla, Federico Becca

Abstract: We show that the measure of the so-called {\it release-energy}, which is an experimentally accessible quantity, makes it possible to assess the value of the Mott gap in the presence of the confinement potential that is unavoidable in the actual experimental setup. Indeed, the curve of the release-energy as a function of the total number of particles shows kinks that are directly related to the exi… ▽ More We show that the measure of the so-called {\it release-energy}, which is an experimentally accessible quantity, makes it possible to assess the value of the Mott gap in the presence of the confinement potential that is unavoidable in the actual experimental setup. Indeed, the curve of the release-energy as a function of the total number of particles shows kinks that are directly related to the existence of excitation gaps. Calculations are presented within the Gutzwiller approach, but the final results go beyond this simple approximation and represent a genuine feature of the real system. In the case of harmonic confinement, the Mott gaps may be renormalized with respect to the uniform case. On the other hand, in the case of the recently proposed off-diagonal confinement, our results show an almost perfect agreement with the homogeneous case. △ Less

Submitted 11 November, 2010; v1 submitted 25 May, 2010; originally announced May 2010.

Comments: 4 pages and 5 figures

Journal ref: Phys. Rev. A 82, 053609 (2010)

Characterization of the Bose-glass phase in low-dimensional lattices

Authors: Juan Carrasquilla, Federico Becca, Andrea Trombettoni, Michele Fabrizio

Abstract: We study by numerical simulation a disordered Bose-Hubbard model in low-dimensional lattices. We show that a proper characterization of the phase diagram on finite disordered clusters requires the knowledge of probability distributions of physical quantities rather than their averages. This holds in particular for determining the stability region of the Bose-glass phase, the compressible but not… ▽ More We study by numerical simulation a disordered Bose-Hubbard model in low-dimensional lattices. We show that a proper characterization of the phase diagram on finite disordered clusters requires the knowledge of probability distributions of physical quantities rather than their averages. This holds in particular for determining the stability region of the Bose-glass phase, the compressible but not superfluid phase that exists whenever disorder is present. This result suggests that a similar statistical analysis should be performed also to interpret experiments on cold gases trapped in disordered lattices, limited as they are to finite sizes. △ Less

Submitted 5 March, 2010; originally announced March 2010.

Comments: 4+ epsilon pages and 4 figures

Journal ref: Physical Review B 81, 195129 (2010)