Light-front computational methods

A systematic approach to discretization of the eigenvalue problem is the DLCQ method originally suggested by Pauli and Brodsky. In essence it is the replacement of integrals by trapezoidal approximations, with equally-spaced intervals in the longitudinal and transverse momenta

:$p^+\backslash rightarrow\backslash frac\{2\backslash pi\}\{L\}n\backslash ,,\backslash ;\backslash ;\; \backslash vec\{p\}\_\backslash perp\backslash rightarrow(\backslash frac\{\backslash pi\}\{L\_\backslash perp\}n\_x,\backslash frac\{\backslash pi\}\{L\_\backslash perp\}n\_y),$

corresponding to periodic boundary conditions on the intervals resolution $K\backslash equiv\backslash frac\{L\}\{2\backslash pi\}P^+$. The combination of momentum components that defines $H\_\{\backslash rm\; LC\}=P^+\backslash mathcal\{P\}^-$ is then independent of $L$. The longitudinal momentum fractions $x\_i\backslash equiv\; p\_i^+/P^+$ become ratios of integers $n\_i/K$. Because the $n\_i$ are all positive, DLCQ automatically limits the number of particles to be no more than $K$. When a limit on transverse momentum is supplied via a chosen cutoff, a finite matrix problem is obtained; however, the matrix may be too large for present numerical techniques. An explicit truncation in particle number, the light-cone equivalent of the Tamm—Dancoff approximation, can then be made. Large basis sizes require special techniques for matrix diagonalization; the one typically used is the Lanczos algorithm. For the case of one space dimension, one can readily solve for the hadron spectrum of QCD for any quark masses and colors.

Most DLCQ calculations are done without zero modes. However, in principle, any DLCQ basis with periodic boundary conditions may include them as constrained modes, dependent on the other modes with nonzero momentum. The constraint comes from the spatial average of the Euler–Lagrange equation for the field. This constraint equation can be difficult to solve, even for the simplest theories. However, an approximate solution can be found, consistent with the underlying approximations of the DLCQ method itself.{{ cite journal | author1=S.S. Chabysheva|author2= J.R. Hiller | title= Zero momentum modes in discrete light-cone quantization | journal=Physical Review D | volume= 79 |issue= 9 | page= 096012 | year=2009 | doi=10.1103/PhysRevD.79.096012|bibcode= 2009PhRvD..79i6012C | arxiv= 0903.1239 |s2cid= 9267391 }} This solution generates the effective zero-mode interactions for the light-front Hamiltonian.

Calculations in the massive sector that are done without zero modes will usually yield the correct answer. The neglect of zero modes merely worsens the convergence. One exception is that of cubic scalar theories, where the spectrum extends to minus infinity. A DLCQ calculation without zero modes will require careful extrapolation to detect this infinity, whereas a calculation that includes zero modes yields the correct result immediately. The zero modes are avoided if one uses antiperiodic boundary conditions.

The supersymmetric form of DLCQ (SDLCQ) is specifically designed to maintain supersymmetry in the discrete approximation. Ordinary DLCQ violates supersymmetry by terms that do not survive the continuum limit. The SDLCQ construction discretizes the supercharge $Q^-$ and defines the Hamiltonian $\backslash mathcal\{P\}^-$ by the superalgebra relation $\backslash mathcal\{P\}^-=\backslash \{Q^-,Q^-\backslash \}/2\backslash sqrt\{2\}$. The range of transverse momentum is limited by a simple cutoff in the momentum value. Effects of zero modes are expected to cancel.

In addition to calculations of spectra, this technique can be used to calculate expectation values. One such quantity, a correlator $\backslash langle\; T^\{++\}(x)\; T^\{++\}(y)\backslash rangle$ of the stress energy tensor, has been computed as a test of a Maldacena conjecture. A very efficient Lanczos-based method was developed for this calculation. The most recent results provide direct evidence for the conjecture.

The transverse lattice method{{ cite journal | author1=W. A. Bardeen|author2= R. B. Pearson | title= Local Gauge Invariance and the Bound State Nature of Hadrons | journal=Physical Review D | volume= 14 |issue= 2 | pages= 547–551 | year=1976 | doi= 10.1103/PhysRevD.14.547|bibcode= 1976PhRvD..14..547B }}{{ cite journal | author1=M. Burkardt|author2= S. Dalley | title=The relativistic bound state problem in QCD: transverse lattice methods | journal=Progress in Particle and Nuclear Physics | volume= 48 |issue= 2 | pages= 317–362 | year=2002 | doi= 10.1016/S0146-6410(02)00140-0|bibcode= 2002PrPNP..48..317B | arxiv= hep-ph/0112007 |s2cid= 15164414 }} brings together two powerful ideas in quantum field theory: light-front Hamiltonian quantization and lattice gauge theory. Lattice gauge theory is a very popular means of regulating for calculation the gauge theories that describe all visible matter in the universe; in particular, it manifestly demonstrates the linear confinement of QCD that holds quarks and gluons inside the protons and neutrons of the atomic nucleus. In general, to obtain solutions of a quantum field theory, with its continuously infinite degrees of freedom, one must put kinematical cutoffs or other restrictions on the space of quantum states. To remove the errors this introduces, one may then extrapolate these cutoffs, provided a continuum limit exists, and/or renormalize observables to account for degrees of freedom above the cutoff. For the purposes of Hamiltonian quantization, one must have a continuous time direction. In the case of light-front Hamiltonian quantization, in addition to continuous light-front time $x^+$, it is necessary to keep the $x^-$ direction continuous if one wants to preserve the manifest Lorentz boost invariance in one direction and to include small light-front energies $p^-$. Therefore, at most one can impose a lattice cutoff on the remaining transverse spatial directions. Such a transverse lattice gauge theory was first suggested by Bardeen and Pearson in 1976.

Most practical calculations performed with transverse lattice gauge theory have utilized one further ingredient: the color-dielectric expansion. A dielectric formulation is one in which the gauge group elements, whose generators are the gluon fields in the case of QCD, are replaced by collective (smeared, blocked, etc.) variables which represent an average over their fluctuations on short distance scales. These dielectric variables are massive, carry color, and form an effective gauge field theory with classical action minimized at zero field, meaning that color flux is expelled from the vacuum at the classical level. This maintains the triviality of the light-front vacuum structure, but arises only for a low momentum cutoff on the effective theory (corresponding to transverse lattice spacings of order 1/2 fm in QCD). As a result, the effective cutoff Hamiltonian is initially poorly constrained. The color-dielectric expansion, together with requirements of Lorentz symmetry restoration, has nevertheless been successfully used to organize the interactions in the Hamiltonian in a way suitable for practical solution. The most accurate spectrum of large-$N\_c$ glueballs has been obtained in this way, and as well as pion light-front wave functions in agreement with a range of experimental data.

