Hamiltonians Without Ultraviolet Divergence for Quantum Field Theories

We propose a way of defining Hamiltonians for quantum field theories without any renormalization procedure. The resulting Hamiltonians, called IBC Hamiltonians, are mathematically well-defined (and in particular, ultraviolet finite) without an ultraviolet cut-off such as smearing out the particles over a nonzero radius; rather, the particles are assigned radius zero. These Hamiltonians agree with those obtained through renormalization whenever both are known to exist. We describe explicit examples of IBC Hamiltonians. Their definition, which is best expressed in the particle-position representation of the wave function, involves a kind of boundary condition on the wave function, which we call an interior-boundary condition (IBC). The relevant configuration space is one of a variable number of particles, and the relevant boundary consists of the configurations with two or more particles at the same location. The IBC relates the value (or derivative) of the wave function at a boundary point to the value of the wave function at an interior point (here, in a sector of configuration space corresponding to a lesser number of particles).


Introduction
In many quantum field theories (QFTs), the formulas that one obtains for the Hamiltonian (by means of quantization or other heuristics) contain terms for the creation and annihilation of particles that are ultraviolet (UV) divergent. The problem can be avoided by a UV cut-off, i.e., by discretizing space or treating the electron (and other particles) not as a point but smearing it out instead over a small positive radius; however, these procedures tend to break the Lorentz invariance, and there is no empirical evidence for either discrete space or a positive electron radius [1,2]. Since QFTs are expected to provide merely effective descriptions and not to be valid on very small scales (see, e.g., [3] for discussion), this problem is usually neither regarded as unexpected nor as threatening. Still, the problem has attracted much interest over the years, and has been investigated with some success, as one can take a limit of removing the cutoff through a renormalization procedure for some QFTs. We report here that the UV divergence problem does not occur if the Hamiltonian is defined in a novel way that proceeds directly without renormalization. 1 We describe these Hamiltonians for some non-relativistic QFTs and report about recent proofs [4,5,6] showing that the Hamiltonians obtained in this way are well-defined and self-adjoint. In this approach, space is continuous (as opposed to a lattice), and the radius of the electron (or other particles) is zero. The key element of the approach is a new type of boundary condition that we call an interior-boundary condition (IBC) because it relates the values of ψ on the boundary of configuration space Q to the values in the interior of Q, as we will explain presently. As such, IBCs have already been considered in the past, in one form or another, starting with Landau and Peierls [7] in 1930, followed by [8,9,10,11,12,13,14] (see below for a brief overview). Our new insight is that IBCs, when used for defining QFTs, can avoid the problem of UV divergence. In the past, IBCs received little attention, were not explored systematically, were not studied rigorously except to a rather limited extent, and were not considered for addressing UV divergence.
IBCs as we use them are formulated in the particle-position representation of the state vector ψ in Hilbert space H . Here, "particle representation" means that H is represented as a Fock space (or, if appropriate, a tensor product of several Fock spaces), and "position representation" that the contribution from the n-particle sector of Fock space is represented (like a wave function in quantum mechanics) as a function of n points in 3-dimensional physical space. Specifically, if H is a single Fock space then ψ ∈ H can be viewed as a function on a configuration space of a variable number of particles, such as see Figure 1, where ∆ n = (x 1 , . . . , x n ) ∈ (R 3 ) n : x i = x j for some i = j is the "diagonal," i.e., the set of collision configurations (i.e., those with two or more particles at the same location). The relevant boundary ∂Q of Q is ∪ ∞ n=0 ∆ n ; the IBC relates the values of ψ on ∂Q n = ∆ n to the values of ψ in the interior of Q n−1 , namely at the configuration with one particle removed (or possibly with more than one particle removed, if more than two particles collide).
The IBC approach allows a nonzero flux of probability out of or into the boundary while a term in the Hamiltonian ensures that all of the |ψ| 2 probability lost in one sector gets added in another sector, so that total probability in Q is conserved. The IBC can thus be regarded as enabling the creation of a contribution to ψ on Q n that flows out of the boundary ∂Q n and depends on the values of ψ on Q n−1 . The approach has been shown in the non-relativistic setting to provide a well-defined Hamiltonian H IBC whose domain consists of functions satisfying the IBC, without the need for a UV cutoff or renormalization.
Moreover, as we will explain, H IBC provides a natural mathematical interpretation of the original (UV divergent) expression H orig for the Hamiltonian. One way of expressing Figure 1: Illustration of a configuration space for a variable number of particles in 1 dimension: (a) The zero-particle sector contains only one configuration, the empty configuration; (b) the one-particle sector is a copy of physical space; (c) the two-particle sector; (d) the three-particle sector. In this example, we have not removed the diagonal ∆ n .
the problem with H orig is to say that, due to divergent creation terms, H orig ψ exists as a distribution (such as the Dirac delta distribution) but not as a square-integrabe function, so that H orig does not define an operator in Hilbert space. The IBC approach actually shares the divergent creation term, while the IBC forces wave functions to diverge at the boundary ∂Q n at a particular rate, with the consequence that the Laplacian of ψ also exists as a distribution but not as a square-integrable function. The IBC ensures that the two distributional contributions to H IBC ψ (one from divergent creation terms and one from the Laplacian acting on a diverging ψ) exactly cancel each other, so that their sum is actually a square-integrable function. In fact, the IBC characterizes the domain of functions on which this cancelation occurs. For ψ diverging at the boundary, the annihilation terms in H orig need a separate definition, which the IBC approach provides as well in a natural manner. After a condition equivalent to an IBC had been considered as early as 1930 [7], this approach was not followed further; instead, much research took for granted that the Hamiltonian of a QFT is the sum of two self-adjoint operators, the free Hamiltonian and the interaction Hamiltonian. This is not so in the IBC approach, where the free Hamiltonian and the interaction Hamiltonian each map square-integrable functions to distributions and thus are not defined as self-adjoint operators in Hilbert space. It is only their sum that is a self-adjoint operator in Hilbert space. In particular, the full Hamiltonian cannot be regarded as a perturbation of the free Hamiltonian. This feature makes a relevant difference, at least for some UV divergence problems.
For some QFTs, it is possible to take a limit, after introducing a UV cutoff, in which the cutoff is removed, possibly after subtracting an infinite constant from the Hamiltonian. It turns out [4,5] that, at least for some models, this limiting Hamiltonian (called the renormalized Hamiltonian) agrees with H IBC ; however, this was not visible from the renormalization procedure before considering the IBC.
In this paper, we describe examples of IBCs and how they help define a Hamiltonian, results about the rigorous existence and self-adjointness of the Hamiltonians, and how these Hamiltonians are related to some known cases in which a UV cut-off can be removed, thus making it plausible that the IBC Hamiltonians are physically relevant and not merely mathematical curiosities. Specifically, we consider two models: In Model 1, x-particles can emit and absorb y-particles, and both kinds of particles are nonrelativistic; we give the full definition in Section 2 below. The x-particles will also be called the "sources." Model 2 is a simplified version of Model 1 in which the x-particles cannot move but are fixed at certain locations; it arises as a limiting case of Model 1 in which the mass m x of the x-particle tends to ∞; for simplicity, we consider only a single x-particle (and call its location the origin). That is, y-particles can be created and annihilated at the origin of R 3 . For Model 2 (the one with fixed x-particles), it is known that, after a UV cut-off is introduced, the cut-off can be removed by means of a suitable limiting procedure (renormalization), and we can show that the limiting Hamiltonian H ∞ thus obtained coincides with H IBC up to addition of a finite constant (see Section 3.6 below). For Model 1, however, such a procedure was not known before the IBC approach succeeded (although a renormalization procedure was known for a similar model, the Nelson model [15]). The IBC approach provides, also for Model 1, a well-defined and self-adjoint Hamiltonian without renormalization [6].
Let us give a brief overview of prior works considering IBCs. Historically, the first equation equivalent to an IBC was written down, as far as we know, by Landau and Peierls in 1930 [7] when trying to formulate quantum electrodynamics in the particleposition representation, although their Hamiltonian was still UV divergent. Moshinsky [8,9] considered, as an effective description of nuclear reactions, a model using the Laplace operator with an IBC for (fixed sources and) a mini-Fock space having only sectors with n = 0 and n = 1 particles. Yafaev [13] independently considered the same model and proved for that Hamiltonian that it is well defined and self-adjoint. Thomas [11] considered a similar model with moving sources, also with only two sectors. Georgii and Tumulka [14,Sec. 6] considered IBCs for boundaries of codimension 1; such IBCs are developed systematically in [16]. Recent and upcoming works exploring various aspects of IBCs include [4,5,6,17,18,19,20,21,22,23,24].
We give a gentle introduction to IBCs and the ideas behind them in [17]. Mathematical proofs of our main claims are provided in [4,5,6]; the case of fixed sources in 3d is treated in [4], the case of moving sources in 2d in [5], and that of moving sources in 3d in [6]. In [5], it is also shown how to generalize the IBC approach to a relativistic dispersion relation E = m 2 c 4 + p 2 c 2 (for scalar, not Dirac wave functions) instead of E = p 2 /2m; the condition replacing the IBC is then no longer literally a boundary condition and therefore called an abstract boundary condition. Keppeler and Sieber [18] study IBCs in one space dimension. Galvan [19] has proposed an approach similar to IBCs. A goal for the future is to apply the IBC approach to quantum electrodynamics and other serious theories; work on implementing IBCs for the Dirac equation is in progress [20,24]. For the time being, we report results for non-relativistic model QFTs, based on the Laplace operator.
This paper is organized as follows. In Section 2, we give a more detailed description of Model 1 and Model 2, two simple models of non-relativistic QFT. In Section 3, we describe for each of the two models how to set up an IBC and the corresponding Hamiltonian. In Section 4, we conclude.

