Misplaced Pages

Hamiltonian truncation is a numerical method used to study quantum field theories (QFTs) in ${\displaystyle d\geq 2}$ spacetime dimensions. Hamiltonian truncation is an adaptation of the Rayleigh–Ritz method from quantum mechanics. It is closely related to the exact diagonalization method used to treat spin systems in condensed matter physics. The method is typically used to study QFTs on spacetimes of the form ${\displaystyle \mathbb {R} \times M}$, specifically to compute the spectrum of the Hamiltonian along ${\displaystyle \mathbb {R} }$. A key feature of Hamiltonian truncation is that an explicit ultraviolet cutoff ${\displaystyle \Lambda }$ is introduced, akin to the lattice spacing a in lattice Monte Carlo methods. Since Hamiltonian truncation is a nonperturbative method, it can be used to study strong-coupling phenomena like spontaneous symmetry breaking.

Principles

Energy cutoff

Local quantum field theories can be defined on any manifold. Often, the spacetime of interest includes a copy of ${\displaystyle \mathbb {R} }$, like ${\displaystyle \mathbb {R} ^{d}}$ (flat space), ${\displaystyle \mathbb {R} \times S^{d-1}}$ (an infinite hollow cylinder), ${\displaystyle \mathbb {R} \times \mathbf {T} ^{d-1}}$ (space is taken to be a torus) or even Anti-de Sitter space in global coordinates. On such a manifold we can take time to run along ${\displaystyle \mathbb {R} }$, such that energies are conserved. Solving such a QFT amounts to finding the spectrum and eigenstates of the Hamiltonian H, which is difficult or impossible to do analytically. Hamiltonian truncation provides a strategy to compute the spectrum of H to arbitrary precision. The idea is that many QFT Hamiltonians can be written as the sum of "free" part ${\displaystyle H_{0}}$ and an "interacting" part that describes interactions (for example a ${\displaystyle \phi ^{4}}$ term or a Yukawa coupling), schematically

${\displaystyle H=H_{0}+gV}$

where V can be written as the integral of a local operator ${\displaystyle {\mathcal {V}}}$ over M. There may be multiple interaction terms ${\displaystyle g_{1}V_{1}+g_{2}V_{2}+\ldots }$, but that case generalizes straightforwardly from the case with a single interaction ${\displaystyle gV}$. Hamiltonian truncation amounts to the following recipe:

1. Fix a UV cutoff ${\displaystyle \Lambda }$, and find all eigenstates ${\displaystyle |i\rangle }$ of ${\displaystyle H_{0}}$ with energy ${\displaystyle e_{i}\leq \Lambda }$. Normalize these eigenstates such that ${\displaystyle \langle i|j\rangle =\delta _{ij}}$. Let ${\displaystyle N(\Lambda )}$ be the number of low-energy states.
2. Compute the Hamiltonian explicitly restricted to these low-energy states. The result will be a matrix of size ${\displaystyle N(\Lambda )\times N(\Lambda )}$, explicitly ${\displaystyle H(\Lambda )_{ij}=e_{i}\delta _{ij}+gV_{ij}}$ with ${\displaystyle V_{ij}=\langle i|V|j\rangle .}$
3. Compute the energies and eigenstates of the finite matrix ${\displaystyle H(\Lambda )}$, obeying ${\displaystyle H(\Lambda )|\psi _{\alpha }\rangle =E_{\alpha }(\Lambda )|\psi _{\alpha }\rangle }$.

In a UV-finite quantum field theory, the resulting energies ${\displaystyle E_{\alpha }(\Lambda )}$ have a finite limit as the cutoff ${\displaystyle \Lambda }$ is taken to infinity, so at least in principle the exact spectrum of the Hamiltonian can be recovered. In practice the cutoff ${\displaystyle \Lambda }$ is always finite, and the procedure is performed on a computer.

Range of validity