The basis light-front quantization (BLFQ) approach{{ cite journal | author1= J. P. Vary|author2= H. Honkanen|author3= J. Li|author4= P. Maris|author5= S. J. Brodsky|author6= A. Harindranath|author7= G. F. de Teramond|author8= P. Sternberg | title= Hamiltonian light-front field theory in a basis function approach | journal=Physical Review C | volume= 81 |issue= 3| page= 035205 | year=2010 | doi= 10.1103/PhysRevC.81.035205 | bibcode=2010PhRvC..81c5205V|arxiv= 0905.1411|s2cid= 33206182}} uses expansions in products of single-particle basis functions to represent the Fock-state wave functions. Typically, the longitudinal ($x^-$) dependence is represented in the DLCQ basis of plane waves, and the transverse dependence is represented by two-dimensional harmonic oscillator functions. The latter are ideal for applications to confining cavities and are consistent with light-front holographic QCD.{{ cite journal | author1= G. F. de Teramond|author2= S. J. Brodsky | title= Hadronic spectrum of a holographic dual of QCD | journal=Physical Review Letters | volume= 94 |issue= 20 | page= 201601 | year=2005 | doi= 10.1103/PhysRevLett.94.201601 | bibcode=2005PhRvL..94t1601D | pmid=16090235| arxiv= hep-th/0501022 |s2cid= 11006078 }}{{ cite journal | author1= G. F. de Teramond|author2= S. J. Brodsky | title= Light-Front Holography: A First Approximation to QCD | journal=Physical Review Letters | volume= 102 |issue= 8 | page= 081601 | year=2009 | doi= 10.1103/PhysRevLett.102.081601 | pmid=19257731|bibcode= 2009PhRvL.102h1601D | arxiv= 0809.4899 |s2cid= 33855116 }}{{ cite journal | author1= S. J. Brodsky|author2= F. -G. Cao|author3= G. F. de Teramond | title= AdS/QCD and Applications of Light-Front Holography | journal= Communications in Theoretical Physics | volume= 57 |issue= 4| pages= 641–664 | year=2012 | doi= 10.1088/0253-6102/57/4/21 | bibcode=2012CoTPh..57..641S| arxiv= 1108.5718 |s2cid= 73629251}}{{ cite journal | author1= T. Gutsche|author2= V. E. Lyubovitskij|author3= I. Schmidt|author4= A. Vega | title= Nucleon resonances in AdS/QCD | journal=Physical Review D | volume= 87 | page= 016017 | year=2013 |issue= 1| doi= 10.1103/PhysRevD.87.016017 | bibcode=2013PhRvD..87a6017G| arxiv= 1212.6252 |s2cid= 118685470}}{{ cite journal | author1= T. Gutsche|author2= V. E. Lyubovitskij|author3= I. Schmidt|author4= A. Vega | title= Chiral Symmetry Breaking and Meson Wave Functions in Soft-Wall AdS/QCD | journal=Physical Review D | volume= 87 |issue= 5| page= 056001 | year=2013 | doi= 10.1103/PhysRevD.87.056001 | bibcode=2013PhRvD..87e6001G| arxiv= 1212.5196 |s2cid= 118377538}} The use of products of single particle basis functions is also convenient for incorporation of boson and fermion statistics, because the products are readily (anti)symmetrized. By employing two-dimensional basis functions with rotational symmetry about the longitudinal direction (where the harmonic oscillator functions serve as an example), one preserves the total angular momentum projection quantum number which facilitates determination of the total angular momentum of the mass eigenstates. For applications without an external cavity, where transverse momentum is conserved, a Lagrange multiplier method is used to separate the relative transverse motion from the total system's motion.

The first application of BLFQ to QED solved for the electron in a two-dimensional transverse confining cavity and showed how the anomalous magnetic moment behaved as a function of the strength of the cavity.{{ cite journal | author1= H. Honkanen|author2= P. Maris|author3= J. P. Vary|author4= S. J. Brodsky | title= Electron in a transverse harmonic cavity | journal=Physical Review Letters | volume= 106 |issue= 6| page= 061603 | year=2011 | doi= 10.1103/PhysRevLett.106.061603 |pmid= 21405457|arxiv= 1008.0068| bibcode=2011PhRvL.106f1603H| doi-access= free }} The second application of BLFQ to QED solved for the electron's anomalous magnetic moment in free space{{ cite journal | author1=X. Zhao|author2= H. Honkanen|author3= P. Maris|author4= J.P. Vary|author5= S.J. Brodsky | title=Electron Anomalous Magnetic Moment in Basis Light-Front Quantization Approach | journal=Few Body Systems | volume= 52 |issue= 3–4| pages= 339–344 | year=2012 | doi= 10.1007/s00601-011-0273-2 | bibcode=2012FBS....52..339Z| arxiv= 1110.0553 |s2cid= 14813759}}{{ cite journal | author1=X. Zhao|author2= H. Honkanen|author3= P. Maris|author4= J. P. Vary|author5= S. J. Brodsky | title=Electron g-2 in Light-Front Quantization |journal= Physics Letters B|volume= 737|pages= 65–69| arxiv=1402.4195|bibcode= 2014PhLB..737...65Z|year= 2014|doi= 10.1016/j.physletb.2014.08.020|s2cid= 44229174}} and demonstrated agreement with the Schwinger moment in the appropriate limit.