Two Models of Non-Relativistic QFT
In this section, we describe Model 1 and Model 2, two simple non-relativistic QFTs. We set up their Hilbert spaces H , describe their original Hamiltonians H orig , and explain why they are UV divergent. Model 1 is a natural, physically reasonable model, while Model 2 is an artificially simplified version that will be useful for an easy discussion of IBCs. In Section 3, we will describe, for each of the two models, our alternative definition of the Hamiltonian involving an interior-boundary condition.

Model 1: x-Particles Emit and Absorb y-Particles
Model 1 is a QFT related to the Lee model [25], Schweber's scalar field model [26, p. 339], and the Nelson model [15]. It involves two species of particles, x and y; the x-particles can emit and absorb y-particles. Both species are spinless and assume the non-relativistic dispersion relation E = p 2 /2m with masses m x , m y > 0. The Hilbert space is a tensor product of Fock spaces, with where S − is the anti-symmetrization operator, S + is the symmetrization operator, and S ± L 2 (· · · ) are their ranges (i.e., the spaces of (anti-)symmetric functions on (R 3 ) n ).
Here, we take x-particles to be fermions and y-particles to be bosons. 2 R 3 is understood as physical space, i.e., ordinary position space. In the following, we simply write L 2 (R 3 ) for L 2 (R 3 , C). A vector Ψ ∈ H can be regarded as a function with Q x = Q y the configuration space of a variable number of particles, where the union is understood as a disjoint union and (R 3 ) 0 = {∅}. We will discuss removing the collision configurations later. We call the function ψ the particle-position representation of the vector Ψ ∈ H . A generic element of Q x × Q y can be written as (x, y) = (x 1 , . . . , x m , y 1 , . . . , y n ), where bold-face symbols denote vectors in 3-space, while x denotes a configuration of x-particles and y one of y-particles; we will often write x m instead of x to convey that the configuration consists of m x-particles, and likewise y n instead of y. We call (R 3 ) m the m-particle sector of Q x and (R 3 ) m × (R 3 ) n the (m, n)-particle sector of Q x × Q y . Likewise, we say (m, n)-particle sector of H (or of Ψ, or of ψ) and write H (m,n) (or Ψ (m,n) or ψ (m,n) ) for S − L 2 (R 3 , C) ⊗m ⊗ S + L 2 (R 3 , C) ⊗n , respectively for the projection of Ψ to that subspace, and for the restriction of ψ to (R 3 ) m × (R 3 ) n . The ψ (m,n) function is anti-symmetric in the x variables and symmetric in the y variables.
The spaces Q x , Q y are equipped with the volume measure The inner product in H is then given by

