# Asymptotic Principal Values and Regularization Methods for Correlation Functions with Reflective Boundary Conditions

Authors: Masafumi Seriu
Journal: Advances in Optical Technologies
 (2010)

---

## Abstract

We introduce a concept of asymptotic principal values which enables us to handle rigorously singular integrals of higher-order poles encountered in the computation of various quantities based
on correlation functions of a vacuum. Several theorems on asymptotic principal values are proved, and they are expected to become bases for investigating and developing some classes of regularization methods for singular integrals.
We make use of these theorems for analyzing mutual relations between some regularization methods,
including a method naturally derived from asymptotic principal values. It turns out that the
concept of asymptotic principal values and the theorems for them are quite useful in this type of
analysis, providing a suitable language to describe what is discarded and what is retained in each
regularization method.


---

## Body

## 1. Introduction

Physics of quantum vacuum fluctuations is one of the intriguing research topics expected to be developed through the interplay between theories, experiments, and practical applications.Investigations of quantum vacuum fluctuations even stimulate the border area between physics and mathematics. As a typical example of this sort, we often encounter singular integrals in computing several quantities based on correlation functions of a vacuum in question. The occurrence of singularity or divergence is often a signal of surpassing the border of validity of a model by too much extrapolation. Furthermore, it could originate from deeper physical processes for which satisfactory consistent mathematics is still unavailable. How to handle singular integrals can be then a challenging topic, requiring both mathematical analysis and physical considerations.Faced with singular integrals, we need to resort to some regularization method to get a finite result. The aim of this paper is to give an organized mathematical basis underlying some typical regularization methods and make clear their mutual relations. We introduce below a concept ofasymptotic principal values which can be a key tool to analyze some classes of regularization methods. We then prove several theorems on the asymptotic principal values useful for studying regularization procedures.There are still various uncertainties to clear up in regularization methods, reflecting our lack of mathematical basis for handling infinities. In this situation, we cannot expect any universal regularization method, but we need to customize the method by try and error depending on the problem in question. It is far from the aim of this paper to judge which method is better than the others. It just tries to present a concrete mathematical basis for further considerations and developments of better regularization methods.The organization of the paper is as follows. In Section2, we present one simple example where singular integrals of a higher-order pole emerge. The origin of singularities in this example is physically clear and we can get some idea on how these integrals should be regularized. In Section 3, we introduce asymptotic principal values which describe precisely how singular integrals behave, which can be useful for investigating and developing several regularization methods. We then prove several useful theorems on the asymptotic principal values. In Section 4, we analyze some typical regularization methods by means of theorems prepared in the previous section. Section 5 is devoted for a summary and several discussions.


## 2. Typical Example of Singular Integrals

