Abstract

We study the action of the free quantum Hamiltonian in the Schrödinger representation on functional spaces generated by Gaussian states with boundary conditions set by localized external sources. It turns out that inspite of very general initial assumptions the proposed scheme requires at least two sources at a distance limited by non-zero constant of the model. In this way we provide example of quntum field evolved by the action of the free Hamiltonian confining external sources which set boundary conditions.

Keywords: extensions of closed semi-bounded quadratic forms, square root of operator, quantum field theory in the Schrödinger representation, Gaussian functional.

Introduction

The Hamilton operator of the free scalar quantum field in the Schrödinger representation is the sum of the second variational derivative and quadratic form of the Laplace operator [1]

$$\mathscr{H}=\int_{\mathbb{R}^{3}}d^{3}x\,\bigl{(}-\frac{\delta^{2}}{\delta A(\boldsymbol{x})^{2}}+(\partial_{m}A)^{2}\bigr{)}.$$ (1)

This operator acts on the set of two times differentiable functionals defined on the domain of its potential term. In order to construct the ground state of $\mathscr{H}$ one needs to deploy spectral representation of the quadratic form by fixing the corresponding self-adjoint operator

$$L:\quad A(\boldsymbol{x})\to-\frac{\partial^{2}A(\boldsymbol{x})}{\partial x_{m}^{2}}$$

and the scalar product

$$(f(\boldsymbol{x}),g(\boldsymbol{x}))=\int_{\mathbb{R}^{3}}d^{3}x\,\overline{f(\boldsymbol{x})}g(\boldsymbol{x}),$$ (2)

so that

$$\mathfrak{l}(A)=\int_{\mathbb{R}^{3}}(\partial_{m}A(\boldsymbol{x}))^{2}\,d^{3}x=(A,LA),\quad A\in\mathcal{D}(L).$$

Then the functional of the ground state of (1) is given by the Gaussian of the quadratic form $\omega(A)$ corresponding to the square root of $L$

$$\Omega_{0}(A)=\exp\{-\frac{1}{2}\omega(A)\}=\exp\{-\frac{1}{2}(A,L^{1/2}A)\}.$$

Due to the equation

$$\frac{1}{4}\int_{\mathbb{R}^{3}}d^{3}x\,\bigl{(}\frac{\delta}{\delta A(\boldsymbol{x})}\omega(A)\bigr{)}^{2}=(L^{1/2}A,L^{1/2}A)=\mathfrak{l}(A),\quad A\in\mathcal{D}(\mathfrak{l}),$$ (3)

the action of $\mathscr{H}$ on the functional $\Omega_{0}(A)$ is reduced to multiplication by the scalar $\operatorname{Tr}L^{1/2}$

$$\mathscr{H}\,\Omega_{0}(A)=\int_{\mathbb{R}^{3}}d^{3}x\,\bigl{(}-\frac{\delta^{2}}{\delta A(\boldsymbol{x})^{2}}+(\partial_{m}A)^{2}\bigr{)}\,\exp\{-\frac{1}{2}\omega(A)\}=\\ =\frac{1}{2}\int_{\mathbb{R}^{3}}d^{3}x\,\bigl{(}\frac{\delta^{2}}{\delta A(\boldsymbol{x})^{2}}\omega(A)\bigr{)}\,\exp\{-\frac{1}{2}\omega(A)\}=(\operatorname{Tr}L^{1/2})\,\Omega_{0}(A).$$

One may note that the flat integral of the product of variational derivatives in the first term of $\mathscr{H}$ leads to the flat scalar product in the quadratic form $\omega(A)$ and so there is no other option in the definition of (2).

Hamiltonian $\mathscr{H}$ has singularities outside of the domain of the quadratic form $\mathcal{D}(\mathfrak{l})$, that is in the region where potential term is infinite. There is no natural measure on the space of functions to describe $\mathcal{D}(\mathfrak{l})$ as an open set which does not include boundary, but we may say that functions with behaviour

$$A(\boldsymbol{x})=\frac{C_{n}}{|\boldsymbol{x}-\boldsymbol{x}_{n}|^{1/2}},\quad\boldsymbol{x}\to\boldsymbol{x}_{n},\quad n=1,\ldots N$$ (4)

reside very close to $\mathcal{D}(\mathfrak{l})$. Indeed, for any $\epsilon>0$, if $A(\boldsymbol{x})$ expands around $\boldsymbol{x}_{n}$ as

$$A(\boldsymbol{x})=\frac{C_{n}}{|\boldsymbol{x}-\boldsymbol{x}_{n}|^{1/2-\epsilon}},\quad\boldsymbol{x}\to\boldsymbol{x}_{n},$$

then the integrand in $\mathfrak{l}(A)$ has singularity of the order

$$(\partial_{m}A(\boldsymbol{x}))^{2}=\frac{C_{n}^{2}}{|\boldsymbol{x}-\boldsymbol{x}_{n}|^{3-2\epsilon}},\quad\boldsymbol{x}\to\boldsymbol{x}_{n},$$

and so is locally integrable in $\mathbb{R}^{3}$. Meanwhile in the case $\epsilon=0$ the integral in the quadratic form diverges logarithmically around $\boldsymbol{x}_{n}$ and the Hamiltonian $\mathscr{H}$ acquires singularity in the vicinity of functions with behaviour (4). As long as the domain of $\omega$ is wider than $\mathcal{D}(\mathfrak{l})$, the functional $\Omega_{0}(A)$ in the vicinity of (4) is well defined and has finite boundary values.

Now let us assume that there exists positive extension $\mathfrak{l}_{\kappa}(A)$ of the quadratic form $\mathfrak{l}(A)$ to some function with the behaviour (4)

$$\mathfrak{l}_{\kappa}(A) =\mathfrak{l}(A),\quad A\in\mathcal{D}(\mathfrak{l}),$$ (5)

$$\mathfrak{l}_{\kappa}(B) <\infty,\quad B\stackrel{{\scriptstyle\boldsymbol{x}\to\boldsymbol{x}_{n}}}{{=}}\frac{C_{n}}{|\boldsymbol{x}-\boldsymbol{x}_{n}|^{1/2}}.$$ (6)

Provided that $\mathfrak{l}_{\kappa}$ is closed w. r. to the scalar product (2) one may construct corresponding self-adjoint operator $L_{\kappa}$ and define its spectral decomposition (on the extensions of closed semi-bounded quadratic forms see [2]). Then one may substitute $\omega(A)$ with the form of the square root of $L_{\kappa}$ and observe that the Gaussian

$$\Omega_{\kappa}(A)=\exp\{-\frac{1}{2}\omega_{\kappa}(A)\}=\exp\{-\frac{1}{2}(A,L_{\kappa}^{1/2}A)\},$$ (7)

also satisfies the eigenstate equation for $\mathscr{H}$

$$\mathscr{H}\,\Omega_{\kappa}(A)=\frac{1}{2}\int_{\mathbb{R}^{3}}d^{3}x\,\bigl{(}\frac{\delta^{2}}{\delta A(\boldsymbol{x})^{2}}\omega_{\kappa}(A)\bigr{)}\,\exp\{-\frac{1}{2}\omega_{\kappa}(A)\}=(\operatorname{Tr}L_{\kappa}^{1/2})\,\Omega_{\kappa}(A),$$ (8)

provided that the following equation holds

$$\frac{1}{4}\int_{\mathbb{R}^{3}}d^{3}x\,\bigl{(}\frac{\delta}{\delta A(\boldsymbol{x})}\omega_{\kappa}(A)\bigr{)}^{2}=(L_{\kappa}^{1/2}A,L_{\kappa}^{1/2}A)=\mathfrak{l}_{\kappa}(A)=\mathfrak{l}(A),\quad A\in\mathcal{D}(\mathfrak{l}).$$ (9)

In order to verify equations (8) and (9) one needs to fix the set w. r. to which the variational derivative is operating. The form $\omega(A)$ in (3) is defined as Fourier preimage of quadratic form of multiplication operator and can be reconstructed by variation w. r. to the functions from the domain of the quadratic form $\mathfrak{l}(A)$. This condition can be written as

$$\overline{\omega|_{\mathcal{D}(\mathfrak{l})}}=\omega,$$

where LHS is the closure of restriction of the form $\omega(A)$ to the set $\mathcal{D}(\mathfrak{l})$. That is for any $A\in\mathcal{D}(\omega)$ there exists a sequence $A_{n}\in\mathcal{D}(\mathfrak{l})$ such that