Original Hamiltonian
As our example of a (non-relativistic) Hamiltonian H orig , we take using the notation δ 3 for the 3-dimensional Dirac delta function and for the configuration of n − 1 y-particles with the j-th particle removed; E 0 is the energy that must be expended for creating a y-particle (the "y rest energy"), and g ∈ R is a coupling constant (i.e., the "charge" of an x-particle). Instead of ψ (m,n) (x m , y n ), we can also simply write ψ(x m , y n ), as the argument uniquely determines which sector of ψ must be used; we sometimes find it useful to use the more explicit notation to make the relations to other sectors more easily visible. The first term in (10) is the free fermion Hamiltonian H x ψ, the second and third terms are the free boson Hamiltonian H y ψ, and the third and fourth line together are the interaction Hamiltonian H inter ψ responsible for the creation and annihilation of y-particles. In terms of creation and annihilation operators, with † denoting the adjoint operator, and a x (q), a y (q) the annihilation operators for an x, y-particle at location q in position space, formally defined by The combinatorial factors √ m + 1 and m −1/2 arise from the fact that every configuration of m x-particles occurs in m! different permutations. Since the operators a x,y (q) contain the point evaluation of functions, they can actually be defined rigorously on a dense subspace of H , while their adjoints are ill defined for any nonzero ψ.
Model 1 is not Galilean covariant but a modification that is can be set up [27] and can be combined with IBCs [21, Section 4.2].

UV Divergence of the Original Hamiltonian
The original Hamiltonian H orig as defined in (10) above is UV divergent and thus ill defined. The source of the difficulty is the delta function in the last line of (10): Since a delta function is not a square-integrable function, the right-hand side of (10) does not lie in H for any choice of ψ; thus, (10) does not define an operator in H . It can also be pointed out that the delta function δ 3 (which plays the role in (10) of the wave function of a newly created y-particle) is a state of infinite energy, The standard procedure for obtaining a well-defined Hamiltonian (UV cut-off) is to replace the delta function by a square-integrable function ϕ : which amounts to saying that an x-particle has an extended charge distribution with density function ϕ (if ϕ is real-valued). Equivalently, the interaction Hamiltonian is replaced by where

