Functional determinant

For a positive selfadjoint operator S on a finite-dimensional Euclidean space V, the formula

${\displaystyle {\frac {1}{\sqrt {\det S}}}=\int _{V}e^{-\pi \langle x,Sx\rangle }\,dx}$

holds.

The problem is to find a way to make sense of the determinant of an operator S on an infinite dimensional function space. One approach, favored in quantum field theory, in which the function space consists of continuous paths on a closed interval, is to formally attempt to calculate the integral

${\displaystyle \int _{V}e^{-\pi \langle \phi ,S\phi \rangle }\,{\mathcal {D}}\phi }$

where V is the function space and ${\displaystyle \langle -,-\rangle }$ the L² inner product, and ${\displaystyle {\mathcal {D}}\phi }$ the Wiener measure. The basic assumption on S is that it should be selfadjoint, and have discrete spectrum λ₁, λ₂, λ₃… with a corresponding set of eigenfunctions f₁, f₂, f₃… which are complete in L² (as would, for example, be the case for the second derivative operator on a compact interval Ω). This roughly means all functions φ can be written as linear combinations of the functions f_i:

${\displaystyle |\phi \rangle =\sum _{i}c_{i}|f_{i}\rangle \quad {\text{with }}c_{i}=\langle f_{i}|\phi \rangle .\,}$

Hence the inner product in the exponential can be written as

${\displaystyle \langle \phi |S|\phi \rangle =\sum _{i,j}c_{i}^{*}c_{j}\langle f_{i}|S|f_{j}\rangle =\sum _{i,j}c_{i}^{*}c_{j}\delta _{ij}\lambda _{i}=\sum _{i}|c_{i}|^{2}\lambda _{i}.}$

In the basis of the functions f_i, the functional integration reduces to an integration over all basisfunctions. Formally, assuming our intuition from the finite dimensional case carries over into the infinite dimensional setting, the measure should then be equal to

${\displaystyle {\mathcal {D}}\phi =\prod _{i}{\frac {dc_{i}}{2\pi }}.}$

This makes the functional integral a product of Gaussian integrals:

${\displaystyle \int _{V}{\mathcal {D}}\phi \;e^{-\langle \phi |S|\phi \rangle }=\prod _{i}\int _{-\infty }^{+\infty }{\frac {dc_{i}}{2\pi }}e^{-\lambda _{i}c_{i}^{2}}.}$

The integrals can then be evaluated, giving

${\displaystyle \int _{V}{\mathcal {D}}\phi \;e^{-\langle \phi |S|\phi \rangle }=\prod _{i}{\frac {1}{2{\sqrt {\pi \lambda _{i}}}}}={\frac {N}{\sqrt {\prod _{i}\lambda _{i}}}}}$

where N is an infinite constant that needs to be dealt with by some regularization procedure. The product of all eigenvalues is equal to the determinant for finite-dimensional spaces, and we formally define this to be the case in our infinite-dimensional case also. This results in the formula

${\displaystyle \int _{V}{\mathcal {D}}\phi \;e^{-\langle \phi |S|\phi \rangle }\propto {\frac {1}{\sqrt {\det S}}}.}$

If all quantities converge in an appropriate sense, then the functional determinant can be described as a classical limit (Watson and Whittaker). Otherwise, it is necessary to perform some kind of regularization. The most popular of which for computing functional determinants is the zeta function regularization.[1] For instance, this allows for the computation of the determinant of the Laplace and Dirac operators on a Riemannian manifold, using the Minakshisundaram–Pleijel zeta function. Otherwise, it is also possible to consider the quotient of two determinants, making the divergent constants cancel.

Let S be an elliptic differential operator with smooth coefficients which is positive on functions of compact support. That is, there exists a constant c > 0 such that

${\displaystyle \langle \phi ,S\phi \rangle \geq c\langle \phi ,\phi \rangle }$

for all compactly supported smooth functions φ. Then S has a self-adjoint extension to an operator on L² with lower bound c. The eigenvalues of S can be arranged in a sequence

${\displaystyle 0<\lambda _{1}\leq \lambda _{2}\leq \cdots ,\qquad \lambda _{n}\to \infty .}$

Then the zeta function of S is defined by the series:[2]

${\displaystyle \zeta _{S}(s)=\sum _{n=1}^{\infty }{\frac {1}{\lambda _{n}^{s}}}.}$

It is known that ζ_S has a meromorphic extension to the entire plane.[3] Moreover, although one can define the zeta function in more general situations, the zeta function of an elliptic differential operator (or pseudodifferential operator) is regular at ${\displaystyle s=0}$.

Formally, differentiating this series term-by-term gives