$$||A_{n}-A||_{L_{2}(\mathbb{R}^{3})}\stackrel{{\scriptstyle n\to\infty}}{{\to}}0,\quad\omega(A_{n}-A)\stackrel{{\scriptstyle n\to\infty}}{{\to}}0.$$

Thus the same set $\mathcal{D}(\mathfrak{l})$ should be used for the variation in the eigenstate equation (8), that is we need to verify that the above closure also holds for the quadratic form $\omega_{\kappa}(A)$

$$\overline{\omega_{\kappa}|_{\mathcal{D}(\mathfrak{l})}}=\omega_{\kappa}.$$ (10)

As we will see in what follows, proof of this equation requires some effort and is not possible for every closable extension $\mathfrak{l}_{\kappa}$ of the quadratic form of the Laplace operator.

Another essential requirement of the presented construction is the positivity of the form $\mathfrak{l}_{\kappa}$. Otherwise the operator $L_{\kappa}$ possesses negative eigenspace and the quadratic form $\omega_{\kappa}$ goes imaginary thereon. This in turn contradicts the self-adjointness of $\mathscr{H}$ and does not allow to introduce a scalar product on the set of functionals with finite normalization for $\Omega_{\kappa}$.

The construction of solutions to the equation of the ground state of the free quantum field by means of the extensions of quadratic form in the potential term of the Hamiltonian was proposed in [3]. There extensions generated by $\delta$-potential were used, that is extensions to functions with the behaviour

$$A(\boldsymbol{x})=\frac{C_{n}}{|\boldsymbol{x}-\boldsymbol{x}_{n}|},\quad\boldsymbol{x}\to\boldsymbol{x_{n}}.$$

However no proof of the possibility to reconstruct the quadratic form of the square root by its action on $\mathcal{D}(\mathfrak{l})$ for that case is still not provided.

In this work we develop the scheme of construction of eigenfunctionals of the free Hamiltonian by extending the quadratic form of Laplacian to functions with behaviour (4). We will show that positive extensions require that behaviour at least in 2 points and will derive the admissible distance between the latter as the function of the extension parameter. Then we will verify that extensions in question satisfy equation (10), discuss normalization of $\Omega_{\kappa}$ and possibility to introduce eigenvalues in (8) with continuous dependence on the position of sources. That is the model can be interpreted as dynamical system, which includes quantum field evolved by the action of the free Hamiltonian and external sources – subject of the further quantization. Thus we expect that presented construction can be employed in description of renormalized asymptotically free quatum field with self-interaction and interaction with external sources.

1 Extensions of the quadratic form

Let us turn to detailed description of extensions of the quadratic form of the scalar Laplace operator oulined in (5), (6). We will work in the momentum representation (in the space of Fourier transforms) of fields

$$\hat{A}(\boldsymbol{p})=\frac{1}{(2\pi)^{3/2}}\int_{\mathbb{R}^{3}}A(\boldsymbol{x})e^{-i\boldsymbol{p}\cdot{\boldsymbol{x}}}\,d^{3}x$$

and for simplicity will omit the Fourier hat in the notation of $A(\boldsymbol{p})$ and other functions. Operator $\mathscr{H}$ in this representation acts in the following way

$$\mathscr{H}\,\Omega_{0}(A)=\int_{\mathbb{R}^{3}}d^{3}p\,\bigl{(}-\frac{\delta^{2}}{\delta A(\boldsymbol{p})^{2}}+p^{2}A^{2}(\boldsymbol{p})\bigr{)}\,\Omega_{0}(A).$$

Quadratic forms $\mathfrak{l}(A)$ and $\omega(A)$ are diagonalized

$$\mathfrak{l}(A) =\int_{\mathbb{R}^{3}}p^{2}|A(\boldsymbol{p})|^{2}\,d^{3}p, L: \,A(\boldsymbol{p})\to p^{2}A(\boldsymbol{p}),$$

$$\omega(A) =\int_{\mathbb{R}^{3}}p|A(\boldsymbol{p})|^{2}\,d^{3}p, L^{1/2}: \,A(\boldsymbol{p})\to pA(\boldsymbol{p}),$$

while the scalar product connecting quadratic forms with their operators is still flat

$$(f(\boldsymbol{p}),g(\boldsymbol{p}))=\int_{\mathbb{R}^{3}}d^{3}p\,\overline{f(\boldsymbol{p})}g(\boldsymbol{p}).$$ (11)

Key objects for construction of extensions of the quadratic form $\mathfrak{l}(A)$ are the resolvent of the self-adjoint operator $L$

$$R_{\mu}=(L-\mu)^{-1}:\quad A(\boldsymbol{p})\to\frac{1}{p^{2}-\mu}A(\boldsymbol{p}),\quad\mu\in\mathbb{C}\setminus\mathbb{R}^{+}\cup 0$$

and singular potential $v(\boldsymbol{p})$ – some generalized function obeying conditions

$$R_{\mu}v\in L_{2}(\mathbb{R}^{3})\text{ or }(R_{\mu}v,R_{\nu}v)<\infty$$ (12)

$$(v,R(\mu)v)=\int_{\mathbb{R}^{3}}\frac{|v(p)|^{2}}{p^{2}-\mu}d^{3}p=\infty.$$ (13)

As candidates for $v(\boldsymbol{p})$ we will consider powers of modulus of $\boldsymbol{p}$

$$v_{\alpha}(\boldsymbol{p})=p^{\alpha}.$$

Such potentials have distinctive singularities in the coordinate representation

$$\frac{1}{(2\pi)^{3/2}}\int_{\mathbb{R}^{3}}v_{\alpha}(\boldsymbol{p})e^{i\boldsymbol{p}\cdot\boldsymbol{x}}\,d^{3}p=\begin{cases}(2\pi)^{3/2}\delta(\boldsymbol{x}),&\alpha=0,\\ \frac{2^{3/2+\alpha}\Gamma(\frac{3+\alpha}{2})}{\Gamma(-\frac{\alpha}{2})}\frac{1}{x^{3+\alpha}},&\alpha\neq 0,\end{cases}$$

so it makes sense to include into consideration their shifts in coordinate space

$$v_{\alpha,\boldsymbol{x}}(\boldsymbol{p})=e^{i\boldsymbol{p}\cdot\boldsymbol{x}}v_{\alpha}(\boldsymbol{p}).$$

The Birman-Krein-Vishik theory [4] describes the domain $\mathcal{D}(\mathfrak{l}_{\kappa})$ of extensions of the quadratic form $\mathfrak{l}$ as the sum of $\mathcal{D}(\mathfrak{l})$ and linear span of the vector $R_{\mu}v$

$$\mathcal{D}(\mathfrak{l}_{\kappa})=\mathcal{D}(\mathfrak{l})\dotplus\{\xi R_{\mu}v\},\quad\mu\in\mathbb{C}\setminus\mathbb{R}^{+}\cup 0,\quad\xi\in\mathbb{C}.$$

Note that particular choice of $\mu$ is irrelevant here since

$$R_{\mu_{1}}v-R_{\mu_{2}}v=(\mu_{1}-\mu_{2})R_{\mu_{1}}R_{\mu_{2}}v\in\mathcal{D}(L),\quad\mathcal{D}(L)\subset\mathcal{D}(\mathfrak{l}).$$

As shown in the work [5] if the potential $v(\boldsymbol{p})$ fails to satisfy (13) then the vector $R_{\mu}v$ belongs to the domain $\mathcal{D}(\mathfrak{l})$, that is $\mathfrak{l}_{\kappa}$ changes the action of $\mathfrak{l}$ and so changes the action of $\mathscr{H}$ which is not our goal. In the terms of $\alpha$ the divergence of the integral in (13) is expressed in the inequality

$$-\frac{1}{2}\leq\alpha.$$

It is not hard to see that condition (12) imposes restriction on $\alpha$ from the other side, namely the strict inequality

$$\alpha<\frac{1}{2}.$$

Domain of the self-adjoint operator $L_{\kappa}$ corresponding to the quadratic form $\mathfrak{l}_{\kappa}$ includes linear combinations of vectors $R_{\rho}v$, $R_{\rho}^{2}v$ of the following type

$$R_{\rho}v+\xi(\kappa,\rho)R_{\rho}^{2}v,\quad\rho<0,\quad\xi\in\mathbb{C},$$