Simplified Version: Model 2
In Model 2, there is only one x-particle, and it is fixed at the origin. The only dynamical (and quantized) degrees of freedom reside in the y-particles; such models were considered in particular by van Hove [28,29] and Lee [25], and they tend to arise in the limit m x → ∞. Specifically, we take the configuration space of Model 2 to be the set where 0 denotes the origin in R 3 ; a configuration is denoted by y n = (y 1 , . . . , y n ). Note that we exclude the possibility that any y-particle can be at the location of the x-particle (i.e., the origin), but we do not exclude the possibility that two y-particles can be at the same location because for the purposes of this model, in which the y-particles do not interact, there is no need to exclude such configurations. Correspondingly, we take the Hilbert space to be the bosonic Fock space Since for the definition of the Hilbert space it plays no role whether the origin is excluded or not, elements ψ ∈ H can be regarded as complex-valued functions on Q that are permutation-symmetric in every sector; thus, H ⊂ L 2 (Q), where Q is thought of as equipped with the measure µ defined in the same way as in (7) (L 2 (Q) contains also non-symmetric functions). The "original" Hamiltonian is now a simplified version of (10): Also this Hamiltonian is UV divergent.

Interior-Boundary Condition and Corresponding Hamiltonian for the Two Models
For Model 1 and Model 2, we describe the IBC and the Hamiltonian H IBC . We also explain why H IBC is a reasonable interpretation of the formula for H orig . We begin with the simpler scenario of Model 2.

IBC for Model 2
Let S 2 denote the unit sphere in R 3 . The IBC demands the following: For every ω ∈ S 2 , The IBC is a condition on the wave function ψ at or near the "boundary" of configuration space; the relevant boundary ∂Q That is, the boundary consists of those configurations at which one of the y-particles collides with the x-particle. Due to the permutation symmetry of ψ (n+1) , we can assume without loss of generality that it is the n + 1-st variable, y n+1 , that approaches 0. The name "interior-boundary condition" reflects the fact that (27) relates the values (or limits) of ψ on the boundary ∂Q y to values of ψ in the interior of configuration space Q (namely, in the interior of a different sector corresponding to a lesser number of y-particles).
The IBC (27) allows in particular that ψ(y n , rω) diverges like 1/r as r → 0. In fact, it requires that ψ so diverges whenever the right-hand side of (27) is nonzero. 3 By permutation symmetry, ψ(y n ) diverges as any y j approaches 0.
On wave functions ψ satisfying the IBC (27), the Hamiltonian H = H IBC is defined by This equation differs from the expression (26) for H orig only in the second line, where ψ(y n , 0) has been replaced by a more complicated expression involving the behavior of ψ near the configuration (y n , 0); after all, ψ diverges at this configuration by virtue of the IBC, so the expression ψ(y n , 0) does not make sense. The great similarity between H IBC and H orig adds to the suggestion that H IBC is the right Hamiltonian because it is a mathematical interpretation of the formal expression H orig .
In view of the delta function appearing in (29), it may seem unlikely that such a Hamiltonian can be well defined. However, a theorem from [4], repeated here as Theorem 1 in Section 3.2 below, shows that in fact it is well defined. It turns out that the last line in (29) always gets canceled by contributions to the first line, which may contain delta functions because we now allow wave functions that diverge like 1/r = 1/|y j |, and the Laplacian of 1/r is −4πδ 3 (as readers may recall from electrostatics, where 1/r occurs as the Coulomb potential φ generated by a point charge, satisfying the Poisson equation ∆φ = −4πρ with charge density ρ). As a consequence of this cancelation, H IBC ψ is a square-integrable function for ψ satisfying the IBC. Conversely, if we want the Laplacian to produce a contribution that is a Dirac delta contribution, we are led to allowing ψ that diverges like 1/r, and then the necessary and sufficient condition for the cancelation is just the IBC. See also Remark 6 in Section 3.3 below for further discussion.
Let us point out a connection between the second line of (29), the line that differed from the original Hamiltonian (26), and the corresponding line in (26). Think of ψ y n , rω as a function of r; as r → 0, a ψ from the domain of H IBC can be expanded in the form with complex coefficients c −1 (y n ), c 0 (y n ); moreover, ∂ ∂r rψ(y n , rω) = c 0 (y n ) + O(r) , which yields c 0 in the limit r → 0. If ψ did not diverge as r → 0, this would be exactly the value lim r 0 ψ(y n , rω) = ψ y n , 0 occurring in (26) in the second line. So, (29) is really quite similar to (26). Indeed, we can regard c 0 as the natural interpretation of the expression ψ(y n , 0) in a situation in which ψ also contains a term c −1 /r diverging at 0; put differently, it is a way of extending the annihilation operator a y (0) to functions containing a term c −1 /r. Thus, in total, H IBC appears as a natural interpretation of H orig , with the IBC characterizing the right domain on which the distribution-valued terms in H IBC ψ cancel each other, leaving a square-integrable function. Keppeler and Sieber [18], apart from developing a 1-dimensional version, give further reasons why the IBC model represents a reasonable interpretation of the original Hamiltonian (26). They argue [18, App. A] that for any eigenfunction ψ of (26), when integrating (26) in y n over a ball of radius r around the source, using the Gauss integral theorem on the Laplacian term, and taking the limit r → 0, one should obtain that Assuming that ψ (n) diverges as r → 0 no worse than 1/r, this relation is equivalent to the IBC (27). Finally, if every eigenfunction satisfies the IBC (27), then so does any function in the domain of the Hamiltonian. The corresponding reasoning in the 1-dimensional case yields that ψ should have a jump in the first derivative at the origin, with the magnitude of the jump proportional to the wave function in the next lower sector.
The Hamiltonian H IBC is not the sum of two self-adjoint operators, the free Hamilonian and an interaction Hamiltonian. That is because the free Hamiltonian is defined on a different domain than H IBC , containing wave functions that do not satisfy the IBC and do not diverge on the diagonal. It is the action of the Laplacian on functions that diverge on the diagonal that leads to delta functions and thus makes the last line of (29) possible, and it is only in conjunction with the IBC that the Hamiltonian (29) leads to conservation of probability.
Indeed, in order to understand why the IBC (27) was chosen this way, and how it works together with the formula (29) for the Hamiltonian, it is illuminating to calculate the balance equation for |ψ| 2 and check that |ψ| 2 is conserved. We go through such a calculation in the next section. The relevance of this theorem is that it makes clear that H IBC is not afflicted by UV divergence. We give here a calculation checking on a non-rigorous level that probability is conserved. Using the symbol