Let us consider as an example the measurement of the electromagnetic vacuum fluctuations in a half-space bounded by a perfectly reflecting infinite mirror. Recently the switching effect [1] and the smearing effect due to the quantum spread of a probe particle [2] have been analyzed by studying the measurement process of the Brownian particle released in this environment.We thus take a model introduced in [3] and reanalyzed in [1, 2]. Suppose that a flat, infinitely spreading mirror of perfect reflectivity is placed on the xy-plane (z=0). Then let us investigate the quantum vacuum fluctuations of the electromagnetic field inside the half-space z>0 by releasing a classical charged probe-particle with mass m and charge e in the environment. We can estimate the quantum fluctuations of the vacuum through the velocity dispersions of the probe-particle released in the environment.When the velocity of the probe-particle is much smaller than the light velocityc, the motion for the particle is described by(1)mdv⃗dt=eE⃗(x⃗,t),
where E⃗(x⃗,t) is the electric field.Within the time-period when the particle does not move so much, (1) along with the initial condition v⃗(0)=v⃗0 is solved approximately as(2)v⃗(t)≃v⃗0+em∫0tE⃗(x⃗,t′)dt′.For simplicity, let us consider only the “sudden-switching” case; the measurement is switched on abruptly, stably continued forτ [sec] before switched off abruptly. It is mathematically described by a step-like switching function without any switching tails. The velocity dispersions of the particle, 〈Δvi2〉(i=x,y,z), are then given by(3)〈Δvi2(x⃗,τ)〉=e2m2∫0τdt′∫0τdt′′〈Ei(x⃗,t′)Ei(x⃗,t′′)〉R,
by noting that 〈Ei(x⃗,t)〉R = 0. Here 〈Ei(x⃗,t′) Ei(x⃗,t′′)〉R  (i=x,y,z) are the renormalized two-point correlation functions of the electric field (the suffix “R” is for “renormalized”). Now 〈Ei(x⃗,t′)Ei(x⃗,t′′)〉R(i=x,y,z) are computed [4] as(4)〈Ez(x⃗,t′)Ez(x⃗,t′′)〉R=1π21(T2-(2z)2)2,〈Ex(x⃗,t′)Ex(x⃗,t′′)〉R=〈Ey(x⃗,t′)Ey(x⃗,t′′)〉R=-1π2T2+4z2(T2-(2z)2)3,
where T:=t′-t′′. (We set c=ℏ=1 hereafter throughout the paper.)It is obvious that the integral in (3) is regular when τ<2z, but singular when τ>2z, reflecting the singularity at |T|=2z inherent in the correlation functions 〈Ei(x⃗,t′)Ei(x⃗,t′′)〉R given in (4).For the present purpose, it suffices to show only the result of〈Δvz2〉 for τ>2z [1](5)〈Δvz2〉=e232π2m2{τz3ln(τ+2zτ-2z)2+8(1-2z/τ)z2ρ+O(ρ)}(6)~e24π2m2z2(1+1ρ)(forτ≫2z),
where ρ(>0) is a dimension-free asymptotic parameter for handling the singular integral properly (see Section 3 for details). Accordingly the above expression should be understood as an asymptotic expression as ρ~0.This result is derived by a formula for anasymptotic principal value, the rigorous definition of which shall be given in the next section,(7)℘(ρ)∫01dx1-x(x2-σ2)2=18σ3ln(1+σ1-σ)2+1-σ2σ2ρ+O(ρ),
for 0<σ<1. One can derive (7) with the help of  Theorem 1; a more direct derivation is also found in AppendixC of [1].Leaving the rigorous treatment of singular integrals for the next section, we here focus on the physical reason why the singularity of correlation functions occurs at|T|=2z. Due to the mirror reflections of signals with the light velocity, the values of the electric field at the two world-points (t′,x,y,z) and (t′′,x,y,z) are expected to be strongly correlated when |t′-t′′|=2z. When the measuring time τ is short enough (shorter than the travel-time of the signal 2z), then it always follows |t′-t′′|<2z, so that these correlations are not captured by the probe. When the measuring time is long enough (τ>2z), however, these strong correlations accumulate in the velocity fluctuations of the particle at z. Therefore it is expected that the resulting singular term of the form A/ρ (A>0) contains information on the reflecting boundary.On the other hand, typical regularization procedures [5] correspond to discarding such a singular term (e.g., the 1/ρ term in (6)) in effect. It should be clarified when this type of regularization is valid and when not. We shall discuss on this point in more detail in Section 4.It turns out that the model given here is too simplified and should be modified taking into account the switching effect [1] and the smearing effect due to the quantum spread of the probe-particle [2]. However, it suffices for the present purpose of giving some example of singular integrals.


## 3. Basic Formulas for Handling Singular Integrals