with some $\xi$ defined by the boundary conditions of $L_{\kappa}$. But we will not study this domain, neither we will not study the action of $\mathfrak{l}_{\kappa}$, the only useful object for the construction of functions of $L_{\kappa}$ is its resolvent given by the Krein formula

$$R_{\mu}^{\kappa}=(L_{\kappa}-\mu)^{-1}=R_{\mu}+\frac{R_{\mu}v(R_{\bar{\mu}}v,\,\cdot\,)}{\kappa-\sigma(\mu)}.$$ (14)

Here $\sigma(\mu)$ is the Nevanlinna (Herglotz) function [6], defined by the equation

$$\sigma(\mu)-\sigma(\nu)=(\mu-\nu)(R_{\bar{\mu}}v,R_{\nu}v)=(v,(R_{\mu}-R_{\nu})v)$$ (15)

and $\kappa$ is the real separation constant of the particular solution which serves in (14) as the extension parameter. Note that in the case $\kappa=\infty$ the second term in the RHS of (14) vanishes and we arrive at the resolvent $R_{\mu}$ of $L$

$$R_{\mu}^{\infty}=R_{\mu}.$$

In what follows we will restrict our consideration to potentials

$$v(\boldsymbol{p})=v_{-1/2}(\boldsymbol{p})=\frac{1}{p^{1/2}},$$

such that vectors $R_{\mu}v$ reside in the region described by (4). Indeed, $R_{\mu}v$ in coordinate representation is given by the integral

$$(R_{\mu}v)(\boldsymbol{x})=\frac{1}{(2\pi)^{3/2}}\int_{\mathbb{R}^{3}}\frac{e^{i\boldsymbol{p}\cdot\boldsymbol{x}}\,d^{3}p}{(p^{2}-\mu)p^{1/2}}=\frac{4\pi}{(2\pi)^{3/2}x}\int_{0}^{\infty}\frac{\sin px}{p^{2}-\mu}p^{1/2}\,dp,$$

which after the change of variable $s=px$ turns into

$$(R_{\mu}v)(\boldsymbol{x})=\frac{4\pi}{(2\pi)^{3/2}x^{1/2}}\int_{0}^{\infty}\frac{\sin s}{s^{2}-\mu x^{2}}s^{1/2}\,ds.$$ (16)

The integral by $s$ converges absolutely, it is bounded for $x\to 0$ and it is not equal to zero in this limit. Then one may conclude that singularity of $(R_{\mu}v)(\boldsymbol{x})$ is defined by coefficient $x^{-1/2}$ and so is described by (4).

The Nevanlinna function for $v(\boldsymbol{p})=p^{-1/2}$ is given by the following equation

$$\sigma(\mu)-\sigma(\lambda)=\int_{\mathbb{R}^{3}}\bigl{(}\frac{1}{p^{2}-\mu}-\frac{1}{p^{2}-\lambda}\bigr{)}\frac{d^{3}p}{p}=2\pi(\ln(-\lambda)-\ln(-\mu)).$$ (17)

Here we choose the main branch of the logarithm with the cut along the negative semi-axis in such a way that RHS is defined for all complex $\mu$ and $\lambda$ outside of the positive semi-axis and zero. This yields that denominator in (14) can be chosen as

$$\kappa-\sigma(\mu)=2\pi\ln(-\frac{\mu}{\tilde{\kappa}}),\quad\tilde{\kappa}>0,$$ (18)

where $\tilde{\kappa}$ is the positive value combining dimensional constant to correctly define the logarithm and the extension parameter $\kappa$.

It is evident that RHS of (18) decreases when $\mu$ goes from $-\infty$ to zero and changes its sign for any value of $\tilde{\kappa}$. Hence resolvent (14) has a pole in the point $\mu=-\tilde{\kappa}$ and corresponding self-adjoint operator and quadratic form possess negative eigenspace. This contradicts the requirement of positivity of the extension $\mathfrak{l}_{\kappa}$ and shows that nontrivial solution of the eigenstate equation of $\mathscr{H}$ produced by a single localized external source does not exist.

It is worth to note that logarithmic Nevanlinna function (17) and analysis of extensions of the quadratic form coincide with the corresponding objects of the theory of singular perturbations of 2-dimensional Laplace operator by the $\delta$-potential described in [7].

2 Interaction with two sources

In order to fix obstacles encountered in the previous section let us turn to a more complex singular potential. In the original work [8] Krein equation was derived for the interaction with multiple sources $v_{m}(\boldsymbol{p})$

$$R^{\kappa}_{\mu}=R_{\mu}+R_{\mu}v_{m}(\kappa-\sigma(\mu))_{ml}^{-1}(R_{\bar{\mu}}v_{l},\,\cdot\,),$$ (19)

where $\kappa$ is a real symmetric matrix and the matrix $\sigma$ is defined by the equation

$$\sigma_{jk}(\mu)-\sigma_{jk}(\nu)=(v_{j},(R_{\mu}-R_{\nu})v_{k}).$$ (20)

In our case it is natural to add to the resolvent the potential

$$v_{\boldsymbol{x}}(\boldsymbol{p})=e^{i\boldsymbol{p}\cdot\boldsymbol{x}}p^{-1/2}$$

shifted by the vector $\boldsymbol{x}$ in the coordinate space. In order to avoid the treatment of zeros of the matrix $\kappa_{jk}-\sigma_{jk}(\lambda)$ one may turn to the basis

$$v_{-}(\boldsymbol{p}) =\frac{1}{\sqrt{2}}(v(\boldsymbol{p})-v_{\boldsymbol{x}}(\boldsymbol{p}))=\frac{1}{\sqrt{2}}p^{-1/2}(1-e^{i\boldsymbol{p}\cdot\boldsymbol{x}}),$$

$$v_{+}(\boldsymbol{p}) =\frac{1}{\sqrt{2}}(v(\boldsymbol{p})+v_{\boldsymbol{x}}(\boldsymbol{p}))=\frac{1}{\sqrt{2}}p^{-1/2}(1+e^{i\boldsymbol{p}\cdot\boldsymbol{x}})$$

and restrict to Krein resolvents with nontrivial second term generated by the antisymmetric potential $v_{-}(\boldsymbol{p})$

$$R_{\mu}^{\kappa}=R_{\mu}+\frac{R_{\mu}v_{-}(R_{\bar{\mu}}v_{-},\,\cdot\,)}{\kappa-\sigma(\mu)}.$$ (21)

In this case by substituting $v_{-}(\boldsymbol{p})$ into (15) we have the following equation for $\sigma(\mu)$

$$\sigma(\mu)-\sigma(\lambda)=\int_{\mathbb{R}^{3}}\frac{1-\cos(\boldsymbol{p}\cdot\boldsymbol{x})}{p}\bigl{(}\frac{1}{p^{2}-\mu}-\frac{1}{p^{2}-\lambda}\bigr{)}\,d^{3}p=\\ =4\pi\int_{0}^{\infty}\bigl{(}p-\frac{\sin px}{x}\bigr{)}\bigl{(}\frac{1}{p^{2}-\mu}-\frac{1}{p^{2}-\lambda}\bigr{)}dp.$$ (22)

Similarly to (18) let us introduce the dimensional parameter $\tilde{\kappa}$ which includes $\kappa$ and define denominator in (21) by means of the integral

$$\sigma(\mu)-\kappa=4\pi\int_{0}^{\infty}\bigl{(}\frac{p}{p^{2}-\mu}-\frac{\sin px}{x(p^{2}-\mu)}-\frac{p}{p^{2}+\tilde{\kappa}}\bigr{)}dp=\\ =2\pi\ln\frac{\tilde{\kappa}}{-\mu}-4\pi\int_{0}^{\infty}\frac{\sin px}{x}\frac{dp}{p^{2}-\mu}.$$ (23)

This expression is valid for all $\mu\in\mathbb{C}\setminus\mathbb{R}^{+}\cup 0$ and satisfies equation (22).

Let us study the RHS of (23) on the negative semi-axis. Its second term is a function of the product $\mu x^{2}$ and for $\mu\to-\infty$ it has the following expansion

$$\mathrm{J}(\mu x^{2})=4\pi\int_{0}^{\infty}\frac{\sin px}{x}\frac{dp}{p^{2}-\mu}\stackrel{{\scriptstyle\mu x^{2}\to-\infty}}{{=}}-\frac{4\pi}{\mu x^{2}}+\frac{8\pi}{\mu^{2}x^{4}}+\mathcal{O}(\frac{1}{\mu^{3}x^{6}}).$$ (24)

Thus the whole sum (23) goes logarithmically to $-\infty$