Self-Adjointness and Conservation of Probability
for the usual probability current, we obtain from (29) that, at any configuration y n (without any y-particle at the origin), If g = 0, then the last line can be re-written, using the IBC (27), as − g √ n + 1 2π Im 2π 2 √ n + 1 g m y lim r 0 rψ * (y n , rω) using that ∂ r (rψ) = ψ + r∂ r ψ and Im(r|ψ| 2 ) = 0. Thus, This equation possesses a simple interpretation: −r 2 times the integral is just the flux of probability current toward the origin across the sphere of radius r in the coordinate space of y n+1 ; the limit of that as r → 0 is the current of probability into the origin in the coordinate space of y n+1 ; summing over all y-particles would yield, due to the bosonic symmetry of ψ, n + 1 equal terms; thus, the second summand on the right-hand side of (37) is the total flux of probability into the boundary configurations obtained from y n by adding one y-particle at the origin. That is, the balance equation (37) asserts that the probability density |ψ| 2 changes in two ways, due to transport of probability in the n-particle sector and by increasing at just the rate at which probability disappears on the (n + 1)-particle sector by flowing into the boundary ∂Q (n+1) y . Therefore, total probability is conserved.
If g is set to 0, then the amount of probability exchanged between different sectors vanishes, and the norm of each sector, R 3n dy n |ψ(y n )| 2 , is conserved separately, corresponding to the fact that the Hamiltonian H IBC of (29) commutes with the y-particle number operator in this case. In fact, H IBC reduces to the free Hamiltonian (i.e., the "second quantization" of the free 1-particle operator −( 2 /2m y )∇ 2 y + E 0 ), and the IBC (27) asserts that, in the expansion ψ(y n , rω) = c −1 /r + c 0 + O(r) considered in (31), c −1 = 0; thus, for g = 0 the IBC demands no more than that ψ is non-singular at r = 0, as functions in the domain of the free Hamiltonian must be.