In view of the example in the previous section, it is clear that we sometimes need to estimate a singular integral whose integrand possesses a higher-order pole. In order to investigate various regularization methods later, we first need some concrete quantity corresponding to a singular integral for which all the information is retained and nothing is discarded. Then, the following asymptotic definition of a singular integral may be relevant.Definition 1.
Letf(x) be an arbitrary real function defined around an interval [a,b], differentiable at x=c (a<c<b) sufficiently many times. For a positive integer n, then, let us introduce an asymptotic principal value of order ρ defined by
(8)℘(c,ρ)(f,n):=℘(ρ)∫abf(x)(x-c)ndx:={∫ac-ρ+∫c+ρb}f(x)(x-c)ndx,
where ρ is a sufficiently small positive parameter.The asymptotic principal value is a generalization of the standard Cauchy principal value, corresponding tolimρ↓0℘(c,ρ)(f,1), in two ways. First, the order of singularity n can be greater than 1. Second, only the asymptotic behavior as ρ~0 is concerned and the convergence for the limit ρ→0 is not necessarily required. In other words, we focus on how the integral behaves near ρ~0 rather than the ρ→0 limit itself. In this sense, all the information is retained and no infinities are discarded in defining the asymptotic principal value.Let us now introduce another asymptotic quantity.Definition 2.
With the same premises as inDefinition1, we define
(9)℘(ρ)[f(x)(x-c)n]ab:=[f(x)(x-c)n]ac-ρ+[f(x)(x-c)n]c+ρb.It is easily shown that(10)℘(ρ)[f(x)(x-c)n]ab=[f(x)(x-c)n]ab-1ρn{f}(n)(c,ρ)
with
(11){f}(n)(c,ρ):=f(c+ρ)-(-)nf(c-ρ).By a Taylor-expansion inρ, it is obvious that {f}(n)(c,ρ) is O(ρ) for an even n (when f′(c)≠0), so that the term (1/ρn){f}(n)(c,ρ) in (10) is O(1/ρn-1) and singular in the ρ→0 limit. Similarly, {f}(n)(c,ρ) is O(1) for an odd n (when f(c)≠0), so that (1/ρn){f}(n)(c,ρ) is O(1/ρn) and singular as ρ→0.The following definition is just for making formulas below concise.Definition 3.
We have(12)ζ(c,ρ)(f,n):=1n℘(ρ)[f(x)(x-c)n]ab.With these preparations, let us start with the following lemma.Lemma 1.
For any functionf(x) differentiable at x=c and for any integer n(≥2), it follows that
(13)℘(c,ρ)(f,n)=1n-1℘(c,ρ)(f′,n-1)-ζ(c,ρ)(f,n-1).Proof.
Noting that(14)1(x-c)n=-1n-1(1(x-c)n-1)′,
we have
(15)℘(c,ρ)(f,n):=-1n-1℘(ρ)∫ab(1(x-c)n-1)′f(x)dx=-1n-1{℘(ρ)[f(x)(x-c)n-1]ab-℘(ρ)∫abf′(x)(x-c)n-1dx},
where the partial-integral has been performed to get the last line. Then the equality follows.We now prove a formula which relates a multipole integral with a simple-pole integral.Theorem 1.
For any functionf(x) differentiable sufficiently many times at x=c and for any integer n (≥2), it follows that
(16)℘(c,ρ)(f,n)=1(n-1)!℘(c,ρ)(f(n-1),1)-∑k=1n-1(n-k)!(n-1)!ζ(c,ρ)(f(k-1),n-k).Proof.
(1∘) For n=2, the claimed equality reduces to
(17)℘(c,ρ)(f,2)=℘(c,ρ)(f′,1)-ζ(c,ρ)(f,1),
which clearly holds due to  Lemma 1.
(2∘) Let us assume that the equality holds for a function F(x) and for n=m (m≥2), that is,
(18)℘(c,ρ)(F,m)=1(m-1)!℘(c,ρ)(F(m-1),1)-∑k=1m-1(m-k)!(m-1)!ζ(c,ρ)(F(k-1),m-k).
Now applying  Lemma 1 for n=m+1, we have
(19)℘(c,ρ)(f,m+1)=1m℘(c,ρ)(f′,m)-ζ(c,ρ)(f,m)=1m{∑k=1m-11(m-1)!℘(c,ρ)(f(m),1)-∑k=1m-1(m-k)!(m-1)!ζ(c,ρ)(f(k),m-k)}-ζ(c,ρ)(f,m),
where the assumed equation (18) for F(x)=f′(x) has been used to get the second equality. Rearranging the summation, the last equality reduces to
(20)℘(c,ρ)(f,m+1)=1m!℘(c,ρ)(f(m),1)-∑k=1m((m+1)-k)!m!ζ(c,ρ)(f(k-1),(m+1)-k).
Thus the claimed equation (16) holds for n=m+1.
(3∘) By the mathematical induction, (16) holds for any integer n (n≥2).Based on  Theorem1, it is natural to introduce a quantity ℘̃(c,ρ)(f,n), which is a simple-pole part plus a regular part of ℘(c,ρ)(f,n), putting aside singular contributions from higher-order poles.Definition 4.
We define(21)℘̃(c,ρ)(f,n):=1(n-1)!℘(c,ρ)(f(n-1),1)-∑k=1n-1(n-k-1)!(n-1)![f(k-1)(x)(x-c)n-k]ab.The “mild part”℘̃(c,ρ)(f,n) of ℘(c,ρ)(f,n) shall be important in the discussion of regularization methods in Section 4.We now have a formula which enables us to separate singular contributions from a multipole integral.Theorem 2.
For any functionf(x) differentiable sufficiently many times at x=c and for any integer n (≥2), it follows that
(22)℘(c,ρ)(f,n)=℘̃(c,ρ)(f,n)+∑k=1n-1(n-k-1)!(n-1)!1ρn-k{f(k-1)}(n-k)(c,ρ).Proof.
It is straightforward to show this formula due to  Theorem1 along with (12) and (10).Lemma 2.
For any functionf(x) differentiable at x=c and for a positive integer n, it follows that
(23)℘(c,ρ)(f,n+1)=1n∂c℘(c,ρ)(f,n)+1nρn{f}n(c,ρ).Proof.
We compute∂c℘(c,ρ)(f,n) directly as
(24)∂c℘(c,ρ)(f,n)=∂c({∫ac-ρ+∫c+ρb}f(x)(x-c)ndx)=n℘(c,ρ)(f,n+1)-1ρn{f}n(c,ρ),
where the second term in the last line comes from the c-derivative applied to the upper and the lower limit of the integral region. Thus the claimed equation follows.Theorem 3.
For any functionf(x) differentiable sufficiently many times at x=c and for an integer n (≥2), it follows that
(25)℘̃(c,ρ)(f,n)=1(n-1)!∂cn-1℘(c,ρ)(f,1).Proof.
Due to  Theorem2, the claimed equation (25) is equivalent to
(26)℘(c,ρ)(f,n)=1(n-1)!∂cn-1℘(c,ρ)(f,1)+∑k=1n-1(n-k-1)!(n-1)!1ρn-k{f(k-1)}(n-k)(c,ρ).
Thus it suffices to show (26).
(1∘) Let us consider the case n=2, where the R.H.S. (right-hand side) of (26) becomes
(27)∂∂c({∫ac-ρ+∫c+ρb}f(x)x-cdx)+1ρ{f(c-ρ)+f(c+ρ)}.
In this expression, the c-derivative applied to the upper and the lower limit of the integral region yields a term which exactly cancels the second term. As a result, the above expression reduces to ℘(c,ρ)(f,2), that is, the L.H.S. (left-hand side) of (26). Thus (26) holds for n=2.
(2∘) Let us now assume that (26) holds for n=m (m≥2), that is,
(28)℘(c,ρ)(f,m)=1(m-1)!∂cm-1℘(c,ρ)(f,1)+∑k=1m-1(m-k-1)!(m-1)!1ρm-k{f(k-1)}(m-k)(c,ρ).
Due to Lemma 2, then, it becomes
(29)℘(c,ρ)(f,m+1)=1m∂c℘(c,ρ)(f,m)+1mρm{f}m(c,ρ)=1m∂c{∑k=1m-11(m-1)!∂cm-1℘(c,ρ)(f,1)+∑k=1m-1(m-k-1)!(m-1)!1ρm-k{f(k-1)}(m-k)(c,ρ)}+1mρm{f}m(c,ρ),
where (28) has been used to get the last line. Noting that the relation
(30)∂c{f}m(c,ρ)={f′}m(c,ρ),
which obviously holds from (11), it reduces to
(31)℘(c,ρ)(f,m+1)=1m!∂cm℘(c,ρ)(f,1)+∑k=1m(m-k)!m!1ρm-k+1{f(k-1)}(m-k+1)(c,ρ).
Thus (26) holds for n=m+1.
(3∘) By the mathematical induction, (26) holds for any integer n (n≥2). Thus the claimed formula (25) has been shown.


