Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2406.17844

Numerical exploration of the bootstrap in spin chain models

David E Berenstein, George Hulsey, P. N. Thomas Lloyd

Published Jun 25, 2024• 0 citations
arXivHTMLPDF

Abstract

We analyze the bootstrap approach (a dual optimization method to the variational approach) to one-dimensional spin chains, leveraging semidefinite programming to extract numerical results. We study how correlation functions in the ground state converge to their true values at and away from criticality and at relaxed optimality. We consider the transverse Ising model, the three state Potts model, and other non-integrable spin chains and investigate to what extent semidefinite methods can reliably extract numerical emergent physical data, including conformal central charges, correlation lengths and scaling dimensions. We demonstrate procedures to extract these data and show preliminary results in the various models considered. We compare to exact analytical results and to exact diagonalization when the system volume is small enough. When we attempt to go to the thermodynamic limit, the semidefinite numerical method with translation invariance imposed as a constraint finds the solution with periodic boundary conditions even if these have not been specified. This implies that the determination of all conformal data in correlators has to be handled at finite volume. Our investigation reveals that the approach has practical challenges. In particular, the correlation functions extracted from the optimal solution, which function as slack variables in the optimization, have convergence issues that suggest an underlying exponential complexity in the system size.

DOI
10.48550/arxiv.2406.17844

Related papers

Successes and challenges of using Semidefinite Programming for the study of Spin Chain Hamiltonians

David E Berenstein, P. N. Thomas Lloyd

We study semidefinite programming (SDP) methods to analyze spin chain Hamiltonians. We examine the ground state energy, the first excited charged states and ground state correlators in two simple models: the Ising model in a transverse magnetic field and the closely related 3-state Potts model. Our goal is to understand precisely what the SDP program is doing and when it works well, why it does so. We focus on the following novel ingredients: using charge constraints to obtain excited states and to see if additional constraints from integrable models are effective at improving the method. At criticality we also explore to what extent we can use approximate Virasoro correlators to extract conformal data: the central charge and some critical exponents of charged states. We also use these to identify the location of the phase transition. In the special case where the system is made of free fermions we prove that the SDP finds the exact energy of the ground state and produces the correct two point functions of the fermions. Away from free fermion theories, the SDP gets progressively worse at estimating data beyond the value of the ground state energy (like correlation functions), although it qualitatively matches these. In order to be effective, the SDP seems to run into scaling issues where the amount of input needed scales poorly with the lattice volume.

arXiv:2606.02013•Jun 1, 2026•0 citations

Bootstrapping ground state properties of classical frustrated magnets

Nisarga Paul, Gil Refael

We introduce a method based on semidefinite programming that produces rigorous two-sided bounds on ground state energy densities and correlation functions of translation-invariant classical spin models on infinite lattices. In this method, the challenge of non-convex optimization on an infinite lattice is replaced with a hierarchy of finite-size convex optimizations arising from positivity conditions that any probability distribution over spin configurations must satisfy. This adapts the Lasserre hierarchy in the theory of polynomial optimization to the context of frustrated magnetism, and we prove convergence of this hierarchy in the thermodynamic limit. Our method subsumes the Luttinger--Tisza method and further applies to non-quadratic Hamiltonians and non-Bravais lattices, thus addressing limitations of prior analytical methods. We apply the method to various two-dimensional frustrated spin models, where it brackets the energy densities and observables accurately across large parameter ranges with typical run times of seconds per parameter point.

arXiv:2605.06784•May 7, 2026•0 citations

Tensor networks, quantum spin chains, and quantum field theory

Yijian Zou