Remarks
1. Comparison to Bethe-Peierls boundary condition. The IBC (48) has some parallels to the Bethe-Peierls boundary condition [30], which reads for a wave function ψ ∈ L 2 (R 3 , C): with given constant α ∈ R. This condition is used for giving precise meaning to a Schrödinger equation for ψ : R 3 → C with a Dirac delta function as the potential, a kind of interaction known as a point interaction [31]. Like our IBCs (27), (45), (48), the Bethe-Peierls boundary condition (38) concerns the "boundary" at r = 0. However, in contrast to the IBC, which connects two sectors of Q, the Bethe-Peierls boundary condition involves only one sector, as the wave function ψ is defined on R 3 . Also, the Bethe-Peierls boundary condition implies zero current into r = 0 (see, e.g., [17] for more detail), whereas the IBC leads to a nonzero current into the boundary.
2. Comparison to vertex conditions in networks. An example of a known boundary condition leading to nonzero current into the boundary is provided by the vertex conditions for quantum mechanics on networks (see, e.g., [32] and references therein). A network, or graph, consists of several 1-dimensional spaces (i.e., intervals) called edges, glued together at their end points called vertices. The wave function is a (say, complex-valued) function on the network (i.e., on the union of the edges), and the Hamiltonian is given by the Laplace operator on each edge, along with boundary conditions for the end points of the edges, also called vertex conditions. The vertex conditions are related to the fact that, since no positive amount of probability can be located at a vertex, all of the probability current into the vertex must be compensated by a current out of the vertex; that is, for each vertex v, the sum of the currents along all edges e connected to v (pointing away from v) must vanish, i.e., that ψ is continuous at vertices, and where α ∈ R is a given constant and ∂ e ψ(v) means the derivative of ψ along the edge e, taken at v.
To draw parallels between (42) and IBCs, we may compare the edges of a network to the sectors of the configuration spaces considered in this paper; let us call both the "blocks" of the space. Of course, the edges of a graph have equal dimension while the sectors of our Q have different dimension. We may note that in both situations the probability loss in one block is compensated by a probability gain in others. Also, the condition of continuity (41) is a relation between the boundary value in one block and a value in another, but in contrast to an IBC the latter is also a boundary value; that is, in analogy to the terminology "interior-boundary condition," (41) is a "boundary-boundary condition." 3. Combinatorial factors. As remarked at the end of Section 2.1.1, the combinatorial factors √ n and √ n + 1 that appear in the Hamiltonian arise from the fact that we use ordered configurations (y 1 , . . . , y n ) although in nature configurations are unordered, as in {y 1 , . . . , y n }. This is also true of the IBC (27) and the associated Hamiltonian (29). If we used unordered configurations, which leads to topologically non-trivial configuration spaces (see, e.g., [33,34]), these factors would not appear.
4. Neumann-type IBC. Two well known types of boundary conditions are Dirichlet boundary conditions, ψ and Neumann boundary conditions, with n the normal vector on the boundary hypersurface. While the IBC (27) resembles the Dirichlet type in that it involves the values of ψ (or rather, rψ) on the boundary, one can as well set up a different IBC that resembles the Neumann type in that it involves the normal derivative of ψ (or rather, rψ). Here, r = 0 plays the role of the boundary, and the radial direction plays the role of the normal to the boundary. The Neumann-type IBC asserts that for any n ∈ {0, 1, 2, . . .}, any configuration y n ∈ (R 3 \ {0}) n , and any ω ∈ S 2 , lim r 0 ∂ ∂r rψ (n+1) (y n , rω) = g m y 2π 2 √ n + 1 ψ (n) (y n ) .
Also this IBC typically leads to ψ (n+1) (y n , rω) diverging like 1/r as r → 0. The name "Neumann-type IBC" should be understood with care since, as we have observed already in (31), the expression lim r→0 ∂ r (rψ) yields, for ψ = c −1 /r + c 0 + O(r), just c 0 ; as a consequence, in the case c −1 = 0 in which ψ can actually be evaluated at r = 0, the expression lim r→0 ∂ r (rψ) just yields ψ(r = 0), so that the IBC (45) has quite the character of a Dirichlet boundary condition.
On wave functions ψ satisfying the IBC (45), the Hamiltonian H = H IBC is defined by (Hψ) (n) (y n ) = − 2 2m y n j=1 ∇ 2 y j ψ (n) (y n ) + nE 0 ψ (n) (y n ) This equation differs from (29) (for the Dirichlet-type case) in the second line, which no longer involves a derivative, and in the last line, which we will discuss further in Remark 6 below.
While (45) and (46) together define a different time evolution for ψ than (27) and (29), Equation (37) is still true and guarantees that the amount of probability lost on the n + 1-sector of Q due to probability flux into the boundary is added on the n-sector. In fact, as we show in [4], also (45) and (46) define a self-adjoint operator, provided E 0 > 0.
For g = 0, it is again true (as in the Dirichlet-type case) that the sectors decouple, and that the y-particle number operator is conserved. However, the Hamiltonian does not become the free Hamiltonian. Rather, the IBC (45) becomes a Bethe-Peierls boundary condition (38) with α = 0, and the Hamiltonian becomes the second quantization of the negative Laplacian with a point interaction at the origin.
We remark that the family of IBCs just described actually depends only on 4 real parameters, although α, β, γ, δ would seem at first to be 4 complex (and thus 8 real) parameters. That is so because (50) requires α and γ to have equal phases, and (51) requires β and δ to have equal phases; since the phase of α * δ is the phase difference between α and δ, which must be equal to that between γ and β and thus by (52) to the phase of −1, α and δ must have equal phases. In other words, all of α, β, γ, δ must have the same phase, while their absolute values are constrained by (52), so that only 4 real parameters are independent.
On a non-rigorous level, the conservation of probability (i.e., self-adjointness of the Hamiltonian (49)) can be checked by means of a calculation similar to the one in Section 3.2 above (and to the one in [16]). The self-adjointness breaks down (already on the non-rigorous level) if we relax the conditions (50)-(52). If, however, (50)-(52) are satisfied, then (48) and (49) define a self-adjoint Hamiltonian, provided E 0 > 0 [4].