## 4. Typical Regularization Methods and Their Mutual Relations

Based on the results in the previous section, let us now come back to the problem of regularization methods for singular integrals.Let us consider a typical singular integral(32)I=∫abf(x)(x-c)ndx
for any function f(x) differentiable sufficiently many times at x=c (a<c<b) and for a positive integer n.


### 4.1. Regularization Method with Partial Integrals

The first method of regularization we consider is amethod of partial integrals which is sometimes made used of. We insert an identity(33)1(x-c)n=-1n-1(1(x-c)n-1)′,
into (32) and formally perform a partial integral(34)I=-1n-1∫ab(1(x-c)n-1)′f(x)dx=1n-1{∫abf′(x)(x-c)n-1dx-[f(x)(x-c)n-1]ab}.
In this way, the order of singularity is reduced by one. Repeating the similar procedure, the integral I is reduced to the n=1 case for which the prescription of the Cauchy principal value may be applied.Due to  Theorem2, however, it is obvious that singular terms should exist and should have been discarded by hand in the above procedure. Indeed, compared with (34) with the rigorous expression (15), it is obvious that the singularities which should reside in the second term on the R.H.S. of (34) are simply discarded by hand. Thus, in view of  Theorem 2, the above method is equivalent to the simple replacement of I as(35)I↦℘̃(c,ρ)(f,n).There is still room, however, to regard the method of partial integrals as a shorthand prescription of what we here call themethod of infinitesimal imaginary part, which is much more of theoretical grounds [5]. We shall consider this method in the next subsection.