The extension of BLFQ to the time-dependent regime, namely, time-dependent BLFQ (tBLFQ) is straightforward and is currently under active development. The goal of tBLFQ is to solve light-front field theory in real-time (with or without time-dependent background fields). The typical application areas include intense lasers (see Light-front quantization#Intense lasers}) and relativistic heavy-ion collisions.

The light-front coupled cluster (LFCC) method{{ cite journal | author1= S. S. Chabysheva|author2= J. R. Hiller | title= A Light-Front Coupled-Cluster Method for the Nonperturbative Solution of Quantum Field Theories | journal=Physics Letters B | volume= 711 |issue= 5 | pages= 417–422 | year=2012 | doi= 10.1016/j.physletb.2012.04.032 |arxiv= 1103.0037 | bibcode=2012PhLB..711..417C| doi-access= free }} is a particular form of truncation for the infinite coupled system of integral equations for light-front wave functions. The system of equations that comes from the field-theoretic Schrödinger equation also requires regularization, to make the integral operators finite. The traditional Fock-space truncation of the system, where the allowed number of particles is limited, typically disrupts the regularization by removing infinite parts that would otherwise cancel against parts that are retained. Although there are ways to circumvent this, they are not completely satisfactory.

The LFCC method avoids these difficulties by truncating the set of equations in a very different way. Instead of truncating the number of particles, it truncates the way in which wave functions are related to each other; the wave functions of higher Fock states are determined by the lower-state wave functions and the exponentiation of an operator $T$. Specifically, the eigenstate is written in the form $\backslash sqrt\{Z\}e^T|\backslash phi\backslash rangle$, where $\backslash sqrt\{Z\}$ is a normalization factor and $|\backslash phi\backslash rangle$ is a state with the minimal number of constituents. The operator $T$ increases particle number and conserves all relevant quantum numbers, including light-front momentum. This is in principle exact but also still infinite, because $T$ can have an infinite number of terms. Zero modes can be included by inclusion of their creation as terms in $T$; this generates a nontrivial vacuum as a generalized coherent state of zero modes.

The truncation made is a truncation of $T$. The original eigenvalue problem becomes a finite-sized eigenvalue problem for the valence state $|\backslash phi\backslash rangle$, combined with auxiliary equations for the terms retained in $T$:

:$P\_v\backslash overline\{\backslash mathcal\{P\}^-\}|\backslash phi\backslash rangle=\backslash frac\{M^2+P\_\backslash perp^2\}\{P^+\}|\backslash phi\backslash rangle,\; \backslash ;\backslash ;\backslash ;\backslash ;\; (1-P\_v)\backslash overline\{\backslash mathcal\{P\}^-\}|\backslash phi\backslash rangle=0.$

Here $P\_v$ is a projection onto the valence sector, and $\backslash overline\{\backslash mathcal\{P\}^-\}\backslash equiv\; e^\{-T\}\backslash mathcal\{P\}^-\; e^T$ is the LFCC effective Hamiltonian. The projection $1-P\_v$ is truncated to provide just enough auxiliary equations to determine the functions in the truncated $T$ operator. The effective Hamiltonian is computed from its Baker–Hausdorff expansion $\backslash overline\{\backslash mathcal\{P\}^-\}=\backslash mathcal\{P\}^-+[\backslash mathcal\{P\}^-,T]+\backslash frac12[[\backslash mathcal\{P\}^-,T],T]+\backslash cdots$, which can be terminated at the point where more particles are being created than are kept by the truncated projection $1-P\_v$. The use of the exponential of $T$ rather than some other function is convenient, not only because of the Baker–Hausdorff expansion but more generally because it is invertible; in principle, other functions could be used and would also provide an exact representation until a truncation is made.

The truncation of $T$ can be handled systematically. Terms can be classified by the number of annihilated constituents and the net increase in particle number. For example, in QCD the lowest-order contributions annihilate one particle and increase the total by one. These are one-gluon emission from a quark, quark pair creation from one gluon, and gluon pair creation from one gluon. Each involves a function of relative momentum for the transition from one to two particles. Higher order terms annihilate more particles and/or increase the total by more than one. These provide additional contributions to higher-order wave functions and even to low-order wave functions for more complicated valence states. For example, the wave function for the $|q\backslash bar\{q\}g\backslash rangle$ Fock state of a meson can have a contribution from a term in $T$ that annihilates a $q\backslash bar\{q\}$ pair and creates a pair plus a gluon, when this acts on the meson valence state $|q\backslash bar\{q\}\backslash rangle$.

The mathematics of the LFCC method has its origin in the many-body coupled cluster method used in nuclear physics and quantum chemistry.{{ cite journal | author1= R.J. Bartlett|author2= M. Musial | title=Coupled-cluster theory in quantum chemistry | journal=Reviews of Modern Physics | volume= 79 | pages= 291–352 | year=2007 |issue= 1 | doi= 10.1103/RevModPhys.79.291 | bibcode=2007RvMP...79..291B}} The physics is, however, quite different. The many-body method works with a state of a large number of particles and uses the exponentiation of $T$ to build in correlations of excitations to higher single-particle states; the particle number does not change. The LFCC method starts from a small number of constituents in a valence state and uses $e^T$ to build states with more particles; the method of solution of the valence-state eigenvalue problem is left unspecified.

The computation of physical observables from matrix elements of operators requires some care. Direct computation would require an infinite sum over Fock space. One can instead borrow from the many-body coupled cluster method a construction that computes expectation values from right and left eigenstates. This construction can be extended to include off-diagonal matrix elements and gauge projections. Physical quantities can then be computed from the right and left LFCC eigenstates.

Renormalization concepts, especially the renormalization group methods in quantum theories and statistical mechanics, have a long history and a very broad scope. The concepts of renormalization that appear useful in theories quantized in the front form of dynamics are essentially of two types, as in other areas of theoretical physics. The two types of concepts are associated with two types of theoretical tasks involved in applications of a theory. One task is to calculate observables (values of operationally defined quantities) in a theory that is unambiguously defined. The other task is to define a theory unambiguously. This is explained below.

Since the front form of dynamics aims at explaining hadrons as bound states of quarks and gluons, and the binding mechanism is not describable using perturbation theory, the definition of a theory needed in this case cannot be limited to perturbative expansions. For example, it is not sufficient to construct a theory using regularization of loop integrals order-by-order and correspondingly redefining the masses, coupling constants, and field normalization constants also order-by-order. In other words, one needs to design the Minkowski space-time formulation of a relativistic theory that is not based on any a priori perturbative scheme. The front form of Hamiltonian dynamics is perceived by many researchers as the most suitable framework for this purpose among the known options.

The desired definition of a relativistic theory involves calculations of as many observables as one must use in order to fix all the parameters that appear in the theory. The relationship between the parameters and observables may depend on the number of degrees of freedom that are included in the theory.

For example, consider virtual particles in a candidate formulation of the theory. Formally, special relativity requires that the range of momenta of the particles is infinite because one can change the momentum of a particle by an arbitrary amount through a change of frame of reference. If the formulation is not to distinguish any inertial frame of reference, the particles must be allowed to carry any value of momentum. Since the quantum field modes corresponding to particles with different momenta form different degrees of freedom, the requirement of including infinitely many values of momentum means that one requires the theory to involve infinitely many degrees of freedom. But for mathematical reasons, being forced to use computers for sufficiently precise calculations, one has to work with a finite number of degrees of freedom. One must limit the momentum range by some cutoff.

Setting up a theory with a finite cutoff for mathematical reasons, one hopes that the cutoff can be made sufficiently large to avoid its appearance in observables of physical interest, but in local quantum field theories that are of interest in hadronic physics the situation is not that simple. Namely, particles of different momenta are coupled through the dynamics in a nontrivial way, and the calculations aiming at predicting observables yield results that depend on the cutoffs. Moreover, they do so in a diverging fashion.

There may be more cutoff parameters than just for momentum. For example, one may assume that the volume of space is limited, which would interfere with translation invariance of a theory, or assume that the number of virtual particles is limited, which would interfere with the assumption that every virtual particle may split into more virtual particles. All such restrictions lead to a set of cutoffs that becomes a part of a definition of a theory.

Consequently, every result of a calculation for any observable $X\_\{\backslash rm\; observable\}(\backslash mu)$ characterized by its physical scale $\backslash mu$ has the form of a function of the set of parameters of the theory, $p$, the set of cutoffs, say $\backslash Lambda$, and the scale $\backslash mu$. Thus, the results take the form

:$X\_\{\backslash rm\; observable\}(\backslash mu)\; =\; X\_\{\backslash rm\; theory\}(p,\; \backslash Lambda,\; \backslash mu)\; \backslash ,\; .$

However, experiments provide values of observables that characterize natural processes irrespective of the cutoffs in a theory used to explain them. If the cutoffs do not describe properties of nature and are introduced merely for making a theory computable, one needs to understand how the dependence on $\backslash Lambda$ may drop out from $X\_\{\backslash rm\; theory\}(p,\backslash Lambda,\; \backslash mu)$. The cutoffs may also reflect some natural features of a physical system at hand, such as in the model case of an ultraviolet cutoff on the wave vectors of sound waves in a crystal due to the spacing of atoms in the crystal lattice. The natural cutoffs may be of enormous size in comparison to the scale $\backslash mu$. Then, one faces the question of how it happens in the theory that its results for observables at scale $\backslash mu$ are not also of the enormous size of the cutoff and, if they are not, then how they depend on the scale $\backslash mu$.

The two types of concepts of renormalization mentioned above are associated with the following two questions:

• How should the parameters $p$ depend on the cutoffs $\backslash Lambda$ so that all observables $X(p,\backslash Lambda,\backslash mu)$ of physical interest do not depend on $\backslash Lambda$, including the case where one removes the cutoffs by sending them formally to infinity?
• What is the required set of parameters $p$?

The renormalization group concept associated with the first question{{ cite journal | author1= E. C. G. Stueckelberg|author2= A. Petermann | title=Normalization of constants in the quanta theory | journal=Helvetica Physica Acta | volume= 26 | page= 499 | year=1953 }}{{ cite journal | author1=M. Gell-Mann|author2= F. E. Low | title=Quantum electrodynamics at small distances | journal=Physical Review | volume= 95 |issue= 5 | pages= 1300–1312 | year=1954 | doi= 10.1103/PhysRev.95.1300|bibcode= 1954PhRv...95.1300G |url= https://authors.library.caltech.edu/60469/1/PhysRev.95.1300.pdf }} predates the concept associated with the second question.{{ cite journal | author=K. G. Wilson | title=Model Hamiltonians for Local Quantum Field Theory | journal=Physical Review | volume= 140 | issue=2B | page= B445 | year=1965 | doi= 10.1103/PhysRev.140.B445 | bibcode=1965PhRv..140..445W }}{{ cite journal | author=K. G. Wilson | title=A Model Of Coupling Constant Renormalization | journal=Physical Review D | volume= 2 | issue=8 | pages= 1438–1472 | year=1970 | doi= 10.1103/PhysRevD.2.1438 | bibcode=1970PhRvD...2.1438W | osti=1444813 }}{{ cite journal | author=G. 't Hooft | title=Dimensional regularization and the renormalization group | journal=Nuclear Physics B | volume= 61 | pages= 455–468 | year=1973 | doi= 10.1016/0550-3213(73)90376-3| bibcode=1973NuPhB..61..455T | url=https://cds.cern.ch/record/880603 }}{{ cite journal | author1=K. G. Wilson|author2= J. B. Kogut | title=The Renormalization group and the epsilon expansion | journal=Physics Reports | volume= 12 |issue= 2 | pages= 75–199 | year=1974 | doi= 10.1016/0370-1573(74)90023-4|bibcode= 1974PhR....12...75W }} Certainly, if one were in possession of a good answer to the second question, the first question could also be answered. In the absence of a good answer to the second question, one may wonder why any specific choice of parameters and their cutoff dependence could secure cutoff independence of all observables $X(p,\backslash Lambda,\backslash mu)$ with finite scales $\backslash mu$.

The renormalization group concept associated with the first question above relies on the circumstance that some finite set $p(\backslash Lambda)$ yields the desired result,

:$X\_\{\backslash rm\; observable\}(\backslash mu)\; =\; \backslash lim\_\{\backslash Lambda\; \backslash rightarrow\; \backslash infty\}\; X\_\{\backslash rm\; theory\}[p(\backslash Lambda),\; \backslash Lambda,\; \backslash mu]\; \backslash ,\; .$

In this way of thinking, one can expect that in a theory with $n$ parameters a calculation of $n$ observables at some scale $\backslash mu$ is sufficient to fix all parameters as functions of $\backslash Lambda$. So, one may hope that there exists a collection of $n$ effective parameters at scale $\backslash mu$, corresponding to $n$ observables at scale $\backslash mu$, that are sufficient to parametrize the theory in such a way that predictions expressed in terms of these parameters are free from dependence on $\backslash Lambda$. Since the scale $\backslash mu$ is arbitrary, a whole family of such $n$-parameter sets labeled by $\backslash mu$ should exist, and every member of that family corresponds to the same physics. Moving from one such family to another by changing one value of $\backslash mu$ to another is described as action of the renormalization group. The word group is justified because the group axioms are satisfied: two such changes form another such change, one can invert a change, etc.

The question remains, however, why fixing the cutoff dependence of $n$ parameters $p$ on $\backslash Lambda$, using $n$ conditions that $n$ selected observables do not depend on $\backslash Lambda$, is good enough to make all observables in the physical range of $\backslash mu$ not depend on $\backslash Lambda$. In some theories such a miracle may happen but in others it may not. The ones where it happens are called renormalizable, because one can normalize the parameters properly to obtain cutoff independent results.

Typically, the set $p(\backslash Lambda)$ is established using perturbative calculations that are combined with models for description of nonperturbative effects. For example, perturbative QCD diagrams for quarks and gluons are combined with the parton models for description of binding of quarks and gluons into hadrons. The set of parameters $p(\backslash Lambda)$ includes cutoff dependent masses, charges and field normalization constants. The predictive power of a theory set up this way relies on the circumstance that the required set of parameters is relatively small. The regularization is designed order-by-order so that as many formal symmetries as possible of a local theory are preserved and employed in calculations, as in the dimensional regularization of Feynman diagrams. The claim that the set of parameters $p(\backslash Lambda)$ leads to finite, cutoff independent limits for all observables is qualified by the need to use some form of perturbation theory and inclusion of model assumptions concerning bound states.

The renormalization group concept associated with the second question above is conceived to explain how it may be so that the concept of renormalization group associated with the first question can make sense, instead of being at best a successful recipe to deal with divergences in perturbative calculations.{{ cite journal | author=P. A. M. Dirac | title=Quantum Electrodynamics without Dead Wood | journal=Physical Review | volume= 139 | issue=3B | page= B684 | year=1965 | doi= 10.1103/PhysRev.139.B684| bibcode=1965PhRv..139..684D }} Namely, to answer the second question, one designs a calculation (see below) that identifies the required set of parameters to define the theory, the starting point being some specific initial assumption, such as some local Lagrangian density which is a function of field variables and needs to be modified by including all the required parameters. Once the required set of parameters is known, one can establish a set of observables that are sufficient to define the cutoff dependence of the required set. The observables can have any finite scale $\backslash mu$, and one can use any scale $\backslash mu$ to define the parameters $p(\backslash Lambda)$, up to their finite parts that must be fitted to experiment, including features such as the observed symmetries.

Thus, not only the possibility that a renormalization group of the first type may exist can be understood, but also the alternative situations are found where the set of required cutoff dependent parameters does not have to be finite. Predictive power of latter theories results from known relationships among the required parameters and options to establish all the relevant ones.{{ cite journal | author1=R. J. Perry|author2= K. G. Wilson | title=Perturbative renormalizability with an infinite number of relevant and marginal operators | journal=Nuclear Physics B | volume= 403 |issue= 3 | pages= 587–601 | year=1993 | doi= 10.1016/0550-3213(93)90363-T|bibcode= 1993NuPhB.403..587P }}

The renormalization group concept of the second kind is associated with the nature of the mathematical computation used to discover the set of parameters $p$. In its essence, the calculation starts with some specific form of a theory with cutoff $\backslash Lambda$ and derives a corresponding theory with a smaller cutoff, in the sense of more restrictive, say $\backslash Lambda/2$. After re-parameterization using the cutoff as a unit, one obtains a new theory of similar type but with new terms. This means that the starting theory with cutoff $\backslash Lambda$ should also contain such new terms for its form to be consistent with the presence of a cutoff. Eventually, one can find a set of terms that reproduces itself up to changes in the coefficients of the required terms. These coefficients evolve with the number of steps one makes, in each and every step reducing the cutoff by factor of two and rescaling variables. One could use other factors than two, but two is convenient.

In summary, one obtains a trajectory of a point in a space of dimension equal to the number of required parameters and motion along the trajectory is described by transformations that form new kind of a group. Different initial points might lead to different trajectories, but if the steps are self-similar and reduce to a multiple action of one and the same transformation, say $T$, one may describe what happens in terms of the features of $T$, called the renormalization group transformation. The transformation $T$ may transform points in the parameter space making some of the parameters decrease, some grow, and some stay unchanged. It may have fixed points, limit cycles, or even lead to chaotic motion.

Suppose that $T$ has a fixed point. If one starts the procedure at this point, an infinitely long sequence of reductions of the cutoff by factors of two changes nothing in the structure of the theory, except the scale of its cutoff. This means that the initial cutoff can be arbitrarily large. Such a theory may possess the symmetries of special relativity, since there is no price to pay for extending the cutoff as required when one wishes to make the Lorentz transformation that yields momenta which exceed the cutoff.

Both concepts of the renormalization group can be considered in quantum theories constructed using the front form of dynamics. The first concept allows one to play with a small set of parameters and seek consistency, which is a useful strategy in perturbation theory if one knows from other approaches what to expect. In particular, one may study new perturbative features that appear in the front form of dynamics, since it differs from the instant form. The main difference is that the front variables $x^-$ (or $p^+$) are considerably different from the transverse variables $x^\backslash perp$ (or $p^\backslash perp$), so that there is no simple rotational symmetry among them.

One can also study sufficiently simplified models for which computers can be used to carry out calculations and see if a procedure suggested by perturbation theory may work beyond it. The second concept allows one to address the issue of defining a relativistic theory ab initio without limiting the definition to perturbative expansions. This option is particularly relevant to the issue of describing bound states in QCD. However, to address this issue one needs to overcome certain difficulties that the renormalization group procedures based on the idea of reduction of cutoffs are not capable of easily resolving. To avoid the difficulties, one can employ the similarity renormalization group procedure. Both the difficulties and similarity are explained in the next section.

A glimpse of the difficulties of the procedure of reducing a cutoff $\backslash Lambda$ to cutoff $\backslash Lambda/2$ in the front form of Hamiltonian dynamics of strong interactions can be gained by considering the eigenvalue problem for the Hamiltonian $H$,

:$H\; \backslash psi\; =\; E\; \backslash psi\; ,$

where $H\; =\; H\_0\; +\; H\_I$, $H\_0$ has a known spectrum and $H\_I$ describes the interactions. Let us assume that the eigenstate $\backslash psi$ can be written as a superposition of eigenstates of $H\_0$ and let us introduce two projection operators, $P$ and $Q$, such that $P$ projects on eigenstates of $H\_0$ with eigenvalues smaller than $\backslash Lambda/2$ and $Q$ projects on eigenstates of $H\_0$ with eigenvalues between $\backslash Lambda/2$ and $\backslash Lambda$. The result of projecting the eigenvalue problem for $H$ using $P$ and $Q$ is a set of two coupled equations

:$H\_0\; Q\; \backslash psi\; +\; Q\; H\_I\; Q\; \backslash psi\; +\; Q\; H\_I\; P\; \backslash psi\; =\; E\; Q\; \backslash psi\; \backslash ,\; ,$

:$H\_0\; P\; \backslash psi\; +\; P\; H\_I\; Q\; \backslash psi\; +\; P\; H\_I\; P\; \backslash psi\; =\; E\; P\; \backslash psi\; \backslash ,\; .$

The first equation can be used to evaluate $Q\; \backslash psi$ in terms of $P\; \backslash psi$,

:$Q\; \backslash psi\; =\; \backslash frac\{1\}\{E\; -\; H\_0\; -\; Q\; H\_I\; Q\; \}\; \backslash \; Q\; H\_I\; P\; \backslash psi\; \backslash ,\; .$

This expression allows one to write an equation for $P\; \backslash psi$ in the form

:$H\_\{\backslash rm\; eff\}\; P\; \backslash psi\; =\; E\; P\backslash psi\; \backslash ,\; ,$

where

:$H\_\{\backslash rm\; eff\}\; =\; H\_0\; +\; P\; H\_I\; P\; +\; P\; H\_I\; Q\; \backslash frac\{\; 1\}\{\; E\; -\; H\_0\; -\; Q\; H\_I\; Q\; \}\; Q\; H\_I\; P\; .$

The equation for $P\backslash psi$ appears to resemble an eigenvalue problem for $H\_\{\backslash rm\; eff\}$. It is valid in a theory with cutoff $\backslash Lambda/2$, but its effective Hamiltonian $H\_\{\backslash rm\; eff\}$ depends on the unknown eigenvalue $E$. However, if $\backslash Lambda/2$ is much greater than $E$ of interest, one can neglect $E$ in comparison to $QH\_0Q$ provided that $QH\_IQ$ is small in comparison to $Q\; H\_0\; Q$.

In QCD, which is asymptotically free, one indeed has $H\_0$ as the dominant term in the energy denominator in $H\_\{\backslash rm\; eff\}$ for small eigenvalues $E$. In practice, this happens for cutoffs $\backslash Lambda$ so much larger than the smallest eigenvalues $E$ of physical interest that the corresponding eigenvalue problems are too complex for solving them with required precision. Namely, there are still too many degrees of freedom. One needs to reduce cutoffs considerably further. This issue appears in all approaches to the bound state problem in QCD, not only in the front form of the dynamics.

Even if interactions are sufficiently small, one faces an additional difficulty with eliminating $Q$-states. Namely, for small interactions one can eliminate the eigenvalue $E$ from a proper effective Hamiltonian in $P$-subspace in favor of eigenvalues of $H\_0$. Consequently, the denominators analogous to the one that appears above in $H\_\{\backslash rm\; eff\}$ only contain differences of eigenvalues of $H\_0$, one above $\backslash Lambda/2$ and one below. Unfortunately, such differences can become arbitrarily small near the cutoff $\backslash Lambda/2$, and they generate strong interactions in the effective theory due to the coupling between the states just below and just above the cutoff $\backslash Lambda/2$. This is particularly bothersome when the eigenstates of $H\_0$ near the cutoff are highly degenerate and splitting of the bound state problem into parts below and above the cutoff cannot be accomplished through any simple expansion in powers of the coupling constant.

In any case, when one reduces the cutoff $\backslash Lambda$ to $\backslash Lambda/2$, and then $\backslash Lambda/2$ to $\backslash Lambda/4$ and so on, the strength of interaction in QCD Hamiltonians increases and, especially if the interaction is attractive, $Q\; H\_I\; Q$ can cancel $H\_0$ and $E$ cannot be ignored no matter how small it is in comparison to the reduced cutoff. In particular, this difficulty concerns bound states, where interactions must prevent free relative motion of constituents from dominating the scene and a spatially compact systems have to be formed. So far, it appears not possible to precisely eliminate the eigenvalue $E$ from the effective dynamics obtained by projecting on sufficiently low energy eigenstates of $H\_0$ to facilitate reliable calculations.

Fortunately, one can use instead a change of basis.{{ cite journal | author1=S. D. Glazek|author2= K. G. Wilson | title=Renormalization of Hamiltonians | journal=Physical Review D | volume= 48 |issue= 12 | pages= 5863–5872 | year=1993 | doi= 10.1103/PhysRevD.48.5863 |pmid= 10016252 | bibcode=1993PhRvD..48.5863G| arxiv=hep-th/9706149 |s2cid= 39086918 }} Namely, it is possible to define a procedure in which the basis states are rotated in such a way that the matrix elements of $H\_I$ vanish between basis states that according to $H\_0$ differ in energy by more than a running cutoff, say $\backslash lambda$. The running cutoff is called the energy bandwidth. The name comes from the band-diagonal form of the Hamiltonian matrix in the new basis ordered in energy using $H\_0$. Different values of the running cutoff $\backslash lambda$ correspond to using differently rotated basis states. The rotation is designed not to depend at all on the eigenvalues $E$ one wants to compute.

As a result, one obtains in the rotated basis an effective Hamiltonian matrix eigenvalue problem in which the dependence on cutoff $\backslash Lambda$ may manifest itself only in the explicit dependence of matrix elements of the new $H\_\{\backslash rm\; eff\}$. The two features of similarity that (1) the $\backslash Lambda$-dependence becomes explicit before one tackles the problem of solving the eigenvalue problem for $H\_\{\backslash rm\; eff\}$ and (2) the effective Hamiltonian with small energy bandwidth may not depend on the eigenvalues one tries to find, allow one to discover in advance the required counterterms to the diverging cutoff dependence. A complete set of counterterms defines the set of parameters required for defining the theory which has a finite energy bandwidth $\backslash lambda$ and no cutoff dependence in the band. In the course of discovering the counterterms and corresponding parameters, one keeps changing the initial Hamiltonian. Eventually, the complete Hamiltonian may have cutoff independent eigenvalues, including bound states.

In the case of the front-form Hamiltonian for QCD, a perturbative version of the similarity renormalization group procedure is outlined by Wilson et al.{{ cite journal | author1=K. G. Wilson|author2= T. S. Walhout|author3= A. Harindranath|author4= W.-M. Zhang|author5= R. J. Perry|author6= S. D. Glazek | title=Nonperturbative QCD: A Weak coupling treatment on the light front | journal=Physical Review D | volume= 49 |issue= 12| pages= 6720–6766 | year=1994 | doi= 10.1103/PhysRevD.49.6720 |pmid= 10016996| bibcode=1994PhRvD..49.6720W| arxiv= hep-th/9401153|s2cid= 119422380}} Further discussion of computational methods stemming from the similarity renormalization group concept is provided in the next section.

The similarity renormalization group procedure, discussed in #Similarity transformations, can be applied to the problem of describing bound states of quarks and gluons using QCD according to the general computational scheme outlined by Wilson et al. and illustrated in a numerically soluble model by Glazek and Wilson.{{ cite journal | author1=S. D. Glazek|author2= K. G. Wilson | journal=Physical Review D | volume= 57 |issue= 6 | pages= 3558–3566 | year=1998 | doi= 10.1103/PhysRevD.57.3558 | title=Asymptotic freedom and bound states in Hamiltonian dynamics | bibcode=1998PhRvD..57.3558G| arxiv=hep-th/9707028 |s2cid= 16805417 }} Since these works were completed, the method has been applied to various physical systems using a weak-coupling expansion. More recently, similarity has evolved into a computational tool called the renormalization group procedure for effective particles, or RGPEP. In principle, the RGPEP is now defined without a need to refer to some perturbative expansion. The most recent explanation of the RGPEP is given by Glazek in terms of an elementary and exactly solvable model for relativistic fermions that interact through a mass mixing term of arbitrary strength in their Hamiltonian.{{ cite journal | author=S. D. Glazek | title=Perturbative Formulae for Relativistic Interactions of Effective Particles | journal=Acta Physica Polonica B | volume= 43 | issue=9 | page= 1843 | year=2012 | doi= 10.5506/APhysPolB.43.1843| doi-access=free }}{{ cite journal | author=S. D. Glazek | title=Fermion mass mixing and vacuum triviality in the renormalization group procedure for effective particles | journal=Physical Review D | volume= 87 | issue=12 | page= 125032 | year=2013 | doi= 10.1103/PhysRevD.87.125032 | bibcode=2013PhRvD..87l5032G| arxiv=1305.3702 | s2cid=119222650 }}

The effective particles can be seen as resulting from a dynamical transformation akin to the Melosh transformation from current to constituent quarks.{{ cite journal | author=H. J. Melosh | title=Quarks: Currents and constituents | journal=Physical Review D | volume= 9 | issue=4 | pages= 1095–1112 | year=1974 | doi= 10.1103/PhysRevD.9.1095 | bibcode=1974PhRvD...9.1095M| url=https://thesis.library.caltech.edu/4807/3/melosh-hj_1973.pdf }} Namely, the RGPEP transformation changes the bare quanta in a canonical theory to the effective quanta in an equivalent effective theory with a Hamiltonian that has the energy bandwidth $\backslash lambda$; see #Similarity transformations and references therein for an explanation of the band. The transformations that change $\backslash lambda$ form a group.

The effective particles are introduced through a transformation

:$\backslash psi\_s\; =\; U\_s\; \backslash ,\; \backslash psi\_0\; \backslash ,\; U\_s^\backslash dagger\; \backslash ,\; ,$

where $\backslash psi\_s$ is a quantum field operator built from creation and annihilation operators for effective particles of size $s\; \backslash sim\; 1/\backslash lambda$ and $\backslash psi\_0$ is the original quantum field operator built from creation and annihilation operators for point-like bare quanta of a canonical theory. In great brevity, a canonical Hamiltonian density is built from fields $\backslash psi\_0$ and the effective Hamiltonian at scale $s$ is built from fields $\backslash psi\_s$, but without actually changing the Hamiltonian. Thus,

:$H\_s(\backslash psi\_s)\; =\; H\_0(\backslash psi\_0)\; \backslash ,\; ,$

which means that the same dynamics is expressed in terms of different operators for different values of $s$. The coefficients $c\_s$ in the expansion of a Hamiltonian in powers of the field operators $\backslash psi\_s$ depend on $s$ and the field operators depend on $s$, but the Hamiltonian is not changing with $s$. The RGPEP provides an equation for the coefficients $c\_s$ as functions of $s$.

In principle, if one had solved the RGPEP equation for the front form Hamiltonian of QCD exactly, the eigenvalue problem could be written using effective quarks and gluons corresponding to any $s$. In particular, for $s$ very small, the eigenvalue problem would involve very large numbers of virtual constituents capable of interacting with large momentum transfers up to about the bandwidth $\backslash lambda\; \backslash sim\; 1/s$. In contrast, the same eigenvalue problem written in terms of quanta corresponding to a large $s$, comparable with the size of hadrons, is hoped to take the form of a simple equation that resembles the constituent quark models. To demonstrate mathematically that this is precisely what happens in the RGPEP in QCD is a serious challenge.

The Bethe–Salpeter amplitude, which satisfies the Bethe–Salpeter equation{{ cite journal | author1=E.E. Salpeter|author2= H.A. Bethe | title=A Relativistic Equation for Bound-State Problems | journal=Physical Review | volume= 84 |issue= 6 | pages= 1232–1242 | year=1951 | doi= 10.1103/PhysRev.84.1232|bibcode= 1951PhRv...84.1232S }}{{ cite journal | author1=J. H. O. Sales|author2= T. Frederico|author3= B. V. Carlson|author4= P. U. Sauer | title=Light-front Bethe-Salpeter equation | journal=Physical Review C | volume= 61 |issue= 4| page= 044003 | year=2000 | doi= 10.1103/PhysRevC.61.044003|bibcode= 2000PhRvC..61d4003S| arxiv= nucl-th/9909029|s2cid= 17341344}}{{ cite journal | author1=J. H. Sales|author2= T. Frederico|author3= B. V. Carlson|author4= P. U. Sauer | title=Renormalization of the ladder light-front Bethe-Salpeter equation in the Yukawa model | volume= 63 |issue= 6| page= 064003 | year=2001 | doi= 10.1103/PhysRevC.63.064003 |journal=Physical Review C|bibcode= 2001PhRvC..63f4003S|hdl= 11449/66525|hdl-access= free}} (see the reviews by Nakanishi{{ cite journal | author=N. Nakanishi | title=A General survey of the theory of the Bethe-Salpeter equation | journal=Progress of Theoretical Physics Supplement | volume= 43 | pages= 1–81 | year=1969 | doi= 10.1143/PTPS.43.1| bibcode=1969PThPS..43....1N | doi-access=free }}{{ cite journal | author=N. Nakanishi | title=Review of the Wick-cutkosky Model | journal=Progress of Theoretical Physics Supplement | volume= 95 | pages= 1–24 | year=1988 | doi= 10.1143/PTPS.95.1| bibcode=1988PThPS..95....1N | doi-access=free }} ), when projected on the light-front plane, results in the light-front wave function. The meaning of the ``light-front projection" is the following. In the coordinate space, the Bethe–Salpeter amplitude is a function of two four-dimensional coordinates $x\_\{1,2\}=(ct\_\{1,2\},\backslash vec\{x\}\_\{1,2\})$, namely: $\backslash Phi=\backslash Phi(x\_1,x\_2;p)$, where $p$ is the total four-momentum of the system. In momentum space, it is given by the Fourier transform:

:$\backslash Phi(k\_1,k\_2;p)=\backslash int\; d^4x\_1\; d^4x\_2\backslash Phi(x\_1,x\_2;p)\backslash exp(ik\_1\; x\_1+ik\_2\; x\_2)$

(the momentum space Bethe–Salpeter amplitude $\backslash Phi(k\_1,k\_2;p)$ defined in this way includes in itself the delta-function responsible for the momenta conservation $k\_1+k\_2=p$). The light-front projection means that the arguments $x\_1,x\_2$ are on the light-front plane, i.e., they are constrained by the condition (in the covariant formulation): $\backslash omega\backslash cdot\; x\_1=\backslash omega\backslash cdot\; x\_2=0$. This is achieved by inserting in the Fourier transform the corresponding delta functions $\backslash delta(\backslash omega\backslash cdot\; x\_\{1,2\})$:

:$\backslash psi\_\{LF\}\backslash propto\backslash int\; d^4x\_1\; d^4x\_2\backslash delta(\backslash omega\backslash cdot\; x\_1)\; \backslash delta(\backslash omega\backslash cdot\; x\_2)\backslash Phi(x\_1,x\_2;p)\backslash exp(ik\_1\; x\_1+ik\_2\; x\_2).$

In this way, we can find the light-front wave function $\backslash psi\_\{LF\}$. Applying this formula to the Bethe–Salpeter amplitude with a given total angular momentum, one reproduces the angular momentum structure of the light-front wave function described in Light front quantization#Angular momentum. In particular, projecting the Bethe–Salpeter amplitude corresponding to a system of two spinless particles with the angular momentum $l$, one reproduces the light-front wave function

:$\backslash psi\_\{lm\}(\backslash vec\{k\},\backslash hat\{n\})=f\_1(k,\; \backslash vec\{k\}\backslash cdot\backslash hat\{n\})Y\_\{lm\}(\backslash hat\{k\})+\; f\_2(k,\backslash vec\{k\}\backslash cdot\backslash hat\{n\})Y\_\{lm\}(\backslash hat\{n\}),$

given in Light front quantization#Angular momentum.

The Bethe–Salpeter amplitude includes the propagators of the external particles, and, therefore, it is singular. It can be represented in the form of the Nakanishi integral{{ cite book | first=N. | last=Nakanishi | title=Graph Theory and Feynman Integrals | publisher=Gordon and Breach | location= New York | date=1971 }} through a non-singular function $g(\backslash gamma,z)$:

{{NumBlk|:|$\backslash Phi(k,p)=\backslash frac\{1\}\{\backslash sqrt\{4\backslash pi\}\}\backslash int\_\{-1\}^1dz\text{'}\backslash int\_0^\{\backslash infty\}d\backslash gamma\text{'}\; \backslash frac\{g(\backslash gamma\text{'},z\text{'})\}\{\backslash left[\backslash gamma\text{'}+m^2\; -\backslash frac\{1\}\{4\}M^2-k^2-p\backslash cdot\; k\backslash ;\; z\text{'}-i\backslash epsilon\backslash right]^3\},$|{{EquationRef|1}}}}

where $k=(k\_1-k\_2)/2$ is the relative four-momentum. The Nakanishi weight function $g(\backslash gamma,z)$ is found from an equation and has the properties: $g(\backslash gamma,z=\backslash pm\; 1)=0$, $g(\backslash gamma\backslash to\backslash infty,z)\backslash to\; 0$. Projecting the Bethe–Salpeter amplitude ({{EquationNote|1}}) on the light-front plane, we get the following useful representation for the light-front wave function (see the review by Carbonell and Karmanov{{ cite journal | author1=J. Carbonell|author2= V.A. Karmanov | title=Solutions of the Bethe-Salpeter equation in Minkowski space and applications to electromagnetic form factors | journal=Few Body Systems | volume= 49 | pages= 205–222 | year=2011 |issue= 1–4 | doi= 10.1007/s00601-010-0133-5|bibcode= 2011FBS....49..205C | arxiv= 1012.0246 |s2cid= 118922610 }}):

:$\backslash psi\_\{LF\}(k\_\backslash perp,x)=\backslash frac\{1\}\{\backslash sqrt\{4\backslash pi\}\}\backslash int\_0^\{\backslash infty\}\backslash frac\{x(1-x)g(\backslash gamma\text{'},1-2x)d\backslash gamma\text{'}\}\; \{\backslash Bigl[\backslash gamma\text{'}+k\_\backslash perp^2\; +m^2-x(1-x)M^2\backslash Bigr]^2\}.$

It turns out that the masses of a two-body system, found from the Bethe–Salpeter equation for $\backslash Phi(k,p)$ and from the light-front equation for $\backslash psi\_\{LF\}(k\_\backslash perp,x)$ with the kernel corresponding to the same physical content, say, one-boson exchange (which, however, in the both approaches have very different analytical forms) are very close to each other. The same is true for the electromagnetic form factors{{ cite journal | author1=V. A. Karmanov|author2= P. Maris | title=Manifestation of three-body forces in three-body Bethe-Salpeter and light-front equations | journal=Few Body Systems | volume= 46 |issue= 2 | pages= 95–113 | year=2009 | doi= 10.1007/s00601-009-0054-3|bibcode= 2009FBS....46...95K | arxiv= 0811.1100 |s2cid= 15070868 }} This undoubtedly proves the existence of three-body forces, though the contribution of relativistic origin does not exhaust, of course, all the contributions. The same relativistic dynamics should generate four-body forces, etc. Since in nuclei the small binding energies (relative to the nucleon mass) result from cancellations between the kinetic and potentials energies (which are comparable with nucleon mass, and, hence relativistic), the relativistic effects in nuclei are noticeable. Therefore, many-body forces should be taken into account for fine tuning to experimental data.

One of the advantages of light-front quantization is that the empty state, the so-called perturbative vacuum, is the physical vacuum.{{ cite journal | author1=Y. Nambu|author2= G. Jona-Lasinio | title=Dynamical model of elementary particles based on an analogy with auperconductivity | journal=Physical Review | volume= 122 | pages= 345–358 | year=1961 |issue= 1 | doi= 10.1103/PhysRev.122.345 | bibcode=1961PhRv..122..345N| doi-access=free }}{{ cite journal | author1=M. Gell-Mann|author2= R. J. Oakes|author3= B. Renner | title=Behavior of current divergences under SU(3) x SU(3) | journal=Physical Review | volume= 175 |issue= 5| pages= 2195–2199 | year=1968 | doi= 10.1103/PhysRev.175.2195 | bibcode=1968PhRv..175.2195G|url= https://authors.library.caltech.edu/3634/1/GELpr68.pdf}}{{ cite journal | author1=G. 't Hooft|author2= M. Veltman | title=Regularization and renormalization of gauge fields | journal=Nuclear Physics B | volume= 44 | pages= 189–213 | year=1972 |issue= 1 | doi= 10.1016/0550-3213(72)90279-9 | bibcode=1972NuPhB..44..189T|hdl= 1874/4845 |url= https://repositorio.unal.edu.co/handle/unal/81144 | hdl-access= free }}{{ cite journal | author1=M. A. Shifman|author2= A.I. Vainshtein|author3= V. I. Zakharov | title=QCD and Resonance Physics: Applications | journal=Nuclear Physics B | volume= 147 |issue= 5| pages= 448–518 | year=1979 | doi= 10.1016/0550-3213(79)90023-3 | bibcode=1979NuPhB.147..448S}}{{ cite journal | author=R. P. Feynman | title=The Qualitative Behavior of Yang-Mills Theory in (2+1)-Dimensions | journal=Nuclear Physics B | volume= 188 | issue=3 | pages= 479–512 | year=1981 | doi= 10.1016/0550-3213(81)90005-5 | bibcode=1981NuPhB.188..479F}}{{ cite journal | author=E. Witten | title=Dynamical Breaking of Supersymmetry | journal=Nuclear Physics B | volume= 188 | issue=3 | pages= 513–554 | year=1981 | doi= 10.1016/0550-3213(81)90006-7 | bibcode=1981NuPhB.188..513W}}{{ cite journal | author1=J. Gasser|author2= H. Leutwyler | title=Chiral Perturbation Theory to One Loop | journal=Annals of Physics | volume= 158 | pages= 142–210 | year=1984 |issue= 1 | doi= 10.1016/0003-4916(84)90242-2 | bibcode=1984AnPhy.158..142G|url= https://cds.cern.ch/record/147992 | url-access= subscription }}{{ cite journal | author=S. D. Glazek | title=Light Front QCD in the Vacuum Background | journal=Physical Review D | volume= 38 | issue=10 | pages= 3277–3286 | year=1988 | doi= 10.1103/PhysRevD.38.3277 | pmid=9959077 | bibcode=1988PhRvD..38.3277G}}{{ cite journal | author1=P. Maris|author2= C. D. Roberts|author3= P. C. Tandy | title=Pion mass and decay constant | journal=Physics Letters B | volume= 420 |issue= 3–4| pages= 267–273 | year=1998 | doi= 10.1016/S0370-2693(97)01535-9 | bibcode=1998PhLB..420..267M| arxiv= nucl-th/9707003 |s2cid= 16778465}}{{ cite journal | author1= S. J. Brodsky|author2= C. D. Roberts|author3= R. Shrock|author4= P. C. Tandy | title= Confinement contains condensates | journal=Physical Review C | volume= 85 |issue= 6| page= 065202 | year=2012 | doi= 10.1103/PhysRevC.85.065202 |arxiv= 1202.2376| bibcode=2012PhRvC..85f5202B| doi-access= free }}{{ cite journal | author1= A. Casher|author2= L. Susskind | title=Chiral magnetism (or magnetohadrochironics) | journal=Physical Review D | volume= 9 |issue= 2 | pages= 436–460 | year=1974 | doi= 10.1103/PhysRevD.9.436 | bibcode=1974PhRvD...9..436C}} The massive states of a theory can then be built on this lowest state without having any contributions from vacuum structure, and the wave functions for these massive states do not contain vacuum contributions. This occurs because each $p\_i^+$ is positive, and the interactions of the theory cannot produce particles from the zero-momentum vacuum without violating momentum conservation. There is no need to normal-order the light-front vacuum.

However, certain aspects of some theories are associated with vacuum structure. For example, the Higgs mechanism of the Standard Model relies on spontaneous symmetry breaking in the vacuum of the theory.{{ cite journal | author1=C. M. Bender|author2= S. S. Pinsky|author3= B. van de Sande | title=Spontaneous symmetry breaking of $\backslash phi^4$ in (1+1)-dimensions in light front field theory | journal=Physical Review D | volume= 48 |issue= 2| pages= 816–821 | year=1993 | doi= 10.1103/PhysRevD.48.816 |pmid= 10016310| bibcode=1993PhRvD..48..816B| arxiv= hep-th/9212009 |s2cid= 14265514}}{{ cite journal | author1=S. S. Pinsky|author2= B. van de Sande | title=Spontaneous symmetry breaking of (1+1)-dimensional $\backslash phi^4$ theory in light front field theory. 2 | journal=Physical Review D | volume= 49 |issue= 4 | pages= 2001–2013 | year=1994 | doi= 10.1103/PhysRevD.49.2001 |pmid= 10017185 | bibcode=1994PhRvD..49.2001P| arxiv=hep-ph/9310330 |s2cid= 17165941 }}{{ cite journal | author1=S. S. Pinsky|author2= B. van de Sande|author3= J.R. Hiller | title=Spontaneous symmetry breaking of (1+1)-dimensional $\backslash phi^4$ theory in light front field theory. 3 | journal=Physical Review D | volume= 51 |issue= 2| pages= 726–733 | year=1995 | doi= 10.1103/PhysRevD.51.726 |pmid= 10018525| bibcode=1995PhRvD..51..726P| arxiv= hep-th/9409019 |s2cid= 15291034}}{{ cite journal | author1= J. S. Rozowsky|author2= C. B. Thorn | title= Spontaneous symmetry breaking at infinite momentum without P+ zero modes | journal=Physical Review Letters | volume= 85 |issue= 8 | pages= 1614–1617 | year=2000 | doi= 10.1103/PhysRevLett.85.1614 | bibcode=2000PhRvL..85.1614R | pmid=10970571| arxiv= hep-th/0003301 |s2cid= 17968437 }}{{ cite journal | author1= D. Chakrabarti|author2= A. Harindranath|author3= L. Martinovic|author4= G. B. Pivovarov|author5= J. P. Vary | title= Ab initio results for the broken phase of scalar light front field theory | journal=Physics Letters B | volume= 617 | pages= 92–98 | year=2005 |issue= 1–2| doi= 10.1016/j.physletb.2005.05.012 |arxiv= hep-th/0310290| bibcode=2005PhLB..617...92C| doi-access= free }}{{ cite journal | author1= V. T. Kim|author2= G. B. Pivovarov|author3= J. P. Vary | title= Phase transition in light front $\backslash phi^4\_\{1+1\}$ | journal=Physical Review D | volume= 69 |issue= 8| page= 085008 | year=2004 | doi= 10.1103/PhysRevD.69.085008 | bibcode=2004PhRvD..69h5008K| arxiv= hep-th/0310216 |s2cid= 119524638}} The usual Higgs vacuum expectation value in the instant form is replaced by $k^+=0$ zero mode analogous to a constant Stark field when one quantizes the Standard model using the front form.{{ cite journal | author1= P. P. Srivastava|author2= S. J. Brodsky | title= A Unitary and renormalizable theory of the standard model in ghost free light cone gauge | journal=Physical Review D | volume= 66 |issue= 4 | page= 045019 | year=2002 | doi= 10.1103/PhysRevD.66.045019|arxiv= hep-ph/0202141 |bibcode= 2002PhRvD..66d5019S }} Chiral symmetry breaking of quantum chromodynamics is often associated in the instant form with quark and gluon condensates in the QCD vacuum. However, these effects become properties of the hadron wave functions themselves using the front form.{{ cite journal | author1= S. J. Brodsky|author2= R. Shrock | title= Condensates in Quantum Chromodynamics and the Cosmological Constant | journal=Proceedings of the National Academy of Sciences | volume= 108 | pages= 45–50 | year=2011 |issue= 1 | doi= 10.1073/pnas.1010113107|arxiv= 0803.2554 |bibcode= 2011PNAS..108...45B | doi-access= free }}{{ cite journal | author1= S. J. Brodsky|author2= C. D. Roberts|author3= R. Shrock|author4= P. C. Tandy | title= Essence of the vacuum quark condensate | journal=Physical Review C | volume= 82 |issue= 2| page= 022201 | year=2010 | doi= 10.1103/PhysRevC.82.022201|bibcode= 2010PhRvC..82b2201B| arxiv= 1005.4610}} This also eliminates the many orders of magnitude conflict between the measured cosmological constant and quantum field theory.

Some aspects of vacuum structure in light-front quantization can be analyzed by studying properties of massive states. In particular, by studying the appearance of degeneracies among the lowest massive states, one can determine the critical coupling strength associated with spontaneous symmetry breaking. One can also use a limiting process, where the analysis begins in equal-time quantization but arrives in light-front coordinates as the limit of some chosen parameter.{{ cite journal | author=K. Hornbostel | title=Nontrivial vacua from equal time to the light cone | journal=Physical Review D | volume= 45 | issue=10 | pages= 3781–3801 | year=1992 | doi= 10.1103/PhysRevD.45.3781| pmid=10014271 | bibcode=1992PhRvD..45.3781H }}{{ cite journal | author1=C.-R. Ji|author2= A. Suzuki | title=Interpolating scattering amplitudes between the instant form and the front form of relativistic dynamics | journal=Physical Review D | volume= 87|issue= 6 |page= 065015 | year=2013 | doi= 10.1103/PhysRevD.87.065015|arxiv= 1212.2265 |bibcode= 2013PhRvD..87f5015J | doi-access= free }} A much more direct approach is to include modes of zero longitudinal momentum (zero modes) in a calculation of a nontrivial light-front vacuum built from these modes; the Hamiltonian then contains effective interactions that determine the vacuum structure and provide for zero-mode exchange interactions between constituents of massive states.

• Light front quantization
• Light-front quantization applications
• Quantum field theories
• Quantum chromodynamics
• Quantum electrodynamics
• Light-front holography

{{reflist}}

• [http://www.ilcacinc.org ILCAC, Inc.], the International Light-Cone Advisory Committee.
• [http://www.saha.ac.in/theory/a.harindranath/light/light.html Publications on light-front dynamics], maintained by A. Harindranath.

Category:Quantum chromodynamics