Dirac delta function terms in H.
Here is a reason for thinking that, among the many different IBCs that are mathematically possible corresponding to different choices of the constants α, β, γ, δ, only the Dirichlet-type IBC (27), corresponding to α = −4π/g, β = 0 = γ, and δ = g/4π, is physically relevant as a replacement of the original (UV divergent) Hamiltonian (26): It is the only choice that leads to a term in H IBC that reproduces the Dirac delta function terms in H orig , i.e., for which the last line of (49) agrees with the last line of (26). That is because the Dirichlet case is the only case in which the last line of (49) can be expressed in terms of ψ (n−1) .
And the last line of (49) is dictated by the condition that Hψ has to be an L 2 function and thus cannot contain contributions that are Dirac delta functions. Indeed, since for any α, β, γ, δ, the wave function ψ (n) diverges at the boundary like 1/r, the Laplacian in the Hamiltonian always yields a distribution of the form with some functions f and g; so the Dirac delta contributions need to be canceled, which leads to the last line of (49).
In the 1-dimensional case, it seems that the physically reasonable IBC is the one given in [18], which is of the Neumann type and involves the normal derivative on both sides of the "diagonal" (i.e., the collision configurations). This is the conclusion one reaches when demanding that the Laplacian term in the Hamiltonian cancels the Dirac delta function terms, starting from the 1-dimensional version of (26): Then, the derivative of ψ needs to have jumps of the appropriate magnitude.

7.
Positivity. Another reason for thinking that the Dirichlet-type IBC (rather than, say, Neumann-type) is the physically relevant choice may be that H IBC for a Dirichlet-type IBC is positive, as mentioned in Theorem 1 in Section 3.2. Generally speaking, the IBC approach neither requires nor guarantees that Hamiltonians are bounded from below. Presently, we do not know for which other choices of α, β, γ, δ the Hamiltonian will be positive, but we see reason to believe that, in the variant of the equations appropriate for Model 1, the Neumann-type IBC leads to a Hamiltonian that is not bounded from below. Further considerations about physical reasonableness of H IBC can be found in Sections 3.5 and 3.6 below.
8. Bohmian trajectories. There is a natural way of defining Bohmian trajectories for the models described in this paper; we describe this in detail elsewhere [21]. The Bohmian configuration Q t follows a Markov jump process in configuration space that is |ψ t | 2 -distributed at every time t. The process has finitely many jumps in every finite time interval. The pieces between the jumps are solutions to Bohm's equation of motion; in particular, they are deterministic, in contrast for example to Nelson's trajectories [35], which follow a diffusion process. In our Models 1 and 2, the jumps correspond to the creation or annihilation of a particle. The jumps to a lower sector (particle annihilation) occur whenever a y-particle hits an x-particle; in that event, the y-particle gets deleted from Q t . While the jumps to a lower sector are deterministic, the jumps to a higher sector (particle creation) are stochastic. They can occur at any configuration Q t , with a rate depending on Q t and ψ t , and lead to a configuration with a new y-particle created at the location of an x-particle; the y-particle then moves in a random direction that is uniformly distributed over the sphere. Compared to previous models of particle creation and annihilation in Bohmian mechanics (see [14] and references therein) that involved a UV cut-off, the difference is that in Models 1 and 2, the y-particle gets created at (rather than near ) an x-particle, and that annihilation is deterministic. The process is time-reversal invariant (notwithstanding that annihilation is deterministic and creation is stochastic).

IBC for Model 1
We now describe an IBC and the corresponding Hamiltonian for Model 1. The IBC demands that for any m, n ∈ {1, 2, . . .}, any configuration x m = (x 1 , . . . , x m ) ∈ R 3m of x-particles, any configuration y n = (y 1 , . . . , y n ) ∈ R 3n of y-particles with x m ∩ y n = ∅ (i.e., x i = y j for all i, j), any i = 1, . . . , m, and any j = 1, . . . , n, where denotes omission, and g ∈ R is the same coupling constant as before.
The IBC is a condition on the wave function ψ near the diagonal ∆ in configuration space Q x × Q y , i.e., the set of "collision configurations," If we regard the collision configurations in ∆ as not admissible configurations then the configuration space is the set difference Q = (Q x × Q y ) \ ∆, and its "boundary" is ∂Q = ∆. The Hamiltonian H IBC is defined [6] on a domain D IBC ⊂ H = F − ⊗ F + consisting of wave functions ψ that satisfy the IBC (54) and obey the following asymptotics, replacing (30), near the boundary surface {x i = y j } in Q (m,n) = R 3m x × R 3n y : where r ij = |x i − y j |, x = (x 1 , . . . , x m ), and y = (y 1 , . . . , y n ) with x k = y r for all k, r. Moreover, in order to make H IBC self-adjoint, the coefficients are related according to [6] with fixed real proportionality factor Note that as m x → ∞, the constant α tends to zero, and by (57) the asymptotics (56) reduces to (30), the one of Model 2. In terms of the c coefficients, the IBC (54) can be reformulated as The Hamiltonian H IBC is given by As shown in [6], H IBC defines a self-adjoint operator on D IBC . Of the five terms in (60), only the nE 0 ψ term is square integrable by itself. All the other terms are distributions that add up to a square integrable function. In particular, in contrast to the Hamiltonian (29) of Model 2, the annihilation term c 0,i (x m , y n ) is not necessarily square integrable but compensates derivatives of the logarithmic contributions in (56). These cancelations are rather subtle, which could be a reason why no renormalization scheme for Model 1 was known before the IBC study in [6].