### 4.2. Regularization Method with Infinitesimal Imaginary Part

The method of infinitesimal imaginary part is based on well-known Dirac’s formula [6] for an integral kernel
(36)1(x-c)±iρ=℘(c,ρ)1x-c∓iπδ(x-c),
which is most easily shown by estimating an integral ∫abdxf(x)/((x-c)±iρ) by means of an appropriate contour-integral for a suitable function f(x).By differentiating the both-sides of (36) n-1 times with respect to c, and by applying  Theorem 3, we get(37)1((x-c)±iρ)n=℘̃(ρ,c)(·,n)∓iπ(n-1)!δ(n-1)(x-c)
in the sense of an integral kernel. Recalling Definition4, however, we see that the R.H.S. is reduced to the n=1 case, for which the prescription of the Cauchy principal value may be applied.It is notable that just the introduction of some infinitesimal imaginary part results in a tamable quantity such as℘̃(ρ,c)(·,n) at the cost of the imaginary contribution of the second term on the R.H.S. of (37). Thus along with some causality arguments [5], it is often argued that the singular integral I should be interpreted as the real part of ∫(f(x)/((x-c)±iρ)n)dx, that is,(38)I↦ℜ∫f(x)((x-c)±iρ)ndx.As far as one is evaluatingreal quantities, one may further argue that the second term on the R.H.S. of (37) shall not contribute. If so, the procedure is in effect equivalent to the replacement (35). In this sense, the method of partial integrals discussed in the previous subsection may be justified provided that it is regarded as a shorthand prescription of the method of infinitesimal imaginary part.Another way of looking at this method is to pay attention to the L.H.S. (rather than the R.H.S.) of (37). As far as computations of real quantities are concerned, then, this method is equivalent to the replacement 1/(x-c)n with ((x-c)/((x-c)2+ρ2))n along with taking the limit ρ→0 after evaluating the integral(39)I↦limρ→0∫abf(x)(x-c(x-c)2+ρ2)ndx.There is some subtle points in this method. One of them is to discard the imaginary part of the R.H.S. of (37) on the grounds that one is evaluating real quantities. Considering that the regularization has been achieved at the cost of introducing the imaginary part though tiny, the imaginary part should carry important information and some concern naturally arises whether one can discard it so freely.Indeed, a simple example can be presented for which this kind of procedure fails. Let us consider an integralI1=∫-11dx which is purposefully regarded as(40)I1=∫-11x·1xdx.
It is obvious that I1=2. The analysis by the asymptotic principal value (see the next subsection) also results in I1(ρ)→2 in the limit ρ→0. This is because all the information is retained in the prescription of the asymptotic principal value.On the other hand, the above-mentioned scheme makes a replacement(41)1x↦xx2+ρ2=12(ln(x2+ρ2))′,
so that
(42)I1→I1(ρ):=∫-11x·12(ln(x2+ρ2))′dx=ln(1+ρ2)-∫01ln(x2+ρ2)dx,
where a partial integral has been performed to get the last line. However, it is clear that I1(ρ)→∞ as ρ→0, contradicting with the obvious result I1=2.Quite interestingly, no contradiction occurs forIm=∫-11xm·(1/x)dx with m≥2 since the second term in (42) becomes -m∫01xm-1ln(x2+ρ2)dx so that no singularity occurs around x~0 for m≥2. More generally, the integral of the form ∫-11(f(x)/x)dx, if treated by the above prescription, gives rise to the dominant contribution -f′(0)∫01ln(x2+ρ2)dx which diverges as ρ→0.With these caveats in mind, let us now move to a new regularization method based on the asymptotic principal values.


