{'fraction_non_alphanumeric': 0.040769185785852126, 'fraction_numerical': 0.007904329350305063, 'mean_word_length': 4.685321660109714, 'pattern_counts': {'":': 0, '<': 15, '<?xml version=': 0, '>': 9, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 19, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A Euclidean Signature Semi-Classical Program

In this article we discuss our ongoing program to extend the scope of certain, well-developed microlocal methods for the asymptotic solution of Schr\"{o}dinger's equation (for suitable `nonlinear oscillatory' quantum mechanical systems) to the treatment of several physically significant, interacting quantum field theories. Our main focus is on applying these `Euclidean-signature semi-classical' methods to self-interacting (real) scalar fields of renormalizable type in 2, 3 and 4 spacetime dimensions and to Yang-Mills fields in 3 and 4 spacetime dimensions. A central argument in favor of our program is that the asymptotic methods for Schr\"{o}dinger operators developed in the microlocal literature are far superior, for the quantum mechanical systems to which they naturally apply, to the conventional WKB methods of the physics literature and that these methods can be modified, by techniques drawn from the calculus of variations and the analysis of elliptic boundary value problems, to apply to certain (bosonic) quantum field theories. Unlike conventional (Rayleigh/ Schr\"{o}dinger) perturbation theory these methods allow one to avoid the artificial decomposition of an interacting system into an approximating `unperturbed' system and its perturbation and instead to keep the nonlinearities (and, if present gauge invariances) of an interacting system intact at every level of the analysis.

which they naturally apply, to the conventional WKB methods of the physics literature and that these methods can be modified, by techniques drawn from the calculus of variations and the analysis of elliptic boundary value problems, to apply to certain (bosonic) quantum field theories. Unlike conventional (Rayleigh/ Schrödinger) perturbation theory these methods allow one to avoid the artificial decomposition of an interacting system into an approximating 'unperturbed' system and its perturbation and instead to keep the nonlinearities (and, if present gauge invariances) of an interacting system intact at every level of the analysis.

I. INTRODUCTION
Perturbation theory, in its various guises, is one of the cornerstone calculational techniques of modern quantum field theory. Its application to quantum electrodynamics, through its Feynman diagrammatic and path integral incarnations among other approaches, was one of the principal success stories of twentieth century physics and gave rise to numerous, spectacularly confirmed experimental predictions. The quantitative success of this particular application of perturbation methods hinged, in large measure, on the feasibility of decomposing the relevant quantum system (of coupled Dirac and Maxwell fields) into a relativistically invariant and explicitly solvable 'unperturbed' system (the corresponding free fields) and a complementary, nonlinear interaction that could realistically be treated as a 'small' perturbation thereto.
When one turns however to the quantization of non-abelian gauge systems (i.e., to Yang-Mills fields either in isolation or coupled to other quantized objects), some of the advantages of the perturbative approach become seriously compromised. First of all, in order to have an explicitly solvable unperturbed system with which to initiate the analysis, one must dramatically modify the Yang-Mills system itself by discarding the very interaction terms that distinguish it from (several copies of) the Maxwell field. In doing so one necessarily disturbs the non-abelian character of the corresponding gauge group, 'truncating' this essential element to the abelian gauge group of the approximating, Maxwell system. Naturally the aim of the associated perturbation theory, at least in part, is to restore these essential features.
But it can only realistically hope to accomplish this through formal, asymptotic expansions in the corresponding Yang-Mills 'coupling constant'. Elegant techniques for doing this are at the forefront of modern research in quantum field theory of which, for example, the discovery of asymptotic freedom represented a crucial advance.
The Feynman path integral approach to field quantization is not intrinsically perturbative in nature (though it is often used in practice to derive elegant, relativistically invariant, perturbative expansions) and can, in its most mathematically coherent, Euclidean-signature formulation be used as a basis for rigorous attacks on quantization problems. Making mathematical sense of the needed (Euclidean-signature) functional integration measures in sufficient generality to handle Yang-Mills theories is however a daunting task, especially in the most physically interesting case of 4 spacetime dimensions wherein even more elementary, interacting quantized fields have continued to defy rigorous construction. Even in ordinary quantum mechanics there are known examples, which gave Feynman himself pause, of explicitly solvable quantum systems that forcefully resisted solution by purely path integral methods alone. Sometimes a direct, 'old-fashioned' attack on Schrödinger's equation is still the most effective approach.
Of course the expression of such misgivings about the intrinsic limitations of perturbation theory and even path integral methods must have a hollow ring to it to experts in the field who must continually deal with the inherent difficulties of quantizing interacting systemsall the more so when this expression is coming from outsiders to this field who may have no concrete alternative to put forward. The present authors however have in fact developed an alternative approximation method for solving the Schrödinger equations that arise in certain interacting quantum field theories and are confident that their method will shed some nonperturbative light on the nature of the associated quantum states. In particular our method does not require that one artificially decompose an interacting quantum system into 'free' and 'interacting' components and thus, in the case of gauge theories, does not require that one, equally artificially, truncate the corresponding gauge group to the abelian gauge group of an associated 'free', approximating system. While we do not expect our approach to have the calculational power of conventional perturbation theory we do anticipate that it will allow an attack on some fundamental questions that more conventional techniques seem inadequate to handle.
To explain this approach let us backtrack slightly and mention that microlocal analysts have long since developed elegant, semi-classical approximation methods for the analysis of certain 'nonlinear oscillatory' purely quantum mechanical systems that are vastly superior to the textbook Wentzel, Brillouin Kramers (or WKB) methods of the physics literature.
The latter, which date back to the early days of quantum mechanics (though with precursors that are much older still), posit a formal expansion in Planck's constant of which the leading order term is precisely a solution of the classical Hamilton-Jacobi equation for the system being quantized. But for Schrödinger operators of the conventional kinetic+potential type this amounts, among other things, to beginning with a complex oscillatory ansatz for a ground state wave function that is known, a priori, to be effectively real and monotonically decreasing in a suitable sense. That one arrives at the classical Hamilton-Jacobi equation at first approximation was an important step in the early days for convincing one that the then new quantum mechanics incorporated classical mechanics in an appropriate (formal) 'correspondence principle' limit. However important this step was conceptually, mathematically it was catastrophic.
First of all the Hamilton-Jacobi equation, even for the simplest, nontrivial, classical systems virtually never has global smooth solutions, not least of all because of the presence of caustics in the associated families of gradient-flow solution curves on the corresponding configuration manifolds. Even for the more or less tractable one dimensional systems (or those reducible to such through a separation of variables) one is obligated to effect a matching of solutions across the boundary between 'classically allowed' and 'classically forbidden' regions in the associated configuration space. It is hardly surprising that, even for such basic problems as the harmonic oscillator one only gets rather crude approximations to its actual, explicitly known wave functions. The lesser known Einstein, Brillouin Keller (or EBK) approximation method does in fact have the capability of treating higher dimensional, interacting quantum systems but only those that are completely integrable at the classical level. Maslov and others have developed ingenious ways of handling caustics when they occur but their methods are aimed at treating the finite time propagation of high frequency 'wave packets' instead of solving the fundamental eigenvalue problems that form the core of quantum mechanics.
In contrast to the above, microlocal analysts proposed a much more natural ansatz for the ground state wave function of a (nonlinear oscillatory) quantum system than that of conventional WKB theory and have shown that the leading order approximation in their approach is determined by a certain (nonlinear) 'eikonal' equation which can, for many such systems, be proven to have a global smooth solution. Higher order corrections to this leading order term can then be computed through the systematic integration of a sequence of linear, 'transport' equations along the flow defined by the previously determined solution to the eikonal problem. Excited state wave functions and their corresponding eigenvalues and higher order corrections thereto can in turn be computed from an equally systematic, somewhat similar-in-spirit transport equations analysis. The eikonal equation alluded to above is nothing more than an 'imaginary time' (or, in field theoretic language, 'Euclidean signature') variant of the classical Hamilton-Jacobi equation for the system actually under study. While one may lose a certain 'physical intuition' in dealing with this equation instead of the original, 'real time' Hamilton-Jacobi one, the mathematical gain is enormous.
In the very special case of harmonic oscillators, for example, this approach reproduces, exactly, all of the correct eigenfunctions and eigenvalues for this fundamental problem and for anharmonic oscillatory systems provides, at the very least, provably asymptotic expansions for the solutions thereto. Most importantly, by avoiding the conventional decomposition into an explicitly solvable 'unperturbed' system and its perturbation this microlocal approach secures, even at leading order, a much more accurate approximation to the actual wave functions of the interacting system than that provided by standard (Rayleigh/Schrödinger) perturbation theory. On the other hand, for reasons that we shall clarify below, its methods for technically carrying out this elegant analysis seem to be essentially limited to finite dimensional, ordinary quantum mechanical systems and, in their current form, not to allow for an extension to quantum field theories.
The authors, however, not being initially aware of these pre-existing microlocal results, embarked some years ago upon a program to develop semi-classical approximation methods that would in fact be applicable to certain interesting classes of interacting quantum field theories. In our approach the leading order term is a solution to the Hamilton-Jacobi equation for the Euclidean-signature variant of the field equations that one is quantizing. This is an enormously more favorable setting than what would be provided by the corresponding Lorentzian-signature Hamilton-Jacobi equation of conventional WKB theory and lends itself to resolution via standard methods drawn from the calculus of variations. For the special cases of massive, polynomially interacting scalar fields of renormalizable type in 2, 3 and 4 spacetime dimensions, for example, we have been able to prove the existence of globally, Fréchet smooth 'fundamental solutions' to the relevant Euclidean-signature Hamilton-Jacobi equations and to substantially characterize their asymptotic behaviors. These solutions provide the bases for the computation, via linear transport equations, of their higher order 'loop' corrections as well as that for the systematic computation of corresponding excited states. These latter steps require a detailed regularization and renormalization procedure that has not yet been carried out in the present context but, being all essentially linear in character, seem analytically less problematic than the ones we have already completed.
For the special case of non-interacting (or 'free') scalar fields our calculations truncate at finite order reproducing the functional representations of the conventional Fock-space exact solutions for these fundamental problems. The current status of this program is reviewed below in section III A.
To illustrate in some detail the general features of our Euclidean-signature semi-classical program, including its approach to the handling of higher order loop corrections and excited state computations we recall, in section II, its 'fallback' application to 'nonlinear oscillatory' problems in ordinary quantum mechanics and emphasize how our approach thereto differs from that of the microlocal literature-the key distinction that allows us to contemplate its application to quantum field theories. We also discuss therein some specific applications of our version of this microlocal, semi-classical program to certain quantum anharmonic oscillators, comparing the results with those derivable from ordinary (Rayleigh/Schrödinger) perturbation theory. A particularly interesting case is that of the quartic oscillator for which we found strong evidence that the asymptotic expansion generated for the (logarithm of its) ground state wave function may itself be Borel summable to what is conjecturally the exact solution of this paradigmatic model problem. We also briefly discuss the applicability of our methods to the problems of solving the so-called Wheeler-DeWitt equation for Bianchi IX models of (minisuperspace) quantum cosmology and to the Schrödinger operators arising in supersymmetric quantum mechanics. The latter represents our first excursion into the issue of whether the methods under study can be generalized to include fermionic degrees of freedom.
The ultimately most ambitious application of our program so far however, is to the case of quantized Yang-Mills fields in (the renormalizable instances of) 3 and 4 spacetime dimensions. In section III B we discuss the current status of this project. In particular we present the essential features of a proof that, in a suitable function space setting, the Euclidean sig-nature Yang-Mills action functional, defined over the Euclidean 'half-space' R − × R 3 , admits an absolute minimizing connection, weakly satisfying the Yang-Mills Dirichlet problem, for 'arbitrary' boundary data specified on the hypersurface {0} × R 3 . This extremized action, regarded as a functional of the boundary data (essentially a 'spatial' connection on R 3 ), provides us with the leading order approximation to (the logarithm of) our proposed ground state wave functional. We also give a preliminary discussion of the differentiability of the action functional so defined and present a sketch of how we hope to extend this argument to a full Fréchet smoothness result. A potential application of our anticipated results to the mass gap problem for quantum Yang-Mills theory is sketched in section III C.
We also include, in section III D, a brief discussion of how our Euclidean signature semiclassical program might be applied to the (currently purely formal ) functional Wheeler-DeWitt equation in a non-symmetric, field theoretic setting. A key open question is whether one can develop and analyze a genuine elliptic boundary value problem for this highly nonlinear system and, if so, how the gravitational positive action theorem can be successfully applied thereto.

MECHANICAL SYSTEMS
Elegant 'microlocal analysis' methods have long since been developed for the study of Schrödinger operators of the formĤ in the special cases for which the configuration manifold M ≈ R n , the metric g is flat and for which the potential energy function V : M → R is of a suitable 'non-linear oscillatory' type [1][2][3][4]. These methods 1 begin with an ansatz for the ground state wave function of the and proceed to derive asymptotic expansions for the logarithm, S : R n → R, expressed formally as a power series in Planck's reduced constant ( := h/2π), together with the associated ground state energy eigenvalue N is a corresponding (for us inessential) normalization constant which one could always evaluate at any (finite) level of the calculation.
When the above ansätze are substituted into the time-independent Schrödinger equation and the latter is required to hold order-by-order in powers of the leading order term in the expansion (II.3) is found to satisfy an inverted-potential-vanishing-energy 'Hamilton-Jacobi' equation given by 1 2m For a large class of (non-linear oscillatory) potential energy functions and when g is flat (with g = n i=1 dx i ⊗ dx i ) this equation can be proven to have a globally-defined, smooth, positive 'fundamental solution' that is unique up to a (trivial) additive constant. In particular this is true whenever 1. V is smooth, non-negative and has a unique global minimum attained at the origin of R n where V vanishes, 2. V can be expressed as where each of the 'frequencies' ω i > 0 for i ∈ {1, . . . , n} and wherein the smooth and the coercivity condition Since only the sufficiency of these conditions was actually established in [1] it is quite conceivable that a satisfactory fundamental solution to Eq. (II.5) exists under weaker hypotheses on the potential energy.
Our approach to proving the existence of a global, smooth fundamental solution to the (inverted-potential-vanishing-energy) Hamilton-Jacobi equation is quite different from that developed previously in the microlocal literature but has the advantage of being applicable to certain field theoretic problems whereas it seems the latter does not. 2 To establish the existence of S (0) we began by proving that the (inverted potential) action functional defined on an appropriate Sobolev space of curves γ : (−∞, 0] → R n , has a unique minimizer, γ x , for any choice of boundary data and that this minimizer always obeys lim tց−∞ γ x (t) = (0, . . . , 0). (II. 13) We then showed that every such minimizing curve is smooth and satisfies the (inverted potential ) Euler-Lagrange equation with vanishing (inverted potential) energy for each x ∈ R n we proceeded to prove, using the (Banach space) implicit function theorem, that the S (0) : R n → R, so-defined, satisfies the Hamilton-Jacobi globally on R n and regenerates the minimizers γ x as the integral curves of its gradient (semi-)flow in the sense that ∀ t ∈ I := (−∞, 0] and ∀ x ∈ R n (II.17) Actually each such integral curve γ x : I → R n extends to a larger interval, (−∞, t * (γ x )) with 0 < t * (γ x ) ≤ ∞ ∀ x ∈ R n but since, in general, t * (γ x ) < ∞ we only have a semi-flow rather than a complete flow generated by 1 m ∇S (0) . Purely harmonic oscillations on the other hand (for which A(x 1 , . . . , x n ) = 0) are an exception, having t * (γ Among the additional properties established for S (0) were the Taylor expansion formulas and where here (exceptionally) no sum on the repeated index is to be taken, and the global lower bound . . , n}. Note especially that this last inequality guarantees that, in particular, e −S (0) / will always be normalizable on The higher order 'quantum corrections' to S (0) (i.e., the functions S (k) for k = 1, 2, . . .) can now be computed through the systematic integration of a sequence of (first order, linear) 'transport equations', derived from Schrödinger's equation, along the integral curves of the gradient (semi-)flow generated by S (0) . The natural demand for global smoothness of these quantum 'loop corrections' forces the (heretofore undetermined) energy coefficients E (2) , . . .} all to take on specific, computable values.
Excited states can now be analyzed by substituting the ansatz while retaining the 'universal' factor e −S (x)/ determined by the ground state calculations.
From the leading order analysis one finds that these excited state expansions naturally allow themselves to be labelled by an n-tuple m = (m 1 , m 2 , . . . , m n ) of non-negative integer 'quantum numbers', m i , so that the foregoing notation can be refined to E (k) , k = 0, 1, 2 . . .} through the solution of a sequence of linear, first order transport equations integrated along the semi-flow generated by S (0) .
A key feature of this program, when applied to an n-dimensional harmonic oscillator, is that it regenerates all the well-known, exact results for both ground and excited states, correctly capturing not only the eigenvalues but the exact eigenfunctions as well [1][2][3].
One finds for example that the fundamental solution to the relevant (inverted-potentialvanishing-energy) Hamilton-Jacobi equation, for an n-dimensional oscillator (with mass m and (strictly positive) oscillation frequencies {ω i }) is given by and that all higher order corrections to the logarithm of the ground state wave function vanish identically leaving the familiar gaussian N is a normalization constant.
where m = (m 1 , m 2 , . . . , m n ) is an n-tuple of non-negative integers with |m| := n i=1 m i > 0, and proceeds after a finite number of unequivocal steps, to assemble the exact excited eigenstate prefactor where H k is the Hermite polynomial of order k (and While there is nothing especially astonishing about being able to rederive such wellknown, exact results in a different way, we invite the reader to compare them with those obtainable via the textbook WKB methods of the physics literature [5,6]. Even for purely harmonic oscillators conventional WKB methods yield only rather rough approximations to the wave functions and are, in any case, practically limited to one-dimensional problems and to those reducible to such through a separation of variables. The lesser known Einstein Brillouin Keller (or EBK) extension of the traditional semi-classical methods does apply to higher (finite-)dimensional systems but only to those that are completely integrable at the classical level [7]. In sharp contrast to these well-established approximation methods the (Euclidean signature 3 ) semi-classical program that we are advocating here requires neither classical integrability nor (as we shall see) finite dimensionality for its implementation.
As was discussed in the concluding section of Ref. [1] our fundamental solution, to the (inverted-potential-vanishing-energy) Hamilton-Jacobi equation for a coupled system of nonlinear oscillators has a natural geometric interpretation. The graphs, in the associated phase space T * R n , of its positive and negative gradients correspond precisely to the stable (W s (p) ⊂ T * R n ) and unstable (W u (p) ⊂ T * R n ) Lagrangian submanifolds of the assumed, isolated equilibrium point p ∈ T * R n : Another result established for the aforementioned nonlinear oscillators of Ref. [1] is that the first quantum 'loop correction', S (1) (x 1 , . . . , x n ), to the ('tree level') fundamental solution, S (0) (x 1 , . . . , x n ), also has a natural geometric interpretation in terms of 'Sternberg coordinates' for the gradient (semi-)flow generated by this fundamental solution. Sternberg coordinates, by construction, linearize the Hamilton-Jacobi flow equation to the form through, as was proven in Ref. [1], the application of a global diffeomorphism µ : R n → µ(R n ) ⊂ R n = (y 1 , . . . , y n ) , (II. 35) x → µ(x) = y 1 (x), . . . , y n (x) (II.36) that maps R n to a star-shaped domain Though not strictly needed for the constructions of Ref. [1], Sternberg coordinates have the natural feature of generating a Jacobian determinant for the Hilbert-space integration measure that exactly cancels the contribution of the first quantum 'loop correction', S (1) (x), to inner product calculations, taking, for example, where, in the last integral, the contribution of S (1) • µ −1 (y) to the wave function has precisely cancelled the non-Cartesian measure factor det g * * (y), leaving the constant (Euclidean) factor det g * * (0) in its place. Roughly speaking therefore, this role of S (1) is to 'flatten out' the Sternberg coordinate volume element, reducing it to ordinary Lebesgue measure (albeit only over the star-shaped domain µ(R n )), by exactly cancelling the Jacobian determinant that arises from the coordinate transformation.
For the nonlinear oscillators discussed in Ref. [1], Sternberg coordinates also have the remarkable property of allowing the leading order transport equation for excited states to be solved in closed form. Indeed, the regular solutions to this equation are comprised of the monomials (m) ϕ (0) (y) = (y 1 ) m 1 (y 2 ) m 2 · · · (y n ) mn (II. 40) wherein, precisely as for the harmonic case, the m i are non-negative integers with |m| := n i=1 m i > 0. On the other hand the higher order corrections, ϕ (k) (y); k = 1, 2, . . . , to these excited state prefactors will not in general terminate at a finite order as they do for strictly harmonic oscillators but they are nevertheless systematically computable through the sequential integration of a set of well-understood linear transport equations [1,2]. Formal expansions (in powers of ) for the corresponding (ground and excited state) energy eigenvalues are uniquely determined by the demand for global regularity of the associated eigenfunction expansions. More precisely one finds, upon integrating the relevant transport equation at a given order, that the only potential breakdown of smoothness for the solution would necessarily occur at the 'origin' x = 0 (chosen here to coincide with the global minimum of the potential energy) but that this loss of regularity can always be uniquely avoided by an appropriate choice of eigenvalue coefficient at the corresponding order.
A number of explicit calculations of the eigenfunctions and eigenvalues for a family of 1-dimensional anharmonic oscillators of quartic, sectic, octic, and dectic types were carried out in Ref. [1] and compared with the corresponding results from conventional Rayleigh/Schrödinger perturbation theory. To the orders considered (and, conjecturally, to all orders) our eigenvalue expansions agreed with those of Rayleigh/Schrödinger theory whereas our wave functions, even at leading order, more accurately captured the more-rapidthan-gaussian decay known rigorously to hold for the exact solutions to these problems. For the quartic oscillator in particular our results strongly suggested that both the ground state energy eigenvalue expansion and its associated wave function expansion are Borel summable to yield natural candidates for the actual exact ground state solution and its energy.
Remarkably all of the integrals involved in computing the quantum corrections S (1) , S (2) , S (3) , · · · to S (0) (up to the highest order computed in [1], namely S (25) ) were expressible explicitly in terms of elementary functions for the quartic and sectic oscillators whereas for the octic and dectic cases some (but not all) of the quantum corrections required, in addition, hypergeometric functions for their evaluation. It seems plausible to conjecture that these patterns persist to all orders in and thus, for the quartic and sectic 4 cases in particular, lead to formal expansions for S in terms of elementary functions. The evidence supporting the conjectured Borel summability of this formal expansion in the quartic case is discussed in detail in Section V.A. of [1]. By contrast, the more conventional approach (in the physics literature) to semiclassical methods leads instead to a standard (non-inverted-potential-non-vanishing-energy) Hamilton-Jacobi equation for which, especially in higher dimensions, caustics are virtually unavoidable and for which, even in their absence, a nontrivial matching of solutions across the boundary separating classically allowed and classically forbidden regions must be performed. While Maslov and others have developed elegant methods for dealing with these complications [9] their techniques are more appropriate in the short wavelength limit wherein wave packets of highly excited states are evolved for finite time intervals. On the other hand our approach is aimed at the ground and lower excited states though, in principle, it is not limited thereto.
As we have already mentioned though, our approach is a natural variation of one that has been extensively developed in the microlocal analysis literature but it also differs from this innovative work in fundamental ways that are crucial for our ultimate, intended application to field theoretic problems. In the microlocal approach [2][3][4] one begins by analyzing the (classical, inverted potential) dynamics locally, near an equilibrium, by appealing to the stable manifold theorem of mechanics [10]. One then shows, by a separate argument, that, for an equilibrium p (lying in some neighborhood U ⊂ R n ) the corresponding stable (W s (p) ⊂ T * U) and unstable (W u (p) ⊂ T * U) submanifolds of the associated phase space T * U are in fact Lagrangian submanifolds that can be characterized as graphs of the (positive and negative) gradients of a smooth function φ : U → R: This function is shown to satisfy a certain 'eikonal' equation (equivalent to our invertedpotential-vanishing-energy Hamilton-Jacobi equation restricted to U ) and φ(x) itself is, of course, nothing but the (locally defined) analogue of our action function S (0) (x). A further argument is then needed to extend φ(x) to a solution globally defined on R n .
The potential energies, V (x), dealt with in the microlocal literature often entail multiple local minima, or "wells", for which our global convexity and coercivity hypotheses are not appropriate. Much of the detailed analysis therein involves a careful matching of locally defined approximate solutions (constructed on suitable neighborhoods of each well) to yield global asymptotic approximations to the eigenvalues and eigenfunctions for such problems.
Since, however, we are focussed primarily on potential energies having single wells (corresponding to unique classical "vacuum states"), many of the technical features of this elegant analysis are not directly relevant to the issues of interest herein.
For the case of a single well, however, we have essentially unified and globalized several of the, aforementioned, local arguments, replacing them with the integrated study of the properties of the (inverted potential) action functional (II.11). When one turns from finite dimensional problems to field theoretic ones [11,12] this change of analytical strategy will be seen to play an absolutely crucial role. For the typical (relativistic, bosonic) field theories of interest to us in this context, the Euler Lagrange equations for the corresponding, inverted potential action functionals that now arise are the Euclidean signature, elliptic analogues of the Lorentzian signature, hyperbolic field equations that one is endeavoring to quantize.
While generalizations of the aforementioned stable manifold theorem do exist for certain types of infinite dimensional dynamical systems, the elliptic field equations of interest to us do not correspond to well-defined dynamical systems at all. In particular their associated Cauchy initial value problems are never well-posed. This is the main reason, in our opinion, why the traditional microlocal methods have not heretofore been applicable to quantum field theories.
On the other hand the direct method of the calculus of variations is applicable to the Euclidean signature action functionals of interest to us here and allows one to generalize the principle arguments discussed above to a natural infinite dimensional setting.
Before turning to such field theoretic generalizations however, we wish to draw attention to two particular, purely quantum mechanical extensions of the methods sketched above.
The first of these entails the application of (a suitably extended version of) the methods in question to so-called SUSY QM or 'supersymmetric quantum mechanics'. In SUSY QM a set of commuting, 'bosonic' quantum degrees of freedom {x 1 , . . . , x n } is introduced and matched by an equal number of anticommuting, 'fermionic' degrees of freedom {Ψ 1 , . . . , Ψ n }.
In terms of these a pair (in the simplest case of so-called N = 2 supersymmetry) of selfadjoint supersymmetry 'charge' operators {Q 1 , Q 2 } is defined through the introduction of a 'superpotential' function W (x 1 , . . . , x n ) and the Hamiltonian H for the system expressed algebraically as a certain anticommutator of these charges. Depending upon the choice of the superpotential such systems may exhibit either unbroken or broken supersymmetry and thus provide simplified, finite dimensional models for corresponding field theoretic systems.
The Schrödinger equation for such a SUSY QM system can be concretely represented in terms of an ordinary matrix partial differential operator acting upon a multicomponent, 'spinorial' wave function. In Ref. [13] we have begun to investigate the applicability of the 'microlocal' methods sketched above to such equations. Since fermion number is exactly conserved for such systems the corresponding Schrödinger operator assumes a block-diagonal form wherein each block corresponds to a fixed number of fermions which, thanks to the Pauli exclusion principle, varies from 0 to n. The associated Hilbert space thus includes a completely empty and a fully filled femionic sector upon each of which the corresponding Schrödinger operator reduces to an elementary, single component form to which the methods discussed above are, with only slight modifications, immediately applicable. Indeed, one can readily apply them to analyze, in depth, both the ground and excited states for the cases of both unbroken and broken supersymmetry.
The (anti-) commutation relations satisfied by the charge operators {Q 1 , Q 2 } and the Hamiltonian allow one to relate, to some extent, the eigenstates and eigenvalues corresponding to different fermionic subspaces of the full Hilbert space by, roughly speaking, using Q † = 1/2(Q 1 − iQ 2 ) and Q = 1/2(Q 1 + iQ 2 ) as (fermion number) raising and lowering operators. However, the nilpotency of these SUSY charges prevents one, when n > 2, from generating all of the eigenvalues and eigenstates through the actions of Q † and Q on the fermion vacuum and fermion-filled eigenspaces. Thus for such higher dimensional problems one must address directly the solution of nontrivial matrix Schrödinger eigenvalue problems. Since our analysis of these multicomponent problems is not yet complete we shall only remark here that, in the corresponding, 'microlocal' approach to their solution, the funda- The second quantum mechanical extension of the microlocal semi-classical methods that we wish to mention is that to a problem in 'quantum cosmology' -namely the problem of solving the relevant Wheeler-DeWitt equation for spatially homogeneous, Bianchi type IX (or 'Mixmaster') universes. Though the (partial differential) Wheeler-DeWitt equation for this model problem was first formulated nearly a half century ago, techniques for its solution that bring to light the discrete, quantized character naturally to be expected for its solution, have only recently been developed. In particular the microlocal analytical methods that we have already sketched for the study of conventional Schrödinger operators can be modified in such a way as to apply to this equation [14,15].
That some essential modification of the microlocal methods is needed for the analysis of Refs. [14,15], the scope of microlocal methods can indeed, in spite of this apparent impasse, be broadened to provide creditable, aesthetically appealing answers to the questions raised above.
While it is far from clear whether such methods could ever be further generalized to apply to the field theoretic Wheeler-DeWitt operator of (formal) canonical quantum gravity, a sketch of how that might be carried out (in a Euclidean-signature semi-classical expansion) by appealing to the 'positive action theorem' is given in the concluding section of Ref. [14] and reviewed in Section III D below. In particular, we draw attention there to several remarkably attractive features of such an approach and show how it avoids some of the serious complications that obstructed progress on the, somewhat similar-in-spirit, Euclidean-pathintegral approach to quantum gravity.

III. MODIFIED SEMI-CLASSICAL METHODS FOR BOSONIC FIELDS
A. An application to polynomial scalar field theories For technical reasons, the elegant 'microlocal' methods developed in the past for the analysis of conventional Schrödinger eigenvalue problems, have not heretofore been applicable to quantum field theories. In this section we describe a 'Euclidean signature semi-classical' program to extend the scope of these analytical techniques to encompass the study of selfinteracting (massive) scalar fields with polynomial renormalizable self-interaction in 1 + 1, In this section, we apply our method to suitable polynomial field theories and establish existence, uniqueness and global regularity for a globally defined 'fundamental solution' of the Hamilton-Jacobi equation mentioned above, via an application of the direct method of the calculus of variations, elliptic theory and a Banach space version of the implicit function theorem.
In particular, such polynomial field theories include the cases of Φ 4 scalar fields in 4dimensional Minkowski spacetime, Φ p fields with exponents p = 4, and 6 in 3-dimensional Minkowski space, and with p an arbitrary, positive even integer greater than 2 in 2dimensional Minkowski space.
For each case, global existence and smoothness (in a suitable function space setting) of the relevant Hamilton-Jacobi functional, S (0) [ϕ] is proven and it is discussed how our tree approximation for these ground state wave functionals captures the more-rapid-thangaussian decay that should hold for the exact solutions. Since, in our setup, the squared modulus of the ground state wave functional provides the natural integration measure for the associated Hilbert space of quantum states, it seems encouraging that our approach exhibits this non-Fock-like behavior already at leading order. By contrast note that, to any finite order, conventional Rayleigh/Schrödinger theory would generate instead an approximate wave functional that decays more slowly than the corresponding gaussian. Because our techniques, even when applied to quantum mechanics problems, do not simply reproduce the standard results of Rayleigh/Schrödinger perturbation theory, we expect that one should be able to generate much more accurate approximations to the Hilbert space of quantum states for certain quantized fields.
Although in this section we only deal with the aforementioned fundamental solutions, the fields considered are rigorously proven to be non-trivial in lower dimensions; [16][17][18][19][20].
On the contrary, Φ 4 fields in 4 dimensions are often believed (though still not rigorously proven [21,22]) to renormalize to (trivial) free fields. If such a conclusion were proven to be correct, it would only emerge at the level of the higher order quantum, 'loop' corrections, because the necessity to regularize otherwise ill-defined functional Laplacians only arises in our approach at the level of the transport equations for these higher order corrections. The study of the higher order, 'loop' corrections, for both ground and excited states, is work currently in progress.
The formal Schrödinger operator we consider for a scalar field in n + 1 dimensions, with n = 1, 2, or 3, iŝ in which x ′ ∈ R n , ∇ ′ is the 'spatial' gradient on R n , and the self-interaction polynomial P(·) ≡ k j=2 a j (·) j is assumed to be convex, to possess no linear term, to include a 'mass' term, a 2 > 0, and to be of even degree k, with in which d ≡ n + 1 is the dimension of the domain. The number 2d/(d − 2) above is the critical exponent for the elliptic theory applied to the Euclidean signature action functional which is guaranteed to be coercive by the hypotheses imposed upon P. (These assumptions are satisfied by the 'massive' Φ 4 theory, for which P The functional Laplacian, formally defined as 'trace of the Hessian', will require some 'regularization' in order to be well-defined on the (wave) functionals of interest; in fact, although the Hessians of the latter functionals will be smooth, they need not be of trace class.
that arises at lowest order from substituting our ansatz for the ground state wave functional into the Schrödinger equation and demanding that the latter hold order-by-order in powers of relative to the formal In the formulas above ϕ(·) denotes a real-valued distribution on R n , boundary data for a real spacetime scalar field Φ defined on (−∞, 0] × R n . As for the finite-dimensional problem described in Section II, a 'fundamental solution' to Eq. (III.4) is constructed by proving the existence of unique minimizers Φ ϕ for the Euclidean signature action functional (III.3) for arbitrary boundary data ϕ specified at t = 0, and The minimization procedure is carried out for Φ in the space of distributions The latter isomorphism of spaces holds because the trace maps, tr : Theorem 7.58 in [28]. In particular, one may choose an extensionφ which is smooth in the interior; cf. Lemma 5.1 in [11]. 5 Coerciveness, weak lower sequential semicontinuity, and strict convexity of I es on A(ϕ) (cf. Theorem 1.1.3 of Ref. [23], and Chapter 3 of [11]) guarantee existence and uniqueness for the absolute minimizer Φ ϕ ; [11]. Uniqueness of the absolute minimizer, although not necessary to guarantee that S (0) [ϕ] be well-defined, is required for 'everywhere' differentiability of S (0) with respect to the boundary data to hold. 6 5 The above choice of the space of boundary data B, although not the natural one for a minimization procedure taking place in H 1 (R − × R n ), is required for the Hamilton-Jacobi equation to hold strongly, besides being instrumental to implementing the first step in the regularity theory, and constitutes a nonsubstantial restriction. Weakening the requirement on the boundary space B, for example by imposing the boundary data to be in H 1/2 ({0} × R n ), would not in substance yield a stronger theorem. In that case in fact, the absolute minimizer for I es would then be in H 1 (R − × R n ) rather than in H 3/2 (R − × R n ), and the Hamilton-Jacobi equations would not hold strongly. An analogous remark can be made for the case of the Dirichlet problem for Yang-Mills connections analyzed in the next section. 6 In the next section we will argue for the case of Yang-Mills connections that non-uniqueness is an obstruction to everywhere differentiability for the suitable analogue of S (0) ; with regard to this, see also the discussion under the title "Smoothness of the Hamilton-Jacobi functional" further below in the present section.
A candidate for the tree approximation to the ground state wave functional for the P(Φ) theory is the functional in which N is a normalization constant.
The weak Dirichlet problem satisfied by the absolute minimizer. With regard to existence and regularity for the elliptic theory associated to the functional I es via the calculus of variations, in [11] the authors follow the essential lines of the classical references [24,25], with some perhaps minor modifications of those. Here, we highlight the general outline and the main differences in comparison to the classical literature. We start by providing fundamental definitions based on the standard theory of the calculus of variations. Once these well-founded definitions are established and acquired, one can develop the resulting theory, with no need for repetitive reproduction of those in each specific context. What perhaps may well be simply a slight shift in the philosophical approach to boundary value problems, affects the coherence and elegance of the theory and, in some instances, its content.
With that in mind, we define a weakly differentiable function Φ to be a weak solution to the Dirichlet Problem , the topological dual space of A(0), and ϕ prescribed in (III.14) In comparison with Ref. [24], notice that here we do not assume the boundary value and its extension to the interior to be C ∞ , in replacement, we elect to , when (III.13) is satisfied, it applies also to all ω ∈ A(0). Definition (III.12) implicitly defines the weak Dirichlet Laplacian on A(ϕ) through the identification ∆Φ ≃ R − ×R n −∇Φ · ∇(·) dx (that is, we regard ∆Φ, for each Φ ∈ A(ϕ), as an In analogy, we define a weakly differentiable function Φ to be a weak solution to the nonlinear Dirichlet Problem in which Q is a given nonlinear function of a single variable, such that Q(Φ) ∈ In the definition above, Λ ≡ −∆ + Q(·) is a nonlinear differential operator with image in In the context of the polynomial theory under study, we set Q ≡ P ′ . In fact, the conditions imposed on P described in the current section guarantee that for some constant C independent of Φ, thus ΛΦ is a continuous linear operator on A(0).
In wide generality, given a differentiable function of a real variable P, satisfying P ′ (Φ) ∈ possesses a stationary point Φ ∈ A(ϕ), such Φ would also be a weak solution to the nonlinear Dirichlet problem (III.15), and vice versa (by a density argument).
In fact, for Φ ∈ A(ϕ), ω ∈ A(0), one has In the context of the polynomial theory under study, P(·) is a polynomial satisfying conditions which guarantee the existence of a unique (absolute) minimizer of J = I es on A(ϕ), with prescribed boundary value ϕ ∈ B, and such a minimizer is also its unique stationary point, that is, the (unique) weak solution to the Euler-Lagrange equations contained in (III.15), with Q substituted by P ′ .
Elliptic regularity and global control. Having ascertained that the unique absolute minimizer of I es over A(ϕ), with prescribed ϕ ∈ B, satisfies weakly the nonlinear Dirichlet with f ≡ −P ′ . In the critical cases instead, because no improvement would be obtained at the first step of the bootstrapping for (III.12) with f ≡ −P ′ , one needs to implement a special initial step which, by exploiting the existence theory, replaces system (III.20) by the 20) satisfies by construction the linear Dirichlet-type problem (III.21), the latter having been derived by freezing the coefficient g(Φ). The estimate g(Φ) ∈ System (III.21) still presents itself as 'borderline', thus requires a non-straightforward technique. Its regularity is accomplished in [11] by applying a 'regularity lifting method by means of a contracting operator' 7 . For the method of contracting operator and more details on its adaptation to the context at hand, see [11,27].
At various stages of the proof of regularity, additional difficulties arise, due to the nontriviality of the boundary data ϕ prescribed at {0} × R n , and to the unboundedness of the domain R − × R n . In fact, the regularity lifting method, as described in [27], directly applies only to linear Dirichlet-type problems with prescribed vanishing boundary data on bounded domains Ω. Further work is also required because, in the context of the Euclidean-signature semi-classical program under study, we are interested in establishing global control of the solution over the unbounded domain R − × R n (and over the slices {t} × R n ).
Boundary issues are dealt with by reflection of operators and associated function spaces across the boundary, and introduction of cut-off functions, in a suitable way, thus transforming Dirichlet-type problems on boundary-type neighborhoods like in which a = (0, a 1 , a 2 , a 3 ) is a fixed boundary point, and R a positive number, into Dirichlettype problems on balls with homogeneous boundary data. The doubling technique may involve technicalities. In our context, because R − × R n is endowed with a Euclidean metric, the resulting doubled boundary-type neighborhoods are flat, while in more general cases, the latter would possess a Lipschitz-bounded metric; cf. for example, [12,26].
To implement the procedure, we first extend the boundary value ϕ ∈ B to the unique Regularity of (III.21) could also have been accomplished by an alternate method, perhaps unknown or scarcely utilized outside the field of geometric analysis, based on the regularity theory as developed several years earlier by one of the authors in the context of Yang-Mills connections; cf. [26] and the 'Acknowledgements' section in the present article. This is one of some instances in which techniques developed in gauge theories anticipate methods later established and formally codified in classical analysis contexts.
Suchφ is proven to satisfy the global estimateφ ∈ H 3/2 (R − × R n ), besides the regularity , that is, smoothness up to and including the boundary; for the technical details, we refer the reader to Lemma 5.1 in [11].
Due to technicalities pertaining to unbounded domains 8 , to the purpose of establishing a global estimate for Φ, we also consider the extension Φ L , unique solution to The same results hold for Φ L as obtained forφ, thus both those functions are smooth in the interior and satisfy the best possible global estimate on R − × R n ; cf. Lemma 5.3 in [11].
After performing the various steps of what presents itself as a somewhat laborious and at times technical procedure, one ultimately obtains the regularity theorem below (Theorem 5.1 in [11]).
Φ is smooth all the way up to and including the boundary.
Thus, also Φ, unique solution to the non-linear problem (III. 20), and unique minimizer of I es over the space A(ϕ), satisfies interior regularity and the best possible global estimate In order to implement the initial step to interior regularity, with focus on the critical cases, one shows that αΦ, with α a suitable cut-off function centered at a point a ∈ R − × R n , is a 8 These arise due to the additional constraints on the Sobolev exponents required for Sobolev embeddings to hold on unbounded domains; namely, on unbounded domains Ω u having the cone property, in dimension d > hq, the embeddings L q h (Ω u ) ⊂ L p (Ω u ) require the additional hypothesis q ≤ p (besides the hypotheses h/d + 1/p − 1/q ≥ 0, p, q ≥ 1, h ≥ 0, required on bounded domains); cf. Theorems 7.57, 7.58 in [28]. In the proof of Theorem III.1, we use repeatedly the embeddings Then, bootstrapping on systems such as The initial step to boundary regularity, still with focus on the critical cases, is not sostraightforward. One first observes that the function Φ−φ is a weak solution to the Dirichlet problem with homogeneous boundary data, (III. 28) in which g(Φ) is the even extension of g(Φ) andř is the odd extension of −g(Φ)φ. Because even and odd extensions both preserve membership in The functionΦ, defined as the odd extension of Φ −φ, is then a weak solution in H 1 (R d ) to Eq. (III.28).
Let then a = (0, a 1 , a 2 , a 3 ) be a fixed boundary point, U R a boundary-type neighborhood as defined in (III.22), B a;R its double, and α a smooth cut-off function, compactly supported in B a;R , with value 1 on B a;R 1 , with 0 < R 1 < R. Then, ∆α is a lower order linear differential operator with smooth coefficients determined by the cut-off α, applied toΦ.
Unfortunately, in dimensions d > 3, the termř, odd extension across {0} × R n of Φ −φ, accounts for the lesser regularity of the right hand side in (III.30), which now no longer satisfies the hypothesis satisfied by its counterpart f 1 , namely, f 1 ∈ L 2 (B a;R ) ⊂ L d/2 (B a;R ) in dimensions 3 and 4, in System (III.25), utilized for the proof of interior regularity. The latter hypothesis is explicitly formulated, and required, in Theorem 3.3.2 of Ref. [27]. The need to work with lesser regularity would often arise when studying boundary regularity for Dirichlet-type problems such as (III.21) in which non-vanishing boundary data in a borderline case are prescribed, and it seems worthwhile to extend the regularity lifting theorem to include those boundary cases. That is the purpose of the following lemma.
The corresponding hypotheses in Theorem 3.3.2 of Ref. [27] are a(x) ∈ L in which X A , X B are the characteristic functions of Ω A , and Ω B . As in Ref. [27], one can show that, if A is a sufficiently large constant, the operator T (aX A ·) : L p (Ω) → L p (Ω) is a contracting operator ∀p ∈ ( d d−2 , ∞), making use of the asymptotic behavior (III.32). There are only minor changes, in comparison with Ref. [27], in the procedure to estimate the term is not present in Ref. [27]. Nevertheless, despite this limitation, a first improvement is all that is needed in a variety of cases, analogous to ours, in which the non-homogeneous term b depends on the solution Φ to a non-linear system related to (III.27), and one can feed that initial improvement obtained for u into the term b and iterate. Thus, Lemma III.1 may be instrumental to a somewhat general extent, if one wishes to achieve boundary regularity by our doubling technique, that is, by transforming a Dirichlet-type problem with non-homogeneous boundary data on a boundary-type neighborhood, into an interior neighborhood with a (possibly) Lipschitz-bounded metric and homogeneous boundary data. (Ω), on bounded domains Ω; cf. Ref. [29].
Having established Lemma III.1, a repeated application of doubling techniques, and procedures conceptually similar to those applied so far, yield the results stated in Theorem III.1; cf. [11] for technical details, in particular with regard to the proof of higher order regularity.
In the latter reference, we show at an early step thatΦ ≡ Φ−φ satisfiesΦ ∈ C 1 ((−∞, 0]×R 2 ) in 3 dimensions, whileΦ ∈ C 0 ((−∞, 0] × R 3 ) in 4 dimensions, for general prescribed boundary data in ϕ ∈ B. If the boundary value satisfies ϕ ∈ C ∞ ({0} × R 2 ), by differentiating the auxiliary equations satisfied byΦ, obtained in [11], at first with respect to tangential derivatives, then relating those to normal derivatives, we eventually obtainΦ and Φ in To the purpose of establishing global control of Φ over the unbounded domain R − × R n , having already proved, besides interior regularity, the global estimatesφ, The coefficients appearing in System (III.27) belong only locally to certain Sobolev spaces; that is because their various summands fail to lie in L p (R − × R n ) for any shared value p.
That issue can be easily circumvented in our case. In fact, replacing −∆ with Λ L , allows one to rewrite System (III.20) as (III. 33) in which the polynomial R no longer contains the lower degree term in P ′ . Replacingφ by , and a cutoff function should not be introduced for the case at hand. To the purpose of establishing global control, the cases d ≡ n + 1 = 4, d = 3, and d = 2 are treated separately in [11]. One ultimately obtains the following The proof relies on an extension to unbounded domains of the regularity lifting method previously shown. Such extension relies on constructing contracting operators here, X A is the characteristic function on Ω A , the definition of which, adapted to the unbounded domain case, is for A a positive fixed constant, and T 1 is the solution operator for , and that is a contracting operator, yielding the Smoothness of the Hamilton-Jacobi functional S. There are in general two obstructions to smoothness of Hamilton-Jacobi functionals G defined, in analogy to the way S has been defined in the present section, as in which G is a given action functional, ϕ are specified boundary data in a function space In the context of the Yang-Mills theory, discussed in the next section, the lack of uniqueness for the absolute minimizer for the Euclidean action corresponds to a lack of 'everywhere' differentiability of the analogue to the functional S. In that context, the relevant Riemannian space will be the moduli space of spatial (i.e., boundary data) connections modulo gauge transformations with its (weak) Riemannian structure induced from the Yang-Mills kinetic energy in a well-known way; cf. , for example, [31] and [32]. In that context, wave functionals exhibit scars localized along lower dimensional subsets of the (moduli-space) configuration manifold, to some extent analogous to the scars that arise, in a totally different setting, in the study of quantum chaos (see, for example, [6]). In the present section, we rule out those obstructions for the polynomial theories, thus showing smoothness of the functional (III.9) as an application of the implicit function theorem between Banach spaces. More precisely, the approach used in [11] entails proving Theorem III.2 below, under the additional strict convexity assumption on P, Theorem III.2 The solution Φ ϕ to the nonlinear Dirichlet problem (III.20) depends smoothly on the boundary data; more precisely, is a smooth map.
Theorem III.2, combined with smoothness of the functional I es : The proof of Theorem III.2 is essentially an application of the Implicit Function Theorem to the functional 42) in which Λ ≡ −∆ + P ′ (·), and Φ L ϕ is the unique solution to the linear boundary value problem (III.24) previously defined. Such functional is smooth because the map ϕ ∈ B → We linearize E(ϕ, h) at the point (ϕ 0 , h 0 ) along the direction of the second variable h, thus deriving the linear operator and H −1 (R − × R n ). That is achieved by showing that the standard inner product in H 1 (R − × R n ) and the inner product induced by P given by determine equivalent norms, thus induce the same natural topology. Such equivalence, yields in particular that D 2 E 0 (ξ) is a bounded linear operator on A(0). Note that, for any g ∈ A(0), one has by definition D 2 E 0 (ξ)(g) = ξ, g P , and D 2 E 0 is a one-to-one, self-adjoint operator. The implicit function theorem between Banach spaces then states that there exist neighborhoods I ⊂ B, J ⊂ A(0) of ϕ 0 , h 0 respectively, such that ∀ϕ ∈ I there exists a unique h(ϕ) ∈ J for which E(ϕ, h(ϕ)) = 0 and that such map, ϕ ∈ B → h(ϕ) ∈ A(0), is , this ultimately yields smoothness for the functional S. The Hamilton-Jacobi equations for the absolute minimizer of the Euclidean signature action functional I es . For any given distribution Φ ∈ A(ϕ), finiteness of the integral R − ×R n P(Φ) dt dx ′ , together with finiteness of the Euclidean action, guarantee that the energy functional having energy density be finite. Here and throughout the text, ∇ ′ denotes the gradient with respect to the spatial variables only, and we identify x 0 ≡ t, ∂ 0 Φ ≡ ∂Φ/∂x 0 ≡Φ .
Thus, by Fubini's theorem, the energy of a distribution Φ ∈ A(ϕ), is finite for almost all t's in R − and integrable over R − .
Further, for Φ a minimizer of I es over the affine space A(ϕ), a simple calculation, using the Euler-Lagrange equations satisfied by Φ as well as the interior regularity proved for such equations in this context, shows that As a result, by interior regularity, if Φ is a minimizer, the Fundamental Theorem of Calculus applied to the smooth functions of time If the boundary value ϕ ∈ B is in addition prescribed to be smooth, ensuring smoothness of the minimizer Φ all the way up to and including the boundary, the equalities above hold for all t 1 < t 2 ≤ 0. Integrating (III.51) over R n , one obtains as long as the quantities on the left hand side are finite, that is, outside possibly a subset of R − of zero measure.
A computation, involving approximation of the integral on the right hand side of (III. 52) via integrals over the regions { x ′ = R} × [t 1 , t 2 ] and an application of Green's theorem, in which dσ is the surface element on the cylinder { x ′ = R} × [t 1 , t 2 ], and t 1 , t 2 can be taken arbitrarily outside possibly a set of zero measure. At that point, one shows that in which C is a constant, whenever Φ is a minimizer of the action functional I es over the space A(ϕ), with boundary value ϕ ∈ B; cf. [11] for details. Such constant C depends on Φ only and, ultimately, by uniqueness of the minimizer in the case of the polynomial scalar field theory, is uniquely determined by the prescribed boundary value ϕ ∈ B.
Further, the global control established for the action minimizer Φ, namely Φ ∈ H 3/2 (R − × R n ) (the best possible global estimate for prescribed boundary conditions ϕ in Having shown already that the solution Φ ϕ to (III.20), unique minimizer of I es with prescribed boundary value ϕ ∈ B, depends smoothly on the latter, one has, in particular, is differentiable with respect to λ. One also has the relation v ′ (0) |{0}×Rn = ψ.
The variational derivative with respect to ϕ ∈ B, evaluated at ψ, of the functional S defined in (III.9) is therefore Combining the vanishing of the energy at time t = 0, with (III.58), one derives the Hamilton-Jacobi equation Decay of the approximate ground state wave functional. A 'virial' argument.
A first straightforward estimate for Ω 0 , the tree approximation to the ground state functional (cf. Eq. (III.11)), of the type in which C is a constant independent of ϕ, can be proven as a consequence of coerciveness of I es and of the Trace Theorem, which guarantees an estimate of ϕ H 1 2 (R n ) in terms of Φ ϕ H 1 (R − ×R n ) . (This argument does not apply to the massless Φ 4 theory on R − × R 3 because I es is not coercive in that case.) A heuristic argument for a more refined estimate, yielding a better decay, is based on a socalled virial argument, which takes into account the presence of the higher order polynomial term in I es . A rigorous argument would require one taking into account the higher order corrections of the ground state functional, S (1) , S (2) . . . , which we are disregarding at the moment. Such quantum corrections all vanish in the case of (massive) free fields, but do not vanish for the more general polynomial scalar field theory under study.
Having made the coerciveness assumption P(Φ) ≥ CΦ 2 for some positive constant C, by defining m 0 ≡ √ 2C and considering the corresponding massive free field functional (For an explicit calculation of S f ree (0) [ϕ], see [33].) Thus, exp{−S (0) [ϕ]/ } decays at least as rapidly as some specific Gaussian. Here, we are still working under the assumption that the ground state quantum corrections to S (0) [ϕ] may be legitimately disregarded for the purpose of estimating the decay of the true ground state functional.
In [11] we make a conjecture on the behavior of the S (0) functional under a constant rescaling ϕ → ϕ λ = e λ ϕ, in the limit of large ϕ. To that purpose, we consider the ratio and observe that, in the case of the massive free field functionals S (0) [ϕ] f ree described above, this ratio would simply be given by R = 2 for all ϕ.
Assuming that lim ϕ H 1/2 (R n ) →∞ R(ϕ), exists, we evaluate lim t→t * R(ϕ t ) , (III. 62) in which ϕ t is a curve satisfying ϕ t H 1/2 (R n ) → ∞, as t approaches its limiting value t * . Using the fact that both numerator and denominator of (III.61) tend to infinity as ϕ H 1/2 (R n ) → ∞, we appeal to L'Hospital's rule and differentiate with respect to t numerator and denominator in (III.62). To that purpose, we further use the technique of estimating such limit along 'solution curves' of the 'gradient semi-flow' of S (0) [ϕ]. We compute formal time derivatives along the flow by applying the functional differential operator The Hamilton-Jacobi equation then simplifies the formula for this ratio of 'time' derivatives to Luckily the formula above is independent of S (0) (ϕ) and only depends upon the given 'potential energy' from the Hamilton-Jacobi equation. This simplification is the main reason for proposing to compute the flow along the (Hamilton-Jacobi) solution curves instead of along the 'rescaling curves' ϕ λ = e λ ϕ. Note also that in the absence of the higher order terms eq. (III.65) immediately reproduces the free field result T → 2 without the need for taking ϕ 'large'. In the general case, Sobolev estimates show that each polynomial term except the top order ones will tend asymptotically to 0 as ϕ(·) H 1 (R n ) → ∞, leaving the result T → 1 + k/2. This is the 'intuitively expected' result since, from the form of the Hamilton-Jacobi equation satisfied by S (0) [ϕ], we seem to need S (0) [ϕ] scaling like ϕ 1+k/2 for large ϕ in order to match the behavior of the given potential energy for large arguments.
A difficulty in making the above argument rigorous arises through the fact that extension of the formal 'integral curves' of the 'gradient semi-flow' of the functional S (0) [ϕ] to positive t's (for data specified at t = 0), does not necessarily make sense. Using our regularity results there is a clear mathematical sense to such 'curves' for t < 0 but, to extend them in the opposite temporal direction seems problematic in general, especially when the boundary data chosen is as 'rough' as possible. However, such rough data cannot arise at an interior point of such a hypothetically extendible curve. In fact, the global control and interior regularity, established in [11] and outlined in the present section for a solution Φ ϕ to the nonlinear Dirichlet problem (III.20), ensure Φ ϕ ∈ H 3/2 (R − × R n ) ∩ C ∞ (R − × R n ). Thus, , ∀t < 0, and one could regard the smoothed interior data at some t < 0 as new 'initial data' for a curve that is in fact extendible (at least back to the original t = 0 starting point) and make presumably precise sense of the argument for a dense subset of the full space of initial data. To compute such higher order 'loop' corrections for bosonic field theories, one will first need to regularize the formal functional Laplacian that arises in the Schrödinger operator (III.1) and that will reoccur in each of the transport equations which result from substituting ansatz (III.5) into the time independent Schrödinger equation (III.6) and expanding formally in powers of . Solving these transport equations for the 'loop corrections', , · · · , to the ground state wave functional simply amounts to evaluating sequentially computable, smooth functionals on the Euclidean action minimizers, Φ ϕ , for arbitrary chosen boundary data, ϕ.
Solving the transport equations for excited states is somewhat more involved since these equations entail a lower order term in the unknown but the technology for handling this, is well understood [1][2][3]. If, for example, a Sternberg diffeomorphism could be shown to exist, then the leading order, excited state transport equation could be solved in closed form (as already established for nonlinear oscillators).
Otherwise, one could simply fall back on the machinery developed in Refs. [1][2][3], which does not assume the existence of Sternberg coordinates, and solve this and the corresponding higher order excited state equations in a less direct fashion. In either case it is intriguing to note that the excited states for interacting field theories would be naturally labeled by sequences of (integral) 'particle excitation numbers' in much the same way that the Fockspace excited states of a free field are characterized.
One often hears that the fundamental particle interpretation of interacting quantized fields hinges upon their approximation by corresponding free fields. This is unsatisfactory at best since, of course, an elementary particle cannot 'turn off' its self-interactions in order to conveniently behave, even asymptotically, like a Fock-space, free field quantum. As we have already emphasized one of the natural features of this (Euclidean signature-semi-classical) program is that it maintains the dynamical nonlinearities of an interacting quantum system intact at every level of the analysis rather than attempting to reinstate nonlinear effects gradually through a perturbative expansion.
B. An application to Yang-Mills fields As in the foregoing models, construction of a Euclidean-signature-semi-classical wave functional for Yang-Mills theory proceeds by showing the existence of a solution to the Euclidean signature Dirichlet problem for initial data prescribed on {0} × R 3 . The existence proof takes the form of a localizing and counting argument, based on the collective work of Uhlenbeck [34,35], Sedlacek [36], and Marini [26]. In particular, two results of Uhlenbeck are crucial: first, that an L n/2 bound on the curvature of a connection yields a representative in Hodge gauge, the L p 1 norm of which is bounded by the L p norm of its curvature [34], and second, the removability of singularities for Yang-Mills connections in dimension 4 through application of an appropriate gauge transformation [35]. In [26], smooth solutions to the Dirichlet and Neumann problems for Euclidean signature Yang-Mills theory are shown to exist on a compact manifold with smooth boundary. Because our base manifold in the present program is R − × R 3 (Euclidean spacetime), we extend this work to accommodate a noncompact base manifold with smooth boundary. The counting argument in [36] depends on compactness of the base manifold; thus we give a new proof of the relevant result. Explicitly allowing for a noncompact base manifold obviates the need to restrict gauge transformations to approach a consistent value at spatial infinity, as would be necessary if one proceeded by compactifying spacetime and invoking results in [26]. If gauge transformations were required to approach a coherent limit at spatial infinity, the space of physical connections on spacetime would divide into disjoint topological sectors, and a distinction could be introduced between 'large' and 'small' gauge transformations as in [37] according to homotopy class. Under the approach in [37], this distinction gives rise to the 'vacuum angle' of quantum Yang-Mills theory, whereby physical states are only invariant up to a phase under large (i.e. non-null-homotopic) gauge transformations. However as noted in [37] and [38], a vacuum angle can alternatively be introduced at the level of the Lagrangian, leaving this avenue open to our approach if a vacuum angle is physically indicated.
We present our analysis for the case of Yang-Mills theory in 4 spacetime dimensions, but we note here that our results allow analogous construction of a candidate leading-order ground state wave functional in lower dimensions. Because the quantization of Yang-Mills theory in 3 spacetime dimensions has been well addressed by the methods of Karabali, Kim, and Nair (see e.g. [39]), in future work we aim to make contact with this program as a potentially illuminating comparison. We discuss the possible application of the Banach space implicit function theorem to establish Fréchet differentiability of S to all orders.
Although our leading-order ground state wave functional will be constructed on R − × R 3 , we carry out the minimization of the Euclidean signature Yang-Mills action more generally over a smooth n-dimensional Riemannian manifold (M, g), not necessarily compact, with smooth boundary ∂M. The Yang-Mills structure group is taken to be a compact semisimple real Lie group G, with Lie algebra g. We denote by P a principal G-bundle over M. Following [40], associated bundles η having fiber V are constructed by specifying a representation ρ : G → Aut V and forming the twisted product 9 η = P × ρ V . In particular, the automorphism bundle of η is Aut η ≡ P × ad G, where ad : G → Aut G denotes the action of G on itself by conjugation. Gauge transformations are given by sections Γ(Aut η). The adjoint bundle Ad η ≡ P × Ad g, where Ad denotes the adjoint representation of G on g, is used to construct g-valued differential k-forms on M as sections of Ad η ⊗Λ k M. A connection on an associated bundle η is a first-order differential operator D : Γ(η) → Γ(η⊗T * M), which can be expressed in terms of a base connection D 0 as D 0 + gA for A ∈ Γ(Ad η ⊗ T * M). The coefficient g represents the Yang-Mills coupling constant. Connections are gauge transformed through 9 Note that, as in [41], for a topological space X admitting a continuous right action by a topological group G and a topological space Y admitting a continuous left action by G, the twisted product X × G Y is defined as the Cartesian product X × Y modulo the equivalence relation (x.g, y) ∼ (x, g.y). conjugation by elements of Γ(Aut η): (III.66) The curvature of a connection F (D) ∈ Γ(Ad η ⊗ Λ 2 M) is defined by and transforms by conjugation under gauge transformations: Locally, we can write Central to defining the Euclidean signature Yang-Mills action and obtaining its minimizer for given boundary data on M is the construction of Sobolev spaces for g-valued differential forms. An L 2 inner product on sections of Ad η ⊗ Λ k M is given by where for θ = T ⊗ µ ∈ Γ(Ad η ⊗ Λ k M), * θ ≡ T † ⊗ * µ, with T † the adjoint in g and * µ the Hodge dual with respect to the metric g. A G-invariant positive-definite inner product results, since in the adjoint representation of g, tr XY † = −tr XY , the negative of the Killing form, which thanks to compactness and semisimplicity of g is negative definite. The L p norm of a g-valued differential k-form φ can then be defined by Accordingly, the Euclidean signature Yang-Mills action is defined as where F µν , F ρσ ≡ tr F µν F ρσ † , and the index I runs over a basis of the Lie algebra g normalized with respect to (·, ·). Summation with respect to repeated up and down indices is implicit (Einstein summation convention). Gauge invariance of the action functional follows from G-invariance of the form inner product. The Yang-Mills Dirichlet problem results from varying this action with fixed initial data to derive Euler-Lagrange equations: where ∼ denotes a gauge transformation on ∂M which can be extended to the interior.
for all bounded open sets U ⊂ M.
Although the setting is different, our proof is similar in method to Sedlacek's proof of Lemma 3.6 in [36]. Since norms are weakly lower semicontinuous, we show that the is immediate from our hypotheses, while for the second, weak convergence of a suitable subsequence follows by using the Principle of Uniform Boundedness together with continuity of the embedding L 2 1 ֒→ L 4 and the multiplication L p 1 (U) (and for a boundary neighborhood, additionally i * A ≡ A τ | U ∈ L p 1 (∂U ∩ ∂M)), where n/2 ≤ p < n. Then there exists K ≡ K(n)/g > 0 and c ≡ c(n) such that every connection D = d + gA satisfying F L n/2 (U ) < K (and for a boundary neighborhood, additionally F τ L n/2 (∂U ∩∂M ) < K) is gauge equivalent to a connection d + gÂ ,Â ∈ L p 1 (U, Ad η ⊗ T * U), satisfying Moreover, the gauge transformation s used to obtainÂ = 1 g s −1 ds + s −1 A s can be taken in L n/2 2 (U) (s will in fact always have one more degree of regularity than A ; see Lemma 1.2 in Ref. [35]).
As observed in [26], conformal invariance of the condition F n/2 < K allows us to achieve the simultaneous condition F τ L n/2 (∂U ∩∂M ) < K for a half-ball U at the boundary by means of a dilationx = rx, since F L n/2 (∂U ∩∂M ) acquires a factor of r. For the case n = 4 as in the physical problem, the requisite L n/2 (U) bound on F is the same as a bound on the Euclidean signature Yang-Mills action restricted to U. Thus to apply Theorem III.4, we seek a cover of M by balls on the interior and half-balls centered at boundary points such that the Yang-Mills action restricted to each ball U in the cover meets the condition F L n/2 (U ) < K. We prove a general result for a sequence of positive L 1 -bounded functions on R n (or indeed an arbitrary Riemannian manifold, with or without boundary) assuring the existence of a cover (missing at most a finite collection of points) for which the L 1 norm of the functions restricted to balls in the cover is less than an arbitrary fixed bound ε. For the case of Yang-Mills theory, the possibility that our cover may miss finitely many points will be addressed by removability of singularities, as discussed below.
For simplicity, we state and prove our 'good cover' theorem on R n , but the modifications for a general Riemannian manifold are obvious.
Theorem III.5 Consider a sequence of positive functions {f i } i∈N on R n such that for some constant C.
Then for any arbitrarily fixed ǫ > 0, there exists a subsequence of {f i } and an open cover G made up of balls B ⊂ R n , covering R n except at most a finite number of points P 1 , P 2 , . . . , P N and satisfying where the index i B depends on the ball B ∈ G .
For the proof, we construct a countable cover C of R n as the union of countable covers of R n by balls of radius 1 j for each j ∈ N. Fixing an enumeration of the balls in the cover, we extract a subsequence {f converges. From this subsequence, we extract a further subsequence {f (2) i } so that on the second ball B 2 , converges, and so on. We observe that for the diagonal subsequence {f converges on each ball B k in C . If the limit of this sequence for some B k is greater than or equal to the bound ε, we discard the ball B k from C . Finally, we show that the resulting collection G of balls, for each of which converges to a value less than ε, covers all but finitely many points of M. To do so, we use the uniform bound C on the L 1 norms of the f to deduce that the number of pairwise disjoint balls of radius 1 j which have been discarded for each j ∈ N cannot exceed [C/ε]. For each j ∈ N we then select a maximal set of pairwise disjoint discarded balls of radius 1 j . By expanding their radii to 3 j , we ensure that their union S j contains all the discarded balls of radius 1 j . This fact together with the bound on the number of discarded balls for each j ∈ N allows us to argue by contradiction that the set of points ∩ ∞ j=1 S j missed by the cover is limited by the same bound [C/ε] as the number of discarded balls of each radius 1 j . (For full details, refer to [12].) A crucial difference between the above proof and that of the analogous result in [36] (namely Proposition 3.3) is that in our version the underlying manifold M need not be compact. In [36], compactness is used in two ways: first to guarantee that any cover of M by balls of radius 1 j can be assumed finite, and second to argue that the centers of balls B on which B f (i) i dx exceeds ε must be fixed, using the fact that a sequence of centers of such balls must have a convergent subsequence. Because we first pass to a subsequence of {f i } on which balls can be consistently labeled as good or bad, we avoid the need to argue after diagonalizing that bad balls stabilize. (c) A is gauge equivalent to a connection in L 2 Our proof of part (a) proceeds by pulling back a conformal mapping f : For part (b), we construct the appropriate gauge transformation s as the composition σ • f −1 of the maps σ and f −1 from part (a). An estimate analogous to that in part (a) using the transformation law for 1-forms yields the desired bound on the decay ofÃ ≡ A s .
To complete the proof of part (c), we simply compute derivatives of the transformation law to show componentwise that |∂A s /∂y α | ≤ K ′ |y| −2 for y ∈ V. We then define a simple extension of s to all of R − × R 3 by extrapolation into R − × R 3 \ V of its values on ∂V, so that A s is now a global gauge representative for A with decay of the connection and its componentwise partial derivatives bounded by a constant multiple of |y| −2 . Since i * A obeys the same decay condition and since the volume element on {0} × R 3 has weight |y| 2 , . Note that as a corollary, any smooth connection A on {0} × R 3 which serves as boundary data for a Yang-Mills connection approaches pure gauge at spatial infinity. This follows simply because the decay condition implies that the curvature of A approaches zero at spatial infinity, and hence A must approach pure gauge (although it need not approach the same value in all spatial directions).
Almost-everywhere Gâteaux differentiability of S. To show that the functional S[·] is Gâteaux differentiable almost everywhere, we employ the Banach space version of Rademacher's theorem derived in [43], Theorem 6. This requires us to view S [·] as a map from an open subset U of a separable real Banach space X to a real Banach space Y satisfying the Radon-Nikodym property (i.e. every function of bounded variation from [0, 1] to Y must be differentiable almost everywhere), and to show that as such it is locally Lipschitz. Thus our first task is to show that with suitable interpretation of its domain and range, S [·] is locally Lipschitz. The Banach space version of Rademacher's theorem then states that such a locally Lipschitz map T from U ⊂ X to Y is Gâteaux differentiable outside a Gaussian null subset of U. A Gaussian null set in a separable Banach space is a Borel subset which has measure zero under every nondegenerate Gaussian measure (for details of defining Gaussian measures on Banach spaces, see [43]).
To achieve Banach space structure on the domain B of S [·], we view the local (see [42]).
Given a fixed reference connectionÃ on {0}×R 3 , we interpret S [·] as a functional defined on the Banach space We prove that thus considered, S [·] is locally Lipschitz: ) be a reference connection such thatF ≡ FÃ satisfies {0}×R 3F 2 dx ′ < ∞, and such that there exists an extensionÃ ∈ L 2 1;loc (R − × R 3 , Ad η⊗T * (R − × R 3 )) ofÃ into R − × R 3 with finite Euclidean signature Yang-Mills action. Then for all sections ω ∈ L 2 is well defined and locally Lipschitz with respect to the norm · 2,1 on We first verify that S [Ã + ω] is well defined, namely thatÃ + ω has an extension A into R − × R 3 with finite Euclidean signature Yang-Mills action. To construct A , we represent the extensionÃ ofÃ into R − × R 3 in a choice of gauge such that i * Ã =Ã (as opposed to merely i * Ã ∼Ã). We then define A piecewise over the time component of R − × R 3 as a linear interpolation betweenÃ + ω andÃ, followed byÃ : . Straightforward computation of the curvature of A followed by Sobolev estimates confirms that S [Ã + ω] is well defined.
where the last equality is satisfied because A ∈ B * . Arguments from Hodge theory imply that the solution to (III.76) is unique and that the map A → A L , A ∈ B * is smooth. Thus we define the smooth map where A 2 1;loc (0) * is the dual of A 2 1;loc (0) and D * A L +h F A L +h is interpreted in a weak sense as where in the last line we denote A L 0 + h 0 by A 0 and F A L 0 +h 0 by F 0 . In order to apply the Banach space implicit function theorem and conclude (almost everywhere) smoothness of S [·], we must show that except for A 0 belonging to a set of measure zero in B * (with respect to a suitable measure), the kernel of the map D h E (A 0 , h 0 ) is trivial (modulo gauge transformations). To this end, we assume that the absolute minimizer of C. An intended application to the mass gap problem A fundamental question in quantum gauge theory is whether the Schrödinger operator for certain non-abelian Yang-Mills fields admits a spectral gap. Such a gap, if it exists, could represent the energy difference between the actual vacuum state and that of the lowest energy, excited 'glueball' states and confirm the expectation that massless gluons (the fundamental quanta of Yang-Mills dynamics) cannot propagate freely as photons do but must instead exhibit a form of 'color confinement'. It seems to be well understood that this question lies beyond the scope of conventional perturbation theory and will require a more global analytical treatment for its ultimate resolution.
Many years ago I. M. Singer proposed an elegant, geometrical approach to this fundamen-tal problem based on the fact that the classical, reduced configuration space for (Lorentzian signature) Yang-Mills dynamics -namely the 'orbit space' of spatial connections modulo gauge transformations -has a naturally induced, curved Riemannian metric with everywhere non-negative sectional curvature [31]. The classical Hamiltonian for the reduced dynamics -a real-valued functional defined on the cotangent bundle of this orbit spaceconsists of a 'kinetic' term induced from the spatial integral of the square of the vectorial electric component of the full, spacetime Yang-Mills curvature tensor and a 'potential' term induced from the spatial integral of the square of its complementary, vectorial magnetic component. The non-vanishing curvature of the Riemannian metric defined by the kinetic term arises from the implementation of the Gauss-law constraint during the process of reduction to the quotient, orbit space and was independently computed by several investigators [31,32,45]. The classical reduced dynamics is thus that for a system point (namely a gauge equivalence class of spatial connections) moving on a positively curved, infinite dimensional manifold under the influence of a (non-negative) potential energy.
Upon canonical quantization the Schrödinger operator for this (pure Yang-Mills) dynamical system will thus include a kinetic term that, formally at least, encompasses the (negative 10 ) Laplace-Beltrami operator for an infinite-dimensional, curved Riemannian manifold -namely the orbit space alluded to above. Whereas the (covariant) Hessian of sufficiently smooth (wave) functionals can be rigorously defined in such infinite-dimensional contexts, its associated trace need not make sense without some suitable regularization since the Hessian will not, in general, be trace class. Singer, in particular proposed an elegant zeta-function regularization scheme to define the needed Laplacian [31].
A classical result in Riemannian geometry due to A. Lichnerowicz [46] shows that the Laplace operator for a complete, connected (finite-dimensional) Riemannian manifold necessarily exhibits a spectral gap provided that the Ricci tensor of this manifold is bounded, positively, away from zero 11 . Such a result however cannot be expected to extend in any straightforward way at least, to the infinite-dimensional manifolds arising in quantum Yang-Mills theory. First of all, as Singer pointed out, their Ricci tensors, which would result from taking traces of corresponding (rigorously computable) curvature tensors, are not in general 10 We here adopt the usual physicists' sign convention for the definition of a Laplacian. 11 It follows from the Bonnet-Myers theorem that such a manifold is necessarily compact [47].
well-defined -the curvature tensors in question not being trace-class -and would require a suitable regularization for their meaningful formulation. Again Singer proposed a zeta function regularization scheme as an elegant means of accomplishing this. Some such regularization, however, is actually a desirable feature of the quantization procedure, at least in 4 spacetime dimensions, since it allows the introduction of a length scale into the quantum formalism. In the absence of such a scale no hypothetical spectral energy gap could even be expressed in terms of the naturally occurring parameters of the theory (Planck's constant, the speed of light and the Yang-Mills coupling constant).
Another difficulty with attempting to extend the Lichnerowicz argument to the infinitedimensional setting of interest here is that, thanks to the Bonnet-Myers theorem, one knows that a complete, finite dimensional Riemannian manifolds with positive Ricci curvature bounded away from zero is necessarily compact [47]. For a connected such manifold the lowest eigenvalue of its associated (negative) Laplacian always vanishes and corresponds to a globally constant eigenfunction. That such an eigenfunction is nevertheless always normalizable follows from the manifold's compactness. The spectral gap referred to in Lichnerowicz's theorem is thus simply the lowest non-vanishing eigenvalue of the manifold's (negative) Laplacian which, in view of compactness, necessarily has a discrete spectrum.
Generalizations of Lichnerowicz's theorem have been established under less stringent conditions on the Ricci tensors provided that the manifolds under study have finite diameters [48,49]. L. Andersson has proven that Riemannian Hilbert manifolds have finite diameters whenever their full sectional curvatures are positively bounded away from zero [50] but this result does not apply to the orbit space sectional curvatures of interest here since these latter admit (infinite dimensional) families of 2-planes on which they actually vanish. In any case the diameters of these Yang-Mills orbit spaces are known to be infinite [51].
The true, normalizable ground state wave functional must necessarily reflect the presence of the potential energy term in the Schrödinger operator. In a recent paper [52] the authors showed how to modify the original Lichnerowicz argument (in a finite dimensional setting) to allow for the occurrence of such a potential energy term and show that a corresponding gap estimate follows therefrom provided that a suitably defined 'Bakry-Emery Ricci tensor' is bounded positively away from zero. This Bakry-Emery Ricci tensor differs from the actual Ricci tensor by a term in the (covariant) Hessian of the logarithm of the true ground state wave function. Its positivity could hold on a flat or even negatively curved space and thus its applicability is not limited to manifolds of finite diameter.
Furthermore the natural integration measure arising in this (generalized Lichnerowicz) analysis includes the squared modulus of the ground state wave function itself so that the total space, even if it has infinite diameter, now has finite measure simply by virtue of the normalizability of the vacuum state. This should prove to be especially significant for any potential extensions to infinite dimensional problems wherein formal Lebesgue measures no longer make sense but for which normalizable vacuum state wave functional are nevertheless expected to exist.
Our Euclidean signature semi-classical program, when applied to Yang-Mills fields, has the significant advantage over conventional, Rayleigh-Schrödinger perturbation theory of keeping the non-linearities and non-abelian gauge invariances fully intact at every level of the analysis. Our expectation is that it should yield an asymptotic expansion for the needed, fully gauge invariant, logarithm of the ground state wave functional that is far superior to any attainable by conventional perturbation methods. The latter, by requiring an expansion in the Yang-Mills coupling constant, disturb both the nonlinear structure and the closely associated (non-abelian) gauge invariance of the Yang-Mills dynamical system at the outset and attempt to reinstate those vital features only gradually, order-by-order in the expansion.
Though our main focus in Ref. [52] was on the Yang-Mills system we proved therein that (non-vanishing) orbit space curvature also arises naturally through the (minimal) coupling of a Maxwell field to a charged scalar field. In this case curvature arises only for the scalar factor of the (product) orbit space and not for the Maxwell factor which remains flat. We were thus led to conjecture that orbit space curvature could even serve as an independent source of mass for matter fields themselves provided that they are (minimally) coupled to (abelian or non-abelian) gauge fields. In the standard model of elementary particle physics such masses arise exclusively (in view of the necessity of maintaining gauge invariance) through the interaction of matter fields with the Higgs field. It is intriguing to speculate whether orbit space curvature effects could supply an alternative source of such masses. into the action, the resulting reduced kinetic term (a quadratic form in the 'velocity' of the spatial connection) is found to be degenerate along gauge orbit directions but smooth, gauge invariant and positive definite in the transversal directions [31,32,45]. It thus follows that this kinetic term defines a smooth, Riemannian metric on the natural 'orbit space' of spatial connections modulo gauge transformations. This orbit space is (at least almost everywhere) itself a smooth, infinite dimensional manifold and provides the geometrically natural (reduced) configuration space for (classical) Yang-Mills dynamics.
A corresponding smooth potential energy function is induced on this orbit space by the integral over R 3 (at fixed t) of the square of the curvature of the spatial connection fieldthe 'magnetic' component of the curvature of the full spacetime connection field. A Legendre transformation leads in turn to the Hamiltonian functional for the classical dynamics which takes the 'standard' form of a sum of (curved space) kinetic and potential energies defined on the cotangent bundle of the aforementioned orbit space.
The sectional curvature of this reduced configuration space was independently computed in [31,32,45] and shown to be everywhere non-negative but almost everywhere non-vanishing whenever the gauge group is non-abelian. Though Singer discussed the need for a suitable regularization scheme to make sense of the formally (positively) divergent Ricci tensor of the orbit space metric, the actual form of such a regulated Ricci tensor seems still to be unknown. It would be most interesting if a suitably defined Ricci tensor could be shown to be bounded, positively away from zero on this orbit space, especially inasmuch as we think it quite unlikely that our proposed Bakry-Emery 'enhancement' of this tensor would nullify its (hypothetical) positivity properties but perhaps, more likely, complement them 12 .
Furthermore, as we discussed near the end of Section II of Ref. [52], it is quite plausible that strict positivity of the Bakry-Emery Ricci tensor, though seemingly sufficient for the implication of a spectral gap, is not absolutely necessary for this conclusion to hold.

D. Euclidean-Signature Asymptotic Methods and the Wheeler-DeWitt Equation
Globally hyperbolic spacetimes, { (4) V, (4) g}, are definable over manifolds with the product structure, (4) V ≈ M ×R. We shall focus here on the 'cosmological' case for which the spatial factor M is a compact, connected, orientable 3-manifold without boundary. The Lorentzian metric, (4) g, of such a spacetime is expressible, relative to a time function x 0 = t, in the 3+1-dimensional form (4) wherein, for each fixed t, the Riemannian metric is the first fundamental form induced by (4) g on the corresponding t = constant, spacelike hypersurface. The unit, future pointing, timelike normal field to the chosen slicing (defined by the level surfaces of t) is expressible in terms of the (strictly positive) 'lapse' function N and 'shift vector' field Y i ∂ ∂x i as or, in covariant form, as The canonical spacetime volume element of (4) g, µ(4) g := − det (4) g, takes the 3+1dimensional form where µ γ := √ det γ is the volume element of γ.
In view of the compactness of M the Hilbert and ADM action functionals, evaluated on domains of the product form, Ω = M × I, with I = [t 0 , t 1 ] ⊂ R, simplify somewhat to and tr γ K := γ ij K ij (III.84) designate the second fundamental form and mean curvature induced by (4) g on the constant t slices. In these formulas spatial coordinate indices, i, j, k, . . . , are raised and lowered with γ and the vertical bar, '|', signifies covariant differentiation with respect to this metric so that, for example, Y i|j = ∇ j (γ)γ iℓ Y ℓ . When the variations of (4) g are appropriately restricted, the boundary term distinguishing I Hilbert from I ADM makes no contribution to the field equations and so can be discarded. Writing with Lagrangian density one defines the momentum conjugate to γ via the Legendre transformation so that p = p ij ∂ ∂x i ⊗ ∂ ∂x j is a symmetric tensor density induced on each t = constant slice. In terms of the variables {γ ij , p ij , N, Y i } the ADM action takes the Hamiltonian form Insofar as a consistent factor ordering for the Hamiltonian constraint operator, H ⊥ (γ, i δ δγ ), also exists, one will be motivated to propose the (Euclidean-signature, semiclassical) ansatz The essential question that now comes to light is thus the following: Is there a well-defined mathematical method for establishing the existence of a In view of the field theoretic examples discussed above one's first thought might be to seek to minimize an appropriate Euclidean-signature action functional subject to suitable boundary and asymptotic conditions. But, as is well-known from the Euclidean-signature path integral program [58], the natural functional to use for this purpose is unbounded from below within any given conformal class -one can make the functional arbitrarily large and negative by deforming any metric (4) g Eucl with a suitable conformal factor [58,59].
But the real point of the constructions above was not to minimize action functionals but rather to generate certain 'fundamental sets' of solutions to the associated Euler-Lagrange equations upon which the relevant action functionals could then be evaluated. But the Einstein equations, in vacuum or even allowing for the coupling to conformally invariant matter sources, encompass, as a special case, the vanishing of the 4-dimensional scalar curvature, (4) R( (4) g Eucl ). Thus there is no essential loss in generality, and indeed a partial simplification of the task at hand to be gained, by first restricting the relevant, Euclideansignature action functional to the 'manifold' of Riemannian metrics satisfying (in the vacuum case) (4) R( (4) g Eucl ) = 0 and then seeking to carry out a constrained minimization of this functional.
Setting (4) R( (4) g Eucl ) = 0 freezes out the conformal degree of freedom that caused such consternation for the Euclidean path integral program [58,59], wherein one felt obligated to integrate over all possible Riemannian metrics having the prescribed boundary behavior, but is perfectly natural in the present context and opens the door to appealing to the positive action theorem which asserts that the relevant functional is indeed positive when evaluated on arbitrary, asymptotically Euclidean metrics that satisfy (4) R( (4) g Eucl ) ≥ 0 [60][61][62][63].
Another complication of the Euclidean path integral program was the apparent necessity to invert, by some still obscure means, something in the nature of a 'Wick rotation' that had presumably been exploited to justify integrating over Riemannian, as opposed to Lorentziansignature, metrics. Without this last step the formal 'propagator' being constructed would presumably be that for the Euclidean-signature variant of the Wheeler-DeWitt equation and not the actual Lorentzian-signature version that one wishes to solve. In ordinary quantum mechanics the corresponding, well-understood step is needed to convert the Feynman-Kac propagator, derivable by rigorous path-integral methods, back to one for the actual Schrödinger equation.
But in the present setting no such hypothetical 'Wick rotation' would ever have been performed in the first place so there is none to invert. Our focus throughout is on constructing asymptotic solutions to the original, Lorentz-signature Wheeler-DeWitt equation and not to its Euclidean-signature counterpart. That a Euclidean-signature Einstein-Hamilton-Jacobi equation emerges in this approach has the very distinct advantage of leading one to specific problems in Riemannian geometry that may well be resolvable by established mathematical methods. By contrast, path integral methods, even for the significantly more accessible gauge theories discussed above, would seem to require innovative new advances in measure theory for their rigorous implementation. Even the simpler scalar field theories, when formulated in the most interesting case of four spacetime dimensions, seem still to defy realization by path integral means. It is conceivable, as was suggested in the concluding section of [1], that focusing predominantly on path integral methods to provide a 'royal road' to quantization may, inadvertently, render some problems more difficult to solve rather than actually facilitating their resolution.
The well-known 'instanton' solutions to the Euclidean-signature Yang-Mills equations present a certain complication for the semi-classical program that we are advocating in that they allow one to establish the existence of non-unique minimizers for the Yang-Mills action functional for certain special choices of boundary data [12]. This in turn can obstruct the global smoothness of the corresponding solution to the Euclidean-signature Hamilton-Jacobi equation. While it is conceivable that the resulting, apparent need to repair the associated 'scars' in the semi-classical wave functionals may have non-perturbative implications for the Yang-Mills energy spectrum -of potential relevance to the 'mass-gap' problem -no such corrections to the spectrum are expected or desired for the gravitational case. Thus it is reassuring to note that analogous 'gravitational instanton' solutions to the Euclideansignature Einstein equations have been proven not to exist [59].
We conclude by noting that other interesting, generally covariant systems of field equations exist to which our ('Euclidean-signature semi-classical') quantization methods could also be applied. Classical relativistic 'membranes', for example, can be viewed as the evo- A solution to the operator ordering problem for these quantized constraints, when the ambient spacetime is Minkowskian, was proposed by one of us in [64]. For the compact, codimension one case (i.e., when M is compact and k = 1) it is not difficult to show that the relevant Euclidean-signature Hamilton-Jacobi equation has a fundamental solution given by the volume functional of the maximal, spacelike hypersurface that uniquely spans, a la Plateau, the arbitrarily chosen embedding [65]. It would be especially interesting to see whether higher-order quantum corrections and excited state wave functionals can be computed for this system in a way that realizes a quantum analogue of general covariance. The reader will not have failed to note the very strong sense in which our work builds on the fundamental contributions of Karen Uhlenbeck.