$$\sigma(\mu)-\kappa=2\pi\ln\frac{\tilde{\kappa}}{-\mu}+\frac{4\pi}{\mu x^{2}}+\mathcal{O}(\mu^{-2})\to-\infty,\quad\mu\to-\infty.$$

Derivative of the function $\sigma(\mu)$ for $\mu<0$ is strictly positive

$$\sigma^{\prime}(\mu)=4\pi\int_{0}^{\infty}\bigl{(}p-\frac{\sin px}{x}\bigr{)}\frac{dp}{(p^{2}-\mu)^{2}}>0,$$ (25)

which means that $\sigma(\mu)$ monotonously increases with the increase of $\mu$. In order to calculate the value $\sigma(0)$ let us write the integral $\mathrm{J}(\mu x^{2})$ for finite $\mu$, $x$ in the terms of the integral sine and cosine

$$4\pi\int_{0}^{\infty}\frac{\sin px}{x}\frac{dp}{p^{2}-\mu}=\frac{4\pi}{\sqrt{\mu}x} \bigl{(}\cos\sqrt{\mu}x\,\operatorname{Si}\sqrt{\mu}x-\sin\sqrt{\mu}x\,\operatorname{Ci}\sqrt{\mu}x\bigr{)}=$$ (26)

$$\stackrel{{\scriptstyle\mu x^{2}\to 0}}{{=}}4\pi(1-\boldsymbol{\gamma}-\frac{1}{2}\ln(-\mu x^{2})+\mathcal{O}(\mu x^{2})).$$ (27)

Here $\operatorname{Si}(\sqrt{\mu}x)$ and $\operatorname{Ci}(\sqrt{\mu}x)$ with non-conventional integration constant are defined as

$$\operatorname{Si}(\sqrt{\mu}x)=\int_{0}^{\sqrt{\mu}x}\frac{\sin t}{t}dt,$$ (28)

$$\operatorname{Ci}(\sqrt{\mu}x)=\boldsymbol{\gamma}+\frac{1}{2}\ln(-\mu x^{2})-\int_{0}^{\sqrt{\mu}x}\frac{1-\cos t}{t}dt,$$ (29)

and $\boldsymbol{\gamma}$ is the Euler constant. One may note that RHS in (26) includes only Taylor-expandable even functions and the logarithm. That is branching by $\mu$ in (26) is defined solely by the logarithm from the integral cosine (29). Substituting (27) into (23) one derives the following expansion for $\sigma(\mu)-\kappa$

$$\sigma(\mu)-\kappa=2\pi\ln\tilde{\kappa}x^{2}+4\pi(\boldsymbol{\gamma}-1)+\ln(-\mu x^{2})\mathcal{O}(\mu x^{2}),\quad\mu\to 0.$$ (30)

This equation shows that we can always choose $\tilde{\kappa}$ in such a way that RHS of (30) is negative when $\mu$ tends to zero. The latter means, due to the positivity of derivative (25), that $\sigma(\mu)-\kappa$ does not turn to zero on the negative semi-axis and that for such $\tilde{\kappa}$ resolvent (21) defines positive operator $L_{\kappa}$ and thus real-valued quadratic form $\omega_{\kappa}(A)$.

Figure 1: Functions $\sigma(\mu,x)$, $\sigma(\mu,x_{\mathrm{b}})$ compared to $2\pi\ln\frac{\tilde{\kappa}}{-\mu}$

Let us look at (30) from another side. Suppose that for chosen values of $\tilde{\kappa}$ and $x$ the limit of RHS in (30) is negative and we can correctly define the functional $\Omega_{\kappa}(A)$. When increasing the distance $x$ the limit of RHS of (30) is increasing as well and for some boundary $x=x_{\mathrm{b}}$ it turns to zero (see Fig. 1)

$$\sigma(\mu,x_{\mathrm{b}})|_{\mu=0}-\kappa=0.$$

For $x>x_{\mathrm{b}}$ the function $\sigma(\mu,x)-\kappa$ changes its sign through zero in the interval $-\infty<\mu<0$ and the corresponding self-adjoint operator $L_{\kappa}$ possesses negative eigenvalue. That is for any chosen extension parameter $\tilde{\kappa}$ admissible relative positions of sources always exist and are always limited by the sphere with radius proportional to inverse square root of $\tilde{\kappa}$.

3 Electromagnetic field

Quantum Hamiltonian of the free electromagnetic field in the Coulomb gauge acts as the following operator [1]

$$\mathscr{H}=\int_{\mathbb{R}^{3}}d^{3}p\,\bigl{(}-\frac{\delta}{\delta A_{l}(\boldsymbol{p})^{2}}P_{lj}^{\bot}\frac{\delta}{\delta A_{j}(\boldsymbol{p})^{2}}+p^{2}A_{l}^{2}(\boldsymbol{p})\bigr{)},$$ (31)

where $P_{lj}^{\bot}$ is the projector to the transverse (solenoidal) component

$$P_{lj}^{\bot}=\delta_{lj}-\frac{p_{l}p_{j}}{p^{2}}.$$

Corresponding functional of the ground state differs from its scalae counterpart only by summation in the quadratic form $\omega(A)$

$$\Omega_{0}(A)=\exp\{-\frac{1}{2}\omega(A)\},\quad\omega(A)=\int_{\mathbb{R}^{3}}pA_{l}^{2}(\boldsymbol{p})\,d^{3}p,$$

provided that $A_{l}(p)$ is transverse

$$p_{l}A_{l}(\boldsymbol{p})=0.$$

The latter condition ensures that linear functional (1-form)

$$\frac{\delta\omega(A)}{\delta A_{l}(\boldsymbol{p})}$$

is transverse w. r. to variables $(\boldsymbol{p},l)$ and that the projector $P_{lj}^{\bot}$ in the middle of (31) acts as unity when $\mathscr{H}$ is applied to $\Omega_{0}(A)$.

In order to construct closed positive extensions of the quadratic form

$$\mathfrak{l}(A)=\int_{\mathbb{R}^{3}}p^{2}A_{l}^{2}(\boldsymbol{p})\,d^{3}p$$

from the second term of (31) let us consider the following set (in the index $j$) of transverse (in the index $l$) vector singular potentials

$$v_{l,\boldsymbol{x}_{n}}^{j}(\boldsymbol{p})=\frac{e^{i\boldsymbol{p}\cdot\boldsymbol{x}_{n}}}{p^{1/2}}\bigl{(}\delta_{lj}-\frac{p_{l}p_{j}}{p^{2}}\bigr{)}=\frac{e^{i\boldsymbol{p}\cdot\boldsymbol{x}_{n}}}{p^{1/2}}P_{lj}^{\bot}.$$

There are 3 potentials $j=1,2,3$ attached to each point which, according to [9], generate $6\times 6$ matrix of the denominator of resolvent (19)

$$\sigma_{nj}^{mk}(\mu)=\begin{pmatrix}\frac{4\pi}{3}\delta_{jk}\ln\frac{\tilde{\kappa}}{-\mu}&\delta_{jk}(\mathrm{J}-2\dot{\sigma})-\frac{x_{j}x_{k}}{x^{2}}(\mathrm{J}-6\dot{\sigma})\\ \delta_{jk}(\mathrm{J}-2\dot{\sigma})-\frac{x_{j}x_{k}}{x^{2}}(\mathrm{J}-6\dot{\sigma})&\frac{4\pi}{3}\delta_{jk}\ln\frac{\tilde{\kappa}}{-\mu}\end{pmatrix},$$ (32)

where $m,n=1,2$; $x_{j}$ and $x_{k}$ are $j$-th and $k$-th coordinates of the vector $\boldsymbol{x}=\boldsymbol{x}_{2}-\boldsymbol{x}_{1}$ and we use the notation $\dot{\sigma}$ for dimensionless derivative of (23)

$$\dot{\sigma}(\mu x^{2})=\frac{1}{x^{2}}\frac{\partial}{\partial\mu}\sigma(\mu,x)=-\frac{2\pi}{\mu x^{2}}-\mathrm{J}^{\prime}(\mu x^{2}).$$

Matrix $\sigma_{nj}^{mk}(\mu)$ satisfies equation (20) with off-diagonal separation constants set to zero