### 4.3. Regularization Method with Asymptotic Principal Values

Let us finally introduce a new regularization method based on the asymptotic principal values.For the simple example in Section2, there has been a definite physical interpretation of the singularity in the correlation function. Furthermore, the system considered there has been a combination of quantum objects with a macroscopic mirror. Therefore it might be also probable that the deepest cause of the singularity resides in the validity issue of the model originating from too much extrapolation from the quantum side to the macroscopic situation. Indeed there is an investigation showing that the quantum fluctuations of the mirror boundary drastically decrease the singular behavior near the mirror [7]. Therefore it is reasonable to take the origin of the singularity more realistically (rather than just mathematical phenomenon), expecting that some physical processes suppress the order of singularity.Going back to the example of the integralI in (32), then, it is possible to interpret I in the sense of an asymptotic principal value,
(43)I↦℘(c,ρ)(f,n)
with the dimension-free parameter ρ being provided by the ratio of some natural cut-off scale with the system-size in question. (E.g., the ratio of the plasma wave-length of the mirror with 2z for the example in Section 2). The advantage of this regularization scheme is that one can explicitly analyze the ρ-dependence of the integral. For instance, one may study the influence of the quantum fluctuations of the mirror by treating ρ as a fluctuation parameter. The result for 〈Δvz2〉 given in (6) along with (7) is an example of the computation by the method of asymptotic principal values.We see that  Theorem2 is the basis for understanding the relation between the regularization methods discussed so far. The difference between the method of asymptotic principal value (℘(c,ρ)(f,n)) and the method of infinitesimal imaginary part (℘̃(c,ρ)(f,n)) is given by the second term on the R.H.S. of (22), which is of O(1/ρn-1).


## 4.1. Regularization Method with Partial Integrals