${\displaystyle \zeta _{S}'(s)=\sum _{n=1}^{\infty }{\frac {-\log \lambda _{n}}{\lambda _{n}^{s}}},}$

and so if the functional determinant is well-defined, then it should be given by

${\displaystyle \det S=\exp \left(-\zeta _{S}'(0)\right).}$

Since the analytic continuation of the zeta function is regular at zero, this can be rigorously adopted as a definition of the determinant.

This kind of Zeta-regularized functional determinant also appears when evaluating sums of the form ${\displaystyle \sum _{n=0}^{\infty }{\frac {1}{(n+a)}}}$, integration over 'a' gives ${\displaystyle \sum _{n=0}^{\infty }\log(n+a)}$ which it just can be considered as the logarithm of the determinant for a Harmonic oscillator this last value is just equal to ${\displaystyle -\partial _{s}\zeta _{H}(0,a)}$, where ${\displaystyle \zeta _{H}(s,a)}$ is the Hurwitz Zeta function.

We will compute the determinant of the following operator describing the motion of a quantum mechanical particle in an infinite potential well:

${\displaystyle \det \left(-{\frac {d^{2}}{dx^{2}}}+A\right)\qquad (x\in [0,L]),}$

where A is the depth of the potential and L is the length of the well. We will compute this determinant by diagonalizing the operator and multiplying the eigenvalues. So as not to have to bother with the uninteresting divergent constant, we will compute the quotient between the determinants of the operator with depth A and the operator with depth A = 0. The eigenvalues of this potential are equal to

${\displaystyle \lambda _{n}={\frac {n^{2}\pi ^{2}}{L^{2}}}+A\qquad (n\in \mathbb {N} \setminus \{0\}).}$

This means that

${\displaystyle {\frac {\det \left(-{\frac {d^{2}}{dx^{2}}}+A\right)}{\det \left(-{\frac {d^{2}}{dx^{2}}}\right)}}=\prod _{n=1}^{+\infty }{\frac {{\frac {n^{2}\pi ^{2}}{L^{2}}}+A}{\frac {n^{2}\pi ^{2}}{L^{2}}}}=\prod _{n=1}^{+\infty }\left(1+{\frac {L^{2}A}{n^{2}\pi ^{2}}}\right).}$

Now we can use Euler's infinite product representation for the sine function:

${\displaystyle \sin z=z\prod _{n=1}^{\infty }\left(1-{\frac {z^{2}}{n^{2}\pi ^{2}}}\right)}$

from which a similar formula for the hyperbolic sine function can be derived:

${\displaystyle \sinh z=-i\sin iz=z\prod _{n=1}^{\infty }\left(1+{\frac {z^{2}}{n^{2}\pi ^{2}}}\right).}$

Applying this, we find that

${\displaystyle {\frac {\det \left(-{\frac {d^{2}}{dx^{2}}}+A\right)}{\det \left(-{\frac {d^{2}}{dx^{2}}}\right)}}=\prod _{n=1}^{+\infty }\left(1+{\frac {L^{2}A}{n^{2}\pi ^{2}}}\right)={\frac {\sinh L{\sqrt {A}}}{L{\sqrt {A}}}}.}$

For one-dimensional potentials, a short-cut yielding the functional determinant exists.[4] It is based on consideration of the following expression:

${\displaystyle {\frac {\det \left(-{\frac {d^{2}}{dx^{2}}}+V_{1}(x)-m\right)}{\det \left(-{\frac {d^{2}}{dx^{2}}}+V_{2}(x)-m\right)}}}$

where m is a complex constant. This expression is a meromorphic function of m, having zeros when m equals an eigenvalue of the operator with potential V₁(x) and a pole when m is an eigenvalue of the operator with potential V₂(x). We now consider the functions ψ^m₁ and ψ^m₂ with

${\displaystyle \left(-{\frac {d^{2}}{dx^{2}}}+V_{i}(x)-m\right)\psi _{i}^{m}(x)=0}$

obeying the boundary conditions

${\displaystyle \psi _{i}^{m}(0)=0,\quad \qquad {\frac {d\psi _{i}^{m}}{dx}}(0)=1.}$

If we construct the function

${\displaystyle \Delta (m)={\frac {\psi _{1}^{m}(L)}{\psi _{2}^{m}(L)}},}$

which is also a meromorphic function of m, we see that it has exactly the same poles and zeroes as the quotient of determinants we are trying to compute: if m is an eigenvalue of the operator number one, then ψ^m₁(x) will be an eigenfunction thereof, meaning ψ^m₁(L) = 0; and analogously for the denominator. By Liouville's theorem, two meromorphic functions with the same zeros and poles must be proportional to one another. In our case, the proportionality constant turns out to be one, and we get

${\displaystyle {\frac {\det \left(-{\frac {d^{2}}{dx^{2}}}+V_{1}(x)-m\right)}{\det \left(-{\frac {d^{2}}{dx^{2}}}+V_{2}(x)-m\right)}}={\frac {\psi _{1}^{m}(L)}{\psi _{2}^{m}(L)}}}$

for all values of m. For m = 0 we get

${\displaystyle {\frac {\det \left(-{\frac {d^{2}}{dx^{2}}}+V_{1}(x)\right)}{\det \left(-{\frac {d^{2}}{dx^{2}}}+V_{2}(x)\right)}}={\frac {\psi _{1}^{0}(L)}{\psi _{2}^{0}(L)}}.}$

The problem in the previous section can be solved more easily with this formalism. The functions ψ⁰_i(x) obey

${\displaystyle {\begin{aligned}&\left(-{\frac {d^{2}}{dx^{2}}}+A\right)\psi _{1}^{0}=0,\qquad \psi _{1}^{0}(0)=0\quad ,\qquad {\frac {d\psi _{1}^{0}}{dx}}(0)=1,\\&-{\frac {d^{2}}{dx^{2}}}\psi _{2}^{0}=0,\qquad \psi _{2}^{0}(0)=0,\qquad {\frac {d\psi _{2}^{0}}{dx}}(0)=1,\end{aligned}}}$

yielding the following solutions:

${\displaystyle {\begin{aligned}&\psi _{1}^{0}(x)={\frac {1}{\sqrt {A}}}\sinh x{\sqrt {A}},\\&\psi _{2}^{0}(x)=x.\end{aligned}}}$

This gives the final expression

${\displaystyle {\frac {\det \left(-{\frac {d^{2}}{dx^{2}}}+A\right)}{\det \left(-{\frac {d^{2}}{dx^{2}}}\right)}}={\frac {\sinh L{\sqrt {A}}}{L{\sqrt {A}}}}.}$