$$\sigma_{nj}^{mk}(\mu)-\sigma_{nj}^{mk}(\nu)=\int_{\mathbb{R}^{3}}\overline{v_{l,\boldsymbol{x}_{n}}^{j}(\boldsymbol{p})}\bigl{(}\frac{1}{p^{2}-\mu}-\frac{1}{p^{2}-\lambda}\bigr{)}v_{l,\boldsymbol{x}_{m}}^{k}(\boldsymbol{p})\,d^{3}p.$$

Provided that all diagonal constants $\tilde{\kappa}$ are equal, this matrix can be diagonalized simultaneously for all $\mu$ by changing the basis to the set of normalized vectors proportional to

$$\begin{pmatrix}\boldsymbol{x}\\ \boldsymbol{x}\end{pmatrix},\quad\begin{pmatrix}\boldsymbol{x}_{1}^{\bot}\\ \boldsymbol{x}_{1}^{\bot}\end{pmatrix},\quad\begin{pmatrix}\boldsymbol{x}_{2}^{\bot}\\ \boldsymbol{x}_{2}^{\bot}\end{pmatrix},\quad\begin{pmatrix}\boldsymbol{x}\\ -\boldsymbol{x}\end{pmatrix},\quad\begin{pmatrix}\boldsymbol{x}_{1}^{\bot}\\ -\boldsymbol{x}_{1}^{\bot}\end{pmatrix},\quad\begin{pmatrix}\boldsymbol{x}_{2}^{\bot}\\ -\boldsymbol{x}_{2}^{\bot}\end{pmatrix},$$ (33)

where 3-dimensional vectors $\boldsymbol{x}_{1}^{\bot}$ and $\boldsymbol{x}_{2}^{\bot}$ are orthogonal to $\boldsymbol{x}$ and to each other. Indeed, vector $\boldsymbol{x}$ is eigenvector for all elements of (32), vectors $\boldsymbol{x}_{1,2}^{\bot}$ are eigenvectors for diagonal elements and for the first terms in off-diagonal ones, but they are annihilated by the second terms. Taking the quadratic form of (32) on normalized vectors (33) we have the following Nevanlinna functions

$$\sigma_{\shortparallel}^{+} =\frac{4\pi}{3}\ln\frac{\tilde{\kappa}}{-\mu}+4\dot{\sigma}(\mu x^{2}), \sigma_{\bot_{1,2}}^{+} =\frac{4\pi}{3}\ln\frac{\tilde{\kappa}}{-\mu}+\mathrm{J}(\mu x^{2})-2\dot{\sigma}(\mu x^{2}),$$ (34)

$$\sigma_{\shortparallel}^{-} =\frac{4\pi}{3}\ln\frac{\tilde{\kappa}}{-\mu}-4\dot{\sigma}(\mu x^{2}), \sigma_{\bot_{1,2}}^{-} =\frac{4\pi}{3}\ln\frac{\tilde{\kappa}}{-\mu}-\mathrm{J}(\mu x^{2})+2\dot{\sigma}(\mu x^{2}).$$ (35)

Calculation of higher terms in the expansion (27) gives us the expansion of $\dot{\sigma}(\mu x^{2})$ for small argument

$$\dot{\sigma}(\mu x^{2})=-2\pi(\frac{1}{6}\ln(-\mu x^{2})+\frac{\boldsymbol{\gamma}}{6}-\frac{2}{9}+\ln(-\mu x^{2})\mathcal{O}(\mu x^{2})).$$

As long as for $\mu\to-\infty$ the value of $\dot{\sigma}(\mu x^{2})$ tends to zero, one may conclude that out of 6 functions (34), (35) only antisymmetric ones, namely $\sigma_{\shortparallel}^{-}$ and 2 functions $\sigma_{\bot}^{-}$ are upper-bounded for $-\infty<\mu<0$ and become zeroless for corresponding values of $\tilde{\kappa}$.

We would like to note that matrix similar to (32) corresponding to 3 or more localized potentials can hardly be simultaneously diagonalized for all $\mu$ for general set of distances even in the scalar case. Hence separation of the subspace of positive extensions of the form $\mathfrak{l}(A)$ may become an untrivial problem for multiple-potential interaction.

4 Square root of the operator $L_{\kappa}$

The positive square root of the self-adjoined operator $L_{\kappa}$ can be written in the terms of the resolvent $R_{\lambda}^{\kappa}$ by means of the spectral decomposition in the following way

$$L_{\kappa}^{1/2}=\lim_{\epsilon\to 0}\frac{1}{2\pi i}\int_{0}^{\infty}(R_{\lambda+i\epsilon}^{\kappa}-R_{\lambda-i\epsilon}^{\kappa})\lambda^{1/2}\,d\lambda,\quad\lambda^{1/2}>0.$$

Substituting the expression for $R_{\lambda}^{\kappa}$ and using the notation $\pm i0$ for the limit in $\epsilon$ we can write the formal representation

$$L_{\kappa}^{1/2}=\frac{1}{2\pi i}\int_{0}^{\infty}\bigl{(}R_{\lambda+i0}-R_{\lambda-i0}+\\ +\frac{R_{\lambda+i0}v(R_{\lambda-i0}v,\,\cdot\,)}{\kappa-\sigma(\lambda+i0)}-\frac{R_{\lambda-i0}v(R_{\lambda+i0}v,\,\cdot\,)}{\kappa-\sigma(\lambda-i0)}\bigr{)}\lambda^{1/2}\,d\lambda=\\ =p+\frac{1}{2\pi i}\int_{0}^{\infty}\bigl{(}\frac{R_{\lambda+i0}v(R_{\lambda-i0}v,\,\cdot\,)}{\kappa-\sigma(\lambda+i0)}-\frac{R_{\lambda-i0}v(R_{\lambda+i0}v,\,\cdot\,)}{\kappa-\sigma(\lambda-i0)}\bigr{)}\lambda^{1/2}\,d\lambda,$$ (36)

where $p$ is the operation of multiplication by the modulus of momentum variable $\boldsymbol{p}$. In particluar

$$L^{1/2}=p,$$ (37)

but we cannot substitute this in (36) as long as $L_{\kappa}^{1/2}$ and $L^{1/2}$ are possibly defined on different domains. Nevertheless the expression (36), in general case with some regularization, can be used as the kernel of the quadratic form of $L_{\kappa}^{1/2}$. For this purpose we will also need the gap of the function $\sigma(\lambda)$ across the positive semi-axis, which in the scalar case is equal to

$$\sigma(\lambda+i0)-\sigma(\lambda-i0)=2\pi\ln\frac{\tilde{\kappa}}{-\lambda-i0}-2\pi\ln\frac{\tilde{\kappa}}{-\lambda+i0}-\\ -4\pi\int_{0}^{\infty}\frac{\sin px}{x}(\frac{1}{p^{2}-\lambda-i0}-\frac{1}{p^{2}-\lambda+i0})\,dp=4\pi^{2}i-4\pi^{2}i\frac{\sin\sqrt{\lambda}x}{\sqrt{\lambda}x}.$$ (38)

And for extensions of the quadratic form of Coulomb-gauged vector Laplacian these gaps are

$$\sigma_{\shortparallel}^{-}(\lambda+i0)-\sigma_{\shortparallel}^{-}(\lambda-i0) =\frac{8\pi^{2}i}{3}+8\pi^{2}i\bigl{(}\frac{\cos\sqrt{\lambda}x}{\lambda x^{2}}-\frac{\sin\sqrt{\lambda}x}{(\sqrt{\lambda}x)^{3}}\bigr{)},$$ (39)

$$\sigma_{\bot}^{-}(\lambda+i0)-\sigma_{\bot}^{-}(\lambda-i0) =\frac{8\pi^{2}i}{3}-4\pi^{2}i\bigl{(}\frac{\sin\sqrt{\lambda}x}{\sqrt{\lambda}x}+\frac{\cos\sqrt{\lambda}x}{\lambda x^{2}}-\frac{\sin\sqrt{\lambda}x}{(\sqrt{\lambda}x)^{3}}\bigr{)}.$$ (40)

Integrals similar to (36) can also be used to calculate other functions of $L_{\kappa}$, for example the logarithm.

5 Reconstruction of the quadratic form of $L_{\kappa}^{1/2}$

As was stated in the Introduction in order to construct eigenfunctionals of operation $\mathscr{H}$ one need to prove that the positive quadratic form of $L_{\kappa}^{1/2}$ can be reconstructed by its action on the set $\mathcal{D}(\mathfrak{l})$ everywhere dense in the space $L_{2}(\mathbb{R}^{3})$. Let the function $A(\boldsymbol{p})$ to be from $\mathcal{D}(\omega_{\kappa})$ then there exists a sequence $A_{n}\in\mathcal{D}(L_{\kappa}^{1/2})$ converging to $A$ in the norm induced by the form $\omega_{\kappa}$