The first method of regularization we consider is amethod of partial integrals which is sometimes made used of. We insert an identity(33)1(x-c)n=-1n-1(1(x-c)n-1)′,
into (32) and formally perform a partial integral(34)I=-1n-1∫ab(1(x-c)n-1)′f(x)dx=1n-1{∫abf′(x)(x-c)n-1dx-[f(x)(x-c)n-1]ab}.
In this way, the order of singularity is reduced by one. Repeating the similar procedure, the integral I is reduced to the n=1 case for which the prescription of the Cauchy principal value may be applied.Due to  Theorem2, however, it is obvious that singular terms should exist and should have been discarded by hand in the above procedure. Indeed, compared with (34) with the rigorous expression (15), it is obvious that the singularities which should reside in the second term on the R.H.S. of (34) are simply discarded by hand. Thus, in view of  Theorem 2, the above method is equivalent to the simple replacement of I as(35)I↦℘̃(c,ρ)(f,n).There is still room, however, to regard the method of partial integrals as a shorthand prescription of what we here call themethod of infinitesimal imaginary part, which is much more of theoretical grounds [5]. We shall consider this method in the next subsection.


## 4.2. Regularization Method with Infinitesimal Imaginary Part

The method of infinitesimal imaginary part is based on well-known Dirac’s formula [6] for an integral kernel
(36)1(x-c)±iρ=℘(c,ρ)1x-c∓iπδ(x-c),
which is most easily shown by estimating an integral ∫abdxf(x)/((x-c)±iρ) by means of an appropriate contour-integral for a suitable function f(x).By differentiating the both-sides of (36) n-1 times with respect to c, and by applying  Theorem 3, we get(37)1((x-c)±iρ)n=℘̃(ρ,c)(·,n)∓iπ(n-1)!δ(n-1)(x-c)
in the sense of an integral kernel. Recalling Definition4, however, we see that the R.H.S. is reduced to the n=1 case, for which the prescription of the Cauchy principal value may be applied.It is notable that just the introduction of some infinitesimal imaginary part results in a tamable quantity such as℘̃(ρ,c)(·,n) at the cost of the imaginary contribution of the second term on the R.H.S. of (37). Thus along with some causality arguments [5], it is often argued that the singular integral I should be interpreted as the real part of ∫(f(x)/((x-c)±iρ)n)dx, that is,(38)I↦ℜ∫f(x)((x-c)±iρ)ndx.As far as one is evaluatingreal quantities, one may further argue that the second term on the R.H.S. of (37) shall not contribute. If so, the procedure is in effect equivalent to the replacement (35). In this sense, the method of partial integrals discussed in the previous subsection may be justified provided that it is regarded as a shorthand prescription of the method of infinitesimal imaginary part.Another way of looking at this method is to pay attention to the L.H.S. (rather than the R.H.S.) of (37). As far as computations of real quantities are concerned, then, this method is equivalent to the replacement 1/(x-c)n with ((x-c)/((x-c)2+ρ2))n along with taking the limit ρ→0 after evaluating the integral(39)I↦limρ→0∫abf(x)(x-c(x-c)2+ρ2)ndx.There is some subtle points in this method. One of them is to discard the imaginary part of the R.H.S. of (37) on the grounds that one is evaluating real quantities. Considering that the regularization has been achieved at the cost of introducing the imaginary part though tiny, the imaginary part should carry important information and some concern naturally arises whether one can discard it so freely.Indeed, a simple example can be presented for which this kind of procedure fails. Let us consider an integralI1=∫-11dx which is purposefully regarded as(40)I1=∫-11x·1xdx.
It is obvious that I1=2. The analysis by the asymptotic principal value (see the next subsection) also results in I1(ρ)→2 in the limit ρ→0. This is because all the information is retained in the prescription of the asymptotic principal value.On the other hand, the above-mentioned scheme makes a replacement(41)1x↦xx2+ρ2=12(ln(x2+ρ2))′,
so that
(42)I1→I1(ρ):=∫-11x·12(ln(x2+ρ2))′dx=ln(1+ρ2)-∫01ln(x2+ρ2)dx,
where a partial integral has been performed to get the last line. However, it is clear that I1(ρ)→∞ as ρ→0, contradicting with the obvious result I1=2.Quite interestingly, no contradiction occurs forIm=∫-11xm·(1/x)dx with m≥2 since the second term in (42) becomes -m∫01xm-1ln(x2+ρ2)dx so that no singularity occurs around x~0 for m≥2. More generally, the integral of the form ∫-11(f(x)/x)dx, if treated by the above prescription, gives rise to the dominant contribution -f′(0)∫01ln(x2+ρ2)dx which diverges as ρ→0.With these caveats in mind, let us now move to a new regularization method based on the asymptotic principal values.