For a given cutoff ${\displaystyle \Lambda }$, Hamiltonian truncation has a finite range of validity, meaning that cutoff errors become important when the coupling g is too large. To make this precise, let's take R to be the rough size of the manifold M, that is to say that

${\displaystyle {\text{vol}}(M)\sim R^{d-1}}$

up to some c-number coefficient. If the deformation V is the integral of a local operator of dimension ${\displaystyle \Delta }$, then the coupling g will have mass dimension ${\displaystyle =d-\Delta }$, so the redefined coupling

${\displaystyle {\bar {g}}\equiv gR^{d-\Delta }}$

is dimensionless. Depending on the order of magnitude of ${\displaystyle {\bar {g}}}$, we can distinguish three different regimes:

• ${\displaystyle {\bar {g}}\ll 1}$: perturbation theory is valid. Generically, perturbation theory is either asymptotic or it converges up to some value ${\displaystyle |{\bar {g}}|\leq {\bar {g}}_{\text{max}}=O(1)}$. For infinitesimal values of ${\displaystyle {\bar {g}}}$, quantum effects can be neglected.
• ${\displaystyle {\bar {g}}=O({\text{few}})}$: perturbation theory is no longer reliable, but the truncated energies ${\displaystyle E_{\alpha }(\Lambda )}$ provide a good approximation to their continuum values for reasonable values of the cutoff ${\displaystyle \Lambda }$.
• ${\displaystyle {\bar {g}}\gg 1}$: for very large values of ${\displaystyle {\bar {g}}}$ (or equivalently, when the volume of M becomes very large), Hamiltonian truncation only provides good results when the cutoff ${\displaystyle \Lambda }$ is taken to be astronomically large. In practice, this regime is not accessible. This is an avatar of the orthogonality catastrophe.

Truncation errors and ultraviolet divergences

There are two intrinsic but related issues with Hamiltonian truncation:

1. In some cases, the ${\displaystyle E_{\alpha }(\Lambda )}$ do not have a finite limit as ${\displaystyle \Lambda \to \infty }$.
2. Even when the continuum limit ${\displaystyle \lim _{\Lambda \to \infty }E_{\alpha }(\Lambda )}$ exists, we only have access to cutoff data ${\displaystyle E_{\alpha }(\Lambda )}$ for a range of finite values of the cutoff.

The first case is due to ultraviolet divergences of the quantum field theory in question. In this case, cutoff-dependent counterterms must be added to the Hamiltonian H in order to obtain a physically meaningful result. In order to understand the second problem, one can perform perturbative computations to understand the continuum limit analytically.

Let us spell this out using an example. We have in mind a perturbation of the form gV with

${\displaystyle V=\int _{M}d{\mathbf {x} }\;{\mathcal {V}}(t=0,{\mathbf {x} })}$

where ${\displaystyle {\mathcal {V}}(t,\mathbf {x} )}$ is a local operator. Suppose that we want to compute the first corrections to the vacuum energy due to V. In Rayleigh–Schrödinger perturbation theory, we know that

${\displaystyle E_{\Omega }(\Lambda )=0+g\langle \Omega |V|\Omega \rangle +g^{2}E_{\Omega }^{(2)}(\Lambda )+O(g^{3})}$

where

${\displaystyle E_{\Omega }^{(2)}(\Lambda )=-\int _{0}^{\Lambda }{\frac {dE}{E}}\;\rho _{\Omega }(E),\quad \rho _{\Omega }(E)=\sum _{i}\delta (E-e_{i})|\langle i|V|\Omega \rangle |^{2}\geq 0}$

where the sum runs over all states ${\displaystyle |i\rangle }$ other than the vacuum ${\displaystyle |\Omega \rangle }$ itself. Whether this integral converges or not depends on the large-E behavior of the spectral density ${\displaystyle \rho _{\Omega }(E)}$. In turn, this depends on the short-distance behavior of the two-point correlation function of the operator ${\displaystyle {\mathcal {V}}}$. Indeed, we can write