$$||A-A_{n}||_{L_{2}(\mathbb{R}^{3})}\to 0,\quad\omega_{\kappa}(A-A_{n})\to 0,\quad n\to\infty.$$

Thus it is sufficient to prove that any function $A_{n}(p)$ from $\mathcal{D}(L_{\kappa}^{1/2})$ can be approximated in the norm of $\omega_{\kappa}$ by functions from $\mathcal{D}(\mathfrak{l})$. If $A_{n}$ belongs to $\mathcal{D}(L_{\kappa}^{1/2})$ then it belongs to the domain of quadratic form $\mathfrak{l}_{\kappa}$ of the operator $L_{\kappa}$

$$\mathfrak{l}_{\kappa}(A_{n})=(L_{\kappa}^{1/2}A_{n},L_{\kappa}^{1/2}A_{n})<\infty.$$

As long as $\mathcal{D}(\mathfrak{l}_{\kappa})$ is the direct sum of $\mathcal{D}(\mathfrak{l})$ and the linear span of some vector $R_{\rho}v$, $\rho<0$, then this $A_{n}$ can be uniquely decomposed into the sum

$$A_{n}=a_{n}+\chi_{n}R_{\rho}v,\quad a_{n}\in\mathcal{D}(\mathfrak{l}),\quad\chi_{n}\in\mathbb{C},\quad\rho<0.$$

Thus it is sufficient to show that function $R_{\rho}v$ can be approximated by functions from $\mathcal{D}(\mathfrak{l})$ in the norm of quadratic form $\omega_{\kappa}$.

In order to verify the above proposition let us use expression (36) for the operator $L_{\kappa}^{1/2}$ and write down its quadratic form

$$\omega_{\kappa}(R_{\rho}v)=\omega(R_{\rho}v)+\frac{1}{2\pi i}\int_{0}^{\infty}\bigl{(}\frac{(R_{\rho}v,R_{\lambda+i0}v)(R_{\lambda-i0}v,R_{\rho}v)}{\kappa-\sigma(\lambda+i0)}-\\ -\frac{(R_{\rho}v,R_{\lambda-i0}v)(R_{\lambda+i0}v,R_{\rho}v)}{\kappa-\sigma(\lambda-i0)}\bigr{)}\lambda^{1/2}\,d\lambda.$$ (41)

Now using (15) each quotient here can be transformed into

$$\frac{(R_{\rho}v,R_{\lambda\pm i0}v)(R_{\lambda\mp i0}v,R_{\rho}v)}{\kappa-\sigma(\lambda\pm i0)}=\frac{(\sigma(\lambda\pm i0)-\sigma(\rho))^{2}}{(\lambda\pm i0-\rho)^{2}(\kappa-\sigma(\lambda\pm i0))}=\\ =\frac{2\sigma(\rho)-\kappa-\sigma(\lambda\pm i0)}{(\lambda\pm i0-\rho)^{2}}+\frac{(\kappa-\sigma(\rho))^{2}}{(\lambda\pm i0-\rho)^{2}(\kappa-\sigma(\lambda\pm i0))},$$

and finally, up to some vanishing terms proportional to $i0$, we have

$$\omega_{\kappa}(R_{\rho}v)=4\pi\int_{0}^{\infty}\frac{p^{2}\,dp}{(p^{2}-\rho)^{2}}-\frac{1}{2\pi i}\int_{0}^{\infty}\frac{\sigma(\lambda+i0)-\sigma(\lambda-i0)}{(\lambda-\rho)^{2}}\lambda^{1/2}d\lambda\\ +\frac{1}{2\pi i}(\kappa-\sigma(\rho))^{2}\int_{0}^{\infty}\frac{\sigma(\lambda+i0)-\sigma(\lambda-i0)}{(\lambda-\rho)^{2}(\kappa-\sigma(\lambda+i0))(\kappa-\sigma(\lambda-i0))}\lambda^{1/2}d\lambda.$$ (42)

The modulus of the gap of $\sigma(\lambda)$ calculated in (38) is limited, thus all three integrals in (42) converge absolutely. Then in order to approximate $R_{\rho}v$ by the function $(\Lambda+\rho)R_{-\Lambda}R_{\rho}v\in\mathcal{D}(\mathfrak{l})$ one may use the following limit in the norm $L_{2}(\mathbb{R}^{3})$

$$R_{\rho}v-(\Lambda+\rho)R_{-\Lambda}R_{\rho}v=R_{-\Lambda}v=\frac{1}{(p^{2}+\Lambda)p^{1/2}}\stackrel{{\scriptstyle L_{2}(\mathbb{R}^{3})}}{{\to}}0,\quad\Lambda\to\infty.$$

The same limit is valid in the norm of the form $\omega$

$$\omega(R_{-\Lambda}v)=\int_{\mathbb{R}_{3}}\frac{d^{3}p}{(p^{2}+\Lambda)^{2}}\to 0,\quad\Lambda\to\infty,$$

thus in order to prove that

$$\omega_{\kappa}(R_{-\Lambda}v)\to 0,\quad\Lambda\to\infty,$$

one needs to show that the second and the third terms in (42), after the substitution of $\rho$ with $-\Lambda$, tend to zero for $\Lambda\to\infty$. By the change of variable $\tilde{\lambda}=\lambda/\Lambda$ those terms can be rewritten as

$$-\frac{1}{2\pi i\Lambda^{1/2}}\int_{0}^{\infty}\frac{\sigma(\tilde{\lambda}\Lambda+i0)-\sigma(\tilde{\lambda}\Lambda-i0)}{(\tilde{\lambda}+1)^{2}}\tilde{\lambda}^{1/2}d\tilde{\lambda}$$ (43)

$$\frac{(\kappa-\sigma(-\Lambda))^{2}}{2\pi i\Lambda^{1/2}}\int_{0}^{\infty}\frac{\sigma(\tilde{\lambda}\Lambda+i0)-\sigma(\tilde{\lambda}\Lambda-i0)}{(\tilde{\lambda}+1)^{2}(\kappa-\sigma(\tilde{\lambda}\Lambda+i0))(\kappa-\sigma(\tilde{\lambda}\Lambda-i0))}\tilde{\lambda}^{1/2}d\tilde{\lambda},$$ (44)

correspondingly. Here both integrals converge absolutely and due to the limited gap of $\sigma(\lambda)$ their values are bounded from above for any $\Lambda$. Taking into account coefficients $\Lambda^{-1/2}$ one may conclude that (43), (44) and the whole quadratic form $\omega_{\kappa}(R_{-\Lambda}v)$ tend to zero for $\Lambda\to\infty$.

Summarizing the above we have shown that the quadratic form of the operator $L_{\kappa}^{1/2}$ can be reconstructed by its action on the set $\mathcal{D}(\mathfrak{l})$. The same proof is valid for quadratic forms of square roots of positive extenstions of the transverse Laplace operator, one just needs to use corresponding potential $v$ and equations (39), (40) instead of (38).

We would like to notice that similar construction does not work for extensions of the quadratic form of Laplacian generated by the interaction with $\delta$-potential. In this case corresponding quadratic form (42) and expressions (43), (44) are dimensionless and counterpart of the last integral does not tend to zero for $\Lambda\to\infty$. That is we may assume that the closure of the quadratic form $\omega_{\kappa}$ acting on the set $\mathcal{D}(\mathfrak{l})$ represents a restriction of the form $\omega_{\kappa}$

$$\overline{\omega_{\kappa}|_{\mathcal{D}(\mathfrak{l})}}=\tilde{\omega}_{\kappa}\subset\omega_{\kappa}\quad\tilde{\omega}_{\kappa}\neq\omega_{\kappa}.$$

This restriction corresponds to a self-adjoint operator different from $L_{\kappa}^{1/2}$ such that quadratic form of its square is not an extension of the quadratic form of Laplacian. And this means that equation (9) is not valid for quadratic forms extended by means of the interaction with $\delta$-potential.

6 Eigenstates normalization and eigenvalues

Further development of Gaussians $\Omega_{\kappa}$ involves introduction of the scalar product and in particular investigation of the possibility to normalize the functionals with respect to the ground state of the free theory $\Omega_{0}$. Such a normalization can be performed by means of the functional integration, the result is proportional to