## 4.3. Regularization Method with Asymptotic Principal Values

Let us finally introduce a new regularization method based on the asymptotic principal values.For the simple example in Section2, there has been a definite physical interpretation of the singularity in the correlation function. Furthermore, the system considered there has been a combination of quantum objects with a macroscopic mirror. Therefore it might be also probable that the deepest cause of the singularity resides in the validity issue of the model originating from too much extrapolation from the quantum side to the macroscopic situation. Indeed there is an investigation showing that the quantum fluctuations of the mirror boundary drastically decrease the singular behavior near the mirror [7]. Therefore it is reasonable to take the origin of the singularity more realistically (rather than just mathematical phenomenon), expecting that some physical processes suppress the order of singularity.Going back to the example of the integralI in (32), then, it is possible to interpret I in the sense of an asymptotic principal value,
(43)I↦℘(c,ρ)(f,n)
with the dimension-free parameter ρ being provided by the ratio of some natural cut-off scale with the system-size in question. (E.g., the ratio of the plasma wave-length of the mirror with 2z for the example in Section 2). The advantage of this regularization scheme is that one can explicitly analyze the ρ-dependence of the integral. For instance, one may study the influence of the quantum fluctuations of the mirror by treating ρ as a fluctuation parameter. The result for 〈Δvz2〉 given in (6) along with (7) is an example of the computation by the method of asymptotic principal values.We see that  Theorem2 is the basis for understanding the relation between the regularization methods discussed so far. The difference between the method of asymptotic principal value (℘(c,ρ)(f,n)) and the method of infinitesimal imaginary part (℘̃(c,ρ)(f,n)) is given by the second term on the R.H.S. of (22), which is of O(1/ρn-1).


## 5. Summary and Discussions

In this paper, we have focussed on singular integrals with a higher-order pole which frequently emerge in computing quantities based on two-point correlation functions of a vacuum.To deal with this type of singular integrals, we have introduced the concept ofasymptotic principal values. The asymptotic principal value of order ρ, which is a generalization of the Cauchy principal value, is defined by introducing a cut-off parameter ρ, focussing solely on the asymptotic behavior of the integral as ρ~0. In this sense, it is a rigorous object retaining all the information on the singular integral.We have then proved several theorems on asymptotic principal values which are expected to serve as bases for studying regularization methods for singular integrals.To see how asymptotic principal values can be made use of, we have selected three typical regularization methods and have analyzed their mutual relations with the help of theorems we have prepared. It has turned out that the concept of asymptotic principal values and related theorems are quite useful in this kind of analysis. Indeed, in terms of asymptotic principal values, it has been possible to describe without ambiguity what is discarded and what is retained in each regularization method.No universal regularization method is available so far and we need to carefully select or invent a suitable method depending on the problem in question. For instance, we recall the example in Section2, where velocity dispersion of the probe, 〈Δvz2〉, is sensitive to the regularization method. In particular, the result expected by the method of infinitesimal imaginary part (Section 4.2) (and the method of partial integrals (Section 4.1)) is(44)〈Δvz2〉~e24π2m2z2(forτ≫2z).
On the other hand, the result expected by the method of asymptotic principal values (Section 4.3) is
(45)〈Δvz2〉~e24π2m2z2(1+1ρ)(forτ≫2z),
by choosing ρ in the order of the ratio of plasma wave-length and the typical size 2z. Strictly speaking, the model in Section 2 is a too simplified one and should be modified taking into account the quantum spread of the probe-particle itself. Then the behavior of 〈Δvz2〉 at late time is corrected to a more reasonable one 〈Δvz2〉~1/τ2 rather than ~1/z2 [2].In any case it is significant to compare the results derived by different regularization methods in more detail for approaching to a more satisfactory mathematical theory of regularization procedures.



---

*Source: 798983-2010-04-13.xml*