${\displaystyle \int _{0}^{\infty }dE\,e^{-Et}\rho _{\Omega }(E)=\int _{M}d{\mathbf {x} }\int _{M}d{\mathbf {y} }\,\left}$

where ${\displaystyle {\mathcal {V}}(t,\mathbf {x} )=e^{H_{0}t}{\mathcal {V}}(0,\mathbf {x} )e^{-H_{0}t}}$ evolves in Euclidean time in the interaction picture. Hence the large-E behavior of the spectral density encodes the short-time behavior of the ${\displaystyle \langle {\mathcal {V}}(t,\mathbf {x} ){\mathcal {V}}(0,\mathbf {y} )\rangle }$ vacuum correlator, where both x,y are integrated over space. The large-E scaling can be computed in explicit theories; in general it goes as

${\displaystyle \rho _{\Omega }(E)\mathrel {\mathop {\sim } \limits _{\scriptstyle {E\to \infty }}} c\cdot E^{2\Delta _{\mathcal {V}}-d}}$

where ${\displaystyle \Delta _{\mathcal {V}}}$ is the scaling or mass dimension of the operator ${\displaystyle {\mathcal {V}}}$ and c is some constant. There are now two possibilities, depending on the value of ${\displaystyle \gamma \equiv d-2\Delta _{\mathcal {V}}}$:

• If ${\displaystyle \gamma \leq 0}$, the truncated Casimir energy ${\displaystyle E_{\Omega }(\Lambda )}$ diverges in the continuum limit. In this case, a cutoff-dependent counterterm must be added to H in order to cancel this divergence.
• If ${\displaystyle \gamma >0}$, the truncated Casimir energy converges as ${\displaystyle \Lambda \to \infty }$. The truncation error can be estimated to be

${\displaystyle \qquad E_{\Omega }(\infty )-E_{\Omega }(\Lambda )\approx -g^{2}{\frac {c}{\gamma }}\Lambda ^{-\gamma }+O(g^{3}).}$

A similar analysis applies to cutoff errors in excited states and at higher orders in perturbation theory.

Example of the massive scalar φ theory

Quantization

As an example, we can consider a massive scalar field ${\displaystyle \phi (t,\mathbf {x} )}$ on some spacetime ${\displaystyle \mathbb {R} \times M}$, where M is compact (possibly having a boundary). The total metric can be written as

${\displaystyle ds^{2}=-dt^{2}+h_{ij}({\mathbf {x} })dx^{i}dx^{j}.}$

Let's consider the action

${\displaystyle S=\int _{\mathbb {R} }dt\int _{M}{\sqrt {h}}d\mathbf {x} \,{\mathcal {L}},\quad {\mathcal {L}}=-{\frac {1}{2}}\phi \triangle \phi +{\frac {1}{2}}m^{2}\phi ^{2}+g\phi ^{4}}$

where ${\displaystyle \triangle }$ is the Laplacian on ${\displaystyle \mathbb {R} \times M}$. The g=0 theory can be canonically quantized, which endows the field ${\displaystyle \phi }$ with a mode decomposition

${\displaystyle \phi (t,\mathbf {x} )=\sum _{n\in \mathbb {N} }{\frac {1}{\sqrt {2\omega _{n}}}}\left(a_{n}\,e^{-i\omega _{n}t}f_{n}(\mathbf {x} )+a_{n}^{\dagger }\,e^{i\omega _{n}t}f_{n}(\mathbf {x} )^{*}\right)}$

where the creation and annihilation operators obey canonical commutation relations ${\displaystyle =\delta _{mn}}$. The single-particle energies ${\displaystyle \omega _{n}>0}$ and the mode functions ${\displaystyle f_{n}(\mathbf {x} )}$ depend on the spatial manifold M. The Hamiltonian at t=0 is then given by

${\displaystyle H=H_{0}+gV,\quad H_{0}=\sum _{n\in \mathbb {N} }\omega _{n}a_{n}^{\dagger }a_{n}\quad {\text{and}}\quad V=\int _{M}{\sqrt {h}}d\mathbf {x} \;\phi (t=0,\mathbf {x} )^{4}.}$