$$\frac{\langle\Omega_{\kappa}|\Omega_{\kappa}\rangle}{\langle\Omega_{0}|\Omega_{0}\rangle}=\frac{\int\exp\{-\omega_{\kappa}(A)\}\,\prod_{\boldsymbol{x}}\delta A(\boldsymbol{x})}{\int\exp\{-\omega(A)\}\,\prod_{\boldsymbol{x}}\delta A(\boldsymbol{x})}=\exp\{\frac{1}{8}\operatorname{Tr}(\ln L-\ln L_{\kappa})\}.$$

Using the equation

$$(R_{\bar{\mu}}v,R_{\mu}v)=\frac{\partial}{\partial\mu}(v,R_{\mu}v)=\sigma^{\prime}(\mu),$$

derived by differentiation of (15), it is not hard to see, from the representation of logarithm similar to (36), that the trace of the difference of logarithms can be written as the following expression

$$\operatorname{Tr}(\ln L_{\kappa}-\ln L)=\frac{1}{2\pi i}\int_{0}^{\infty}\bigl{(}\frac{\sigma^{\prime}(\lambda+i0)}{\kappa-\sigma(\lambda+i0)}-\frac{\sigma^{\prime}(\lambda-i0)}{\kappa-\sigma(\lambda-i0)}\bigr{)}\ln\lambda\,d\lambda.$$ (45)

The latter can be further estimated as the integral

$$\operatorname{Tr}(\ln L_{\kappa}-\ln L)\sim\\ \sim\frac{1}{2\pi i}\int_{0}^{\infty}\frac{\sigma(\lambda+i0)-\sigma(\lambda-i0)}{(\kappa-\sigma(\lambda+i0))(\kappa-\sigma(\lambda-i0))}\sigma^{\prime}_{\text{ln}}(\lambda)\ln\lambda\,d\lambda\sim-\int_{C}^{\infty}\frac{d\lambda}{\lambda\ln\lambda},$$ (46)

which diverges as dilogarithm of the upper limit. Here $\sigma^{\prime}_{\text{ln}}(\lambda)$ is the term of $\sigma^{\prime}(\lambda)$ stemming from the logarithm in $\sigma(\lambda)$ – it is the first term in the exression

$$\frac{\partial\sigma}{\partial\lambda}=-\frac{2\pi}{\lambda}+\frac{2\pi}{\lambda^{3/2}x}(\cos\sqrt{\lambda}x\operatorname{Si}\sqrt{\lambda}x-\sin\sqrt{\lambda}x\operatorname{Ci}\sqrt{\lambda}x)+\\ +\frac{2\pi}{\lambda}(\sin\sqrt{\lambda}x\operatorname{Si}\sqrt{\lambda}x-\cos\sqrt{\lambda}x\operatorname{Ci}\sqrt{\lambda}x)$$ (47)

for $\sigma(\lambda)$ from (23) or the first term in the expression

$$\frac{\partial\sigma_{\shortparallel}^{-}}{\partial\lambda}=-\frac{4\pi}{3\lambda}-\frac{12\pi}{\lambda^{2}x^{2}}+\bigl{(}\frac{12\pi}{\lambda^{5/2}x^{3}}-\frac{4\pi}{\lambda^{3/2}x}\bigr{)}\bigl{(}\cos\sqrt{\lambda}x\operatorname{Si}\sqrt{\lambda}x-\sin\sqrt{\lambda}x\operatorname{Ci}\sqrt{\lambda}x\bigr{)}\\ +\frac{12\pi}{\lambda^{2}x^{2}}\bigl{(}\sin\sqrt{\lambda}x\operatorname{Si}\sqrt{\lambda}x+\cos\sqrt{\lambda}x\operatorname{Ci}\sqrt{\lambda}x\bigr{)},$$ (48)

for $\sigma_{\shortparallel}^{-}(\lambda)$ from (35). Further one may note that contributions to (45) from other terms of derivatives (47), (48) are finite and that the infinity of (46) is eliminated when substracting traces of logarithms of $L_{\kappa}$ for sufficiently close pairs $\boldsymbol{x}_{1}$, $\boldsymbol{x}_{2}$ and $\boldsymbol{x}_{3}$, $\boldsymbol{x}_{4}$

$$|\operatorname{Tr}\ln L_{\kappa;\boldsymbol{x}_{1},\boldsymbol{x}_{2}}-\operatorname{Tr}\ln L_{\kappa;\boldsymbol{x}_{3},\boldsymbol{x}_{4}}|<\infty.$$

Infinitesimal version of this statement also holds, that is the derivative

$$\frac{\partial\operatorname{Tr}\ln L_{\kappa;\boldsymbol{x}_{1},\boldsymbol{x}_{2}}}{\partial x}=\frac{1}{2\pi i}\int_{0}^{\infty}\Bigl{(}\bigl{(}\frac{\partial_{x}\partial_{\lambda}\sigma}{\kappa-\sigma}\bigr{|}_{\lambda+i0}-\frac{\partial_{x}\partial_{\lambda}\sigma}{\kappa-\sigma}\bigr{|}_{\lambda-i0}\bigr{)}-\\ -\bigl{(}\frac{\partial_{x}\sigma\partial_{\lambda}\sigma}{(\kappa-\sigma)^{2}}\bigr{|}_{\lambda+i0}-\frac{\partial_{x}\sigma\partial_{\lambda}\sigma}{(\kappa-\sigma)^{2}}|_{\lambda-i0}\bigr{)}\Bigr{)}\ln\lambda\,d\lambda$$ (49)

is finite. This finiteness can be directly verified for $\sigma_{\shortparallel}^{-}(\lambda)$, but for other discsussed Nevanlinna functions some integration by part is required in (45) before applying derivative w. r. to $x$.

One can also note that denominators in (45) and (49)

$$\kappa-\sigma(\lambda\pm i0)$$

have zeros at $\lambda=0$ for the distance $|\boldsymbol{x}_{2}-\boldsymbol{x}_{1}|$ tending to the limit $x_{\mathrm{b}}$ imposed by $\kappa$ (sketch in Fig. 1 is valid for all discussed Nevanlinna functions). This results in the infinite growth of the norm of the eigenstate corresponding to $\boldsymbol{x}_{1}$, $\boldsymbol{x}_{2}$ such that $|\boldsymbol{x}_{2}-\boldsymbol{x}_{1}|=x_{b}$

$$\frac{\langle\Omega_{\kappa,x_{\mathrm{b}}}|\Omega_{\kappa,x_{\mathrm{b}}}\rangle}{\langle\Omega_{\kappa,x}|\Omega_{\kappa,x}\rangle}=\frac{1}{4}\operatorname{Tr}(\ln L_{\kappa,x}-\ln L_{\kappa,x_{\mathrm{b}}})\sim\frac{1}{4}\int_{0}^{C}\frac{d\lambda}{\lambda\ln^{2}\lambda}.$$

The above considerations show that constructed eigenstates of the free Hamiltonian for sufficiently close pairs $\boldsymbol{x}_{1}$, $\boldsymbol{x}_{2}$ and $\boldsymbol{x}_{3}$, $\boldsymbol{x}_{4}$ reside in their own 2-point functional space $H_{2}$ not connected to the space $H_{0}$ generated by excitations of $\Omega_{0}$. The space $H_{2}$ can be represented as Hilbert space with basis of excitations of Gaussian $\Omega_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$ with elements of the form

$$\Omega(A)=\sum_{n,\boldsymbol{p}_{1},\ldots\boldsymbol{p}_{n}}\psi(\boldsymbol{p}_{1},\ldots\boldsymbol{p}_{n})\Psi_{\kappa}^{{\dagger}}(\boldsymbol{p}_{1})\dots\Psi_{\kappa}^{{\dagger}}(\boldsymbol{p}_{n})\,\Omega_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}(A),$$ (50)

where $\Psi_{\kappa}^{{\dagger}}(\boldsymbol{p})$ are creation operators depending on $\boldsymbol{x}_{1}$, $\boldsymbol{x}_{2}$

$$\Psi_{\kappa;\boldsymbol{x}_{1},\boldsymbol{x}_{2}}^{{\dagger}}(\boldsymbol{p})=L_{\kappa;\boldsymbol{x}_{1},\boldsymbol{x}_{2}}^{1/2}A-\frac{i\delta}{\delta A(\boldsymbol{p})}.$$