Ground State Energy and Effective Yukawa Potential
Elsewhere [4], we show that H IBC as in (29) for Model 2 with Dirichlet-type IBC (27) and E 0 > 0 possesses a non-degenerate ground state ψ min . It is given by with normalization constant and eigenvalue That is, the state is a superposition of different numbers of y-particles, and in each sector all y-particles have the same wave function; so the x-particle at the origin is dressed with a cloud of y-particles. The probability distribution of the number n of y-particles is a Poisson distribution with mean value We claim further that in Model 1, in which the x-particles interact by exchanging y-particles, the x-particles effectively interact through a Yukawa potential; this in fact agrees with the result of Yukawa's original reasoning [36], see below. A simple way of computing the effective interaction potential is to consider (as "Model 2b") N x-particles fixed at x 1 , . . . , x N ∈ R 3 and to find the ground state energy. To this end, consider wave functions of y-configurations that simultaneously satisfy N IBCs, for every i = 1, . . . , N , ω ∈ S 2 , n ∈ {0, 1, 2, . . .}, and y n ∈ (R 3 \ {x 1 , . . . , x N }) n . The corresponding Hamiltonian reads For E 0 > 0, the ground state is with suitable factors c n and eigenvalue That is, the ground state energy of the y-particles, given the x-particles at x 1 , . . . , x N , is given by (68). Regarding this energy function of x 1 , . . . , x N as an effective potential for the x-particles (which is appropriate when the x-particles move slowly), we see that x-particles effectively interact through an attractive Yukawa pair potential, with R the distance between two x-particles. If we take the energy needed to create a y-particle to be E 0 = m y c 2 , then which is, up to the factor √ 2, the value originally obtained by Yukawa [36] considering the effective interaction of nucleons by exchange of pions. We expect that the factor √ 2 is owed to the non-relativistic nature of our model.

IBC Hamiltonians as a Limit of Removing the Cut-Off
If one introduces a UV cut-off into the UV divergent original Hamiltonian (26) of Model 2, it becomes the well-defined operator H ϕ given by (H ϕ ψ) (n) (y n ) = − 2 2m y n j=1 ∇ 2 y j ψ (n) (y n ) + nE 0 ψ (n) (y n ) + g √ n + 1 R 3 d 3 y ϕ * (y) ψ (n+1) (y n , y) analogous to (20) for Model 1. Here, the Dirac delta function δ 3 has been replaced by ϕ : R 3 → C, a square-integrable function describing the charge density of the x-particle, and the limit ϕ → δ 3 would correspond to removing the UV cut-off. It is well known [28,29] that, if E 0 > 0, there are numbers E ϕ ∈ R such that the operator H ϕ − E ϕ possesses a limit as ϕ → δ 3 . Since E ϕ → ∞ in this limit, H ϕ tends to infinity in a sense, but this sense is harmless because two Hamiltonians that differ only by a multiple of the identity operator can be regarded as equivalent, as they generate the same time evolution (if we regard wave functions that differ only by a global phase factor as equivalent). So, in a relevant sense, the limit ϕ → δ 3 can indeed be taken, which suggests regarding the limiting Hamiltonian as the "physically correct" Hamiltonian. Obviously, the limiting Hamiltonian is defined only up to addition of a constant, as this constant can be added to each of the E ϕ . However, for the H ϕ of this concrete model, there is a particular natural choice of E ϕ and thus of the limiting Hamiltonian, described, e.g., in [29,4], for which we will write H ∞ .

Theorem 2. [4]
For the Hamiltonian H IBC of the Dirichlet type for Model 2, defined by (27) and (29), with the constant E min as in (63).
On the one hand, this result lends further support to regarding H IBC and the IBC approach as physically reasonable. On the other hand, the result provides a more direct and explicit representation of H ∞ than was available so far.

Conclusions
If configuration space has a boundary (be it of codimension 1 or codimension 3), then an interior-boundary condition on this boundary can serve to ensure that the probability flux into the boundary is compensated by an equally large gain of probability in another place in configuration space, and thus to ensure the overall conservation of probability. At the same time, an interior-boundary condition has the consequence that contributions to the wave function can flow out of the boundary. That is just what is needed for defining an evolution on Fock space representing the emission and absorption of particles. The crucial point here is that this approach is also possible if the sources (the x-particles in Model 1 and Model 2) are point-shaped, without a UV divergence problem arising. IBCs thus seem like a natural way of implementing Hamiltonians with particle creation and annihilation. The mathematical viability of this approach for nonrelativistic QFTs has been established in recent works [4,5,6]. Its applicability to more serious QFTs remains to be explored.

Conflict of interest.
On behalf of all authors, the corresponding author states that there is no conflict of interest.