Hamiltonian truncation

The Hilbert space of the ${\displaystyle g=0}$ theory is the Fock space of the modes ${\displaystyle \{a_{n}^{\dagger }\}}$. That is to say that there exists a vacuum state ${\displaystyle |\Omega \rangle }$ obeying ${\displaystyle a_{n}|\Omega \rangle =0}$ for all n, and on top of that there are single- and multi-particle states. Explicitly, a general eigenstate of ${\displaystyle H_{0}}$ is labeled by a tuple ${\displaystyle \{k_{n}\}}$ of occupation numbers:

${\displaystyle |{\mathbf {k} }\rangle =\prod _{n\in \mathbb {N} }{\frac {1}{\sqrt {k_{n}!}}}(a_{n}^{\dagger })^{k_{n}}|\Omega \rangle }$

where the ${\displaystyle k_{n}}$ can take values in the integers: ${\displaystyle k_{n}\in \{0,1,2,\ldots \}}$. Such a state has energy

${\displaystyle H_{0}|{\mathbf {k} }\rangle =e({\mathbf {k} })|{\mathbf {k} }\rangle ,\quad e({\mathbf {k} })=\sum _{n\in \mathbb {N} }k_{n}\omega _{n}}$

so finding a basis of low-energy states amounts to finding all tuples ${\displaystyle \{k_{n}\}}$ obeying ${\displaystyle e(\mathbf {k} )\leq \Lambda }$. Let's denote all such states schematically as ${\displaystyle |i\rangle }$. Next, the matrix elements ${\displaystyle V_{ij}}$ can be computed explicitly using the canonical commutation relations. Finally, the explicit Hamiltonian ${\displaystyle H(\Lambda )_{ij}=e_{i}\delta _{ij}+gV_{ij}}$ has to be diagonalized.

The resulting spectra can be used to study precision physics. Depending on the values of g and ${\displaystyle m^{2}}$, the above ${\displaystyle \phi ^{4}}$ theory can be in a symmetry-preserving or a symmetry-broken phase, which can be studied explicitly using the above algorithm. The continuous phase transition between these two phases can also be analyzed, in which case the spectrum and eigenstates of H contain information about the conformal field theory of the Ising universality class.

Special cases

Truncated Conformal Space Approach

The truncated conformal space approach (TCSA) is a version of the Hamiltonian truncation that applies to perturbed conformal field theories. This approach was introduced by Yurov and Al. Zamolodchikov in 1990 and has become a standard ingredient used to study two-dimensional QFTs. The d-dimensional version of TCSA was first studied in 2014.

A RG flow emanating from a conformal field theory (CFT) is described by an action

${\displaystyle S=S_{\text{CFT}}+g\int d^{d}x\,{\mathcal {V}}(x)}$

where ${\displaystyle {\mathcal {V}}}$ is a scalar operator in the CFT of scaling dimension ${\displaystyle \Delta _{\mathcal {V}}\leq d}$. At large distances, such theories are strongly coupled. It is convenient to study such RG flows on the cylinder ${\displaystyle \mathbb {R} \times S^{d-1}}$, taking the sphere to have radius R and endowing the full space with coordinates ${\displaystyle (t,\mathbf {n} )}$. The reason is that the unperturbed (g=0) theory admits a simple description owing to radial quantization. Schematically, states ${\displaystyle |i\rangle }$ on the cylinder are in one-to-one correspondence with local operators ${\displaystyle {\mathcal {O}}_{i}}$ inserted at the origin of flat space:

${\displaystyle |i\rangle =\lim _{x\to 0}{\mathcal {O}}_{i}(x)|\Omega \rangle }$

where ${\displaystyle |\Omega \rangle }$ is the CFT vacuum state. The Hamiltonian on the cylinder is precisely the dilatation operator D of the CFT: the unperturbed energies are given by

${\displaystyle H_{0}|i\rangle ={\frac {\Delta _{i}}{R}}|i\rangle }$