Understanding the universality class of continuous phase transition is of central importance in condensed matter physics. In one spatial dimension, the universal properties are encoded in the conformal field theory (CFT), which is in turn specified by the conformal data. In this thesis, we propose a systematic method to extract complete and accurate conformal data from the critical quantum spin chain based on the operator-state correspondence, thus completing the project initiated by Cardy and others in the 80's. \n \nOur method is based on the low-energy eigenstates of the critical quantum spin chain with periodic boundary conditions. First, scaling dimensions and conformal spins are extracted by the energies and momenta. Second, the primary states and conformal towers are identified by using the Koo-Saleur lattice Virasoro generators. Third, we propose a systematic way of identifying lattice operators with CFT operators, which enables us to compute operator product expansion coefficients from the low-energy eigenstates. Finally, the whole approach is generalized to critical quantum spin chains with antiperiodic boundary condition. In order to reduce finite-size corrections, we put forward the periodic uniform matrix product state (puMPS) algorithm which enables us to compute low-energy eigenstates of a critical quantum spin chain up to several hundreds of spins. Our method also enables us to study nonperturbatively the renormalization group flow between two CFTs as well as explore the emergence of extended symmetries beyond conformal symmetry. \n \nWe test our method with the Ising model and its generalization due to O'Brien and Fendley. The latter model is featured by a tricritical point described by the tricritical Ising CFT and a line of critical points interpolating between the Ising CFT and the tricritical Ising CFT. We extract complete conformal data from the two models and find excellent agreement with analytical results. Furthermore, we study the spectral renormalization group flow between the two CFTs nonperturbatively. At the tricritical point, the underlying CFT has an extended symmetry, the superconformal symmetry. We propose lattice operators that correspond to supervirasoro generators and verify their action on low-energy eigenstates. In this way we study the emergence of superconformal symmetry in critical quantum spin chains.

•Jul 7, 2020•0 citations

Finite lattice methods in statistical mechanics

Murray T. Batchelor

This thesis is concerned with the finite lattice study of spin models. The underlying theme in Parts I and II is the exploitation of Bethe ansatz equations to provide finite lattice data well beyond the range normally available. In Part I the Bethe ansatz equations for eigenvalues of the eight-vertex model are solved numerically to yield mass gap data on infinitely long strips of up to 512 sites in width. The finite-size corrections, at criticality, to the free energy per site and polarization gap are found to agree with recent studies of the XXZ spin chain. The leading corrections to the finite-size scaling estimates of the critical line and thermal exponent are also found, providing an explanation of the poor convergence seen in earlier studies. Away from criticality, the linear scaling fields are derived exactly in the full parameter space of the spin system, allowing a thorough test of a recently proposed method of extracting linear scaling fields and related exponents from finite lattice data. In Part II, the numerical solutions of the Bethe ansatz equations for the eigen-energies of the XXZ Hamiltonian on very large chains are used to identify, via conformal invariance, the scaling dimensions of various two-dimensional models. With periodic boundary conditions, eight-vertex and Gaussian model operators are found. The scaling dimensions of the Ashkm-Teller and Potts models are obtained by the exact relating of eigenstates of their quantum Hamiltonians to those of the XXZ chain with modified boundary conditions. Eigenstates of the Ashkin-Teller and Potts models with free boundaries are also obtained, allowing an examination of their critical surface properties. In Part HI the critical behaviour of an Ising model with competing firstand third- nearest neighbour interactions on the square lattice is investigated using the finite lattice method. In the ferromagnetic region, the phase boundary is located with an accuracy at least equal to that of alternative methods. In the antiphase region, distinctive structure in the finite lattice estimators is found over an extended temperature range. However, the nature of the transition remains unclear.

•Jan 1, 1987•2 citations

Emergence of conformal symmetry in quantum spin chains: Antiperiodic boundary conditions and supersymmetry

Yijian Zou, Guifré Vidal

Universal properties of a critical quantum spin chain are encoded in the underlying conformal field theory (CFT). This underlying CFT is fully characterized by its conformal data. We propose a method to extract the conformal data from a critical quantum spin chain with both periodic and antiperiodic boundary conditions (PBC and APBC) based on low-energy eigenstates, generalizing previous work on spin chains with only PBC. First, scaling dimensions and conformal spins are extracted from the energies and momenta of the eigenstates. Second, the Koo-Saleur formula of lattice Virasoro generators is generalized to APBC and used to identify conformal towers. Third, local and string operators on the lattice are identified with CFT operators with PBC and APBC, respectively. Finally, operator product expansion coefficients are extracted by computing matrix elements of lattice primary operators in the low-energy subspaces with PBC and APBC. To go beyond exact diagonalization, tensor network methods based on periodic uniform matrix product states are used. We illustrate our approach with critical and tricritical Ising quantum spin chains. In the latter case, we propose lattice operators that correspond to super-Virasoro generators and verify their action on low-energy eigenstates. In this way, we explore the emergence of superconformal symmetry in the quantum spin chain.

arXiv:1907.10704•Physical review. B./Physical review. B•Jan 28, 2020•25 citations

The Discrete Bootstrap for Quantum Thermodynamics

Sandor, Joshua