At the same time functionals $\Omega_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$ and $\Omega_{\kappa}^{\boldsymbol{x}_{3},\boldsymbol{x}_{4}}$ have different behaviour in the vicinity of functions with expansion (4) in pairs of points $\boldsymbol{x}_{1}$, $\boldsymbol{x}_{2}$ and $\boldsymbol{x}_{3}$, $\boldsymbol{x}_{4}$. Thus they corrspond to different self-adjoint operators in the same Hilbert space $H_{2}$, denote them as $\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$ and $\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$. These operators are defined on different domains for non-coinciding pairs $\boldsymbol{x}_{1}$, $\boldsymbol{x}_{2}$ and $\boldsymbol{x}_{3}$, $\boldsymbol{x}_{4}$, so that the action of $\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$ is not defined on $\Omega_{\kappa}^{\boldsymbol{x}_{3},\boldsymbol{x}_{4}}$ and vice versa. But it turns out that it is possible to construct a model including as dynamical variables not only the field $A(\boldsymbol{p})$ but also centers of external sources $\boldsymbol{x}_{1}$ and $\boldsymbol{x}_{2}$. For this one need to define action of $\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$ at least on functionals generated by excitations of $\Omega_{\kappa}^{\boldsymbol{x}_{3},\boldsymbol{x}_{4}}$ for sufficiently close pairs $\boldsymbol{x}_{1}$, $\boldsymbol{x}_{2}$ and $\boldsymbol{x}_{3}$, $\boldsymbol{x}_{4}$.

We cannot identify the action $\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}\Omega_{\kappa}^{\boldsymbol{x}_{3},\boldsymbol{x}_{4}}$ with an element of the Hilbert space, as long as appearance of “additional” eigenfunctional contradicts the completeness of the set of excitations of $\Omega_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$ from the decomposion (50). But we can expand $\Omega_{\kappa}^{\boldsymbol{x}_{3},\boldsymbol{x}_{4}}$ via ecxitations of $\Omega_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$ and calculate the action of $\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$ in the term by term way

$$\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}\Omega_{\kappa}^{\boldsymbol{x}_{3},\boldsymbol{x}_{4}}=\sum_{n,\boldsymbol{p}_{1},\ldots\boldsymbol{p}_{n}}\psi(\boldsymbol{p}_{1},\ldots\boldsymbol{p}_{n})\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}\Psi_{\kappa}^{{\dagger}}(\boldsymbol{p}_{1})\dots\Psi_{\kappa}^{{\dagger}}(\boldsymbol{p}_{n})\,\Omega_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}(A),$$ (51)

and then use the result to define expectation (quadratic form)

$$\langle\Omega_{\kappa}^{\boldsymbol{x}_{3},\boldsymbol{x}_{4}}|\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}|\Omega_{\kappa}^{\boldsymbol{x}_{3},\boldsymbol{x}_{4}}\rangle,$$ (52)

which should be finite. To justify this assumption one can remind the finite-dimensional theory of singular perturbations of differential operators [10]. There different perturbations of the same operator act in different spaces but the scalar product like (52) is well defined. That is we can fix some functional $\Omega(A)$ and then using different expansions

$$\Omega(A) =\sum_{n,\boldsymbol{p}_{1},\ldots\boldsymbol{p}_{n}}\psi_{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}(\boldsymbol{p}_{1},\ldots\boldsymbol{p}_{n})\Psi_{\kappa;\boldsymbol{x}_{1},\boldsymbol{x}_{2}}^{{\dagger}}(\boldsymbol{p}_{1})\dots\Psi_{\kappa;\boldsymbol{x}_{1},\boldsymbol{x}_{2}}^{{\dagger}}(\boldsymbol{p}_{n})\,\Omega_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}(A)=$$

$$=\sum_{n,\boldsymbol{p}_{1},\ldots\boldsymbol{p}_{n}}\psi_{\boldsymbol{x}_{3},\boldsymbol{x}_{4}}(\boldsymbol{p}_{1},\ldots\boldsymbol{p}_{n})\Psi_{\kappa;\boldsymbol{x}_{3},\boldsymbol{x}_{4}}^{{\dagger}}(\boldsymbol{p}_{1})\dots\Psi_{\kappa;\boldsymbol{x}_{3},\boldsymbol{x}_{4}}^{{\dagger}}(\boldsymbol{p}_{n})\,\Omega_{\kappa}^{\boldsymbol{x}_{3},\boldsymbol{x}_{4}}(A)$$

define on it actions of operators $\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$ with varying position of centers of external sources $\boldsymbol{x}_{1}$, $\boldsymbol{x}_{2}$ and further calculate expectation values $\langle\Omega|\mathscr{H}_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}|\Omega\rangle$. Important condition for such calculations is the possibility to uniformly define eigenvalues $\operatorname{Tr}L_{\kappa;\boldsymbol{x}_{1},\boldsymbol{x}_{2}}^{1/2}$ of functionals $\Omega_{\kappa}^{\boldsymbol{x}_{1},\boldsymbol{x}_{2}}$ from equation (8) for different pairs $\boldsymbol{x}_{1}$, $\boldsymbol{x}_{2}$. Spectral integral for the difference of traces $\operatorname{Tr}L_{\kappa}^{1/2}-\operatorname{Tr}L^{1/2}$ diverges even worse than (45)

$$\operatorname{Tr}(L_{\kappa}^{1/2}-L^{1/2})=\frac{1}{2\pi i}\int_{0}^{\infty}\bigl{(}\frac{\sigma^{\prime}(\lambda+i0)}{\kappa-\sigma(\lambda+i0)}-\frac{\sigma^{\prime}(\lambda-i0)}{\kappa-\sigma(\lambda-i0)}\bigr{)}\lambda^{1/2}\,d\lambda\sim\\ \sim\frac{1}{2\pi i}\int_{0}^{\infty}\frac{(\sigma(\lambda+i0)-\sigma(\lambda-i0))\sigma^{\prime}_{\text{ln}}(\lambda)\lambda^{1/2}}{(\kappa-\sigma(\lambda+i0))(\kappa-\sigma(\lambda-i0))}\,d\lambda\sim-\int_{0}^{\infty}\frac{d\lambda}{\lambda^{1/2}\ln^{2}\lambda}.$$ (53)

However the derivative of this integral w. r. to the distance between $\boldsymbol{x}_{1}$ and $\boldsymbol{x}_{2}$

$$\frac{\partial\operatorname{Tr}L_{\kappa;\boldsymbol{x}_{1},\boldsymbol{x}_{2}}^{1/2}}{\partial x}=\frac{1}{2\pi i}\int_{0}^{\infty}\Bigl{(}\bigl{(}\frac{\partial_{x}\partial_{\lambda}\sigma}{\kappa-\sigma}(\lambda+i0)-\frac{\partial_{x}\partial_{\lambda}\sigma}{\kappa-\sigma}(\lambda-i0)\bigr{)}-\\ -\bigl{(}\frac{\partial_{x}\sigma\partial_{\lambda}\sigma}{(\kappa-\sigma)^{2}}(\lambda+i0)-\frac{\partial_{x}\sigma\partial_{\lambda}\sigma}{(\kappa-\sigma)^{2}}(\lambda-i0)\bigr{)}\Bigr{)}\lambda^{1/2}\,d\lambda$$ (54)

as well as the difference

$$\operatorname{Tr}L_{\kappa;\boldsymbol{x}_{1},\boldsymbol{x}_{2}}^{1/2}-\operatorname{Tr}L_{\kappa;\boldsymbol{x}_{3},\boldsymbol{x}_{4}}^{1/2}$$

are finite. Here considerations that were used for traces of logarithms also work. That is only the first terms from (47) and (48) bring divergency to (53), and the latter is generated by the constant in the gap (38), (39) or (40) and in this way it does not depend on $\boldsymbol{x}_{1}$ and $\boldsymbol{x}_{2}$. Another noticable property of divergency (53) is that it is negative and hence it forms infinitely deep potential well, which separates space of functionals of this model from functionals of the free field.

Conclusion

We have shown that in order to construct alternative solutions to the eigenstate equation of the free scalar quantum theory the external potential incorporating at least two connected singularities has to be used. The distance between the singularities is proved to be limited by the extension parameter. We have shown for the potential in question that the quadratic form of the square root is reconstructed by its action on the domain of quadratic form of Laplacian. We have estimated spectral integrals for traces of logarithms and find that constructed eigenfunctionals reside in one Hilbert space for sufficiently close positions of pairs of singularities in external potentials. We have also shown that corresponding eigenvalues of quantum Hamiltonians continuously depend on the above positions.