where ${\displaystyle \Delta _{i}}$ is the scaling dimension of the operator ${\displaystyle {\mathcal {O}}_{i}}$. Finally, the matrix elements of the deformation V

${\displaystyle V_{ij}=R^{d-1}\int _{S^{d-1}}d\mathbf {n} \;\langle i|{\mathcal {V}}(t=0,\mathbf {n} )|j\rangle }$

are proportional to OPE coefficients ${\displaystyle {\mathcal {V}}\times {\mathcal {O}}_{j}\sim {\mathcal {O}}_{i}}$ in the original CFT.

Lightcone truncation methods

Real-time QFTs are often studied in lightcone coordinates

${\displaystyle ds^{2}=2dx^{+}dx^{-}-(dx^{i})^{2}.}$

Although the spectrum of the lightcone Hamiltonian ${\displaystyle P_{+}=i\partial /\partial x^{+}}$ is continuous, it is still possible to compute certain observables using truncation methods. The most commonly used scheme, used when the UV theory is conformal, is known as lightcone conformal truncation (LCT). Notably, the spatial manifold M is non-compact in this case, unlike the equal-time quantization described previously. See also the page for light-front computational methods, which describes related computational setups.

Numerical implementation

Hamiltonian truncation computations are normally performed using a computer algebra system, or a programming language like Python or C++.

The number of low-energy states ${\displaystyle N(\Lambda )}$ tends to grow rapidly with the UV cutoff, and it is common to perform Hamiltonian truncation computations taking into account several thousand states. Nonetheless, one is often only interested in the first O(10) energies and eigenstates of H. Instead of diagonalizing the full Hamiltonian explicitly (which is numerically very costly), approximation methods like Arnoldi iteration and the Lanczos algorithm are commonly used.

In some cases, it is not possible to orthonormalize the low-energy states ${\displaystyle |i\rangle }$, either because this is numerically expensive or because the underlying Hilbert space is not positive definite. In that case, one has to solve the generalized eigenvalue problem

${\displaystyle H_{ij}(\Lambda )v_{\alpha }^{j}=E_{\alpha }(\Lambda )\,G_{ij}v_{\alpha }^{j}}$

where ${\displaystyle H_{ij}(\Lambda )=\langle i|H(\Lambda )|j\rangle }$ and ${\displaystyle G_{ij}=\langle i|j\rangle }$ is the Gram matrix of the theory. In this formulation, the eigenstates of the truncated Hamiltonian are ${\displaystyle |\psi _{\alpha }\rangle =v_{\alpha }^{i}|i\rangle }$.

In practice, it is important to keep track of the symmetries of the theory, that is to say all generators ${\displaystyle G_{i}}$ that satisfy ${\displaystyle =0}$. There are two types of symmetries in Hamiltonian truncation:

1. Global symmetries, for instance the ${\displaystyle \mathbb {Z} _{2}}$ symmetry ${\displaystyle \phi \mapsto -\phi }$ in ${\displaystyle \phi ^{4}}$ theory.
2. Symmetries of the spatial manifold M, for instance the orthogonal group ${\displaystyle O(d)}$ when ${\displaystyle M=S^{d-1}}$.

When all states are organized in symmetry sectors with respect to the ${\displaystyle G_{i}}$ the Hamiltonian is block diagonal, so the effort required to diagonalize H is reduced.

References