A novel bootstrap method is exemplified by lattice theories and quantum mechanics. Positivity of expectation values combined with sets of non-trivial constraint equations and semidefinite optimization rigorously and precisely bound correlation functions. Starting with the example of a one-dimensional integral, the technique is built upon to extract limit cycles in a two-dimensional classical dynamical system and spin correlation functions in the Statistical Ising Model. Further positivity constraints and simplifications arising in a large rank limit are implemented to bound the energy expectation value as a function of temperature in Matrix Quantum Mechanics. Along the way, the S-matrix of Chern-Simons-matter theory is computed as an example of computational simplification at large gauge group rank.

•Jan 1, 2025•0 citations
We analyze the bootstrap approach (a dual optimization method to the variational approach) to one-dimensional spin chains, leveraging semidefinite programming to extract numerical results. We study how correlation functions in the ground state converge to their true values at and away from criticality and at relaxed optimality. We consider the transverse Ising model, the three state Potts model, and other non-integrable spin chains and investigate to what extent semidefinite methods can reliably extract numerical emergent physical data, including conformal central charges, correlation lengths and scaling dimensions. We demonstrate procedures to extract these data and show preliminary results in the various models considered. We compare to exact analytical results and to exact diagonalization when the system volume is small enough. When we attempt to go to the thermodynamic limit, the semidefinite numerical method with translation invariance imposed as a constraint finds the solution with periodic boundary conditions even if these have not been specified. This implies that the determination of all conformal data in correlators has to be handled at finite volume. Our investigation reveals that the approach has practical challenges. In particular, the correlation functions extracted from the optimal solution, which function as slack variables in the optimization, have convergence issues that suggest an underlying exponential complexity in the system size.
We study semidefinite programming (SDP) methods to analyze spin chain Hamiltonians. We examine the ground state energy, the first excited charged states and ground state correlators in two simple models: the Ising model in a transverse magnetic field and the closely related 3-state Potts model. Our goal is to understand precisely what the SDP program is doing and when it works well, why it does so. We focus on the following novel ingredients: using charge constraints to obtain excited states and to see if additional constraints from integrable models are effective at improving the method. At criticality we also explore to what extent we can use approximate Virasoro correlators to extract conformal data: the central charge and some critical exponents of charged states. We also use these to identify the location of the phase transition. In the special case where the system is made of free fermions we prove that the SDP finds the exact energy of the ground state and produces the correct two point functions of the fermions. Away from free fermion theories, the SDP gets progressively worse at estimating data beyond the value of the ground state energy (like correlation functions), although it qualitatively matches these. In order to be effective, the SDP seems to run into scaling issues where the amount of input needed scales poorly with the lattice volume.
We introduce a method based on semidefinite programming that produces rigorous two-sided bounds on ground state energy densities and correlation functions of translation-invariant classical spin models on infinite lattices. In this method, the challenge of non-convex optimization on an infinite lattice is replaced with a hierarchy of finite-size convex optimizations arising from positivity conditions that any probability distribution over spin configurations must satisfy. This adapts the Lasserre hierarchy in the theory of polynomial optimization to the context of frustrated magnetism, and we prove convergence of this hierarchy in the thermodynamic limit. Our method subsumes the Luttinger--Tisza method and further applies to non-quadratic Hamiltonians and non-Bravais lattices, thus addressing limitations of prior analytical methods. We apply the method to various two-dimensional frustrated spin models, where it brackets the energy densities and observables accurately across large parameter ranges with typical run times of seconds per parameter point.
Understanding the universality class of continuous phase transition is of central importance in condensed matter physics. In one spatial dimension, the universal properties are encoded in the conformal field theory (CFT), which is in turn specified by the conformal data. In this thesis, we propose a systematic method to extract complete and accurate conformal data from the critical quantum spin chain based on the operator-state correspondence, thus completing the project initiated by Cardy and others in the 80's. \n \nOur method is based on the low-energy eigenstates of the critical quantum spin chain with periodic boundary conditions. First, scaling dimensions and conformal spins are extracted by the energies and momenta. Second, the primary states and conformal towers are identified by using the Koo-Saleur lattice Virasoro generators. Third, we propose a systematic way of identifying lattice operators with CFT operators, which enables us to compute operator product expansion coefficients from the low-energy eigenstates. Finally, the whole approach is generalized to critical quantum spin chains with antiperiodic boundary condition. In order to reduce finite-size corrections, we put forward the periodic uniform matrix product state (puMPS) algorithm which enables us to compute low-energy eigenstates of a critical quantum spin chain up to several hundreds of spins. Our method also enables us to study nonperturbatively the renormalization group flow between two CFTs as well as explore the emergence of extended symmetries beyond conformal symmetry. \n \nWe test our method with the Ising model and its generalization due to O'Brien and Fendley. The latter model is featured by a tricritical point described by the tricritical Ising CFT and a line of critical points interpolating between the Ising CFT and the tricritical Ising CFT. We extract complete conformal data from the two models and find excellent agreement with analytical results. Furthermore, we study the spectral renormalization group flow between the two CFTs nonperturbatively. At the tricritical point, the underlying CFT has an extended symmetry, the superconformal symmetry. We propose lattice operators that correspond to supervirasoro generators and verify their action on low-energy eigenstates. In this way we study the emergence of superconformal symmetry in critical quantum spin chains.
This thesis is concerned with the finite lattice study of spin models. The underlying theme in Parts I and II is the exploitation of Bethe ansatz equations to provide finite lattice data well beyond the range normally available. In Part I the Bethe ansatz equations for eigenvalues of the eight-vertex model are solved numerically to yield mass gap data on infinitely long strips of up to 512 sites in width. The finite-size corrections, at criticality, to the free energy per site and polarization gap are found to agree with recent studies of the XXZ spin chain. The leading corrections to the finite-size scaling estimates of the critical line and thermal exponent are also found, providing an explanation of the poor convergence seen in earlier studies. Away from criticality, the linear scaling fields are derived exactly in the full parameter space of the spin system, allowing a thorough test of a recently proposed method of extracting linear scaling fields and related exponents from finite lattice data. In Part II, the numerical solutions of the Bethe ansatz equations for the eigen-energies of the XXZ Hamiltonian on very large chains are used to identify, via conformal invariance, the scaling dimensions of various two-dimensional models. With periodic boundary conditions, eight-vertex and Gaussian model operators are found. The scaling dimensions of the Ashkm-Teller and Potts models are obtained by the exact relating of eigenstates of their quantum Hamiltonians to those of the XXZ chain with modified boundary conditions. Eigenstates of the Ashkin-Teller and Potts models with free boundaries are also obtained, allowing an examination of their critical surface properties. In Part HI the critical behaviour of an Ising model with competing firstand third- nearest neighbour interactions on the square lattice is investigated using the finite lattice method. In the ferromagnetic region, the phase boundary is located with an accuracy at least equal to that of alternative methods. In the antiphase region, distinctive structure in the finite lattice estimators is found over an extended temperature range. However, the nature of the transition remains unclear.
Universal properties of a critical quantum spin chain are encoded in the underlying conformal field theory (CFT). This underlying CFT is fully characterized by its conformal data. We propose a method to extract the conformal data from a critical quantum spin chain with both periodic and antiperiodic boundary conditions (PBC and APBC) based on low-energy eigenstates, generalizing previous work on spin chains with only PBC. First, scaling dimensions and conformal spins are extracted from the energies and momenta of the eigenstates. Second, the Koo-Saleur formula of lattice Virasoro generators is generalized to APBC and used to identify conformal towers. Third, local and string operators on the lattice are identified with CFT operators with PBC and APBC, respectively. Finally, operator product expansion coefficients are extracted by computing matrix elements of lattice primary operators in the low-energy subspaces with PBC and APBC. To go beyond exact diagonalization, tensor network methods based on periodic uniform matrix product states are used. We illustrate our approach with critical and tricritical Ising quantum spin chains. In the latter case, we propose lattice operators that correspond to super-Virasoro generators and verify their action on low-energy eigenstates. In this way, we explore the emergence of superconformal symmetry in the quantum spin chain.
A novel bootstrap method is exemplified by lattice theories and quantum mechanics. Positivity of expectation values combined with sets of non-trivial constraint equations and semidefinite optimization rigorously and precisely bound correlation functions. Starting with the example of a one-dimensional integral, the technique is built upon to extract limit cycles in a two-dimensional classical dynamical system and spin correlation functions in the Statistical Ising Model. Further positivity constraints and simplifications arising in a large rank limit are implemented to bound the energy expectation value as a function of temperature in Matrix Quantum Mechanics. Along the way, the S-matrix of Chern-Simons-matter theory is computed as an example of computational simplification at large gauge group rank.