1. Whitsitt, Seth; Schuler, Michael; Henry, Louis-Paul; Läuchli, Andreas M.; Sachdev, Subir (2017-07-24). "Spectrum of the Wilson-Fisher conformal field theory on the torus". Physical Review B. 96 (3). American Physical Society (APS): 035142. arXiv:1701.03111. Bibcode:2017PhRvB..96c5142W. doi:10.1103/physrevb.96.035142. ISSN 2469-9950. S2CID 119336211.
2. Philip Giokas; Gerard Watts (2011). "The renormalisation group for the truncated conformal space approach on the cylinder". arXiv:1106.2448 .
3. ^ Matthijs Hogervorst, Slava Rychkov and Balt C. van Rees (2015). "Truncated conformal space approach in d dimensions: A cheap alternative to lattice field theory?". Phys. Rev. D. 91 (2): 025005. arXiv:1409.1581. Bibcode:2015PhRvD..91b5005H. doi:10.1103/PhysRevD.91.025005. S2CID 119273053.
4. Joan Elias-Miro, Slava Rychkov and Lorenzo G. Vitale (2017). "NLO Renormalization in the Hamiltonian Truncation". Phys. Rev. D. 96 (6): 065024. arXiv:1706.09929. Bibcode:2017PhRvD..96f5024E. doi:10.1103/PhysRevD.96.065024. S2CID 119367912.
5. Daniel Rutter; Balt C. van Rees (2018). "Counterterms in Truncated Conformal Perturbation Theory". arXiv:1803.05798 .
6. Hogervorst, Matthijs; Meineri, Marco; Penedones, João; Vaziri, Kamran Salehi (2021). "Hamiltonian truncation in Anti-de Sitter spacetime". Journal of High Energy Physics. 2021 (8). Springer: 63. arXiv:2104.10689. Bibcode:2021JHEP...08..063H. doi:10.1007/jhep08(2021)063. ISSN 1029-8479. S2CID 233346724.
7. Slava Rychkov and Lorenzo G. Vitale (2015). "Hamiltonian truncation study of the ${\displaystyle \phi ^{4}}$ theory in two dimensions". Phys. Rev. D. 91: 085011. arXiv:1412.3460. doi:10.1103/PhysRevD.91.085011. S2CID 119343646.
8. Joan Elias-Miro, Slava Rychkov and Lorenzo G. Vitale (2017). "High-Precision Calculations in Strongly Coupled Quantum Field Theory with Next-to-Leading-Order Renormalized Hamiltonian Truncation". JHEP. 2017 (10): 213. arXiv:1706.06121. Bibcode:2017JHEP...10..213E. doi:10.1007/JHEP10(2017)213. S2CID 56094893.
9. Joan Elias-Miro and Edward Hardy (2020). "Exploring Hamiltonian Truncation in d=2+1". Phys. Rev. D. 102 (6): 065001. arXiv:2003.08405. Bibcode:2020PhRvD.102f5001E. doi:10.1103/PhysRevD.102.065001. S2CID 213004924.
10. V.P. Yurov and Al. B. Zamolodchikov (1990). "Truncated Conformal Space Approach to Scaling Lee-Yang Model". Int. J. Mod. Phys. A. 5 (16): 3221–3246. Bibcode:1990IJMPA...5.3221Y. doi:10.1142/S0217751X9000218X.
11. James, Andrew J A; Konik, Robert M; Lecheminant, Philippe; Robinson, Neil J; Tsvelik, Alexei M (2018-02-26). "Non-perturbative methodologies for low-dimensional strongly-correlated systems: From non-Abelian bosonization to truncated spectrum methods". Reports on Progress in Physics. 81 (4). IOP Publishing: 046002. arXiv:1703.08421. Bibcode:2018RPPh...81d6002J. doi:10.1088/1361-6633/aa91ea. ISSN 0034-4885. S2CID 206095591.
12. Emanuel Katz, Zuhair U. Khandker and Matthew T. Walters (2016). "A Conformal Truncation Framework for Infinite-Volume Dynamics". JHEP. 2016 (7): 140. arXiv:1604.01766. Bibcode:2016JHEP...07..140K. doi:10.1007/JHEP07(2016)140. S2CID 92981707.
13. Nikhil Anand; A. Liam Fitzpatrick; Emanuel Katz; Zuhair U. Khandker; Matthew T. Walters; Yuan Xin (2020). "Introduction to Lightcone Conformal Truncation: QFT Dynamics from CFT Data". arXiv:2005.13544 .

• Quantum field theory
• Quantum mechanics