This paper constructs and analyzes infinite-dimensional stochastic differential equations for Coulomb particle systems, proving existence, uniqueness, and connecting finite and infinite particle dynamics.
Contribution
It introduces a new stochastic analysis method for Coulomb systems, extending beyond determinantal point fields, and provides rigorous construction of Coulomb interacting Brownian motions.
Findings
01
Constructed strong solutions for Coulomb ISDEs in all dimensions d ≥ 2.
02
Proved pathwise uniqueness and reversibility of the infinite-particle dynamics.
03
Connected finite-particle systems with infinite-particle limits through approximation schemes.
Abstract
We study the infinite-dimensional stochastic differential equations (ISDEs) of infinite-particle systems associated with Coulomb random point fields. The stochastic dynamics described by these ISDEs are referred to as Coulomb interacting Brownian motions. In all spatial dimensions d≥2 and for all inverse temperatures β>0, we construct the Coulomb interacting Brownian motions. We prove that the ISDEs admit strong solutions and that pathwise uniqueness holds. The resulting labeled dynamics form an \RdN-valued diffusion, possibly without an invariant measure, while the corresponding unlabeled process is a reversible diffusion with respect to the underlying Coulomb random point field. Moreover, we identify the infinite-particle stochastic dynamics as the limit in path space of finite-particle systems driven by stochastic differential equations. This identification…
n→∞limsi∈SR∑Nn(∇Ψ(x−si)−∇xΨNn(x,si))=0 for some x∈SRℓ.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Full text
Infinite-dimensional stochastic differential equations for Coulomb random point fields
Hirofumi Osada and Shota Osada
Abstract.
We study the infinite-dimensional stochastic differential equations (ISDEs)
of infinite-particle systems associated with Coulomb random point fields.
The stochastic dynamics described by these ISDEs are referred to as
Coulomb interacting Brownian motions.
In all spatial dimensions d≥2 and for all inverse temperatures
β>0, we construct the Coulomb interacting Brownian motions.
We prove that the ISDEs admit strong solutions and that pathwise uniqueness holds.
The resulting labeled dynamics form an (Rd)N-valued diffusion, possibly without
an invariant measure, while the corresponding unlabeled process is a reversible
diffusion with respect to the underlying Coulomb random point field.
Moreover, we identify the infinite-particle stochastic dynamics as
the limit in path space of finite-particle systems driven by stochastic differential equations.
This identification is achieved through two approximation schemes: finite-domain
systems with reflecting boundary conditions and N-particle systems.
Although the N-particle approximation is more fundamental,
its justification relies crucially on the finite-domain approximation
together with the uniqueness of solutions to the ISDEs.
Previously, only the case d=2 and β=2, known as the Ginibre interacting Brownian motion,
was understood through random matrix theory and determinantal random point fields.
Extending this result beyond the determinantal setting has remained a major difficulty.
We introduce a new, conceptually clear method based on stochastic analysis of infinite-particle systems
with long-range interactions that yields a rigorous construction of Coulomb interacting Brownian motions.
A key ingredient is an explicit computation of the logarithmic derivatives
of Coulomb random point fields.
keywords: Coulomb random point fields; stochastic dynamics of infinite-particle systems;
logarithmic derivatives of random point fields; infinite-dimensional stochastic differential equations
The aim of this paper is to develop a theory of infinite-dimensional stochastic
differential equations (ISDEs) arising from long-range interacting particle systems.
We focus on the construction and characterization of strong solutions,
their pathwise uniqueness, and the relationship between labeled and unlabeled
dynamics within the framework of stochastic analysis on infinite-particle systems.
Moreover, we derive the limiting stochastic dynamics as the limit of
N-particle systems described by solutions of finite-dimensional SDEs.
The Coulomb potential occupies a central position in mathematics and the natural sciences.
For d≥2, define Ψ:Rd→R∪{∞} by
[TABLE]
Then the gradient of Ψ has the simple form
[TABLE]
We call Ψ the d-dimensional Coulomb potential.
The purpose of this paper is to construct the stochastic dynamics of
infinite-particle systems in Rd interacting through Ψ
at inverse temperature β>0.
These dynamics are expected to satisfy an ISDE, whose prototypical form is
[TABLE]
which we call the Coulomb interacting Brownian motion.
Equilibrium states of the unlabeled dynamics of Coulomb interacting Brownian motion
are given by Coulomb random point fields (RPFs).
Coulomb RPFs are probability measures describing
unlabeled particle systems with interaction potential Ψ.
They are obtained as limits of finite-particle systems whose labeled
distributions μN are given by the density
[TABLE]
where 0<β<∞, ΦN denotes a confining potential,
and ΨN an interaction potential at the N-particle level,
which may coincide with Ψ.
We assume that these potentials converge to limiting potentials:
[TABLE]
We refer to (A1) for the precise meaning of (1.4).
The N-particle dynamics
XN=(XN,i)i=1N
are gradient dynamics associated with the Dirichlet form
EN on L2(μN),
which is defined by
[TABLE]
Here we denote by (⋅,⋅)Rd the Euclidean inner product on Rd.
By integration by parts, we have
[TABLE]
See Subsection 13.2 for the definition of Dirichlet forms used here.
By (1.3),
[TABLE]
Consequently, XN satisfies the following SDE in differential form:
[TABLE]
Taking the limit N→∞, we formally obtain the limiting ISDE:
[TABLE]
If, in addition, μ is translation invariant, this equation reduces to the ISDE (1.2).
Such multiple representations of the stochastic dynamics arise from the
long-range nature of the Coulomb interaction.
A striking feature of this problem is the exact correspondence between the ambient spatial dimension d and the exponent in the Coulomb potential. This coincidence is not accidental: it is the source of both the depth and the difficulty of the problem. On the one hand, it ensures that Coulomb random point fields (RPFs) are natural equilibrium distributions in any dimension d≥2; on the other hand, it makes the analysis of the associated ISDE (1.6) extremely delicate, due to the long-range nature and singularity of the Coulomb force.
To date, the only case where (1.6) has been solved is d=β=2, corresponding to the Ginibre RPF. In this special situation, the determinantal structure of the equilibrium measure provides powerful tools, including sharp hyperuniform estimates for particle number fluctuations. Beyond d=β=2, however, such determinantal methods break down, and no general construction of the dynamics has been available.
From the equilibrium viewpoint, the classical Gibbsian framework excludes Coulomb interactions because the potential is neither superstable nor integrable at infinity. The construction of Coulomb RPFs remained open for decades, until very recently it was resolved in general by Armstrong–Serfaty and Thoma [1, 44]. These works establish the existence of Coulomb RPFs for all d≥2 and β>0. The present paper addresses the complementary problem: the construction of the corresponding stochastic dynamics, i.e., reversible diffusions with respect to these measures.
Our approach is new, conceptually simple, and robust.
For the first time, it constructs a stochastic dynamics
for infinite-particle Coulomb systems in all spatial dimensions d≥2
and for all inverse temperatures β>0 (Theorem 2.2).
Moreover, this dynamics is realized as a
pathwise unique strong solution
of the Coulomb ISDE (1.6)
with explicit coefficients
(Theorem 2.3).
This follows from the explicit representation of the logarithmic derivative of μ
(Theorem 2.1).
The central innovation of the method is the introduction of a new family of
analytic estimates for the logarithmic derivatives of Coulomb RPFs.
These estimates replace the determinantal and random matrix arguments that were
previously available only in the Ginibre case.
The solution X is constructed as an (Rd)N-valued diffusion process
without relying on an invariant or reference measure,
while the associated unlabeled process is shown to be a reversible diffusion
with respect to the Coulomb RPF (Theorem 2.4).
We derive the ISDE under two different thermodynamic limits.
One approach is based on finite-domain approximations, in which
finite-particle systems are connected through a two-step thermodynamic
limit (Theorem 12.1).
For each radius R, we freeze the exterior configuration on
SRc={∣s∣>R} and consider the N-particle
gradient dynamics inside SR={∣s∣<R} with normal reflection on ∂SR.
Denoting these finite-volume systems by XRN, we first take the limit
N→∞ with R fixed to obtain an infinite-particle
dynamics XR in SR, which is reversible with respect to the
corresponding conditional measure. We then let R→∞ and prove
that XR converges in law to X, where
X is the unique strong solution to the Coulomb ISDE (1.6).
The other approach is based on N-particle approximations.
In this approach, solutions XN of the finite-dimensional SDE
(1.5) converge in the path space to the solution X of the ISDE (Theorem 2.5).
More precisely, we regard XN as an (Rd)N-valued process
by introducing frozen dummy particles for each i>N,
which do not affect the dynamics.
Although the second derivation is the most fundamental, it relies
on the first derivation,
on the convergence of a related family of Dirichlet forms,
and on a sandwich-type convergence theorem built upon these results.
The uniqueness of solutions to the ISDE is indispensable
for implementing the sandwich-type convergence.
These results justify the ISDE and provide a simulation scheme for the limiting stochastic dynamics.
In [31], the vanishing of self-diffusion was proved for the Ginibre interacting Brownian motion.
The decisive mechanism behind this phenomenon is the number rigidity of the stationary RPF.
In contrast, in [25], it was proved that the self-diffusion coefficient is strictly positive for any potential in Ruelle’s class with convex hardcore in dimensions d≥2.
Taken together, these results highlight a structural distinction between Coulomb interactions and Ruelle-type interactions.
With the present results in place, it is natural to formulate the following conjecture:
Conjecture.
For translation-invariant Coulomb interacting Brownian motions in two dimensions, the self-diffusion coefficient vanishes at every inverse temperature β>0.
In contrast, for every d≥3, we conjecture that Coulomb interacting Brownian motions are always diffusive, i.e., their self-diffusion coefficient is strictly positive for all β>0.
In [45], Thoma proved number rigidity for Coulomb RPFs in the case
d=2, and established the absence of number rigidity for Coulomb RPFs
in all dimensions d≥3.
This dichotomy strongly suggests that the self-diffusion coefficient
is strictly positive in dimensions d≥3.
If confirmed, this phenomenon would vividly highlight the profound role
played by the long-range nature of the two-dimensional Coulomb potential,
and would firmly place Coulomb systems within the universality paradigm
of singular interacting diffusions.
Organization of the paper.
The remainder of the paper is organized as follows.
In Section 2, we state our main results and discuss related work and positioning.
Section 3 presents examples.
Sections 4 and 5
establish the existence of the logarithmic derivative of the equilibrium measure and provide its explicit representation.
Section 6 proves non-collision estimates for solutions of ISDEs.
Sections 7–9 introduce Dirichlet form techniques and construct weak solutions under both lower and upper schemes.
Section 10 establishes the existence of unique strong solutions, while Section 11 demonstrates the uniqueness of Dirichlet forms and completes the proof of our diffusion theorem.
Section 12 shows the convergence of finite-particle dynamics to the infinite-particle ISDE dynamics.
Finally, the Appendix (Section 13) collects auxiliary results on configuration spaces, Dirichlet forms, and the general theory of ISDEs.
2. Setup and main theorems
We now formulate the problem in detail and state our main results (Theorems 2.1–2.5).
Throughout the paper, we occasionally attach labels directly to constants
to facilitate later references.
These labels are independent of equation numbering and are used consistently
whenever such constants reappear.
Moreover, when an unlabeled configuration
s is written together with points si,
we implicitly use a labeled representation s=∑iδsi.
2.1. Set up and logarithmic derivatives
Let SR={s∈Rd;∣s∣<R}.
Define the configuration space on Rd by
[TABLE]
We equip S with the vague topology, under which it becomes a Polish space, i.e., homeomorphic to a complete and separable metric space.
A function f:S→R is called smooth if fR,sm is smooth on SRm for all R,m∈N and s,
where fR,sm is the SRm-representation of f defined in Definition 13.6.
It is called local if f is σ[πSR]-measurable for some R∈N, where
πA(s)=s(⋅∩A) for A⊂Rd.
We write
[TABLE]
A probability measure ν on (S,B(S)) is called an RPF on Rd.
For an RPF ν with the one-point correlation function ρν1, we define
the one-reduced Campbell measure ν[1] of ν by
[TABLE]
where νx=ν(⋅+δx∣s({x})≥1) is the reduced Palm measure of ν conditioned at x∈Rd.
We set Lloc1(ν[1])=⋂R=1∞L1(ν[1](⋅∩SR×S)). We refer to Definition 13.5 for the definition of correlation functions.
Let D∘b[1]:=C0∞(Rd)⊗D∘b and D∙b[1]:=C0∞(Rd)⊗D∙b.
We now quote the concept of the logarithmic derivative dν of ν [27].
Definition 2.1*.*
i
The logarithmic derivative dν of ν is an Rd-valued function such that
dν∈Lloc1(ν[1]) and that, for all f∈D∘b[1],
[TABLE]
ii
The ∙-logarithmic derivative dν of ν is the Rd-valued function defined similarly to the logarithmic derivative, but with D∘b[1] replaced by D∙b[1].
A ∙-logarithmic derivative is a logarithmic derivative because D∘b[1]⊂D∙b[1]. Both logarithmic derivatives coincide, which follows from Lemma 13.6.
Intuitively, the logarithmic derivative dν(x,s) represents the force received by the tagged particle at x from (infinitely many) other particles s=∑iδsi.
Informally, dν is defined as the differential of the logarithm of the Hamiltonian if the potential gives the equilibrium state ν.
One can easily calculate the logarithmic derivatives of canonical Gibbs measures with Ruelle’s class potentials using the DLR equation [27].
Even if the DLR equation does not hold, dν can be well defined—this is the case for the Coulomb potentials.
The logarithmic derivative dν of ν plays a vital role because, when it exists, it allows us to derive the stochastic dynamics of an infinite-particle system with the reversible equilibrium distribution ν. This connection leads to the formulation of solutions X=(Xi)i∈N through ISDEs, underscoring the power of this concept in understanding complex systems [27, 28, 36]:
[TABLE]
where {Bi} are independent standard Brownian motions.
2.2. Assumptions
We begin by introducing the assumptions used throughout the paper.
For ϵ>0, we define SR,ϵ2={(x,y)∈SR2;∣x−y∣>ϵ}.
(A1)
Let {ON} be an increasing sequence of open sets with ⋃NON=Rd.
i
Suppose Φ∈C2(Rd).
For each N,R∈N, let ΦN∈Cb2(ON) satisfy
[TABLE]
ii
Let ΨN:(Rd)2→R∪{∞} satisfy, for each R∈N,ϵ>0,
[TABLE]
where c2.1N is a positive constant, δy is the delta measure at y, and ψN∈Cb(ON2) satisfies, for some constant c2.2,
[TABLE]
(A2) i
The RPF μN admits the labeled density given in (1.3),
where β>0 is arbitrary and fixed throughout this paper.
ii
There exist an ℓ0>0 and a constant c2.3 such that,
[TABLE]
iii
There exists a weakly convergent subsequence of {μN}, denoted by {μNn}, with limit μ such that
This estimate is used in the construction of the logarithmic derivative of μ.
ii
Assumption (2.12) merely excludes the trivial finite-particle case.
iii
We do not assume a priori that μ admits a density or correlation functions;
these properties follow from the above assumptions (see Proposition 4.1).
For a function f(x), we set f(i)=∂kf/∂xi for i∈I(k), where
[TABLE]
and xi=x1i1⋯xdid for x=(x1…,xd)∈Rd.
Let I(l)=⊔k=0lI(k) and
[TABLE]
We introduce the following assumption.
When ΨN=Ψ, (A3) below holds automatically.
Indeed, (2.16), (2.18) and (2.19) are obvious, and
(2.17) follows from Lemma 5.8.
(A3)
For the subsequence Nn in (A2) iii,
there exists an ℓ∈N satisfying ℓ>ℓ0−d
such that the following hold:
i
For each i∈I(ℓ+1), R∈N, and ϵ>0,
[TABLE]
ii
For each i∈I(ℓ+1), R∈N, and μ-a.s. s,
[TABLE]
These three assumptions are sufficient to obtain an explicit formula for the coefficient of the ISDE arising from μ and to construct the stochastic dynamics as the unique, strong solution of the ISDE.
To obtain the ISDE in translation-invariant form, we further assume the following:
(A4)
The following convergence holds in Lloc2(μ[1]), for μ[1]-a.e. (x,s), or in
C1(SR,s,ϵ) for μ-a.s. s and ϵ>0:
[TABLE]
Here SR,s,ϵ={x∈SR;∣x−si∣≥ϵ for all i}.
Remark 2.2*.*
The assumptions imposed on the limit RPF are mild.
We believe that they allow the construction of many examples of RPFs with Coulomb interactions to which our theorems apply.
**Unlabeled and labeling maps. **
Let
u:(⊔k=0∞(Rd)k)⊔(Rd)N→M(Rd)
be the unlabeled map defined by
[TABLE]
where M(Rd) is the set of all measures on
(Rd,B(Rd)).
If k=0, we set u(x)=0 for x∈(Rd)0 by convention,
where 0 denotes the zero measure.
A label l:S\{0}→(⊔k=1∞(Rd)k)⊔(Rd)N is a measurable map such that
[TABLE]
We write l(s)=(li(s))i.
Throughout this paper, we fix a label such that the absolute values of particles are non-decreasing:
[TABLE]
By (2.12), we have l(s)∈(Rd)N for μ-a.s. s.
For s∈S, we often write l(s)=(si)i.
For y∈S, we define yN=0 if y=0 and, for y=0,
[TABLE]
**Conditional RPFs and couplings. **
Let SRc={∣s∣>R} as before.
Let μR,yNN be the regular conditional probability of μN given by
[TABLE]
From (A1) and (A2), the labeled density mN is continuous on ONN.
Hence, μR,yNN is well defined for all y∈S.
Moreover, μR,yNN admits the k-density functions mR,yN,k.
For μ-a.s. y, the sequence {yN}N∈N defines a coupling for the family of RPFs {μR,yNN}N∈N.
2.3. Explicit formulas for logarithmic derivatives: Theorem 2.1
In this subsection, we present explicit formulas for the logarithmic derivatives.
This result is crucial for proving the uniqueness of solutions
and the existence of strong solutions to the ISDE.
Let I(k) as in (2.14). Let ℓ be as in (A3). For i∈I(ℓ+1), let
[TABLE]
Here we use the integral form of the Taylor expansion of −∇Ψ(x−si) around x=0.
By (2.17), the infinite sum in (2.25) converges absolutely in C1(SR) for μ-a.s. s
(see Proposition 5.1).
We define the function RR,sℓ∈C1(SR) by
[TABLE]
Let ηR,si be the constant in (5.25).
Let I(ℓ)=⊔k=0ℓI(k). Let
[TABLE]
Let μ[1] be the one-reduced Campbell measure of μ.
Theorem 2.1** (Explicit formulas for dμ).**
i*
Assume (A1)–(A3). Then μ admits a ∙-logarithmic
derivative dμ∈Lloc2(μ[1]), which is given by, for each
R∈N and for μ[1]-a.e. (x,s)∈SR×S,*
[TABLE]
ii*
Assume (A1)–(A3). Then, for μ-a.s. s, each
Q∈N, i∈I(ℓ), and ϵ>0,*
[TABLE]
and dμ(x,s) admits a μ[1]-version that is locally Lipschitz in x on R=d(s).
iii*
Assume (A1)–(A4).
Then dμ can be expressed as*
[TABLE]
The convergence in (2.29) is understood in the sense specified by (A4).
Remark 2.3*.*
By construction, ηR,si and RR,sℓ are independent of πSR(s).
The existence of the logarithmic derivative dμ yields a natural μ-reversible diffusion
[36].
Using dμ, the ISDE for the labeled dynamics X=(Xi)i∈N associated
with the S-valued diffusion Xt=∑iδXti can be written as in [27]:
[TABLE]
From (2.28), for each i∈N, the ISDE (2.30) can be rewritten as
Note that the drift coefficient 21dμ(x,s) is not defined on all of Rd×S.
A key point is therefore to ensure that the dynamics X a.s. avoid the singular set,
and to control the long-range divergence coming from the non-integrability of ∇Ψ at infinity.
These issues are central in making the above ISDE formulations rigorous.
In this subsection, we state a generalized ISDE and present the corresponding weak solutions.
Let a:Rd→(R)d2, where a=(akl)k,l=1d, be a uniformly elliptic, bounded, symmetric matrix-valued function such that a∈Cb2(Rd).
Assume that there exists a constant c2.4≥1 satisfying
[TABLE]
Let σ=(σk,l)k,l=1d∈Cb1(Rd) be a matrix such that σtσ=a. Define
[TABLE]
We consider the following elliptic-type ISDE, which generalizes (2.30):
[TABLE]
Here, Bi=(B1i,…,Bdi) is a d-dimensional Brownian motion, and
[TABLE]
Loosely speaking, a pair (X,B) of an (Rd)N-valued continuous process
X=(Xi)i∈N and an (Rd)N-valued {Ft}-Brownian motion
B=(Bi)i∈N, defined on a filtered probability space
(Ω,F,{Ft},P), is called a weak solution of (2.33)
if (2.33) holds P-a.s.
A weak solution (X,B) is called a strong solution if, in addition, X is a measurable function of B.
See Subsection 13.5 for precise definitions.
For X=(Xi)i∈N, the associated m-labeled process X[m] is
[TABLE]
We set Xt[0]=∑i∈NδXti and μ[0]=μ.
We also write Xt=∑i∈NδXti.
Let μ[m] be the m-reduced Campbell measure of μ.
If the m-point correlation function ρm of μ exists, then
[TABLE]
where μx=μ(⋅+∑i=1mδxi∣s({xi})≥1,1≤i≤m) is the m-reduced Palm measure conditioned at x=(x1,…,xm)∈(Rd)m.
Even if μ does not admit an m-point correlation function, the m-reduced Campbell measure can be defined by localization (see (13.6) and [11]).
Theorem 2.2** (Existence of weak solutions).**
Assume (A1)–(A3).
Then ISDE (2.33) has a weak solution X=(Xi)i∈N starting at s=l(s) for μ-a.s. s.
Moreover, for each m∈{0}∪N, the associated m-labeled process X[m] is a μ[m]-symmetric and conservative diffusion.
We next turn to the construction of unique strong solutions to the ISDE (2.33).
In Theorem 2.3, we establish the existence of strong solutions and prove their pathwise uniqueness under suitable conditions.
To this end, we introduce the notion of the tail decomposition of the underlying RPF.
Let T(S) be the tail σ-field defined by
[TABLE]
Let μa=μ(⋅∣T(S))(a)
denote the regular conditional probability of μ
given T(S).
Then
[TABLE]
Definition 2.2* (Tail decomposition).*
An RPF μ is said to admit a tail decomposition if there exists a set
S0∈B(S) such that μ(S0)=1 and, for all a∈S0, the following hold:
[TABLE]
This decomposition provides a natural framework for the construction of strong solutions and the proof of pathwise uniqueness.
Remark 2.4*.*
i
Let S0 be as in Definition 2.2 and a∈S0.
If A∈T(S) satisfies a∈A, then μa(A)=1 by (2.37).
ii
As in [34, Lem. 14.2], one can prove that μ admits a tail decomposition.
Let W(Rd)=C([0,∞);Rd).
We equip W(Rd) with the topology of uniform convergence on compact time intervals and endow W(Rd)N with the corresponding product topology.
With this topology, W(Rd)N is a Polish space.
Let
[TABLE]
For w=(wi)i=1∞∈W(Rd)N,
define the unlabeled path w:=upath(w) by
wt=upath(w)t=∑i=1∞δwti.
We set
[TABLE]
For an (Rd)N-valued process X=(Xi)i∈N, let X=upath(X) denote the associated unlabeled process. If X∈WNE(Ss,i), then IR,T(X) is well defined.
We impose the following conditions, from basic state-space regularity
to absolute continuity with respect to a reference measure ν.
(SIN)P(X∈WNE(Ss,i))=1.
(NBJ)P(IR,T(X)<∞)=1 for each R,T∈N.
(AC)νP∘Xt−1≪ν for all 0<t<∞.
Condition (SIN) ensures that the dynamics remain in the admissible configuration space.
Condition (NBJ) guarantees that, within any finite time interval, only finitely many particles enter a given bounded region.
Condition (AC)ν ensures absolute continuity of the time marginals with respect to the reference measure ν.
In addition, we impose two further conditions.
The condition (MF) is a mild measurability assumption,
whereas (IFC) is a substantive condition whose formulation
requires the finite-dimensional SDE scheme introduced later.
They are introduced in Subsection 13.6 and Subsection 13.5, respectively.
Precise definitions of strong solutions, unique strong solutions,
and pathwise uniqueness are given in Definitions 13.11–13.13.
The following Theorem 2.3 establishes the existence and uniqueness of a
strong solution of (2.33) in the sense of Definition 13.13.
Theorem 2.3** (Unique strong solutions).**
Assume (A1)–(A3). Then the following hold.
i*
For μ-a.s. a and for μa-a.s. s,
the ISDE (2.33) admits a unique strong solution X starting at s=l(s)
under the constraints (SIN), (NBJ), (AC)$${}_{\mu_{\mathsf{a}}}, (MF), and (IFC).*
ii*
For each m∈{0}∪N, the associated m-labeled process X[m] is a μa[m]-symmetric and conservative diffusion.*
By Theorem 13.4, we obtain the following corollary of Theorem 2.3.
Let ξa be an RPF such that ξa≪μa.
Corollary 2.1**.**
Assume (A1)–(A3).
Then, for μ-a.s. a, any weak solution of (2.33) with initial distribution
ξa∘l−1 is unique in law and pathwise unique
under (SIN), (NBJ), (AC)$${}_{\mu_{\mathsf{a}}}, and (IFC).
From Theorems 2.1–2.3, we obtain the following corollary.
Corollary 2.2**.**
Assume (A1)–(A3).
Then the conclusions of Theorem 2.2 and Theorem 2.3 apply to the ISDE (2.41):
The analysis of fully labeled dynamics on the full label space (Rd)N
plays a crucial role in understanding infinite-particle systems.
While the m-labeled dynamics X[m] provide finite-dimensional projections
that admit invariant measures, they capture only partial aspects of
the underlying infinite system.
While the m-labeled dynamics X[m] provide projections to the tame infinite-dimensional spaces
that admit invariant measures, they capture only partial aspects of
the underlying infinite system.
To obtain a complete description, it is necessary to go beyond finite labels
and establish the existence of a genuine diffusion on (Rd)N.
This step is indispensable, as it ensures that the dynamics are well defined
for the entire infinite system rather than merely for restricted subsystems.
Historically, the construction of fully labeled dynamics has been regarded as
a highly challenging problem.
Even for determinantal processes such as the Dyson model, establishing such
dynamics required delicate arguments based on random matrix theory.
For general Coulomb systems, the situation is more severe:
the long-range nature of the interaction and the absence of a classical
Gibbsian structure prevent a direct application of Dirichlet form methods
or invariant-measure techniques.
Against this backdrop, proving the existence of fully labeled conservative
diffusions represents a substantial advance.
In Theorem 2.4, we prove that the fully labeled process X
constructed in Theorem 2.3 is an (Rd)N-valued conservative diffusion.
Unlike the m-labeled processes X[m], the diffusion X
is obtained without relying on the existence of an invariant measure,
reflecting a fundamental difference between finite-label and full-label dynamics.
Nevertheless, the result guarantees that the fully labeled dynamics form a
legitimate stochastic process with infinite lifetime, thereby providing a
rigorous foundation for the dynamical theory of Coulomb systems.
Let S⋆⋆⊂(Rd)N be the subset defined in (11.14).
Let X be the solution of the ISDE (2.33) constructed in Theorem 2.3.
Let u denote the unlabeled map as in (2.21).
Theorem 2.4** (Diffusion processes in (Rd)N).**
Assume (A1)–(A3).
Then X is a conservative diffusion with state space
S⋆⋆ and
μ(u(S⋆⋆))=1.
Here, saying that X is a conservative diffusion with state space
S⋆⋆ means that the regular conditional law
[TABLE]
admits a version such that the family
{Px}x∈S⋆⋆
consists of probability measures on
C([0,∞);S⋆⋆)
satisfying the strong Markov property and having infinite lifetime
[5].
In particular, the diffusion law Px is well defined for every
x∈S⋆⋆.
The state space can, in principle, be extended to the whole (Rd)N
by freezing the dynamics outside S⋆⋆.
We study finite-volume approximations of the infinite-particle dynamics.
These approximations are realized in two complementary ways:
by restricting the system to a bounded domain with reflecting boundary conditions (Proposition 12.1),
and by truncating the system to N particles (Theorem 2.5).
These theorems rigorously establish the convergence of these finite-volume approximations
to the infinite-dimensional dynamics.
Let dN be the logarithmic derivative of μN.
By (1.3) and (A1),
[TABLE]
Set bN=21(diva+dNa).
Then the SDE describing the N-particle dynamics XN=(XN,i)i=1N is given by,
for 1≤i≤N,
[TABLE]
Here 21nNa denotes the inward unit normal vector on ∂ON,
and LNN,i is a nonnegative, nondecreasing continuous process,
called the local time on the boundary ∂ON
(cf. [8, pp. 251–256]).
If ON=Rd, no local time term appears in (2.42). The SDE then becomes
[TABLE]
We impose
[TABLE]
Here l=(li)i denotes the labeling map introduced in
(2.22), and lN=(li)i=1N.
The initial conditions of these solutions are coupled through the labeling map.
Let c2.5 be a positive constant. Let
[TABLE]
Define the RPF φdμ by (φdμ)(A)=∫Aφdμ for
A∈B(S).
Let XN and X be the (unique strong) solutions of SDEs (2.42) and (2.33), respectively.
We assume:
[TABLE]
By setting XtN,i=li(s) for all t≥0
and i>N, we regard XN=(XN,i) as a W(Rd)N-valued random variable.
We begin by recalling earlier developments on logarithmic derivatives,
which constitute a foundational concept in the theory of infinite-particle dynamics.
Logarithmic derivatives.
The concept of the logarithmic derivative for RPFs
was introduced in [27].
In [27], it was proved that if ν admits a logarithmic derivative
dν and satisfies the quasi-Gibbs property, then (2.3)
admits a weak solution, and the associated unlabeled dynamics
Xt=∑iδXti
is a ν-reversible diffusion.
In [36], it was proved that, for general random point fields,
the mere existence of a logarithmic derivative suffices to construct
such dynamics.
However, the mere existence of a logarithmic derivative
and the construction of a weak solution to the ISDE
do not guarantee uniqueness of solutions
or the existence of strong solutions.
Moreover, they do not allow a detailed analysis
of the qualitative behavior of the solutions.
Addressing these issues is therefore essential.
Hyperuniformity.
The logarithmic derivatives for the sineβ, Airyβ(β=1,2,4), Bessel2, and Ginibre RPFs have been computed [27, 28, 35, 9].
The former three are one-dimensional systems
with two-dimensional Coulomb interaction potentials, and hence are not Coulomb RPFs.
These computations require strong hyperuniformity of particle numbers within SR for each R∈N, a property that must be established separately for each RPF.
For example, in the case of the Ginibre RPF,
we use the following estimate
[38, 32]:
[TABLE]
Such an estimate is called hyperuniform because,
if μN is replaced by ΛN,
then the right-hand side becomes
O(R2) [23].
Here ΛN is the conditioned Poisson RPF defined by
[TABLE]
Here Λ is the Poisson RPF on Rd
with Lebesgue intensity dx,
and R(N) is the positive number such that
the volume of SR(N)
coincides with N.
While Leblé [23] established a hyperuniform estimate
for general β>0, it is weaker than (2.46) and is insufficient for our purposes.
The proofs of the strong hyperuniform estimates for the above-mentioned examples
rely on techniques specific to determinantal RPFs,
and, except in the Ginibre case, do not extend to general Coulomb RPFs,
which necessitates a new approach.
Logarithmic derivatives of Coulomb RPFs.
Until now, the logarithmic derivatives of Coulomb RPFs have been computed only in the case of the two-dimensional Coulomb RPF at inverse temperature β=2, namely the Ginibre RPF [27].
In contrast, we present here a new, remarkably simple and transparent method
that yields uniform estimates for the logarithmic derivative of μN
(2.47) and for its conditional analogue (2.48),
valid simultaneously for all spatial dimensions d≥2 and all inverse temperatures β>0.
The key insight is to exploit a certain harmonicity property of the interaction potential, which permits effective control over the infinite-dimensional interactions. Our approach is robust, conceptually clear, and provides a versatile tool for analyzing the stochastic dynamics of d-dimensional Coulomb RPFs, thereby overcoming the limitations of all previously known methods.
**Annealed and quenched estimates for the logarithmic derivative. **
We introduce a robust method for estimating the logarithmic derivative of μN
and refine it to conditional probabilities μR,yN,
exploiting a form of harmonicity of the interaction potential.
As a first step, we obtain the following uniform L2-bound:
[TABLE]
This bound provides the control needed to define logarithmic derivatives
rigorously in the Coulomb setting (Lemma 4.3).
The proof of (2.47) is remarkably simple:
it relies only on (2.9), (2.10), and the basic definition of the logarithmic derivative.
Moreover, the estimate (2.47) may be viewed as a form of hyperuniformity,
since the quantity in (2.47) diverges when μN,[1] is replaced by
the one-Campbell measure ΛN,[1] of the conditioned Poisson point process
ΛN, even after excluding a neighborhood of x.
This hyperuniformity is much weaker than that used in the previous works
[16, 17, 18, 28, 27, 34, 35, 32],
which motivates the development in this paper of a new framework for the
analysis of stochastic dynamics associated with Coulomb random point fields.
We next establish the quenched analogue of (2.47):
Inequality (2.48) plays a decisive role in the proof of Theorem 2.1.
Indeed, it fulfills the role previously played by strong hyperuniform estimates,
but with a crucial difference:
its proof is significantly simpler.
Earlier results relied on bounding averaged fluctuations
of the exterior configuration,
whereas here we directly control the external condition
in a quenched form.
This shift in perspective—replacing averaged estimates
by a quenched analysis—is one of the key conceptual innovations
of the present paper.
Interacting Brownian motions.
Lang initiated the systematic study of interacting Brownian motions, which was subsequently developed by Fritz, Tanemura, and others in the case of finite-range interaction potentials [20, 21, 7, 42].
When the interaction becomes long-ranged, however, the problem of constructing solutions to ISDEs poses significant difficulties, even for potentials within Ruelle’s class—namely, those that are superstable and integrable at infinity.
Later, the first author introduced the Dirichlet form approach, which made it possible to establish the existence of weak solutions to ISDEs for a broad family of interaction potentials [26, 27, 28, 29]. This framework encompasses, in particular, the Ginibre interacting Brownian motion.
Nonetheless, the stronger questions of constructing pathwise unique strong solutions and establishing their uniqueness remained open.
In [46], Tsai succeeded in constructing pathwise unique strong solutions of ISDEs for the sinβ RPFs with all β≥1. These processes correspond to infinite-particle systems in one dimension with logarithmic interactions, that is, with the two-dimensional Coulomb potential. In this sense, the ISDE in [46] can be regarded as involving a Riesz potential. Tsai’s proof relies on a delicate coupling mechanism that exploits special properties of the logarithmic interaction in one dimension. Such an approach, however, has no counterpart in higher dimensions with genuine Coulomb potentials.
Subsequently, in [34, 18], a general method for constructing pathwise unique strong solutions of ISDEs was developed. The central idea is to represent the ISDE as an infinite hierarchy of finite-dimensional SDEs evolving in a random environment. The well-posedness of this hierarchy is then reduced to a problem of tail triviality for the equilibrium measures of the associated unlabeled dynamics.
To apply this framework, one requires both an equilibrium state (RPF)
and suitable regularity properties of the corresponding logarithmic derivative.
The construction of Coulomb RPFs was established in [1, 44].
By exploiting the quenched bound (2.48),
we apply a Taylor expansion method adapted to the Coulomb interaction potential
and prove the following explicit representation of the logarithmic derivative
(Theorem 2.1):
[TABLE]
Using this representation, we further show that the associated ISDE satisfies the key assumption
known as the IFC condition, which is required in the above theory (Theorem 2.3).
Further related works.
Katori and Tanemura proposed an alternative approach
to constructing stochastic dynamics with logarithmic interaction.
Their method, based on space–time correlation functions,
is restricted to the special case of inverse temperature β=2 in one-dimensional systems.
For further details, see [12] and the references therein.
Moreover, [33] investigated in detail
the relationship between the logarithmic derivative
and the space–time correlation approach.
It was rigorously established that the stochastic dynamics
constructed by these two seemingly distinct methods
are, in fact, identical.
This coincidence highlights the fundamental role
of the logarithmic derivative
in unifying different perspectives
on infinite-particle dynamics with logarithmic interactions.
In [2], Assiotis and Mirsajjadi constructed a non-equilibrium dynamics whose equilibrium distribution is the inverse Bessel RPF, and they identified the ISDE governing this dynamics.
Their method is algebraic in nature, relying on intertwining theory and the use of derivatives of the so-called characteristic polynomial.
This construction bears a striking resemblance to the framework of logarithmic derivatives developed in [27].
In [4], Bufetov and Kawamoto studied intertwining relations for Laguerre processes with parameter α across different particle numbers.
They computed the logarithmic derivatives of the corresponding finite-particle systems and pointed out that the construction of the infinite-particle limit process remains an open problem, to be addressed in subsequent work.
In [14], Kawamoto constructed a unique strong solution to the ISDE associated with the generalized sine random point field at inverse temperature β=2, combining random matrix theory with the general ISDE theory [34].
In [6], Esaki and Tanemura constructed a unique strong solution to a jump-type ISDE associated with a jump-type infinite-particle system with long-range interaction, and constructed the interacting Cauchy processes.
In [40], Suzuki established the ergodicity of the unlabeled dynamics of the Ginibre interacting Brownian motion, thereby providing a concrete confirmation of the general ergodicity theory.
Key difficulties and techniques.
A major obstacle underlying this program is the pronounced and extremely strong long-range interactions
characteristic of Coulomb systems.
Previous methods relied on delicate hyperuniform estimates,
which are available only in the unique determinantal case d=2 and β=2.
Our approach replaces this reliance
with new analytic estimates for the logarithmic derivative,
together with a quenched formulation that yields robust L2 bounds.
This combination provides a conceptually simple yet powerful framework,
which enables us to establish strong well-posedness
of Coulomb interacting Brownian motions
in all spatial dimensions d≥2 and for all inverse temperatures β>0.
Structure of the proofs.
The proofs of the main theorems are distributed as follows:
Theorem 2.1 is proved in Section 5,
Theorem 2.2 in Section 9,
Theorem 2.3 in Section 10,
Theorem 2.4 in Section 11, and
Theorem 2.5 in Section 12.
Sketch of the proof strategy.
The main theorems are proved in Sections 4–12,
and the overall strategy is as follows.
In Section 4, we establish the existence of the logarithmic derivative.
A key ingredient is a robust uniform L2-estimate for the logarithmic derivative,
which exploits the harmonicity of the Coulomb potential.
Section 5 strengthens this result by proving a quenched estimate
and deriving an explicit representation of the logarithmic derivative.
In Section 6, we establish non-collision properties of ISDE solutions
by direct SDE techniques.
These arguments do not rely on Dirichlet form theory;
nevertheless, the resulting non-collision property is later used
to prove uniqueness of Dirichlet forms.
Sections 7–9 form the technical core of the paper:
we introduce Dirichlet form schemes for m-labeled dynamics
and construct infinite-volume processes under both lower and upper Dirichlet forms.
In Section 10, using the explicit representation of the logarithmic derivatives,
we verify the IFC condition ensuring consistency of the dynamics
and upgrade the constructions to pathwise unique strong solutions of the ISDEs.
In Section 11, using the uniqueness of solutions of the ISDE,
we establish uniqueness of extensions of the Dirichlet forms
and prove the existence of the fully labeled diffusion process in (Rd)N.
In Section 12, we prove the convergence of finite-particle systems with reflecting boundary conditions
using the convergence of the approximating lower Dirichlet forms.
Finally, combining the uniqueness of extensions of Dirichlet forms with a sandwich argument,
we prove that the N-particle dynamics converge
to the infinite-particle ISDE dynamics.
3. Concrete examples
We present concrete examples of RPFs that satisfy assumptions (A1)–(A4).
Throughout this section, we assume d≥2 and β>0.
Here aN is a convergent sequence in Rd with limit a,
and a belongs to the droplet,
that is, the support of the equilibrium measure associated with the potentials Φ and Ψ.
Take ΨN=Ψ and ON=Rd.
Then (A1) ii and (A3) are automatically satisfied.
Moreover, (A1) i is satisfied in a wide class of examples.
(A2) i is obvious,
and ii–iii follow from [44].
Indeed, in [44],
a convergent subsequence was constructed,
yielding the limiting Coulomb RPF μ,
which possesses bounded correlation functions of all orders,
together with other regularity properties.
Moreover, the assumptions in [44]
are weaker than those stated above;
see [44, 1.25, Th. 6].
If d=2 and Φ(x)=21∣x∣2, then the resulting ISDE is
[TABLE]
In this way, the term −2βat appears, canceling the effect of the macroscopic position a in the droplet.
We now present two further specific examples of RPFs that satisfy assumptions (A1)–(A4).
These examples exhibit a form of “invariance” under certain transformations, which follows from the invariance of the corresponding finite-particle systems.
As a consequence, the one-point correlation function is constant for each N-particle system.
Example 3.2* (Periodic approximation).*
Let T=(−21,21]d denote the unit torus.
Set ΦN=Φ=0, and let ΨT(x,y) denote the Green’s function of
−c\ref;perΔ on T subject to periodic boundary conditions,
where c3.1 is the positive constant determined according to
(2.7).
Set ON=(−2N1/d,2N1/d)d and
ΨN(x,y)=Nd2−dΨT(N1/dx,N1/dy).
By construction, the measure μ is invariant under translations on Rd.
Example 3.3* (Sphere approximation).*
Let
S={x^∈Rd+1:∣x^∣=1}
denote the unit sphere, and let ed+1=(0,…,0,1) denote the north pole.
Define the stereographic projection
ϖ:S∖{ed+1}→Rd
by
[TABLE]
so that ϖ(x^) is the intersection of the line through ed+1
and x^ with the hyperplane Rd×{0}.
The inverse map ϖ−1:Rd→S∖{ed+1} is
[TABLE]
with Jacobian J(x)=2d/(1+∣x∣2)d. Define
ΨS(x^,y^) as the Green’s function of −c\ref;spΔS on S,
where ΔS is the Laplace–Beltrami operator and c3.2 is determined according to (2.7) and (2.9).
Define ON=Rd and
[TABLE]
Note that, by construction, the measure μ is invariant under the full Euclidean group E(d), i.e., under all translations, rotations, and reflections.
The resulting diffusion X is invariant under the full Euclidean group E(d).
Remark 3.1*.*
The construction of the Coulomb RPFs can be extended to homogeneous spaces and,
more generally, to compact or weighted Riemannian manifolds,
including models beyond the Gaussian, periodic, and spherical cases.
These generalizations will be developed in future work.
4. Construction of the logarithmic derivative dμ of μ
In this section we establish the existence of the logarithmic derivative dμ of μ.
The section is organized into four parts.
In Subsection 4.1, we prove a uniform L2-bound for the logarithmic derivatives dN of the N-particle systems.
In Subsection 4.2, we derive a uniform L2-bound for the logarithmic derivatives of μRN.
In Subsection 4.3, we prove the strong convergence of the density functions of μRN as N→∞ for each fixed R∈N.
Finally, in Subsection 4.4, we construct the logarithmic derivative dμ of μ.
All results in Section 4 are proved under Assumptions (A1)–(A3).
4.1. L2-uniform bound of the logarithmic derivative dN of μN
Let μN be as in (A1)–(A3).
Let μN,[1] denote the one-reduced Campbell measure associated with μN, i.e., the specialization of the definition (2.1) (given for a general RPF ν) to the case ν=μN.
Let mN be the labeled density function of μN given by (1.3).
Let
Moreover, by (A1) and (A2), the density mN vanishes on collision configurations,
that is, mN(x)=0 whenever xi=xj for some i=j.
Hence the singular delta terms in (4.3) do not contribute, and we obtain
[TABLE]
Combining this with (4.2), the assertion follows.
∎
For each Q∈N, let fQ:Rd→R be a cut-off function such that
[TABLE]
In equation (4.4), we adopt the convention such that
[TABLE]
Let c4.1, c4.2, and c4.3 be the quantities defined by
[TABLE]
Lemma 4.2**.**
For each Q, the quantity c\ref;22z is finite.
Proof.
By (4.4), c\ref;22b(N) are nonnegative and bounded.
By (2.6) and (2.10), c\ref;22c(N) are non-negative.
By (2.4) and (4.4),
4.2. L2-uniform bound of logarithmic derivative of μRN
We continue to work under Assumptions (A1)–(A3).
Let μRN=μN∘πSR−1.
We regard μRN as a probability measure on SR, which is the configuration space over SR.
Let dRN be the logarithmic derivative of μRN, that is, for each FR[1]-measurable h∈C0∞(SR)⊗D∘b,
FR[1]=B(SR)×FR and FR=σ[πSR]
[TABLE]
Here μRN,[1] is the one-reduced Campbell measure of μRN.
We regard dRN as an FR[1]-measurable function on Rd×S
satisfying
dRN(x,s)=0 for x∈/SR and
dRN(x,s):=dRN(x,πSR(s)).
We regard μRN,[1] as a measure on
(Rd×S,B(Rd)×FR)
with μRN,[1](SRc×S)=0.
Then (4.12) can be written as follows: for each
FR[1]-measurable h∈C0∞(SR)⊗D∘b,
[TABLE]
We often use the latter interpretation of dRN,
as it allows us to regard dRN(x,s) as a function defined on
Rd×S independently of R.
Remark 4.1*.*
If we replace C0∞(SR)⊗D∘b by C∞(SR)⊗D∘b as the space of test functions in (4.12), then a local-time-type singular drift appears in the representation of dRN, indicating the reflection of particles on ∂SR.
Let mN be as in (1.3).
Then μRN has the labeled density mRN defined on
⊔k=0NSRk such that, for 0≤k≤N and k+l=N,
[TABLE]
Here ZRN is the normalization.
Because mN(x) is a symmetric function of x=(xi)i=1N, mRN is well defined.
For (x,s)∈SR×S and s(SR)≤N−1, dRN satifies
[TABLE]
and by dRN(x,s)=0 for s(SR)≥N, where s=u(s).
Definition 4.1*.*
Let 0≤k≤N−1 and k+l=N−1.
Let s=u(s) and t=u(t).
We set, for s∈SR satisfying s(SR)=k,
[TABLE]
Here f:SR×SN−1→R, SN−1={u∈S;u(Rd)=N−1}, and
fˇ:SR×(Rd)N−1→R is the function satisfying that
fˇ(x,s,t)=f(x,s,t) and that fˇ(x,s,t) is symmetric in (s,t).
The following identity plays an important role in our analysis.
Lemma 4.4**.**
Let dRN and dN be as in (4.14) and (4.1), respectively. Then
[TABLE]
Proof.
Let k+l=N−1. Let s, s, and t be as in (4.15).
Then
4.3. Strong convergence of density functions of μRN
We continue to work under Assumptions (A1)–(A3).
Let ΛR be the Poisson RPF with intensity measure 1SRdx.
Let ΛR[1] be the one-reduced Campbell measure of ΛR.
We denote by mRN,[1] the Radon–Nikodym density of μRN,[1] with respect to ΛR[1].
The existence of mRN,[1] is obvious from (1.3).
Let
[TABLE]
We regard SR,k[1] as a subset of SR×SR, where
SR is the configuration space over SR.
Let μR[1] be the one-reduced Campbell measure of μR=μ∘πSR−1.
Lemma 4.6**.**
Let R∈N such that R≥Q+1.
Let k,l∈{0}∪N.
i*
{fQmR,k,lN,[1]}n∈N is relatively compact in L2(ΛR[1]).*
ii*
Let Nn be as in (A2).
There exist functions mR,k,l[1] satisfying*
[TABLE]
iii*
μR[1] has the Radon–Nikodym density
mR[1]=(dμR[1]/dΛR[1]), which is defined by
mR[1]∧l=mR,k,l[1] on SR,k[1].*
Proof.
Note that mR,k,lN,[1](x,s) is symmetric in (x,s) when regarded as a function on SRk+1 by s=u(s). Let
[TABLE]
Using the permutation invariance of mR,k,lN,[1] in (x,s) and noting that k+1 particles exist in SR for (x,s)∈SR,k[1], we obtain
[TABLE]
Hence, combining this with (4.4) and Lemma 4.3 yields the following:
[TABLE]
As before, we regard fQmR,k,lN,[1](x,s) as a function on SRk+1.
Note that fQmR,k,lN,[1]=0 on (∂SR)×SRk from (4.4) and R≥Q+1.
Hence, from the Rellich embedding theorem and (4.21), there exists q>1 such that
{fQmR,k,lN,[1]}N∈N is relatively compact in Lq(SRk+1,dx),
where we take 1<q<n/(n−1) for n=d(k+1).
According to this and (4.4), {fQmR,k,lN,[1]}n∈N is relatively compact in Lq(ΛR[1]).
From (4.18), {fQmR,k,lN,[1]}N∈N is bounded in L∞(ΛR[1]).
Hence, {fQmR,k,lN,[1]}N∈N is relatively compact in L2(ΛR[1]).
This proves i.
From i and a diagonal argument, there exists a subsequence of {Nn},
denoted by the same symbol, and a limit {mQ,R,k,l[1]} such that
[TABLE]
for all R≥Q+1 and k,l∈{0}∪N.
From the weak convergence
limn→∞μNn=μ
given in (A2),
we obtain the uniqueness of limit points of
{fQmR,k,lNn,[1]}n∈N.
Hence the full sequence converges to
mQ,R,k,l[1]
in (4.22).
Using the consistency (4.20), we define the function mR,k[1] by
[TABLE]
By the decomposition SR×SR=⊔k=0∞SR,k[1],
we define the function mR[1] by mR[1](x,s)=mR,k[1](x,s), (x,s)∈SR,k[1].
We define μR[1]=mR[1]dΛR[1].
Then μR[1] is the one-reduced Campbell measure of μR=μ∘πSR−1.
This proves iii.
∎
Proposition 4.1**.**
i*
μ has a k-density function mR,k on SRk for each R,k∈N.*
The restriction of mR[1](x,s) to SRk is symmetric for each k∈N.
Hence, we construct the k-density function
mR,k of μ on SRk from mR[1] by the formula
[TABLE]
The constant c4.4
is determined by the normalization.
This proves i.
The one-point correlation function ρ1 is given by the formula
[TABLE]
By (2.13), ∫Ss(SR)μ(ds)<∞.
Hence, the right-hand side is finite.
Hence, ρ1 is the one-point correlation function of μ.
This proves ii.
∎
4.4. The logarithmic derivative dμ of μ
We continue to work under Assumptions (A1)–(A3).
This subsection constructs the logarithmic derivative dμ of μ.
The construction proceeds in two steps.
First, we define the logarithmic derivative dR of μR as the limit of the logarithmic derivatives of μRNn as Nn→∞.
Next, we define dμ as the limit of dR as R→∞.
Let SR,k[1] and mRN,[1] be as in (4.18).
Let MR,k,lN,[1]={(x,s)∈SR,k[1];0<mRN,[1](x,s)<l}.
We define
[TABLE]
We define the values of dR,k,lN(x,s) on the boundary ∂SR×SR
as limits of the values from the interior SR×SR.
Let mR[1] be as in Lemma 4.6.
Replacing mRN,[1] by mR[1], we define MR,k,l[1] and dR,k,l similarly.
Let {Nn} be as in (A2).
Lemma 4.7**.**
For each k,l∈{0}∪N and R∈N,
[TABLE]
Proof.
From Lemma 4.5, (4.18), and (4.24), we easily deduce that
[TABLE]
is bounded in L2(ΛR[1]).
Hence, we have a subsequence of {Nn}, denoted by the same symbol, and a limit dR,k,l such that
Hence, from (4.27) and (4.28), we obtain for any h∈C0∞(SR)⊗D∘b
[TABLE]
Hence, dR,k,l=dR,k,l(fQmR,k,l[1])1/2.
Together with (4.27), this yields (4.25).
From (4.24), dR,k,lN satisfies dR,k,lN(x,s)=dR,k,l+1N(x,s)
for (x,s)∈MR,k,lN,[1] and for all l∈N.
This property is inherited by dR,k,l
from (4.27)–(4.29).
Hence, we obtain (4.26).
∎
Let MR[1]={(x,s)∈SR×SR;0<mR[1](x,s)<∞}.
We define
[TABLE]
We define the values of dR(x,s) on the boundary ∂SR×SR
as limits of the values from the interior SR×SR.
Lemma 4.8**.**
For all R≥Q+1,
[TABLE]
Proof.
Recall that dμRN,[1]=mRN,[1]dΛR[1]. By Lemma 4.5, we obtain (4.31).
From (4.31), {dRN(fQmRN,[1])1/2}N∈N is relatively compact under the weak convergence in L2(ΛR[1]).
We easily see
[TABLE]
From Lemma 4.7 and (4.34), the limit points are unique and coincide with dR(fQmR[1])1/2. This proves (4.32).
By (A2) iii, we have μ=limn→∞μNn weakly. Hence by (2.11),
A function satisfying (4.37) for all h in D∘b[1] is called a logarithmic derivative of μR.
We readily see that dR is such a function.
We now construct the logarithmic derivative dμ of μ using the consistency of dR.
In the following theorem, we regard dR as a function on Rd×S, denoted by the same symbol, and define
dR(x,s)=0 for x∈/SR, and dR(x,s):=dR(x,πSR(s)).
Proposition 4.2**.**
There exists a limit of dR for μ[1]-a.e. and in Lloc2(μ[1]), which we denote by dμ.
The limit dμ is the ∙-logarithmic derivative of μ.
Proof.
By (4.37) and μR=μ∘πSR−1,
we deduce that, for each FR[1]-measurable h=f⊗g in D∘b[1] with f in C0∞(SR),
[TABLE]
For a non-negative f satisfying ∫Rdfdx=1,
consider the probability measure f⊗1μ[1].
By (4.40), dRf is an {FR[1]}-martingale under f⊗1μ[1].
Hence the first assertion follows from the martingale convergence theorems and (4.35).
By the first assertion and (4.40), for each h in D∘b,
[TABLE]
By Lemma 13.6, (4.41) holds for all h in D∙b.
Hence dμ is the ∙-logarithmic derivative of μ.
∎
In this section, we present explicit formulas for the logarithmic derivative
dμ and complete the proof of Theorem 2.1.
Although Proposition 4.2 established the existence of a logarithmic derivative
dμ, from which the existence of a weak solution to the ISDE can be deduced,
it does not address the uniqueness of the solution
nor the existence of a strong solution.
To resolve these issues, we need an explicit representation of dμ.
Unless stated otherwise, all results in Section 5 are proved under Assumptions (A1)–(A3).
5.1. Construction of the logarithmic derivative dR,y of μR,y
The goal of this subsection is to construct the logarithmic derivative of
μR,y=μ(⋅∣πSRc(s)=πSRc(y)).
Let μN be the distribution of the N-particle system described in Section 1.
For μN-a.s. y∈S, let
[TABLE]
Let ρR,yN,1 be the one-point correlation function of μR,yN (on SR).
Let μR,y,xN be the reduced Palm measure of μR,yN conditioned at x.
Let
[TABLE]
be the one-reduced Campbell measure of μR,yN.
From (A1) and (A2), ρR,yN,1(x) is continuous in x on SR∩ON.
Let a={aq}q∈N be a family of increasing sequences
aq={aq(R)}R∈N of natural numbers satisfying (13.3).
Let K[aq] be as in (13.4).
By (A2), {μN}N∈N is tight.
Hence, for each q∈N, there exists an aq={aq(R)}R∈N satisfying
[TABLE]
Then {K[aq]}q∈N becomes an increasing sequence of compact sets in S.
Let μR,yN,qN,[1] be the one-reduced Campbell measure of μR,yN,qN.
Let fQ∈C0∞(Rd) be as in (4.4), and let Q<R.
Fix y∈∪q∈NK[aq], and choose q∈N such that
y∈K[aq]\K[aq−1], where K[a0]=∅ by convention.
Note that for each y∈∪q∈NK[aq], q is uniquely determined.
We introduce the quantities
c5.1, c5.2, and c5.3 as follows:
[TABLE]
By (2.10), we have ψN≥0.
Hence c\ref;42a+4βc\ref;42b≥0.
Lemma 5.1**.**
For each Q∈N and all y∈∪q∈NK[aq], c\ref;42c(Q,q) is finite.
Proof.
From (13.4) and (5.3), we have, for all Q<R, N∈N, and y∈S,
[TABLE]
Combining this with (5.4), (2.8), and (2.10), we obtain, for
y∈K[aq]\K[aq+1],
[TABLE]
This, together with (4.4) and (5.4), proves Lemma 5.1.
∎
We define, for μN-a.s. y and x∈SR,
[TABLE]
Here πA(s)=∑si∈Aδsi and SR={∣s∣≤R}.
Let yN be as in (2.23).
By (4.2) and (5.5), we obtain, for x∈ON,
[TABLE]
The following result is a local, quench version of Lemmas 4.3 and 4.5.
Lemma 5.2**.**
Let R≥Q+1. Then for y∈∪q∈NK[aq],
[TABLE]
Proof.
(5.4) corresponds to (4.5) used in the proof of Lemma 4.3.
The remainder of the proof of Lemma 5.2 is similar to that of Lemma 4.3, and is therefore omitted.
∎
Let mR,yN,[1] be the Radon–Nikodym density of μR,yN,[1] with respect to ΛR[1].
Such densities mR,yN,[1] exist by (1.3).
Define
[TABLE]
where SR,k[1]={(x,s)∈SR×SR;s(SR)=k}.
By definition, μR,y=μ(⋅∣πSRc(s)=πSRc(y)) is an RPF on Rd concentrated at πSRc(y) on SRc.
We regard μR,y as an RPF on SR as well as that on Rd.
Let μR,y[1] be the one-reduced Campbell measure of μR,y.
Lemma 5.3**.**
Let R≥Q+1, k,l∈{0}∪N.
For μ-a.s. y, the following hold:
i*
{fQmR,y,k,lN,[1]}N∈N is relatively compact in L2(ΛR[1]).*
ii* There exist functions mR,y,k,l[1] satisfying*
[TABLE]
iii*
μR,y[1] has the Radon-Nikodym density mR,y[1]=(dμR,y[1]/dΛR[1]),
which is defined by
mR,y[1]∧l=mR,y,k,l[1] on SR,k[1].*
Proof.
Lemma 5.3 follows from Lemma 5.2 using similar arguments as for the proof of Lemma 4.6.
Hence, the proof is omitted.
∎
We proceed in the same manner for the conditional measures.
Replacing mRN,[1] by mR,yN,[1] (resp. mR,y[1]) in (4.24),
we define MR,y,k,lN,[1] and dR,y,k,lN (resp. MR,y,k,l[1] and dR,y,k,l) similarly.
Lemma 5.4**.**
For R≥Q+1, k,l∈{0}∪N, and μ-a.s. y,
[TABLE]
Proof.
The proof of Lemma 5.4 is identical to that of Lemma 4.7, using Lemmas 5.2 and 5.3.
We omit the details.
∎
Replacing mR[1] by mR,y[1] in (4.30), we define dR,y analogously to dR.
Lemma 5.5**.**
For all R≥Q+1 and μ-a.s. y∈K[aq]\K[aq−1],
[TABLE]
Proof.
Lemma 5.5 corresponds to Lemma 4.8 and follows from Lemmas 5.2–5.4 using the same argument as for the derivation of Lemma 4.8 from Lemmas 4.5–4.7. Hence, the proof is omitted here.
∎
Lemma 5.6**.**
Let R≥Q+1 and q∈N.
For μ-a.s. y∈K[aq]\K[aq−1], the following hold:
[TABLE]
for all FR[1]-measurable h=(fQf^)⊗g∈D∘b[1].
Proof.
Note that Lemmas 5.2–5.6 correspond to Lemmas 4.5–4.9, respectively.
Hence, Lemma 5.6 follows from Lemmas 5.2–5.5 using similar arguments as for the proof of Lemma 4.9. Hence, the proof is omitted here.
∎
5.2. Explicit formulae of the logarithmic derivative dR,y of μR,y
We continue to work under Assumptions (A1)–(A3).
Let I(k) be as in (2.14).
Let f(i)=∂kf/∂xi,
where ∂xi=(∂x1)i1⋯(∂xd)id.
Let I(ℓ)=⊔k=0ℓI(k) as in Theorem 2.1.
From (2.8), −∇xΨN(x,y)=0 for (x,y)∈/ON2.
Using Taylor expansions of
−∇xΨNn(x,y)
and
−∇Ψ(x−y)
in x
at the origin x=0
for each ∣y∣>R and (x,y)∈ON2,
we obtain, for ∣x∣≤R,
[TABLE]
Here Qℓ,n and Qℓ are defined by
[TABLE]
Let RR,yℓ,n,n, RR,yℓ,n, and RR+ϵ,yℓ be such that, for μ-a.s. y=∑iδyi,
[TABLE]
Proposition 5.1**.**
For each R∈N and μ-a.s. y, we obtain
[TABLE]
Proof.
From (5.15), (5.14), and (2.17), we obtain, for each ϵ>0,
[TABLE]
The last term is finite by (2.17). This yields (5.16).
From the above inequality, we obtain, for μ-a.s. y=∑iδyi,
[TABLE]
Since y(SR+ϵ)<∞, this yields, for μ-a.s. y=∑iδyi,
The convergence above holds in C1(SR,ϵ(s)) for each ϵ>0 and all s such that s=∑si∈SRδsi, where SR,ϵ(s)={x∈SR;∣x−si∣≥ϵ for all i}.
By Lemma 5.5 and the above equation, we obtain, weakly in L1(ΛR[1]),
[TABLE]
From μR,y[1](SR×SR)>0, we have ∫SR×SRmR,y[1]dΛR[1]>0. Hence, the family
{xifQ(x,s)mR,y[1](x,s);i∈I(ℓ)}
is linearly independent in L1(ΛR[1]).
From this and the display equation right above, each numerical sequence {ηR,yi,n,n}n∈N has a finite limit ηR,yi satisfying
Let xji for x∈SRℓ be as in (2.15).
For each j∈I(ℓ),
[TABLE]
Let xj=(xji)i∈I(ℓ) for j∈I(ℓ).
Then the vectors {xj}j∈I(ℓ) are linearly independent in RI(ℓ)
for each x∈SRℓ by (2.15).
Hence, using the display equation right above, we find that, for each i∈I(ℓ), the numerical sequence ηR,yi,n−ηR,yi,n,n converges to zero as n→∞ for μ-a.s. y. From this and limn→∞ηR,yi,n,n=ηR,yi<∞ for each i∈I(ℓ), we obtain (5.25).
Clearly, ηR,yi,n and ηR,yi,n,n are independent of πSR(y).
Hence, (5.25) yields (5.26).
From (5.22) and (5.25), we have
[TABLE]
Because the sum in the display equation right above converges, we deduce
We regard f(x,s)∈C1(SR,ϵ[1],k) as a function of the variables
(x,si)i=1k,
symmetric in (si)i=1k∈SRk, and equip C1(SR,ϵ[1],k) with the C1-norm.
Here C1(SR,ϵ[1],k)=C1(SR) for k=0 by convention.
holds for all (x,s)∈SR×S such that x=si for all i.
ii*
For μ-a.s. y and each R,k∈N and each ϵ>0, we have*
[TABLE]
Here Nn and dR,yNN are as in (A2)iii and (5.6), respectively.
Proof.
Let SRm={s∈S;s(SR)=m}.
Let h=f⊗g∈C0∞(SR)⊗D∘b be an FR[1]-measurable function such that
g(s)=0 on ⨆m=M∞SRm for some M∈N, and
h(x,s)=0 whenever ∣x−si∣≤ϵ for some i and some ϵ>0. Then
[TABLE]
Let μR,yNnNn,[1] and yNn be as in (5.1) and (2.23) with N=Nn, respectively.
By (2.4), (2.7), (5.9), and (5.32), we obtain
[TABLE]
Hence, we obtain
[TABLE]
This yields (5.30) for μR,y[1]-a.s. (x,s). The right-hand side of (5.30) is continuous on
{(x,s)∈SR×SR;x=si}.
Hence, there exists a μR,y[1]-version of dR,y such that
(5.30) holds for all {(x,s)∈SR×SR;x=si}.
This yields i.
For μ-a.s. y, by (5.6) and (5.23), we have, x∈ONn,
[TABLE]
Hence, by (2.4) and (2.7),
together with Propositions 5.1 and 5.2,
we obtain
[TABLE]
in C1(SR,ϵ[1],k) for each ϵ>0. From this and i, we obtain ii.
∎
5.3. Explicit formulas for dμy and dμ: Proof of Theorem 2.1
We continue to work under Assumptions (A1)–(A3).
Let dR,y be as in Lemma 5.5. The consistency of dR,y in the following lemma is crucial for demonstrating the explicit expression of the logarithmic derivatives of μ in Theorem 2.1.
Hence, applying (5.31) to (5.35), we obtain (5.34).
∎
Let S0 be as in Definition 2.2.
We introduce the equivalence relation in S0 such that
a∼Tb if and only if (2.37) holds.
Let [y]={z∈S0;y∼Tz}.
Let S[1][y]={(x,s);δx+s∈[y]} and
[TABLE]
where Ss={s∈S;s({x})∈{0,1} for all x∈Rd}.
For μ-a.s. y, we define the function dy(x,s) on S=[1][y] by
[TABLE]
According to (5.34), dy(x,s) is well defined. Indeed, from (5.34), we have
[TABLE]
We take a version of dy(x,s) such that, for any y∈S0,
[TABLE]
We also set dy(x,s)=0 for y∈/S0.
Because the equivalence relation ∼T gives the partition of S0 and μ(S0)=1, such a version of dy is well defined for all (x,s)∈Rd×S.
Thus, for any (x,s)∈Rd×S, we define d by
[TABLE]
Let μy=μ(⋅∣T(S))(y) be the tail decomposition of μ as in Definition 2.2.
Let μy[1] be the one-reduced Campbell measure of μy.
Proposition 5.3**.**
i* For μ-a.s. y, d is
the ∙-logarithmic derivative dμy of μy.*
ii*
d is the ∙-logarithmic derivative dμ of μ and satisfies (2.26).*
Proof.
From (5.12) and (5.37), we obtain that, for all h∈D∘b[1],
[TABLE]
By Lemma 13.6, (5.39) holds for all h∈D∙b[1], and hence
dy is the ∙-logarithmic derivative of μy.
Combining this with (5.38), we obtain i.
Using μ[1]=∫Sμy[1]μ(dy) and (5.39),
we further obtain, for all h∈D∙b[1],
[TABLE]
Hence, d is the ∙-logarithmic derivative of μ.
Finally, from this, (5.37), and (5.38), for μ-a.s. y and
μy[1]-a.e. (x,s),
[TABLE]
Together with (5.30), this yields (2.26) and completes the proof of ii.
∎
**Proof of Theorem 2.1. **
We obtain the ∙-logarithmic derivative dμ∈Lloc2(μ[1]) denoted by (2.26) from Propositions 4.2 and 5.3.
We will prove the remaining claims of i later.
From (5.27), limR→∞ηR,si=0. From (5.17), limR→∞RR,sℓ=0 in C1(SQ).
Hence, we deduce (2.27).
From (2.26) and (2.27), we have, for μ-a.s. s,
[TABLE]
for all Q,R∈N and each ϵ>0.
Hence, we obtain (2.28). Thus, we have obtained ii.
Because (2.28) holds for all ϵ>0, dμ(x,s) has a μ[1]-version such that dμ(x,s) is locally Lipschitz continuous in x on R=d(s):={x∈Rd;x=si}.
Thus, we obtain i.
In addition to (A1)–(A3), we assume (A4) in the proof of iii.
From (2.28) and (2.20), it follows that, for μ[1]-a.e. x,
[TABLE]
This proves iii.
∎
As mentioned before (A3), (A3) is automatically satisfied when ΨN=Ψ.
The following lemma verifies the only nontrivial condition in (A3).
Lemma 5.8**.**
Let ℓ0 be as in (A2).
Let d+ℓ>ℓ0.
Then (2.17) holds.
Proof.
By (2.13), ∫Ss(SR)μ(ds)≤c\ref;A2Rℓ0, R∈N.
Hence,
[TABLE]
which is finite by (2.13).
This together with ΨN=Ψ and (1.1) proves (2.17).
∎
6. Non-collision in ISDEs with infinitely many particles
In this section, we prove the non-collision property of the stochastic dynamics
for infinite-particle systems described by solutions of the ISDE (2.33).
Importantly, we do not rely on properties of the associated quasi-regular
Dirichlet forms, such as capacity.
The main result of this section, Proposition 6.1, applies to weak solutions of the ISDE (2.33).
These solutions are independent of Dirichlet forms.
We apply Proposition 6.1 to the solutions X and Xa constructed in Subsection 8.2 and Subsection 10.1 via lower Dirichlet forms.
All results in Section 6 are proved under (A1)–(A3).
Proposition 6.1**.**
Let X=(Xi)i∈N be a weak solution of (2.33).
We define
Xui=(Xui,∑j=i∞δXuj).
Assume that
[TABLE]
Let Xt=∑i∈NδXti and Ss={s∈S;s({x})∈{0,1} for all x∈Rd}.
Then
[TABLE]
Proof.
Let i=j∈N be fixed.
Without loss of generality, we can assume X0i,X0j∈SR.
We divide the case into two parts: d≥3 and d=2.
Suppose d≥3. Let
[TABLE]
Here, we suppress i,j,R,T∈N from the notation of τϵ.
Let ∑⟨i,j⟩=∑(k,l)=(i,j),(j,i).
From (2.32), (2.33), and Itô’s formula,
[TABLE]
Note that, for X0i,X0j∈SR, we have
0∧(−log∣Xt∧τϵi−Xt∧τϵj∣)≥−log2R.
Combining this with (2.31) and (6.4), we have
[TABLE]
Let c6.1 and c6.2 be the positive constants defined by
[TABLE]
By (6.2) and (6.1), we have c\ref;K5<∞ and c\ref;K5a<∞.
Let
[TABLE]
By (6.5)–(6.7) and the Schwarz inequality, we have
[TABLE]
By (6.8) and d≥3, we have limsupϵ→0xϵ<∞ and limsupϵ→0yϵ<∞.
Hence by Fatou’s lemma and limsupϵ→0xϵ<∞, for each t≥0,
[TABLE]
Hence, from Fatou’s lemma and (6.9), we obtain, for each i=j∈N,
[TABLE]
Let τ~R,Ti,j=inf{0≤t≤T;Xti∈/SR or Xtj∈/SR}.
By (6.10),
We next suppose d=2.
Applying Itô’s formula to log(0∨(−log∣x∣)) instead of −log∣x∣ in (6.4) yields (6.3) for d=2 in a similar fashion.
For the sake of completeness, we present further details.
Let v be such that v(x)=0∨−log∣x∣ for x=0 and v(0)=∞.
Applying Itô’s formula to logv=log(0∨−log∣x∣), we have
[TABLE]
From (2.31) and the display equation right above, we obtain
[TABLE]
Let xϵ,yϵ≥0 be such that
[TABLE]
Let c6.3=E[v(X0i−X0j)].
By (6.2), c\ref;K5c<∞. From (6.11) and (6.12),
[TABLE]
From (6.13), we have limsupϵ→0xϵ<∞ and limsupϵ→0yϵ<∞.
Hence, using Fatou’s lemma,
we obtain that, for each R∈N and t≥0,
[TABLE]
This corresponds to (6.9) for d≥3.
The rest of the proof for d=2 is the same as the case d=3. Thus, we obtain (6.3) for d=2.
∎
Proposition 6.2**.**
Let X=(Xi)i∈N be a weak solution of (2.33). Assume
Hence, applying Proposition 6.1 to X such that X0=l(s), we have
[TABLE]
Integrating the left-hand side with respect to μ∘l−1 yields (6.3).
∎
Let X[m]=(X1,…,Xm,∑j>mδXj) be an m-labeled process.
Suppose that X[m] is a weak solution to the SDE
[TABLE]
Note that this SDE is the restriction of the ISDE (2.33) to the first m components.
Let l[m](s)=(l1(s),…,lm(s),∑i>mδli(s)) for a label l=(li)i∈N.
Proposition 6.3**.**
Let 2≤m<∞. Let X[m] be as above.
Assume that
[TABLE]
Then P(Xtm∈Ss for all 0≤t<∞)=1, where
Xtm=∑i=1mδXti.
Proof.
The proof is the same as Proposition 6.2, so is omitted.
∎
7. Dirichlet forms in finite and infinite volume
A symmetric Dirichlet form (E,D) on L2(ν)
is understood in the standard sense [8, 5] (see Definition 13.2 for details).
To prove Theorem 2.2–2.5, we employ the Dirichlet form approach
[17, 18, 24, 26, 27, 34].
We introduce Dirichlet forms associated with the ISDE (2.33)
and construct infinite-volume Dirichlet forms as limits of
finite-volume ones.
These finite-volume forms naturally split into two classes,
called the lower and upper Dirichlet forms.
Using the family of lower Dirichlet forms, we prove the closability of the
limiting infinite-volume form.
On the other hand, we show that the upper Dirichlet form is quasi-regular,
which ensures the existence of the associated diffusion.
Unless stated otherwise, all results in Section 7 are proved under Assumptions (A1)–(A3).
7.1. Convergence of Dirichlet forms
For μ-a.s. y, let μR,y=μ(⋅∣πSRc(s)=πSRc(y)) as defined previously.
Let μR,y[m] be the m-reduced Campbell measure of μR,y (see (13.6)).
Let SR be the configuration space over SR.
Let SRk={s∈SR;s(SR)=k}.
Lemma 7.1**.**
For μ-a.s. y and each R∈N, the following hold:
i*
For m,k∈{0}∪N such that m+k≥1,
μR,y[m] has a bounded and continuous density mR,y[m] with respect to ΛR[m] on SRm×SRk.*
ii* μR,y has a one-point correlation function ρR,y1.*
Proof.
Let mR,ym+k be the (m+k)-labeled density of μR,y.
From (5.12) and the definition of the one-reduced Campbell measure, we obtain
[TABLE]
in the distributional sense on SRm+k, with test functions in
C0∞(SRm+k).
By (5.30), the logarithmic derivative dR,y admits the representation
[TABLE]
By (5.28), ∑i∈I(ℓ)ηR,yixi+RR,yℓ(x)∈C1(SR).
Combining these yields i.
From (2.13), ∫Ss(SR)μ(ds)<∞.
Then ∫Ss(SR)μR,y(ds)<∞ for μ-a.s. y.
From this, ii follows by the same argument as in Proposition 4.1 ii.
∎
Remark 7.1*.*
Leblé [22, Th. 1] proved the existence of continuous, local density of μ in case of d=2.
For f∈D∙ and s∈S, let
fR,s and xR(s) be as in Definition 13.6.
Let
[TABLE]
We introduce the m-labeled carré du champs on (Rd)m×S.
For m≥1,
[TABLE]
Let μ[m] be the m-reduced Campbell measure of μ defined in (13.6). Let
[TABLE]
Although ER[m], DR∙[m], and DR∘[m] depend on a and μ,
we suppress this dependence in the notation.
Replacing μ[m] with μR,y[m], we write ER,y[m] and DR,y∙[m].
Lemma 7.2**.**
*Let R∈N and m∈{0}∪N.
i(ER,y[m],DR,y∙[m]) is closable on L2(μR,y[m]) for μ-a.s. y.
ii(ER[m],DR∙[m]) is closable on L2(μ[m]).
iii(ER[m],DR∘[m]) is closable on L2(μ[m]).*
Proof.
Let SRk={s∈S;s(SR)=k} and
μR,y,k[m]=μR,y[m](⋅∩SRm×SRk).
Let
[TABLE]
By Lemma 7.1 i, μR,y,k[m] has a bounded and continuous density on SRm×SRk
for each k∈{0}∪N.
Hence, (ER,y,k[m],DR,y∙[m]) is closable on L2(μR,y,k[m]).
Let CnfR,k[m]={(x,s)∈(Rd)m×S;s(SR)=k}.
Then (Rd)m×S can be decomposed as the disjoint union
(Rd)m×S=⨆k∈{0}∪NCnfR,k[m].
Note that
[TABLE]
By (7.3) and (7.4),
(ER,y,k[m],DR,y∙[m]) is closable on L2(μR,y[m]) for each k.
Since (ER,y[m],DR,y∙[m]) is the countable sum of such closable forms,
(ER,y[m],DR,y∙[m]) is closable on L2(μR,y[m]).
For the case m=0, see [24, Lem. 3.2].
By i, we obtain the superposition
(ER[m],DR⋆[m]),
with respect to y under μ, of the closures
(ER,y[m],DR,y[m])
of
(ER,y[m],DR,y∙[m]) on L2(μR,y[m]).
That is, (ER[m],DR⋆[m]) and (ER[m],DR⋆[m])
is the closed form on L2(μ[m]) defined by
[TABLE]
By the superposition theorem for closed forms [3, Prop. 3.1.1],
this superposition yields a closed form on L2(μ[m]), which extends (ER[m],DR∙[m]).
Hence, ii follows from Lemma 13.4.
Finally, (ER[m],DR∙[m]) is clearly an extension of
(ER[m],DR∘[m]).
Therefore, iii follows from ii by Lemma 13.4.
∎
Remark 7.2*.*
In [36], it was proved that the existence of logarithmic derivative implies the closability of bilinear forms in Lemma 7.2. This result is a general theory applicable to all RPFs ν with logarithmic derivative.
Using this, we can prove Lemma 7.2. We present here direct proof using the explicit representation of the logarithmic derivative for μ in Theorem 2.1.
Let (ER,y[m],DR,y[m]), (ER[m],DR[m]), and
(ER[m],DR[m]) be the closures of the closable forms appearing in Lemma 7.2
i, ii, and iii, respectively.
These closed forms are symmetric Dirichlet forms (cf. Definition 13.2).
In the following lemma, we verify that they are strongly local and quasi-regular.
See Subsection 13.2 for the definitions of strong locality and quasi-regularity.
Lemma 7.3**.**
*Let R∈N and m∈{0}∪N.
i(ER,y[m],DR,y[m])
is a strongly local, quasi-regular Dirichlet form on L2(μR,y[m]) for μ-a.s. y, where we regard μR,y[m] as a measure on (Rd)m×S.
ii(ER[m],DR[m]) is a strongly local, quasi-regular Dirichlet form on L2(μ[m]).
iii(ER[m],DR[m]) is a strongly local closed form on L2(μ[m]).*
Proof.
Let (ER,y,k[m],DR,y,k[m]) be the closure of the closable form
(ER,y,k[m],DR,y∙[m]) on L2(μR,y,k[m]),
introduced in the proof of Theorem 7.2.
Then (ER,y,k[m],DR,y,k[m]) is a strongly local, quasi-regular Dirichlet form
on L2(μR,y,k[m]) satisfying the reflecting boundary condition.
The properly associated diffusion is frozen outside CnfR,k[m],
and CnfR,k[m] is its invariant set.
Moreover, the sets CnfR,k[m] are mutually disjoint for k∈{0}∪N and
(Rd)m×S=⊔k=0∞CnfR,k[m].
The Dirichlet form (ER,y[m],DR,y[m]) coincides with the countable sum
of the closed forms (ER,y,k[m],DR,y,k[m]), namely,
[TABLE]
By a gluing argument, the diffusions associated with
(ER,y,k[m],DR,y,k[m]), k∈{0}∪N,
can be glued together to construct on (Rd)m×S a diffusion properly associated with
(ER,y[m],DR,y[m]).
Thus the form (ER,y[m],DR,y[m]) on L2(μR,y[m]) admits a properly associated diffusion.
Hence, by the general theory of quasi-regular Dirichlet forms,
(ER,y[m],DR,y[m])
is strongly local and quasi-regular on L2(μR,y[m]).
This proves i.
Let CnfR,y[m]={(x,s)∈(Rd)m×S;πSRc(s)=πSRc(y)}.
We denote by the same symbol the associated partition of (Rd)m×S.
With CnfR,k[m] replaced by CnfR,y[m] and the countable sum of
(ER,y,k[m],DR,y,k[m]) replaced by the superposition of
(ER,y[m],DR,y[m]),
ii follows in the same way as i.
Since the associated carré du champ consists of differentials,
the form (ER[m],DR[m]) is strongly local,
which yields iii.
∎
We examine the infinite-volume Dirichlet forms.
Clearly, Da,R[f,f](s) is non-decreasing in R for all f∈D∙ and s∈S.
Hence, we define
[TABLE]
We set the carré du champ Da[f,g](s) by polarization for
s such that Da[f,f](s)<∞ and Da[g,g](s)<∞.
For m∈{0}∪N and f,g∈C0∞((Rd)m)⊗D∙, let
[TABLE]
The notion of convergence in the strong resolvent sense
is defined before Lemma 13.5.
Lemma 7.4**.**
*For m∈{0}∪N, the following hold:
i(E[m],Dlwr[m]) is a closed form on L2(μ[m]) and is the limit of
(ER[m],DR[m])
in the strong resolvent sense on L2(μ[m]).
ii(E[m],∪R∈NDR[m]) is closable on L2(μ[m]).
Its closure (E[m],Dupr[m]) is the limit of
(E[m],DR[m]) in the strong resolvent sense on L2(μ[m]).*
Proof.
It is clear that (ER[m],DR[m]) is an extension of (ER+1[m],DR+1[m])
(see Definition 13.1 for the notion of extension).
Hence, by (7.2)–(7.6), (E[m],Dlwr[m]) is the increasing limit of (ER[m],DR[m])
in the sense of Lemma 13.5.
Hence, we obtain i from Lemma 13.5 i.
Since (E[m],Dlwr[m]) is an extension of
(E[m],∪R∈NDR[m]),
the first assertion of ii follows from i and Lemma 13.4.
It is clear that (ER+1[m],DR+1[m]) is an extension of
(ER[m],DR[m]).
Hence, (ER[m],DR[m]) is decreasing with limit
(E[m],∪R∈NDR[m]).
By definition, (E[m],Dupr[m]) is the closure of the maximal closable part of
(E[m],∪R∈NDR[m]) on L2(μ[m]).
Hence ii follows from Lemma 13.5 ii.
∎
We refer to (E[m],Dlwr[m]) and (E[m],Dupr[m]) as the lower and upper Dirichlet forms for μ[m], respectively.
7.2. Dirichlet forms in infinite volume
We now introduce new bilinear forms in addition to the ones we had before.
We define
[TABLE]
Although D∙[m] and D∘[m] depend on a and μ,
we suppress this dependence in the notation.
Lemma 7.5**.**
*For each m∈{0}∪N, the following hold:
i(E[m],D∙[m]) is closable on L2(μ[m]).
ii(E[m],D∘[m]) is closable on L2(μ[m]).*
Proof.
Clearly, (E[m],Dlwr[m]) is an extension of (E[m],D∙[m]).
By Lemma 7.4 i, (E[m],Dlwr[m]) is a closed form on L2(μ[m]).
Hence, i follows from this and Lemma 13.4.
Note that (E[m],D∙[m]) is an extension of (E[m],D∘[m]).
Hence, ii follows from i and Lemma 13.4.
∎
Let (E[m],D∙[m]) and (E[m],D∘[m]) be the closures of (E[m],D∙[m]) and (E[m],D∘[m]) on L2(μ[m]), respectively.
Lemma 7.6**.**
*For each m∈{0}∪N, the following hold.
i(E[m],D∙[m])=(E[m],Dlwr[m]).
ii(E[m],D∘[m])=(E[m],Dupr[m]).*
Proof.
Let DR∙[m] be as in (7.2).
Since D∙[m]⊂DR∙[m] for all R∈N and E[m](f,f)<∞ for all f∈D∙[m], we obtain D∙[m]⊂Dlwr[m].
It is not difficult to show that D∙[m] is dense in Dlwr[m].
This yields i.
From Lemma 7.4 ii, (E[m],Dupr[m]) is the closure of (E[m],∪R∈NDR[m]) on L2(μ[m]).
Note that (E[m],∪R∈NDR[m]) is an extension of (E[m],D∘[m]).
Hence, (E[m],Dupr[m]) is an extension of (E[m],D∘[m]).
From [17, Lem. 2.5], (E[m],D∘[m]) is dense in (E[m],Dupr[m]).
Hence, we have ii.
∎
Lemma 7.7**.**
*For each m∈{0}∪N,
i(E[m],Dlwr[m]) is a strongly local Dirichlet form on L2(μ[m]).
ii(E[m],Dupr[m]) is a strongly local, quasi-regular Dirichlet form on L2(μ[m]).*
Proof.
By Lemma 7.4 i, (E[m],Dlwr[m]) is a closed form on L2(μ[m]).
It is clearly symmetric, and by (7.6),
[TABLE]
Since Da[m] is the carré du champ given by differentials,
(E[m],Dlwr[m]) satisfies (13.2) and is strongly local.
This proves i.
By Theorem 13.2, (E[m],D∘[m]) is a quasi-regular Dirichlet form on L2(μ[m]).
Moreover, by Lemma 7.6 ii,
(E[m],D∘[m])=(E[m],Dupr[m]).
Hence (E[m],Dupr[m]) is quasi-regular on L2(μ[m]).
Its strong locality follows in the same way as in i.
∎
8. Weak solutions of ISDEs for the lower Dirichlet form
The purpose of Section 8 is to construct weak solutions X
to the ISDEs (2.33)
associated with the lower Dirichlet form
(E[m],Dlwr[m]) on L2(μ[m]).
In Section 8.1, we construct the m-labeled process XR[m],
which is properly associated with the lower Dirichlet form
(ER[m],DR[m]) on L2(μ[m]).
Utilizing XR[m], we then construct the fully labeled process
XR
and show that XR
is a weak solution of the ISDEs (8.1) and (8.2) (Lemma 8.1):
[TABLE]
Furthermore, we define the fully labeled process X
as the limit of XR
and show that the m-labeled process
X[m], derived from X,
is associated with the lower Dirichlet form
(E[m],Dlwr[m]) on L2(μ[m])
(Theorem 8.1).
In Section 8.2, we establish that X
indeed solves the ISDEs (2.33)
(Theorem 8.2).
We impose
[TABLE]
Here l=(li)i denotes the labeling map introduced in (2.22).
The initial conditions of these solutions are coupled through the labeling map.
We assume
[TABLE]
All results in Section 8 are proved under Assumptions (A1)–(A3), (2.44), and (8.4).
8.1. The m-labeled process for the lower Dirichlet form
We briefly recall standard notions from Dirichlet form theory.
A Markov process {Xt} is said to be associated with
a Dirichlet form (E,D) on L2(ν)
if {Xt} is associated with the L2-Markovian semigroup
Tt induced by (E,D) on L2(ν),
that is,
Ex[f(Xt)]=Ttf(x)
for ν-a.s. x, all t≥0, and all f∈L2(ν).
We say that {Xt} is properly associated with (E,D) on L2(ν)
if, in addition, (E,D) is quasi-regular and Ex[f(Xt)]
is a quasi-continuous ν-version of Ttf(x).
If (E,D) is strongly local,
then the associated Markov process is a diffusion,
that is, a continuous Markov process with the strong Markov property.
See Subsection 13.2 or [5] for strong locality and quasi-continuity.
Let (ER[m],DR[m]) be the quasi-regular, strongly local Dirichlet form on L2(μ[m]) as in Lemma 7.3.
Let XR[m] be the diffusion properly associated with (ER[m],DR[m]) on L2(μ[m]).
For s(SR)≥m, we take
[TABLE]
Here l(s)=(li(s))i is the label defined as (2.22).
Let W[m]:=C([0,∞);(Rd)m×S)(m≥1) and W[0]:=C([0,∞);S).
For w=(wi)i∈N, let w[m]=(w1…,wm,wm∗)∈W[m], where
wm∗(t)=∑j>mδwj(t).
Let upath[n,m]:W[n]→W[m] for m≤n such that
[TABLE]
Let PR[m] be the distribution of XR[m].
We can easily show that (ER[m],DR[m]), m∈{0}∪N, are consistent in the sense that
[TABLE]
We take the space where XR[m] is defined as W[m].
From (8.5), we can write XR=(XRi)i∈N.
Here, XRi=wi and (wi)i∈N∈W(Rd)N.
We set XR[m]=(XRm,XRm∗), where
XRm=(XRi)i=1m and XRm∗(t)=∑i>m∞δXRi(t).
Let (Bi)i∈N be a sequence of independent d-dimensional standard Brownian motions.
Let rt be the time-reversal operator on the path space on [0,∞) such that
rt(ω)(s)=ω(t−s) if 0≤s≤t and rt(ω)(s)=ω(0) if t≤s.
Let C[m]={(x,s)∈(Rd)m×S;xk=si for some k,i},
where x=(xk)k=1m and s=∑iδsi.
Lemma 8.1**.**
For all R∈N, the following hold:
i*
XR[m] is a weak solution of (8.1) and
XRm∗(t)=XRm∗(0) for all t≥0.*
ii*
XR is a weak solution of ISDE (8.1) and (8.2).*
iii*
Assume (2.45). Then for i∈N,*
[TABLE]
iv*
The set C[m] has zero capacity with respect to (ER[m],DR[m]) on L2(μ[m]).*
Proof.
Without loss of generality, we may assume that m is the number of particles in SR.
Indeed, particles outside SR remain fixed, and the case m smaller than the number of particles in SR reduces to the case of equality.
The Fukushima decomposition is the counterpart of Itô’s formula in the theory of Dirichlet forms [5, Th. 4.2.6]. Because (ER[m],DR[m]) is a quasi-regular Dirichlet form on L2(μ[m]), we can use the Fukushima decomposition.
We regard xi, i=1,…,m, as a function on (Rd)m×S in an obvious manner.
Applying the Fukushima decomposition to xi, we see that XRm=(XRi)i=1m satisfies
[TABLE]
where M[xi] is the martingale additive functional of finite energy and
N[xi] is the continuous additive functional of zero energy of the Fukushima decomposition for the additive functional At[xi]=XRi(t)−XRi(0).
Using a straightforward calculation, we obtain that (cf. [5, pp.164-165])
[TABLE]
and that the zero energy additive functional N[xi] satisfies
[TABLE]
We can prove (8.9) similarly to [8, Example 5.2.2]. In [8, Example 5.2.2], it was assumed that a is the unit matrix; however, generalizing to the current case is straightforward.
From (8.8) and (8.9), XRm satisfies (8.1).
Because (ER[m],DR[m]) has no energy outside SR, we deduce XRm∗(t)=XRm∗(0) for all t≥0. Thus, we obtain i.
We next consider XR=(XRi)i∈N. From (8.7)–(8.9), XRi
satisfy (8.2) for i=1,…,m.
For i>m, SDE (8.2) obviously holds because
XRm∗(t)=∑i>m∞δXRi(t) is frozen outside SR under the dynamics of (ER[m],DR[m]).
This implies ii.
From the Lyons–Zheng decomposition [5, Th. 6.7.2], we write XRi as the sum of
the martingale additive functionals:
[TABLE]
Indeed, from [5, Th. 6.7.2], (8.10) holds Pμ[m]-a.e. up to lifetime, where Pμ[m] is the distribution of XR[m] satisfying XR[m](0)law∼μ[m]. By (2.45), the distribution of XR[m]≪Pμ[m].
Hence, (8.10) holds.
From (8.8) and (8.10), we obtain (8.6). This implies iii.
Clearly, μ[m](C[m])=0.
Because the m-density function of μR,y is bounded on SRm, d≥2, and (ER[m],DR[m]) has no energy outside SR, we obtain iv.
∎
Lemma 8.2**.**
The following hold:
[TABLE]
Here, c8.1>0 is a constant.
Proof.
From (2.44) and (8.4), we can assume φ=1 without loss of generality.
We have
⟨∫0tσ(XRi(u))dBui⟩t=∫0tσtσ(XRi(u))du and σtσ=a. From (2.31), a is uniformly elliptic and bounded.
Hence, from (8.6) and the martingale inequality, we have (8.11)–(8.13).
The proof of (8.13) is not as easy as other claims.
We can prove (8.13) in the same fashion as Lemma 8.7 in [34, pp. 1212–1214] and the detail is omitted here.
∎
For a fully labeled process X=(Xi)i∈N, let X[m]=(Xm,Xm∗) be the m-labeled process defined by Xm=(Xi)i=1m and Xtm∗=∑i>m∞δXti.
Let (E[m],Dlwr[m]) be the Dirichlet form as in Lemma 7.7.
Theorem 8.1**.**
There exists an X satisfying the following:
i*
For each m∈{0}∪N, X[m] is the μ[m]-symmetric, conservative, continuous Markov process associated with (E[m],Dlwr[m]) on L2(μ[m]).*
ii*
For each m∈{0}∪N, X[m] satisfies*
[TABLE]
iii*
X satisfies*
[TABLE]
In particular, each tagged particle Xi of X=(Xi)i∈N does not explode.
Proof.
Let Wm=W(Rd)m and WN=W(Rd)N with the product topology.
From (2.45), XR(0)law∼(φdμ)∘l−1.
Combining this with (8.11), we deduce the tightness of XRi in W(Rd) for each i∈N.
Hence, we obtain the tightness of XRm and XR in Wm and WN, respectively, because
Wm and WN are endowed with the product topology.
Combining the tightness of XR in WN with
(8.13) and using Lemma 13.3, we obtain that XRm∗ is tight in C([0,∞);S).
Hence, XR[m]=(XRm,XRm∗) is tight in W[m].
Let X be any limit point of XR.
According to Lemma 7.4 i, XR[m] converges to X[m] in law in W[m]
and X[m] is a μ[m]-symmetric and continuous Markov process associated with (E[m],Dlwr[m]) on L2(μ[m]).
In particular, X[m] are unique in law and the convergence in (8.14) holds.
Because XR[m]=(XRm,XRm∗) is tight in W[m], the limit continuous Markov process X[m] in (8.14) is conservative.
Because (8.14) holds for all m∈N, the convergence in (8.15) holds and X is unique in law.
The second statement of iii follows from the first since X is a WN-valued random variable by (8.15).
∎
8.2. Weak solutions of ISDEs for the lower Dirichlet form
Let X be as in Theorem 8.1. Subsection 8.2 aims to prove that X is a weak solution of ISDE (2.33) satisfying (SIN), (NBJ), and (AC)μ.
Let a∈Cb2(Rd) and b=21{diva+dμa}
be as in (2.31) and (2.32).
Let σ∈Cb1(Rd) be a matrix-valued function such that σtσ=a,
Lemma 8.3**.**
The following convergence holds in law in W(Rd)N.
[TABLE]
Proof.
Let ⟨M⟩t be the quadratic variation process of a continuous martingale Mt.
By (8.15), σtσ=a, and a∈Cb2(Rd),
we obtain, for all 0≤t<∞,
[TABLE]
From (8.18), we obtain (8.16),
because convergence in law in W(Rd)
of continuous martingales follows from that of their quadratic variations.
From Theorem 2.1, we have
dμ∈Lloc2(μ[1]).
Then bQ:=1SQb∈L2(μ[1])
for any Q∈N.
Hence for any ϵ>0, there exists bQϵ∈Cb(Rd×S) such that
bQϵ=1SQbQϵ and
[TABLE]
Hence for any T,i∈N,
[TABLE]
Here we used (2.44) and (8.4), the fact that XR[1] is a μ[1]-symmetric Markov process, and (8.19).
Similarly, we have, for any T,i∈N,
[TABLE]
Note that bQϵ∈Cb(Rd×S). For any F∈Cb(W(Rd)), we have
[TABLE]
Hence, from (8.15) and (8.22), we obtain, for each i∈N,
[TABLE]
From this, (8.20), and (8.21), we easily obtain (8.17).
∎
Lemma 8.4**.**
For each i,T∈N,
[TABLE]
Proof.
From (8.15), {XQi(⋅)−XQi(0)}Q∈N is tight in W(Rd).
Hence,
[TABLE]
By Lemma 8.3, {∫0⋅σ(XQi(u))dBui}Q∈N and
{∫0⋅b(XQi(u),XQi⋄(u))du}Q∈N are tight in W(Rd).
Similarly as (8.24), we obtain
[TABLE]
From ISDE (8.1) and XRi⋄(u)=∑j=iδXRj(u), we have
[TABLE]
Note that LRi only increases when XRi(u) is on the boundary ∂SR from the second equation of (8.1).
Hence, from this and (8.24)–(8.26), we obtain
Let XRi and Xi be C([0,∞);(Rd)4)-valued random variables such that
[TABLE]
Lemma 8.5**.**
We obtain
[TABLE]
Proof.
We can rewrite (Xi)i∈N as the C([0,∞);(Rd)N)4-valued process:
[TABLE]
We rewrite (XRi)i∈N similarly.
The convergence of XR(⋅)−XR(0) to X(⋅)−X(0) follows from (8.15) and (2.45).
The convergence of the other terms follows from Lemmas 8.3 and 8.4.
From these, we easily obtain (8.28).
∎
Theorem 8.2**.**
X* is a weak solution of (2.33)
satisfying X(0)law∼(φμ)∘l−1, (SIN), (NBJ), and (AC)μ.*
Proof.
Let {hk} be an increasing sequence of positive numbers such that limkhk=∞.
Let Shk={s∈Rd;∣s∣<hk}.
For w=(wi)i=14∈W(Rd)4, let
[TABLE]
Let Vk be the set of continuity points of the map w↦τk(w):
[TABLE]
We can and do take {hk} such that, for all i,k∈N,
[TABLE]
Let Fk∈C0((Rd)4) such that
Fk(x,y,z,u)=x−y−z−u if (x,y,z,u)∈Shk4.
Using Lemma 8.5 and (8.30), we obtain, for each i,k∈N and t∈[0,∞),
[TABLE]
Recall that XR=(XRi)i=1∞ satisfies SDE (8.1).
By the definition of Fk, (8.27), (8.29), and (8.1), we obtain, for each R,i,k∈N and t∈[0,∞),
[TABLE]
From this and (8.31), we obtain, for each i,k∈N and t∈[0,∞),
[TABLE]
which implies that
[TABLE]
By (8.27) and Lemma 8.5, the equation above holds as the identity of continuous processes.
Hence, using limk→∞τk(Xi)=∞ yields ISDE (2.33).
We proceed with (SIN). Let WNE(Ss,i) be as in (2.39).
We will prove X∈WNE(Ss,i) a.s., that is, each tagged particle neither explodes nor collides with any other particle.
From Theorem 8.1 iii, Xi does not explode for each i∈N.
To prove the non-collision property of X, we use Proposition 6.2.
From Theorem 2.1, dμ∈Lloc2(μ[1]).
From Theorem 8.1, the one-labeled process X[1] given by X is a μ[1]-symmetric Markov process.
Hence, the μ[1]-symmetry of X[1] and (2.44) and (8.4) imply that, for all R,T∈N,
[TABLE]
From this, we obtain (6.1).
Note that (6.14) follows from (2.44) and (8.4).
Thus, we have verified the assumptions (6.1) and (6.14) of Proposition 6.2.
Hence, we obtain the non-collision property (6.3).
Thus, we obtain (SIN).
From (8.6) and Theorem 8.1, we have, for each i∈N,
[TABLE]
From (2.31), (8.32), and the martingale inequality, we have (NBJ).
Indeed, we can prove (NBJ) in the same fashion as Lemma 8.7 in [34, pp. 1212–1214].
Assumption 8.81 in [34] follows from (2.13) in the present paper.
The detail is omitted here.
Assumption
(AC)μ is obvious because X is associated with the reversible Dirichlet form
(Eμ,Dμ) on L2(μ) and X0law∼φμ.
∎
9. Weak solutions of ISDEs for the upper Dirichlet form: Proof of Theorem 2.2
The aim of Section 9 is to prove Theorem 2.2.
In Section 9.1, we construct the m-labeled diffusion X[m].
In Section 9.2, we then construct the associated fully labeled process
X
and show that X is a weak solution of (2.33)
satisfying (SIN), (NBJ), and (AC)μ.
These conditions are essential for establishing pathwise uniqueness
and the existence of strong solutions to the ISDE (2.33).
Unlike in Section 8,
the upper Dirichlet form (E[m],Dupr[m])
is quasi-regular.
Therefore, we construct the m-labeled diffusion X[m]
directly from the infinite-volume upper Dirichlet form.
This approach differs markedly from the finite-volume approximation
used in Section 8.
All results in Section 9 are proved under Assumptions (A1)–(A3).
9.1. The m-labeled process for the upper Dirichlet form
Let (E[m],Dupr[m]) be the Dirichlet form on L2(μ[m]) as in Lemma 7.7.
Let Cap[m] denote the capacity with respect to (E[m],Dupr[m])
on L2(μ[m]).
We write Capμ for the case m=0.
Define u[m,0]:(Rd)m×S→S by
u[m,0](x,s)=u(x)+s.
Lemma 9.1**.**
Let A⊂S.
If Capμ(A)=0,
then, for each m∈N,
[TABLE]
Proof.
The proof follows from the argument of [26, Lem. 4.1].
∎
Theorem 9.1**.**
i*
There exists a μ[m]-symmetric, conservative diffusion X[m] properly associated with (E[m],Dupr[m]) on L2(μ[m]) for each m∈{0}∪N.*
ii*
None of the labeled particles Xi, where i=1,…,m, in X[m] explode or collide with other labeled particles:*
[TABLE]
iii*
For each 1≤i≤m,*
[TABLE]
Here rt is the time-reversal operator
defined before Lemma 8.1, and
[TABLE]
Proof.
From Lemma 7.7 ii, (E[m],Dupr[m]) is a quasi-regular, strongly local Dirichlet form on L2(μ[m]).
Hence, we obtain a μ[m]-symmetric diffusion X[m] properly associated with (E[m],Dupr[m]) on L2(μ[m]) by the general theory of Dirichlet forms [5].
Because (E[m],Dupr[m]) is a quasi-regular Dirichlet form, we use the Fukushima decomposition (cf. [5, Th. 4.2.6]).
We regard the coordinate function xi as a function on (Rd)m×S in an obvious manner for 1≤i≤m.
Then xi∈Dupr,loc[m].
Although xi∈/Dupr[m], we can apply the Fukushima decomposition to xi by localization.
Consider the additive functional A[xi](t):=Xti−X0i of X[m].
Let ζ be the lifetime of the diffusion X[m].
Applying the Fukushima decomposition to xi, we have a decomposition of X[m] such that, for 1≤i≤m,
[TABLE]
where M[xi] is the martingale additive functional of finite energy and N[xi] is the continuous additive functional of zero energy [8, 5].
By a standard calculation (cf. [5, pp.164-165]), we can show that M[xi], 1≤i≤m, are
martingales such that
[TABLE]
For any g∈D∘b[m],
[TABLE]
Here, we set s[m]=(s1,…,sm,∑j>mδsj).
Using (9.5) and (9.8) together with localization and [8, Th. 5.2.4], we have
Applying the Lyons–Zheng decomposition [5, p.284] to xi, we see, by (9.10),
[TABLE]
Let
Xm=(Xi)i=1m and
Xtm∗=∑j>m∞δXtj.
Then X[m]=(Xm,Xm∗).
Let
[TABLE]
Then ζ=ζtp∧ζcnf.
Suppose ζ<∞.
Taking t→∞ in (9.11) yields
[TABLE]
The right-hand side is finite. Hence, Xζi is finite.
By limt→ζXt∧ζi=Xζi, Xi does not explode at finite time ζ.
Hence, we obtain ζ=ζcnf.
Because (Eμ,Dμ) is a quasi-regular Dirichlet form on L2(μ), there exists an increasing sequence of compact sets Kq in S such that
[TABLE]
Let a={aq}q∈N and K[aq] be as in (13.3) and (13.4), respectively.
By the property of compact sets of S [11], there exists such an {aq}q∈N satisfying
Hence, using the general theory of Dirichlet forms [8, 5], we obtain
[TABLE]
Here P[m] is the diffusion measure associated with (E[m],Dupr[m]) on L2(μ[m]).
From this, we deduce ζcnf=∞, P[m]-a.s.
We have already proved that ζ=ζcnf. Therefore, ζ=∞, P[m]-a.s.
Thus X[m] is conservative, which implies i.
We obtain (9.1) from ζtp=ζ=∞.
Since X[m] is μ[m]-symmetric and conservative by i, for all R,T∈N,
there exists a constant c9.1 such that
[TABLE]
From this and Proposition 6.3, we obtain (9.2).
This proves ii.
By ζ=∞, iii follows from (9.10) and (9.11). This yields iii.
∎
Remark 9.1*.*
In the proof of Theorem 9.1, we applied the Lyons–Zheng decomposition
directly to the limit process X[m], since the upper Dirichlet form
(E[m],Dupr[m]) is quasi-regular by Lemma 7.7 ii.
At that point, however, the quasi-regularity of the lower Dirichlet form
(E[m],Dlwr[m]) had not yet been established.
For this reason, in the proof of Theorem 8.1 we instead applied the
Lyons–Zheng decomposition to the finite-volume Dirichlet forms (ER[m],DR[m]), which are known to be quasi-regular by Lemma 7.3.
9.2. Weak solutions of ISDEs for the upper Dirichlet form and the proof of Theorem 2.2
In Subsection 9.1, we have constructed the m-labeled process X[m].
The construction of the fully labeled process X for these partially labeled processes have not yet done.
The goal of the subsection is to construct the fully labeled process X and prove that it satisfies ISDE (2.33) (Theorem 9.2) and complete the proof of Theorem 2.2.
The proof in the present section is different from that of Subsection 8.2.
The upper Dirichlet form Dirichlet form (E[m],Dupr[m]) on L2(μ[m]) is quasi-regular.
Hence, we will utilize the Dirichlet form theory and the general theorems developed in [26, 27].
For a label l=(li)i∈N, let l[0](s)=s and
l[m](s)=((li(s))i=1m,∑j>mδlj(s)).
Let lpath:WNE(Ss,i)→W(Rd)N be as in Lemma 13.2. We define
[TABLE]
Here, for w=(wi), we define the measure-valued path w by wt=∑iδwti.
From Lemma 13.2, lpath(w)0=l(w0). Hence, lpath[m](w)0=l[m](w0).
Let Pl[m](s)[m] be the diffusion measure of X[m] staring at l[m](s)
properly associated with (E[m],Dupr[m]) on L2(μ[m]).
By Lemma 7.7, (E[m],Dupr[m]) is a quasi-regular, strongly local symmetric Dirichlet form on L2(μ[m]).
Thus, Lemma 9.2 can be proved by following the proof of [26, Th. 2.4].
∎
For w=(wi)∈W(Rd)N, let upath[m](w)=(w1,…,wm,∑i>mδwi).
Theorem 9.2**.**
i*
There exists a family of probability measures PxN on
W(Rd)N such that, for μ-a.s. s,*
[TABLE]
ii*
Let X=w. Then X under Pφμ∘l−1N is a weak solution of (2.33) with initial distribution φμ∘l−1 that satisfies (SIN), (NBJ), and (AC)μ.*
Proof.
From (9.15), the family of distributions Pl[m](s)[m] is consistent.
Hence, using the Kolmogorov construction theorem, we obtain i.
Let X[m] be as in Theorem 9.1. Then its first m components form a weak solution of (2.33).
By the consistency (9.16), X is a weak solution of (2.33).
From (9.3), for X[m], none of the m-labeled particles explode or collide with the other labeled particles for each m∈N.
Combining this with the consistency (9.16), each tagged particle Xi of X=(Xi)i∈N neither explodes nor collides with the others.
We have thus obtained (SIN).
We obtain (NBJ) from the martingale decomposition (9.4).
The proof is analogous to that of Theorem 8.2 and is omitted.
Finally, (AC)μ is obvious because the unlabeled dynamics X=upath[0](X) are associated with the Dirichlet form on L2(μ), where Xt=∑i∈NδXti.
∎
**Proof of Theorem 2.2. **
From Theorem 9.1, X[m] is a μ[m]-symmetric,
conservative diffusion for each m∈{0}∪N.
From Theorem 9.2, there exists a fully labeled process X
associated with X[m],
and X is a weak solution of (2.33)
with initial distribution
φμ∘l−1.
Here we may take any function φ
satisfying (2.44).
Using Fubini’s theorem,
we obtain X such that
X is a weak solution of (2.33)
starting at l(s)
for μ-a.s. s.
This completes the proof.
∎
10. Unique strong solutions of ISDEs: Proof of Theorem 2.3
To prove Theorem 2.3, we utilize a general result developed in [34, 18].
A key assumption for the uniqueness of solutions is the tail triviality of the associated RPFs.
Starting from the tail decomposition of μ, we introduce the tail decomposition of Dirichlet forms in Section 10.1.
Using this, we complete the proof of Theorem 2.3 in Section 10.2.
Unless stated otherwise, all results in Section 7 are proved under Assumptions (A1)–(A3).
10.1. Tail decomposition of Dirichlet forms and weak solution
Let μa, a∈S0, be the tail decomposition of μ given in Definition 2.2.
Then μ=∫μadμ(a) and μ(S\S0)=0.
For all R∈N, μ-a.s. a, and μa-a.s. y, let
[TABLE]
Let μa,R,y[m] be the m-reduced Campbell measure of μa,R,y.
Let SR be the configuration space over SR and SRk={s∈SR;s(SR)=k}.
Lemma 10.1**.**
For μ-a.s. a, μa-a.s. y, and each R∈N, the following hold:
i*
Let m∈{0}∪N. Then μa,R,y[m] has a density ma,R,y[m] with respect to ΛR[m] such that ma,R,y[m] is bounded and continuous on SRm×SRk for each k∈{0}∪N.*
ii* μa,R,y has a one-point correlation function ρa,R,y1.*
Proof.
Let dR,y and da,R,y be the logarithmic derivatives of μR,y and μa,R,y, respectively.
From Proposition 5.3, dR,y(x,s)=da,R,y(x,s) for μa[1]-a.e. (x,s) for μ-a.s. a.
Hence, from Theorem 5.1 i,
[TABLE]
In Lemma 7.1, the existence of a bounded continuous density
mR,y[m] of μR,y[m] was derived from the explicit representation
(5.30) of dR,y.
Using the above representation of da,R,y,
assertion i follows by the same argument.
Assertion ii can be proved analogously to Lemma 7.1 ii.
∎
Lemmas 10.2–10.7 below correspond to Lemmas 7.2–7.7, respectively.
The proof of Lemma 10.k is similar to that of Lemma 7.k for 1≤k≤6; therefore, we omit those proofs.
Let μa[m] be the m-reduced Campbell measure of μa.
By replacing μ[m] with μa[m] and μa,R,y[m] in (7.2), we define
(Ea,R[m],Da,R∙[m]) and (Ea,R,y[m],Da,R,y∙[m]), respectively.
Lemma 10.2**.**
*For μ-a.s. a, R∈N, and m∈{0}∪N, the following hold:
i(Ea,R,y[m],Da,R,y∙[m])
is closable on L2(μa,R,y[m]) for μa-a.s. y.
ii(Ea,R[m],Da,R∙[m])
is closable on L2(μa[m]).
iii(Ea,R[m],Da,R∘[m]) is closable on L2(μa[m]).*
Let
(Ea,R,y[m],Da,R,y[m]) be the closure of
(Ea,R,y[m],Da,R,y∙[m]) on L2(μa,R,y[m]).
Let
(Ea,R[m],Da,R[m]) and
(Ea,R[m],Da,R[m]) be
the closures of (Ea,R[m],Da,R∙[m]) and (Ea,R[m],Da,R∘[m]) on L2(μa[m]), respectively.
Lemma 10.3**.**
*For μ-a.s. a, R∈N, and m∈{0}∪N, the following hold:
i(Ea,R,y[m],Da,R,y[m])
is a strongly local, quasi-regular Dirichlet form on
L2(μa,R,y[m]) for μ-a.s. y.
ii(Ea,R[m],Da,R[m])
is a strongly local, quasi-regular Dirichlet form on L2(μa[m]).
iii(Ea,R[m],Da,R[m]) is a strongly local closed form on L2(μa[m]).*
Let (Ea[m],Da,lwr[m]) be the increasing limit of (Ea,R[m],Da,R[m]), R∈N.
Lemma 10.4**.**
For μ-a.s. a and m∈{0}∪N, the following hold:
i*
The form (Ea[m],Da,lwr[m]) is a closed form on L2(μa[m]) and is the strong resolvent
limit of (Ea,R[m],Da,R[m]) on L2(μa[m]).*
ii* (Ea[m],∪R∈NDa,R[m]) is closable on L2(μa[m]).
The closure (Ea[m],Da,upr[m]) is the limit of (Ea[m],Da,R[m])
in the strong resolvent sense on L2(μa[m]).*
Let Da∙[m]={f∈C0∞((Rd)m)⊗D∙;Ea[m](f,f)<∞,f∈L2(μa[m])}.
We define Da∘[m] by replacing D∙ with D∘.
Lemma 10.5**.**
*For μ-a.s. a and m∈{0}∪N, the following hold:
i(Ea[m],Da∙[m]) is closable on L2(μa[m]).
ii(Ea[m],Da∘[m]) is closable on L2(μa[m]).*
Let (Ea[m],Da∙[m]) and (Ea[m],Da∘[m]) be the closure of
(Ea[m],Da∙[m]) and (Ea[m],Da∘[m]) on L2(μa[m]), respectively.
Lemma 10.6**.**
*For μ-a.s. a and each m∈{0}∪N, the following hold:
i(Ea[m],Da∙[m])=(Ea[m],Da,lwr[m]).
ii(Ea[m],Da∘[m])=(Ea[m],Da,upr[m]).*
Lemma 10.7**.**
*For μ-a.s. a and m∈{0}∪N, the following hold:
i(Ea[m],Da,lwr[m]) is a strongly local Dirichlet form on L2(μa[m]).
ii(Ea[m],Da,upr[m]) is a strongly local, quasi-regular Dirichlet form on L2(μa[m]).*
Similarly as dR,y for μR,y[1] in Lemma 5.6, we obtain da,R,y for μa,R,y[1].
In Lemma 10.8, we prove that da,R,y coincides with dμ under μa,R,y[1].
Lemma 10.8**.**
For μ-a.s. a, all R∈N, and μa-a.s. y, da,R,y satisfies
[TABLE]
Proof.
By μa,R,y[1]μa(dy)μ(da)=μ[1], we obtain, for h∈C0∞(SR)⊗D∘b,
Let (Ea,R[m],Da,R[m]) be as in Lemma 10.3.
Let Xa,R[m] be the diffusion properly associated with (Ea,R[m],Da,R[m]) on L2(μa[m]).
Proposition 10.1**.**
In addition to (A1)–(A3), we assume (2.45).
Then for μ-a.s. a, there exists a fully labeled process Xa satisfying the following:
i*
For each m∈{0}∪N, the associated m-labeled process Xa[m] is the μa[m]-symmetric and conservative continuous Markov properly associated with
(Ea[m],Da,lwr[m]) on L2(μa[m]).*
ii*
For each m∈{0}∪N, the m-labeled process Xa[m] and
Xa satisfy*
[TABLE]
In particular, each tagged particle Xai of Xa does not explode.
Proof.
The proof is the same as that of Theorem 8.1 and is omitted.
∎
Proposition 10.2**.**
For μ-a.s. a, the following hold.
i*
There exists a μa[m]-symmetric, conservative diffusion Xa[m] properly associated with (Ea[m],Da,upr[m]) on L2(μa[m]) for each m∈{0}∪N.*
ii*
None of the labeled particles Xai, 1≤i≤m, in Xa[m] explode or collide with other labeled particles.*
iii*
We write
Xa[m]=(Xa1,…,Xam,∑j>m∞δXaj).
For each 1≤i≤m,*
[TABLE]
where rt is the time reversal operator defined before Lemma 8.1.
Proof.
The proof is similar to that of Theorem 9.1. Hence, it is omitted.
∎
Let Xa be the fully labeled process constructed in the same fashion as X in Theorem 9.2.
Then Xa[m] in Theorem 10.2 is the m-labeled process of Xa.
Proposition 10.3**.**
For μ-a.s. a,
Xa and Xa are weak solutions of (2.33) with initial distribution φμa∘l−1 satisfying (SIN), (NBJ), and (AC)$${}_{\mu_{\mathsf{a}}}.
Proof.
The proof of Proposition 10.3 is identical to the combination of those for
Theorems 8.2 and 9.2.
Hence, it is omitted.
∎
10.2. Unique strong solution: Proof of Theorem 2.3
From Proposition 10.3, we obtain weak solutions Xa and Xa of (2.33) with initial distribution φμa∘l−1, which satisfy (SIN), (NBJ), and (AC)$${}_{\mu_{\mathsf{a}}}.
Our goal is to prove that Xa[m] and Xa[m]
satisfy (IFC). To this end, we apply Theorem 13.5. More precisely, we
verify assumptions {B1} and {B2} of Theorem 13.5 for
Xa and Xa. These assumptions concern the set
H[a]=⋃k∈N3H[a]k introduced in (13.10)–(13.11).
See Subsection 13.7 for details on {B1} and {B2}.
Once these assumptions are verified, (IFC) follows from Theorem 13.5.
Recall that {B1} is a condition imposed on a general RPF ν.
Lemma 10.9**.**
For μ-a.s. a, the processes Xa and Xa satisfy {B1} with ν replaced by μa.
Proof.
Let K[aq] be the compact set defined in (13.4).
We can easily find a sequence a={aq}q∈N
as in (13.3) satisfying (13.5) and
[TABLE]
Then (10.2) and (13.5) yield μa(S\⋃qK[aq+])=0.
This proves {B1} i.
We take a μa[1]-version of dμa as in Theorem 2.1 and Proposition 5.3.
Here y in Proposition 5.3 is replaced by a.
Then, for μa-a.s. s, dμa(x,s) is locally Lipschitz in x on R=d(s).
Let S=[1][a] be as in (5.36). Let
[TABLE]
where u[1,0](x,s)=δx+s.
Then {B1} ii is satisfied.
In what follows, we denote both Xa and Xa by
the same symbol X, and denote by X the unlabeled dynamics
associated with X.
Let Ss={s∈S;s({x})∈{0,1} for all x∈Rd}.
Note that K[aq] is increasing concerning q from (13.5).
Let
[TABLE]
According to (13.9)–(13.11), condition (B1) iii follows from
[TABLE]
We have already proved (10.4) in Proposition 6.1. Thus, it only remains to prove (10.5).
To verify (10.5), we use Lemma 13.9.
We check all assumptions of Lemma 13.9 for μa:
{UB}, {MF}, {BX}, {Sμa}, and (13.16). See Subsection 13.7.
{MF} is clear since X is associated with a Dirichlet form on L2(μa).
{Sμa} is also clear because the m-labeled process X[m] given by X is associated with a symmetric Dirichlet form on L2(μa[m]).
Let ν=φμa and Qν be the probability such that (X,B) is defined on (Ω,F,Qν,{Ft}).
{BX} makes the following statement: for Qν -a.s.,
[TABLE]
Let M=(Mi)i=1∞ such that Mti=∫0tσ(Xui)dBui.
Because a=σtσ, σ∈Cb1(Rd), and a is uniformly elliptic, we can represent B=(Bi)i=1∞ as a function of M. Thus, B is a σ[Ms;s≤t]-adapted process. Because (X,B) is a solution of ISDE (2.33),
[TABLE]
From (10.7), M is a σ[Xs;s≤t]-adapted process.
Hence, we obtain (10.6).
By Lemma 13.8 and (10.2), we may use (13.16).
Thus all assumptions of Lemma 13.9 are satisfied, and hence (10.5) follows.
This proves {B1} iiii.
∎
Let H[a]k be as in (13.11).
We choose a in the definition of H[a]k as in Lemma 10.9 so that H[a]k satisfies {B1}.
Lemma 10.10**.**
For μ-a.s. a, H[a]k satisfies {B2}.
Proof.
Let (σ,b) be the coefficient of ISDE (2.33).
We define σm and bm by
[TABLE]
for x=(x1,…,xm) and s∈S.
We verify that σm and bm in (10.8)
satisfy (13.13) and (13.14).
By assumption, σ∈Cb1(Rd) and that σ is independent of s.
Hence, σm is Lipschitz continuous, which yields (13.13).
By (5.18), we have RR,sℓ∈C1(SR). From (2.26), we obtain
[TABLE]
Hence, dμ is Lipschitz continuous on H[a]k from (13.10) and (13.11).
Then from a∈Cb2(Rd), b=21{diva+dμa} is Lipschitz continuous on H[a]k.
Hence, by (10.8), bm is Lipschitz continuous on H[a]k, which yields (13.14).
∎
Proposition 10.4**.**
For μ-a.s. a, Xa and Xa satisfy (IFC).
Proof.
By Lemma 10.9, Xa and Xa satisfy {B1}.
By Lemma 10.10, H[a]k satisfies {B2}.
Hence, (IFC) follows from Theorem 13.5.
∎
Theorem 10.1**.**
For μ-a.s. a, the following hold:
i*
For φμa∘l−1-a.s. s,
the processes Xa and Xa are the unique strong solutions of (2.33)
starting at s under the constraints
(SIN), (NBJ), (AC)$${}_{\mu_{\mathsf{a}}}, (MF), and (IFC).
Moreover, the labeled dynamics Xa and Xa are indistinguishable.*
ii*
The m-labeled processes
Xa[m] and Xa[m],
induced by Xa and Xa, respectively,
are conservative, μa[m]-symmetric diffusions properly associated
with the same strongly local, quasi-regular Dirichlet form on L2(μa[m]):*
[TABLE]
Here (Ea[m],Da,lwr[m]) and (Ea[m],Da,upr[m]) are given by Lemma 10.4.
Proof.
Applying Theorem 13.3 to μa, we shall prove Theorem 10.1. Hence, we verify the assumptions of Theorem 13.3 to μa.
Note that φμ has a tail decomposition such that (φμ)a=φμa for μ-a.s. a.
By construction, μa is tail trivial.
Applying Proposition 10.3 to μa, we have weak solutions Xa and Xa to (2.33) with initial distribution φμa∘l−1 that satisfy (SIN), (NBJ), and (AC)$${}_{\mu_{\mathsf{a}}}.
From Proposition 10.4, Xa and Xa satisfy (IFC).
Thus, all assumptions of Theorem 13.3 are fulfilled.
Hence, according to Theorem 13.3, Xa and Xa are unique strong solutions of (2.33) under the constraints (SIN), (NBJ), (AC)$${}_{\mu_{\mathsf{a}}}, (MF), and (IFC).
In particular, Xa and Xa are indistinguishable.
This implies i.
Let Xa[m] and Xa[m] be the m-labeled processes of Xa and Xa, respectively.
Because Xa and Xa are indistinguishable, we obtain that Xa[m] and Xa[m] are indistinguishable.
From Propositions 10.1 and 10.2, Xa[m] and Xa[m] are associated with the m-labeled Dirichlet forms (Ea[m],Da,lwr[m]) and (Ea[m],Da,upr[m]) on L2(μa[m]), respectively.
Because Xa[m] and Xa[m] are indistinguishable, we obtain (10.9).
From Theorem 10.2, (Ea[m],Da,upr[m]) is a strongly local, conservative, and quasi-regular Dirichlet form on L2(μa[m]). Hence, from (10.9), (Ea[m],Da,lwr[m]) is also a strongly local, conservative, and quasi-regular Dirichlet form on L2(μa[m]). This implies ii.
∎
Remark 10.1*.*
The method of using the uniqueness of strong solutions of ISDEs
to prove the uniqueness of extensions of Dirichlet forms
for interacting Brownian motions
is due to Tanemura [43].
This approach was subsequently adopted in [17].
**Proof of Theorem 2.3. **
Let Xa and Xa be as in Theorem 10.1.
By Theorem 10.1, they are indistinguishable.
Hence, we write X:=Xa=Xa.
By Theorem 10.1, for μ-a.s. a and for φμa∘l−1-a.s. s,
(2.33) admits a unique strong solution X starting at s under the constraints
(SIN), (NBJ), (AC)$${}_{\mu_{\mathsf{a}}}, (MF), and (IFC).
Since φ ranges over all functions satisfying (2.44), this remains valid for
s=l(s) for μa-a.s. s.
For each m∈{0}∪N, the m-labeled process X[m] induced by X is properly associated with the strongly local, quasi-regular Dirichlet form in (10.9) on L2(μa[m]).
Hence, X[m] is a μa[m]-symmetric diffusion.
By Proposition 10.1, X[m] is conservative.
This completes the proof of Theorem 2.3.
∎
The aim of this section is to prove the uniqueness of extensions of Dirichlet forms
and to construct the diffusion on (Rd)N.
We first establish the uniqueness of Dirichlet forms and, using this result, show that the fully labeled process
X=(Xi)i∈N constructed in Theorem 2.3 is a conservative diffusion on (Rd)N.
This completes the proof of Theorem 2.4.
Unless stated otherwise, all results in Section 11 are proved under Assumptions (A1)–(A3).
11.1. Uniqueness of extensions of Dirichlet forms
We now analyze the Dirichlet forms directly associated with solutions of the ISDE (2.33)
to prove the uniqueness of extensions of Dirichlet forms asserted in Theorem 11.1.
Let μa be as in Definition 2.2. Let μa[m] be the m-Campbell measure of μa. Recall that
[TABLE]
Let (Ea,R,y[m],Da,R,y[m]) be as in Lemma 10.3.
Let (Ea,R[m],Da,R⋆[m]) be the superposition of (Ea,R,y[m],Da,R,y[m])
with respect to y under the measure μa:
[TABLE]
where [a] is the set in Definition 2.2. We recall that μa([a])=1.
By the general theorem on the superposition of closed forms [3],
(Ea,R[m],Da,R⋆[m]) is a closed form on L2(μa[m]). Let (Ea,R[m],Da,R[m]) be as in Lemma 10.3.
Lemma 11.1**.**
Let m∈{0}∪N. For μ-a.s. a and all R∈N,
[TABLE]
Proof.
Let Xa,R[m] and Xa,R,y[m] be the m-labeled diffusions associated with the Dirichlet forms
(Ea,R[m],Da,R[m]) on L2(μa[m]) and (Ea,R,y[m],Da,R,y[m]) on L2(μa,R,y[m]), respectively.
Let Xa,R and Xa,R,y be the fully labeled processes associated with the diffusions
Xa,R[m] and Xa,R,y[m], respectively.
By the argument preceding Lemma 8.1, such processes exist and are unique.
By Lemma 8.1, both fully labeled processes Xa,R and Xa,R,y
are weak solutions of the pair of SDEs (8.1) and (8.2).
From the representation of the logarithmic derivative dR,y(x,s) given by (5.6), the pair of SDEs (8.1) and (8.2) has a unique strong solution.
Indeed, the coefficients given by the logarithmic derivative are smooth outside
C[1]={(x,s)∈Rd×S;x=si for some i}, where s=∑iδsi, and each particle is subject to the reflecting boundary condition on ∂SR.
Hence, (8.1) admits a unique strong solution until the solution hits C[1], which never occurs by Lemma 8.1 iv.
Here, Lemma 8.1 iv is clearly valid for (Ea,R,y[1],Da,R,y[1]) on
L2(μa,R,y[1]) as well, similarly to (Ea,R[1],Da,R[1]) on
L2(μa[1]).
Thus, (8.1) admits a unique strong solution.
Equation (8.2) means that all particles outside SR are frozen.
Hence, (8.2) clearly admits a trivial unique strong solution.
Let CnfR,y[m] be as in the proof of Lemma 7.3.
Then CnfR,y[m] is an invariant set of the diffusion Xa,R,y[m] properly associated with
(Ea,R,y[m],Da,R,y[m]) on L2(μa,R,y[m]).
Recall that (Ea,R[m],Da,R⋆[m]) is the superposition of (Ea,R,y[m],Da,R,y[m])
with respect to y under μa, and that CnfR,y[m] yields a partition of (Rd)m×S.
Hence the family of diffusions Xa,R,y[m] associated with
(Ea,R,y[m],Da,R,y[m]) on L2(μa,R,y[m])
gives rise to a diffusion associated with (Ea,R[m],Da,R⋆[m]) on L2(μa[m]).
Here we used (11.1).
As observed at the beginning of the proof,
Xa,R,y[m] is the unique solution of (8.1) and (8.2).
Hence, for each y∈S, the diffusion obtained from
(Ea,R[m],Da,R⋆[m]) is this unique solution.
Recall that Xa,R is a weak solution of SDEs (8.1) and (8.2).
Hence, the diffusion Xa,R,y[m], y∈S, obtained from (Ea,R[m],Da,R⋆[m]) and the diffusion Xa,R[m] associated with (Ea,R,y[m],Da,R,y[m]) coincide.
This implies the identity of the Dirichlet forms in (11.2).
∎
We see that (Ea,R[m],Da,R⋆[m]), R∈N, is increasing.
Let (Ea[m],Da,lwr⋆[m]) be the increasing limit of (Ea,R[m],Da,R⋆[m]).
By Lemma 13.5 i, (Ea[m],Da,lwr⋆[m]) is a closed form on L2(μa[m]).
Let (Ea[m],Da,lwr[m]) be as in Lemma 10.4.
Recall that (Ea,R[m],Da,R⋆[m]) is a closed form on L2(μa,R,y[m]) and that
D∘b[m]⊂Da,R,y[m].
Then (Ea,R,y[m],D∘b[m]) is closable on L2(μa,R,y[m]).
Let (Ea,R,y[m],Da,R,y[m]) be its closure.
Let (Ea,R⋆[m],Da,R⋆[m]) be the superposition of
(Ea,R,y[m],Da,R,y[m]) with respect to y under μ.
Let (ER[m],D⋆,R[m]) and
(ER,D⋆,R⋆[m]) be the superpositions of
(Ea,R[m],Da,R[m]) and (Ea,R[m],Da,R⋆[m]), respectively, with respect to a under μ.
Lemma 11.3**.**
Let m∈{0}∪N and R∈N. Then the following hold:
[TABLE]
Proof.
Since the carré du champ of ER[m] vanishes outside SRc, we obtain (11.4).
Equation (11.5) follows immediately from (11.4).
∎
Let (Ea[m],Da,upr[m]) be the closed form given in Lemma 10.4.
Let (Ea[m],Da,upr⋆[m]) be the decreasing limit of (Ea,R[m],Da,R⋆[m]).
Lemma 11.4**.**
Let m∈{0}∪N and R∈N. Then the following hold:
[TABLE]
Proof.
(11.6) follows immediately from Lemma 10.4 and (11.5).
∎
Let (E[m],D⋆,upr[m]),
(E[m],D⋆,upr⋆[m]),
(E[m],D⋆,lwr[m]), and
(E[m],D⋆,lwr⋆[m])
be the superpositions of
(Ea[m],Da,upr[m]), (Ea[m],Da,upr⋆[m]),
(Ea[m],Da,lwr[m]), and (Ea[m],Da,lwr⋆[m]), respectively,
with respect to a under μ.
[TABLE]
Replacing (Ea[m],Da,upr[m]) with (Ea[m],Da,upr∗[m]),
(Ea[m],Da,lwr[m]), and (Ea[m],Da,lwr⋆[m]), we obtain the superpositions
(E[m],D⋆,upr⋆[m]),
(E[m],D⋆,lwr[m]) and
(E[m],D⋆,lwr⋆[m]), respectively.
The following theorem shows that these Dirichlet forms coincide and are strongly local and quasi-regular.
Theorem 11.1** (Uniqueness of extensions of Dirichlet forms).**
i*
Let S0 be the set in Definition 2.2. For μ-a.s. a∈S0,*
[TABLE]
ii*
All Dirichlet forms in (11.7) and (11.8) are strongly local and quasi-regular on L2(μa[m]) and L2(μ[m]), respectively.*
Proof.
From (11.6), (10.9) and (11.3), we obtain (11.7).
Taking the superposition of (11.7) with respect to a under μ yields (11.8).
This proves i.
By Lemma 10.7 ii, (Ea[m],Da,upr[m]) is strongly local and quasi-regular on L2(μa[m]).
Hence, all Dirichlet forms in (11.7) are strongly local and quasi-regular.
Let [a]={b∈S0;μa=μb} be as in Definition 2.2.
We have μa([a])=1, and [a]∩[b]=∅ if [a]=[b].
Note that [a] gives a partition of S0 and μ(S0)=1 by construction.
Let Capa[m] be the capacity given by (Ea[m],Da,upr[m]) on
L2((Rd)m×[a],μa).
There exists an increasing sequence of compact sets Kan satisfying
[TABLE]
Indeed, (11.9) follows from (Q1) in Definition 13.4.
Here, we take, for each n∈N,
[TABLE]
Let Pa[m] be the diffusion properly associated with
the Dirichlet form (Ea[m],Da,upr[m]) on L2(μa[m]),
which is strongly local and quasi-regular by Lemma 10.7 ii.
From (11.9), the state space of Pa[m] can be taken to be disjoint for [a]=[b].
Thus, the collection of diffusions Pa[m], a∈S0, gives the diffusion properly associated with (E[m],D⋆,upr[m]) on L2(μ[m]).
Hence, from [5, Th. 1.5.3], (E[m],D⋆,upr[m]) on L2(μ[m]) is quasi-regular.
The strong locality of (E[m],D⋆,upr[m]) is clear by construction.
Thus, all Dirichlet forms in (11.8) are strongly local and quasi-regular on L2(μ[m]).
This completes the proof of ii.
∎
Theorem 11.1 plays a crucial role in the proof of Theorems 2.4–2.5.
We present further applications of Theorem 11.1 in the following remark.
Remark 11.1*.*
i
The sub-diffusivity of tagged particles in the Ginibre interacting Brownian motion
(d=β=2) follows from (11.7) [31], reflecting the long-range nature
of the two-dimensional Coulomb potential.
In contrast, tagged particles are always diffusive for Ruelle-class potentials with
a convex hardcore in dimensions d≥2 [25].
ii
Suzuki proved that the unlabeled Markov processes associated with a Dirichlet form are ergodic if the reference measure μ is number rigid and the Dirichlet form is locally ergodic [40].
This theorem can be applied to the Ginibre RPF.
In [40], Suzuki considered the Dirichlet form that corresponds to (E[m],D⋆,lwr⋆[m]).
By Theorem 11.1, the continuous Markov process considered in [40] is an ergodic diffusion.
iii Suzuki [41] studied curvature bounds for Dyson’s Brownian motion, an infinite-particle system on R with the two-dimensional Coulomb interaction at inverse temperature β=2, that is, with a Riesz potential. He analyzed Dirichlet forms of the type (E[m],D⋆,lwr⋆[m]). Since our argument for Theorem 11.1 may be valid for Riesz potentials, his results can apply to Dyson’s model studied in [27, 34] by the uniqueness of Dirichlet forms.
Remark 11.2*.*
Kawamoto proved the identity (11.11) for Dyson’s model [13]:
[TABLE]
This identity, together with the uniqueness theorem for solutions of ISDEs in [34],
implies that the tail σ-field of μ is invariant under the dynamics of the associated diffusion.
Recently, Kawamoto has announced an extension of this identity to a broader class of RPFs
(see [15]).
It is therefore plausible that the same identity also holds for the Coulomb RPFs.
Utilizing the uniqueness results of solutions of ISDEs (Theorem 10.1) and Dirichlet forms (Theorem 11.1), we establish Theorem 2.4 in this subsection. The goal is to prove that the fully labeled process X in Theorem 2.3 is an (Rd)N-valued diffusion.
Let Xa and Xa be as in Theorem 10.1.
In the proof of Theorem 2.3, we constructed X as X:=Xa=Xa by proving Xa and Xa are indistinguishable.
Let X[m] be the m-labeled process of X, see (2.34).
In Theorem 2.3, we have already proved that X[m] is a μa[m]-symmetric diffusion
for μ-a.s. a.
However, this does not necessarily imply that the family X[m] indexed by a∈S
is a μ[m]-symmetric diffusion.
Indeed, let Sa[m] be an invariant set of the
μa[m]-symmetric diffusion X[m].
Then
μa[m]((Rd)m×S\Sa[m])=0.
However, it does not necessarily hold that
[TABLE]
where [a] is the equivalence class introduced in (2.37).
If (11.12) holds, then the family X[m] indexed by a is a diffusion
with state space ∪aSa[m],
because each Sa[m] is an invariant set
of the conservative diffusion X[m].
In the following theorem, we prove (11.12) for Ka[m],
defined in the proof, and show that the m-labeled process X[m]
is a diffusion with state space (Rd)m×S.
Proposition 11.1**.**
i* For μ-a.s. a,
X[m] is a μa[m]-symmetric, conservative diffusion properly associated with the Dirichlet form on L2(μa[m]) in (11.7) for each m∈{0}∪N.*
ii*
X[m] is a μ[m]-symmetric, conservative diffusion properly associated with the Dirichlet form on L2(μ[m]) in (11.8) for each m∈{0}∪N.*
Proof.
By Theorem 10.1 ii, X[m] is a μa[m]-symmetric diffusion
properly associated with (Ea[m],Da,lwr[m]) on L2(μa[m]).
Hence by (11.7), we obtain i.
Let Kan be the compact set given in the proof of Theorem 11.1. Let [a] be as in (2.37).
Kan is increasing in n for each a and that Kan⊂[a] and
[TABLE]
Let Ka[m]=(Rd)m×∪n=1∞Kan.
By (11.10), Ka[m]=Kb[m] if and only if [a]=[b].
By (11.13), Ka[m]∩Kb[m]=∅ if [a]=[b].
By the proof of Theorem 11.1, Ka[m] is the state space of the diffusion properly associated with (Ea[m],Da,upr[m]) on L2((Rd)m×[a],μa[m]).
By Theorem 10.1 ii, Xa[m] is properly associated with (Ea[m],Da,upr[m]).
Thus, Ka[m] is an invariant set of the diffusion Xa[m].
Since X[m] coincides with Xa[m] on Ka[m], X[m] is a μ[m]-symmetric, conservative diffusion properly associated with the superposition (E[m],D⋆,upr[m]) of (Ea[m],Da,upr[m]) on L2(μ[m]) with the state space ⊔[a]Ka[m].
Here {Ka[m]}[a] gives a partition of the state space into Ka[m]. Each Ka[m] is an invariant sets of X[m]. We have thus proved ii.
∎
**Proof of Theorem 2.4. **
Redefining S0 in Proposition 11.1 if necessary, we may assume that
its conclusion holds for all a∈S0, with μ(S0)=1.
Let X be the fully labeled process in Theorem 2.3.
Let X=X[0] be the zero-labeled process given by X.
Recall that μ=μ[0].
From Proposition 11.1 ii, X is a μ-reversible diffusion with state space S⋆⋆ properly associated with the Dirichlet form
(Eμ,D⋆,lwr⋆μ) on L2(μ).
S⋆⋆ is given by
[TABLE]
where S0 is as in Definition 2.2 and Kan is given by (11.9) and
∪n=1∞Kan=Ka[0].
From Theorem 11.1 ii, (Eμ,D⋆,lwr⋆μ) on L2(μ) is strongly local and quasi-regular.
The state space S⋆⋆ satisfies Cap(S⋆⋆c)=0,
where Cap is the capacity associated with (E[m],D⋆,lwr⋆[m]) on L2(μ).
This follows from that Ka[0] is an invariant set of the diffusion Xa=Xa[0] for each a∈S0.
Because X=Xa, Ka[0] is the invariant set of the diffusion X for each a∈S0.
We can take a quasi-continuous version of the diffusion measure Px of X such that Px is defined for all x∈S⋆⋆. Note that Px(X∈C([0,∞);S⋆⋆))=1 because S⋆⋆ is the state space of the conservative diffusion X.
Let lpath:WNE(Ss,i)→W(Rd)N be as in Lemma 13.2.
Let upath be the map defined by upath(w)=w, where w(t)=∑iδwi(t) for w=(wi)i.
upath∘lpath is the identity map. Hence, lpath is injective. We note that X=lpath(X).
The value lpath(w)(t)∈(Rd)N is determined by l(w(0)) and w(u), 0≤u≤t.
Hence, lpath(w)(t) is Ft=σ[l(w(0))]∨σ[w(u),0≤u≤t]-measurable.
Thus, the time-shifted path lpath(w)(⋅+a) is determined by lpath(w)(a).
Let S⋆⋆⊂(Rd)N be the set defined by
[TABLE]
We note that u(S⋆⋆)=S⋆⋆.
Let Px=P(X∈⋅∣X0=x) be the distribution of X starting at x∈u−1(S⋆⋆). As before, let Px be the distribution of X starting at x.
We have
[TABLE]
Let σ be an arbitrary {Ft}-stopping time such that σ<∞ a.s.
Let τ=σ∘lpath and Gt=lpath−1(Ft).
We have Gτ=lpath−1(Fσ).
lpath is injective and X is a diffusion with state space S⋆⋆.
For each x∈S⋆⋆ such that x=Xσ(X),
[TABLE]
Here, we used (11.15) for the second line, the strong Markov property of X for the forth line, and the injectivity of l and lpath for the last line. Because X is μ-reversible and satisfies (SIN), X has infinite lifetime. Hence, so is X=lpath(X).
Thus, we conclude that X is a conservative diffusion with state space S⋆⋆ such that
μ(u(S⋆⋆))=1.
∎
12. Finite-volume approximation: Proof of Theorem 2.5
All results in Section 12 are proved under Assumptions (A1)–(A3) and (2.44)–(2.45).
12.1. Finite-domain approximation
We now impose reflecting boundary conditions on ∂SR
in (2.42) and fix the particles outside the region
SR={s∈Rd;∣s∣≤R}.
We suppose SR⊂ON. Then the dynamics are described as follows:
[TABLE]
Here 21nRa denotes the inward unit normal vector on ∂SR,
and LRN,i is the local time on the boundary ∂SR.
Let {Nn}n=1∞ be the sequence in (A2) such that μNn converges to μ.
Taking the limit Nn→∞ in the SDE (12.1)–(12.2), we obtain the ISDE
(8.1)–(8.2).
Remark 12.1*.*
i
The solutions of the above SDE remains fixed outside the region SR.
This follows from (2.22) and (12.2).
ii
Let s=(si)i=1N satisfy si=sj for all i=j.
Then the SDE (12.1)–(12.2) with XRN(0)=s admits a unique strong solution XRN.
Indeed, the coefficients are smooth outside the collision set
N=⋃i=j{si=sj}, and XRN does not hit
N since d≥2; see Proposition 6.1.
Lemma 12.1**.**
i* The sequence XRN(N,R∈N) is tight in W(Rd)N.*
ii*
The sequence XN,[m](N∈N) is tight in W((Rd)m×S) for m∈{0}∪N.*
Proof.
We suppose N≥s(SR).
We only consider the case i≤s(SR) because XRN,i(t) is frozen for i>s(SR).
The unlabeled dynamics XRN are μR,sN-reversible.
Suppose that φ=const.
Then from the Lyons-Zheng decomposition [5, p.284],
[TABLE]
From this, (2.31), (2.44)–(2.45), and the martingale inequality, we obtain
[TABLE]
where c12.1 is a constant independent of N,R∈N.
By the Lyons-Zheng decomposition, XN is tight in W(Rd)N. This implies i.
Moreover, by the Lyons-Zheng decomposition, any limit point satisfies (NBJ).
From this and Lemma 13.3, we obtain ii.
∎
Recall that, for y∈S, yNn=(l1(y),…,lNn(y)), where l=(li)i is the label as in (2.22). Thus, yNn is a function of y.
Lemma 12.2**.**
For μ-a.s. y, the following hold:
i*
The Dirichlet form (ER,yNnNn,[m],DR,yNnNn,[m])
on L2(μR,yNnNn,[m]) converges to the Dirichlet form
(ER,y[m],DR,y[m]) on L2(μR,y[m])
in the strong resolvent sense
after identifying L2(μR,yNnNn,[m]) and L2(μR,y[m])
with their images in L2(ΛR[m])
under the canonical isometric embeddings induced by the Radon–Nikodym derivatives.*
ii*
The Dirichlet form (ER,yNnNn,[m],DR,yNnNn,[m])
on L2(μR,yNnNn,[m]) converges to the Dirichlet form
(ER,y[m],DR,y[m]) on L2(μR,y[m])
in the strong resolvent sense on L2(ΛR[m]) under the identification in i.*
Proof.
Let mR,yNnNn,k and mR,yk
be the k-density functions of μR,yNnNn and μR,y, respectively, on
SRk. By (5.31), mR,yNnNn,k converges to mR,yk
in C(SRk).
Let
pR,yNnNn,k(t,x,y)mR,yNnNn,k(y)dy
and
pR,yk(t,x,y)mR,yk(y)dy
be the transition probabilities of XRNn and XR on SRk, respectively.
For each R∈N and k∈N, these transition densities exist by the classical theory of heat kernels in divergence form.
We set
[TABLE]
By the classical theory of heat kernels in divergence form,
pR,yNnNn,k and pR,yk
are locally Hölder continuous on O, and pR,yNnNn,k converges to
pR,yk uniformly on each compact subset of O.
Moreover, mR,yNnNn,k converges to mR,yk in C(SRk), as noted at the beginning of the proof.
Consequently, the corresponding α-resolvent densities converge, and i follows.
The Dirichlet form in ii is associated with a different boundary condition.
Nevertheless, its heat kernel converges locally uniformly in the same manner as that in i.
This implies the assertion of ii.
∎
Let (ERNn,[m],DR⋆Nn,[m]) and (ERNn,[m],DR⋆Nn,[m]) be the superpositions, with respect to y under μ, of
(ER,yNn,[m],DR,yNn,[m]) and
(ER,yNn,[m],DR,yNn,[m]), respectively.
Lemma 12.3**.**
i*
The Dirichlet form
(ERNn,[m],DR⋆Nn,[m]) on L2(μRNn,[m]) converges to the Dirichlet form
(ER[m],DR⋆[m]) on L2(μR[m])
in the strong resolvent sense on L2(ΛR[m])
after identifying L2(μRNn,[m]) and L2(μR[m])
with their images in L2(ΛR[m])
under the canonical isometric embeddings induced by the Radon–Nikodym derivatives.*
ii*
The Dirichlet form (ERNn,[m],DR⋆Nn,[m]) on L2(μRNn,[m])
converges to the Dirichlet form (ER[m],DR⋆[m])
on L2(μR[m]) in the strong resolvent sense on
L2(ΛR[m]) under the identification in i.*
Proof.
By the definition of superposition together with Fubini’s theorem and Jensen’s inequality,
Lemma 12.3 follows from Lemma 12.2.
∎
We define the particles with indices greater than Nn to remain fixed for all time.
With this convention, the dynamics of the first Nn particles can be regarded as an infinite-particle system.
Under this convention, all solutions start from the same initial condition l(s) almost surely.
Since XR is a solution of SDE (8.1)–(8.2), it follows that XR[m] is associated with (ER[m],DR⋆[m]) on L2(μR[m]).
Hence, by Lemma 12.1 i and Lemma 12.3 i, we obtain (12.3).
By Theorem 8.1 iii, we obtain (12.4), which completes the proof.
∎
12.2. Finite-particle approximation: Proof of Theorem 2.5
**Proof of Theorem 2.5. **
Let XN be the solution of SDE (2.42) with (2.44)–(2.45).
For each m∈{0}∪N with m≤N, let
[TABLE]
Then (EN,[m],D∘b[m]) is closable on L2(μN,[m]).
Let (EN,[m],DN,[m]) be its closure.
The diffusion XN,[m] is properly associated with (EN,[m],DN,[m]) on L2(μN,[m]).
By the Lyons–Zheng decomposition and (2.44)–(2.45),
XN and XN,[m], N∈N, are tight in W(Rd)N and W((Rd)m×S), respectively.
Hence, there exist subsequences of XNn and XNn,[m],
still denoted by the same symbols, and limits X and
X[m] such that X[m] is the
m-labeled process of X. Furthermore, for any nonnegative
f,g∈D∘b[m] and α>0,
[TABLE]
Here Rα,Nn,[m] denotes the α-resolvent of (EN,[m],DN,[m])
on L2(μN,[m]) and
[TABLE]
At this stage, it is not known whether X[m] is a Markov process.
We now prove that this is indeed the case.
Let (ERNn,[m],DR⋆Nn,[m]) and (ERNn,[m],DR⋆Nn,[m]) be as in Lemma 12.3.
For symmetric bilinear forms A and B,
A≻B means that A is an extension of B in the sense of (13.1).
Then it is clear that
[TABLE]
Hence,
[TABLE]
By Lemma 12.3 i, (ERNn,[m],DR⋆Nn,[m]) converges to
(ER[m],DR⋆[m]) in the strong resolvent sense as n→∞.
By Lemma 13.5 i, (ER[m],DR⋆[m]) converges to the increasing limit (E[m],Dlwr⋆[m]) in the strong resolvent sense as R→∞.
Thus, the left closed form in (12.6) converges to (E[m],Dlwr⋆[m])
in the strong resolvent sense in the two-step limits.
Let f∈DRNn,[m].
Then f is B((SR)m)×σ[πSR]-measurable.
Moreover, μRNn,[m] is absolutely continuous with respect to μR[m].
Hence, we regard (ERNn,[m],DR⋆Nn,[m]) as a closed form on L2(μR[m]).
By Lemma 12.3 ii, (ERNn,[m],DR⋆Nn,[m]) converges to (ER[m],DR⋆[m])
in the strong resolvent sense. Moreover, by Lemma 13.5 ii, (ER[m],DR⋆[m]) converges to the decreasing limit (E[m],Dupr⋆[m]) in the strong resolvent sense on L2(μ[m]).
Thus, the right closed form in (12.6) converges to (E[m],Dupr⋆[m])
in the strong resolvent sense in the two-step limits.
Hence, any limit point of
(ENn,[m],DNn,[m]) lies between (E[m],Dlwr⋆[m]) and
(E[m],Dupr⋆[m]) in the sense that, for any such nonnegative f,g∈D∘b,
the quantity in the right-hand side of (12.5) is between
∫SfR⋆,∞α,[m]gdμ[m]
and
∫SfR⋆,∞α,[m]gdμ[m].
Here R⋆,∞α,[m] and R⋆,∞α,[m] are the
α-resolvent of (E[m],Dlwr⋆[m]) and (E[m],Dupr⋆[m]) , respectively.
By Theorem 11.1, these two Dirichlet forms coincide. We write their common form as
(E[m],D⋆[m]),
and denote its α-resolvent on L2(μ[m]) by R⋆α,[m].
Therefore, (ENn,[m],DNn,[m]) converges to
(E[m],D⋆[m]) in the following sense.
For any nonnegative f,g∈D∘b,
[TABLE]
Hence Rα,[m]=R⋆α,[m].
This implies that X[m] is a Markov process with resolvents Rα,[m], α>0, and that X[m]=X[m] in law.
Hence
[TABLE]
This, together with tightness in W(Rd)N, completes the proof.
∎
Remark 12.2*.*
The main difficulty in proving Theorem 2.5 is that the state space of the N-particle system depends on N.
In [16], using a general theorem on the convergence of Dirichlet forms in distinct state spaces [19], a result similar to Theorem 2.5 was proved.
In that work, the N-particle dynamics is considered in domains SRN with reflecting boundary conditions such that RN→∞.
However, this does not yield convergence of the particle systems in the whole space and is therefore unsatisfactory.
The key idea in the proof of Theorem 2.5 is the identification of a relation connecting the N-particle Dirichlet form with finite-domain Dirichlet forms.
13. Appendices
13.1. Topology of the unlabeled and labeled path spaces
Let Ss,i={s∈S;s(Rd)=∞,s({s})∈{0,1} for all s∈Rd}.
Let W(Ss,i) be the set of all Ss,i-valued continuous paths defined on [0,∞).
Remark 13.1*.*
The preprint [30] is a longer version of the published paper
[31], containing additional topological and labeling results
that were omitted from the journal version.
For any label l, there exists a unique map
lpath:WNE(Ss,i)→W(Rd)N such that
l(w0)=lpath(w)0 and upath∘lpath(w)=w.
Set ∥u∥T=sup0≤t≤T∣u(t)∣.
For u=(ui)i∈N and v=(vi)i∈N, define
[TABLE]
Then (W(Rd)N,ϱlpath) is a complete separable metric space.
Let fR∈C0(Rd) satisfy 0≤fR≤1, fR=1 on SR, and
fR+1=0 on SR+1c. For s=∑iδsi, set
fRs=∑ifR(si)δsi.
Let ϱPrh denote the Prohorov metric on S
[11, p. 11].
The vague topology on S is induced by the complete separable metric
[TABLE]
Then C([0,∞);S) is a complete separable metric space under the metric
13.2. Closability and extensions of Dirichlet forms,
and strongly local, quasi-regular Dirichlet Forms
In this subsection, we recall notions related to Dirichlet forms
from [5, 8].
Let S be a Polish space, that is, a topological space homeomorphic to a separable and complete metric space. Let ν be a Radon measure on S.
We set L2(ν)=L2(S,ν).
For non-negative symmetric bilinear forms (E1,D1) and (E2,D2),
we say (E2,D2) is an extension of (E1,D1) if
[TABLE]
A non-negative, closed form (E,D) on L2(ν) is a densely defined bilinear form defined on D that is complete under the inner product E(f,g)+(f,g)L2(ν).
i
A non-negative symmetric bilinear form (E,D0) densely defined on L2(ν) is called closable on L2(ν) if, for any E-Cauchy sequence fn∈D0 such that lim∥fn∥L2(ν)=0 as n→∞,
it holds that limn→∞E(fn,fn)=0.
ii
If (E,D0) is closable on L2(ν), then there exists
a closed extension of (E,D0).
The smallest closed extension (E,D) is called the closure of (E,D0).
Lemma 13.4**.**
Let (E1,D1) and (E2,D2) be non-negative symmetric bilinear forms on L2(ν).
Let (E2,D2) be an extension of (E1,D1). Let (E2,D2) be closable on L2(ν). It holds that (E1,D1) is closable on L2(ν).
Proof.
If {fn} is a Cauchy sequence of (E1,D1), then {fn} is a Cauchy sequence of (E2,D2) by (13.1).
Then limn→∞E2(fn,fn)=0 because of the closability of (E2,D2) on L2(ν).
Hence, limn→∞E1(fn,fn)=0 from (13.1). This completes the proof.
∎
A sequence of closed, nonnegative bilinear forms
(En,Dn) on L2(ν)
is said to converge to (E,D) on L2(ν)
in the strong resolvent sense
if the associated α-resolvents Gαn
converge strongly in L2(ν) to the α-resolvent Gα
of (E,D) for each α>0.
This is equivalent to the same convergence holding for some α>0
(cf. [37]).
In this case, the L2(ν)-semigroups Ttn
associated with (En,Dn)
converge strongly in L2(ν) to the L2(ν)-semigroup Tt
associated with (E,D)
for each t≥0.
For non-negative symmetric bilinear forms (E1,D1) and (E2,D2), we write (E1,D1)≤(E2,D2) if
D1⊃D2 and
E1(f,f)≤E2(f,f) for all f∈D2.
We write (E1,D1)≥(E2,D2) if
(E2,D2)≤(E1,D1).
A sequence of symmetric non-negative bilinear forms {(En,Dn)} on L2(ν) is said to be increasing if (En,Dn)≤(En+1,Dn+1) for all n and decreasing if (En,Dn)≥(En+1,Dn+1) for all n.
i*
Let {(En,Dn)}n∈N be an increasing sequence of non-negative symmetric closed bilinear forms on L2(ν).
Let (E∞,D∞) be the symmetric bilinear form defined by*
[TABLE]
It holds that (E∞,D∞) is a closed symmetric bilinear form on L2(ν).
{(En,Dn)}n∈N converges to (E∞,D∞) in the strong resolvent sense on L2(ν).
ii*
Let {(En,Dn)}n∈N be a decreasing sequence of non-negative symmetric closed bilinear forms on L2(ν).
Let (E∞,D∞) be the symmetric bilinear form defined by*
[TABLE]
It holds that
{(En,Dn)}n∈N converges to (Emax∞,Dmax∞) in the strong resolvent sense on L2(ν).
Here, (Emax∞,Dmax∞) is the maximal closable part of (Emax∞,Dmax∞) on L2(ν).
Let O be the family of all open subsets of S.
For A∈O, let
LA,1={f∈D;f≥1ν-a.e.on A}.
We set O0={A∈O;LA,1=∅}.
Let E1=E+(⋅,∗)L2(ν).
For an open set A∈O, we set
[TABLE]
For any set B⊂S, we set
Cap(B)=inf{Cap(A);A∈O,A⊃B}.
Let
\mathscr{D}(F_{k})=\{f\in\mathscr{D};f=0\text{ \nu-a.e.\,on }F_{k}^{c}\} and
E1=E+(⋅,⋅)L2(ν).
An increasing sequence of closed sets {Fk} is called an E-nest if ∪k≥1D(Fk) is E1-dense.
A function f is called E-quasi-continuous if for any ϵ>0, there exists an open set O with Cap(O)<ϵ such that f∣S\O is finite and continuous.
We call a subset N⊂S an E-polar set if
there exists an E-nest {Fk} such that
N⊂∩k(S\Fk).
Definition 13.4*.*
A Dirichlet form (E,D) on L2(ν) is quasi-regular if:
(Q1) there exists an E-nest {Fk,k≥1} consisting of compact sets;
(Q2)
there exists an E1-dense subset of D
whose elements admit E-quasi-continuous ν-versions;
(Q3)
there exists {fk,k≥1}⊂D having E-quasi-continuous ν-versions {f~k,k≥1}⊂D and an E-polar set N⊂S such that
{f~k,k≥1}⊂D separates the points on S\N.
A Dirichlet form (E,D) on L2(ν) is called strongly local if E(f,g)=0 for any f,g∈D such that f is constant on a neighborhood of the support of g [5].
A diffusion process is a family of Markov processes with continuous sample path and has the strong Markov property. We say a diffusion is conservative if it has an infinite lifetime.
If the Dirichlet form is strongly local and quasi-regular and the state space is homeomorphic to a complete separable metric space, then there exists a properly associate diffusion exists [5].
We refer to [5] for the concept of properly associated.
In general, a Dirichlet form is not necessarily symmetric.
In the present paper, a Dirichlet form means a symmetric Dirichlet form.
13.3. Representation of functions on S and correlation functions
Definition 13.5*.*
A symmetric, locally integrable function ρνm on (Rd)m
is called the m-point correlation function
of an RPF ν (with respect to Lebesgue measure) if
[TABLE]
for any disjoint bounded measurable sets
A1,…,An
and nonnegative integers
k1+⋯+kn=m.
Here, s(A) denotes the number of particles in A
when s=∑iδsi,
and the fraction is interpreted as zero
if s(Ai)−ki<0.
Let a={aq}q∈N be a family of sequences
aq={aq(R)}R∈N of natural numbers.
Assume that, for all q,R∈N,
[TABLE]
where aq+(R)=1+aq(R+1).
For aq={aq(R)}R∈N, let
[TABLE]
Then K[aq] is a compact set in S for each q∈N such that
[TABLE]
Let ν be an RPF and m∈N. Let νq=ν(⋅∩K[aq]).
Suppose that ν has density functions.
Then by (13.4), νq has the m-point correlation function and the m-reduced Campbell measure νq[m] similarly as (2.35).
We define ν[m] as the increasing limit of νq[m]:
[TABLE]
Note that ν[m] is independent of the particular choice of compact sets K[aq].
Let ∥h∥:=∥h∥L1(SQm×S,ν[m])+∥∇xh∥L1(SQm×S,ν[m]).
We note that, in the definition of ∥⋅∥, the differential
is taken only with respect to the x-variable,
and not with respect to s∈S.
This point is very different from the energy norms
used in Theorems 13.1 and 13.2.
Lemma 13.6**.**
The set D∘b[m] is dense in D∙b[m] with respect to ∥⋅∥, Q∈N.
Proof.
Recall that D∘b[m]=C0∞((Rd)m)⊗D∘b and D∙b[m]=C0∞((Rd)m)⊗D∙b.
From the proof of Lemma 2.4 in [24], we easily see that both D∘b and D∙b are dense in L2(ν). Hence, we can prove Lemma 13.6 from this.
∎
For R∈N and m∈{0}∪N, set
SRm={s∈S;s(SR)=m}.
For A∈B(S), define πA(s)=s(⋅∩A).
Let SRm denote the m-fold product of SR.
Definition 13.6*.*
i
For s∈S, a coordinate of s on SR is a tuple
xR(s)=(xRi(s))i∈⊔m=1∞SRm such that πSR(s)=∑iδxRi(s).
ii The restriction xRm(s) of xR(s) on SRm is called the SRm-coordinate.
iii
Let f:S→R and R,m∈N.
We say that fR,sm is the SRm-representation of f if
fR,sm(x) is a function on SRm such that
[TABLE]
and the following conditions hold:
[TABLE]
We define fR,s on ⊔m∈NSRm by fR,s(x)=fR,sm(x) for x∈SRm.
Note that fR,sm is uniquely determined, and for every s∈S,
[TABLE]
By convention, xR(s)=∅ for s∈SR0, and fR,s is constant on SR0.
For a bounded set A, one defines xAm(s), fAm(x) and fA,s(x) analogously by replacing SR with A.
13.4. Quasi-regularity of m-labeled Dirichlet forms
Let the coefficient a satisfy (2.31) and Da[m] be as (7.1).
Let
[TABLE]
Although E[m], DR∙[m], and D∘[m] depend on a and μ,
we suppress this dependence in the notation.
We make assumptions:
A.1
(E[m],D∘[m]) is closable on L2(ν[m]).
A.2
∫Ss(SR)ν(ds)<∞ for each R∈N.
A.3
ν has a bounded k-density function on SRk for each R,k∈N.
Let (E[m],D∘[m]) be the closure of (E[m],D∘[m]) on L2(ν[m]).
We write (Eν,D∘ν) for m=0 since ν[0]=ν.
Assume A.1 for m=0, A.2, and A.3.
Then the Dirichlet form (Eν,D∘ν) on L2(ν) is quasi-regular.
Proof.
The theorem was proved in [24, Th. 1] under the assumption that the diffusion matrix a is the identity.
In the present setting, a is uniformly elliptic and bounded by (2.31),
and hence the same argument applies verbatim.
∎
Lemma 13.7**.**
Assume A.2. Then D∘ν is dense in L2(ν).
Proof.
The assertion was proved in the proof of Lemma 2.5 in [17].
See the last line of the proof of Lemma 2.5 in [17].
∎
In Theorem 13.2, we significantly relax the assumptions introduced in Theorem 13.1.
Specifically, if (E[m],D∘[m]) is found to be closable on L2(ν[m]), then the existence of a density function is no longer necessary.
Theorem 13.2**.**
Assume A.1 and A.2.
Then (E[m],D∘[m]) is a quasi-regular Dirichlet form on L2(ν[m]).
Proof.
Assumptions (Q1) and (Q3) follow from A.1 and A.2
by exactly the same arguments as those used in the proofs of
[24, Th. 1] and [26, Lem. 2.3].
Thus, it only remains to prove (Q2).
For m=0, A.3 is used only in the proof of
[24, Lem. 2.4] to establish Lemma 13.7.
Hence, (Q2) holds for m=0.
For m∈N, the result follows from [26, Lem. 2.3]
by reducing to the case m=0,
and we therefore omit the proof.
∎
13.5. Weak solution and the IFC condition to ISDEs
In Subsection 13.5, we quote the IFC condition from [34].
Let Ss,i={s∈S;s(Rd)=∞,s({s})∈{0,1} for all s∈Rd}.
Let Ssde be a Borel set such that Ssde⊂Ss,i.
Let u[1,0]:Rd×S→S such that u[1,0](x,s)=δx+s.
Let
Ssde=u−1(Ssde) and Ssde[1]=(u[1,0])−1(Ssde).
Let
σ:Ssde[1]→Rd2 and b:Ssde[1]→Rd be Borel measurable functions.
We consider the ISDE of X=(Xi)i∈N with state space Ssde such that
[TABLE]
Here, Xti⋄=∑j=iδXtj and
B=(Bi)i∈N is an (Rd)N-Brownian motion.
Let Lp be the set of all measurable {Ft}t≥0-adapted
processes α such that
E[∫0T∣α(t,ω)∣pdt]<∞ for all T.
Definition 13.7* (weak solution).*
By a weak solution of ISDE (13.7), we mean an (Rd)N×(Rd)N-valued stochastic process (X,B) defined on a probability space (Ω,F,P)
with a reference family {Ft}t≥0 such that
i X=(Xi)i=1∞ is an {Ft}t≥0-adapted, Ssde-valued continuous process.
ii B=(Bi)i=1∞ is an (Rd)N-valued
{Ft}-Brownian motion with B0=0,
iii
the family of measurable {Ft}t≥0-adapted
processes Φ and Ψ defined by
[TABLE]
belong to L2 and L1, respectively.
Here, we can and do take a predictable version of Φi and Ψi (see pp 45-46 in [10]).
iv with probability one, the process (X,B) satisfies for all t
[TABLE]
Let X=(Xi)i∈N be a weak solution to (13.7) starting at s=l(s). Here, l is the label given in (2.22).
Let X such that Xt=∑iδXti.
Define σXm:[0,∞)×(Rd)m→Rd2 and
bXm:[0,∞)×(Rd)m→Rd such that, for (u,v)∈(Rd)m and
v=∑i=1m−1δvi∈S, where
v=(vi)i=1m−1∈(Rd)m−1,
[TABLE]
Here, Xtm∗=∑i=m+1∞δXti.
For s=(si)i∈N, we set sm∗=∑i=m+1∞δsi.
Recall that X0=s=l(s).
We have X0m∗=sm∗ by construction.
The coefficients σXm and bXm depend on
both Xm∗ and the label l.
In particular, X gives a part of σXm and bXm.
Let
[TABLE]
where wtm∗=∑i=m+1∞δwti for
wt=∑i=1∞δwti.
By definition, Ssdem(t,w) is a time-dependent domain in (Rd)m given by wtm∗.
We consider the SDE to Ym=(Ym,i)i=1m with random environment X defined on (Ω,F,Ps,{Ft}) such that
[TABLE]
Here, Ym,i⋄=(Ym,j)j=im and
Ytm,i⋄=∑j=imδYtm,j.
For X=(Xi)i∈N, let Xm∗=(0,…,0,Xm+1,Xm+2,...).
The first m components of Xm∗ are constant paths [math]. A triplet (Ym,Bm,Xm∗) of continuous processes on (Ω,F,Ps,{Ft}) satisfying (13.8) is called a weak solution.
Let
W=W(Rd) and W0m={w∈Wm;w(0)=0}.
Let Btm=σ[w(s);0≤s≤t,w∈Wm].
Let Bt(W0m×WN)=σ[(v(s),w(s));0≤s≤t].
We set Ctm=Bt(W0m×WN), and Cm=B(W0m×WN), where ⋅ denotes the completion with respect to
Psm=Ps∘(Bm,Xm∗)−1.
Definition 13.8*.*
i
Ym is called a strong solution of (13.8) for (X,B) under Ps if
(Ym,Bm,Xm∗) satisfies (13.8) and
there exists a Cm-measurable function Fsm:W0m×WN→Wm
such that Fsm is Ctm/Btm-measurable for each t, and Fsm satisfies Ym=Fsm(Bm,Xm∗)Ps-a.s..
ii
The SDE (13.8) is said to have a unique strong solution for (X,B) under Ps if there exists a function Fsm satisfying the conditions in Definition 13.8 i and,
for any weak solution (Y^m,Bm,Xm∗) of (13.8) under Ps,
Y^m=Fsm(Bm,Xm∗) for Ps-a.s.
By construction, the function Fsm is unique for Psm-a.s.
We introduce the IFC condition of (X,B) defined on (Ω,F,P,{Ft}):
($$\mathbf{IFC}$$) The SDE (13.8) has a unique strong solution Fsm(Bm,Xm∗)
for (X,B) under Ps for P∘X0−1-a.s. s
for all m∈N, where Ps=P(⋅∣X0=s).
13.6. A unique strong solution of ISDEs
Definition 13.9* (uniqueness in law).*
Uniqueness in law of weak solutions for (13.7) with initial distribution ν holds if
the laws of the processes X and X′ in WN coincide for any
weak solutions X and X′ with initial distribution ν.
Definition 13.10* (pathwise uniqueness).*
Pathwise uniqueness of weak solutions of
(13.7) holds if whenever X and X′ are two weak solutions
defined on the same probability space (Ω,F,P) with the same reference family {Ft}t≥0 and the same (Rd)N-valued {Ft}-Brownian motion B such that X0=X0′ a.s.,
then
[TABLE]
Let PBr∞ be the distribution of an (Rd)N-valued Brownian motion
B with B0=0.
Let Bt(PBr∞) and B(PBr∞) be the completion of
σ[ws;0≤s≤t,w∈W0N] and B(W0N)
with respect to PBr∞, respectively.
Let Bt∞=σ[w(s);0≤s≤t,w∈WN].
Definition 13.11* (a strong solution starting at s).*
A weak solution X of (13.7)
with an (Rd)N-valued Ft-Brownian motion
B is called a strong solution starting at s defined on (Ω,F,P,{Ft})
if X0=s a.s. and if there exists a function
Fs:W0N→WN such that
Fs is B(PBr∞)/B(WN)-measurable, and that
Fs is Bt(PBr∞)/Bt∞-measurable for each t, and
that Fs satisfies
[TABLE]
Definition 13.12* (a unique strong solution starting at s).*
The ISDE (13.7) admits a unique strong solution starting at s
if there exists a B(PBr∞)/B(WN)-measurable mapping
Fs:W0N→WN such that the following hold:
i
For any weak solution (X^,B^) of (13.7)
starting at s, it holds that
[TABLE]
ii
For any (Rd)N-valued {Ft}-Brownian motion B
defined on (Ω,F,P,{Ft}) with B0=0,
the continuous process Fs(B)
is a strong solution of (13.7) starting at s.
Definition 13.13* (a unique strong solution under constraints).*
Let ($$\bullet$$) be a condition.
The ISDE (13.7) admits a unique strong solution starting at s
under the constraints ($$\bullet$$) if there exists a
B(PBr∞)/B(WN)-measurable mapping
Fs:W0N→WN such that the following hold:
i
For any weak solution (X^,B^) of (13.7)
starting at s and satisfying ($$\bullet$$),
it holds that
[TABLE]
ii
For any (Rd)N-valued {Ft}-Brownian motion B
defined on (Ω,F,P,{Ft}) with B0=0,
the continuous process Fs(B)
is a strong solution of (13.7) starting at s
and satisfying ($$\bullet$$).
Assume that ν is tail trivial.
Assume that (13.7) has a weak solution (X,B) under P satisfying (SIN), (NBJ), (AC)ν, and (IFC).
It holds that (13.7) has a family of unique strong solutions {Fs} starting at s for P∘X0−1-a.s. s under the constraints (SIN), (NBJ), (AC)ν, (IFC), and (MF).
*Under the same assumptions as Theorem 13.3, the following hold:
i
The uniqueness in law of weak solutions of (13.7) with initial distribution ν holds under the constraints (SIN), (NBJ), (AC)ν, and (IFC).
ii
The pathwise uniqueness of weak solutions of (13.7) with initial distribution ν holds under the constraints (SIN), (NBJ), (AC)ν, and (IFC).*
13.7. Cut-off compact sets and a sufficient condition for (IFC)
In Subsection 13.7, we present a sufficient condition for the IFC condition.
We localize the coefficients of (13.8) to prove the IFC condition.
To do this, we introduce a sequence of compact subsets H[a]k in (Rd)m×S.
Let Ss,i be the set consisting of infinite configurations with no multiple points.
Let x=(x1,…,xm)∈(Rd)m, u(x)=∑i=1mδxi,
and s=∑iδsi. We set
[TABLE]
We set Srm={x∈Rd;∣x∣≤r}m. Let j,k,l=1,…,m. Let
[TABLE]
Then Sp,rm(s) is an open set and Sp,rm(s) is its closure in (Rd)m.
Let a={aq}q∈N and K[aq+] be as in (13.3)–(13.5).
For k=(p,q,r), let
[TABLE]
Although H[a]k=H[a]p,q,r depends on m∈N, we suppress the
dependence on m in the notation. Clearly, the set H[a]p,q,r
is increasing in each parameter p,q,r∈N. Hence, by monotonicity,
the limit of the following exit times ςp,q,r as
p,q,r→∞, taken successively in p, q, and r,
is well defined.
Let Ssde be the set defined before (13.7).
By definition, Ssde is the set such that the coefficient of (13.7) is defined on
Ssde[1]=(u[1,0])−1(Ssde).
Here we extend the domain of u to (Rd)m×S such that u[m,0](x,s)=u(x)+s.
Let (X,B) be a weak solution of (13.7) defined on (Ω,F,P,{Ft}).
Let X[m] be the associated m-labeled process.
Let ν be the RPF in (AC)ν.
{B1} There exists a sequence a={aq}q∈N
of sequences such that the following conditions hold.
i
ν(S\⋃qK[aq+])=0.
ii
The inclusion u[m,0](H[a])⊂Ssde holds.
iii
X[m]=(Xm,Xm∗) does not exit from H[a]=⋃(p,q,r)∈N3H[a]p,q,r:
[TABLE]
Here ςp,q,r(w[m])=inf{t>0;wt[m]∈H[a]p,q,r} is the exit time from H[a]p,q,r.
For k=(p,q,r)∈N3 and (x,s),(y,s)∈Sp,rm(s), we set
(x,s)∼k(y,s)
if x and y are in the same connected component of
Sp,rm(s) and s∈Π2(H[a]k).
Here, Π2 is a projection Π2:(Rd)m×S→S given by Π2(x,s)=s.
With these preparation, we then make the assumptions:
{B2} For each k=(p,q,r)∈N3, there exists a constant Fk such that
[TABLE]
for all (x,s),(y,s)∈H[a]k such that
(x,s)∼k(y,s), where H[a]k is the set defined in (13.11) and ∼k is the equivalence relation defined after (13.12).
Theorem 13.5**.**
Suppose that (13.7) has a weak solution (X,B) satisfying {B1} and that {B2} holds.
Then (X,B) satisfies (IFC).
Proof.
Theorem 13.5 is essentially proved in [34, Prop. 11.1].
Since the formulation there is slightly different, we explain the
correspondence between the assumptions.
(A1)–(A4) in [34] are basic assumptions on
the dynamics, which have already been verified in the present paper in a
form adapted to the present setting.
(B1) in [34] corresponds to
ν(S\⋃qK[aq+])=0, which follows from {B1} i.
(B2) in [34] follows from {B1} and {B2}.
In [34, Prop. 11.1], it is further assumed that the
coefficient σm of the ISDE is constant. This assumption can,
however, be removed. Indeed, it suffices to use
[34, Lem. 11.2]1 in place of
[34, Lem. 11.2]3 in the proof, which can be verified directly.
∎
The next step is to provide a sufficient condition for {B1}.
{UB}a=(akl(x,s))k,l=1d is uniformly elliptic with upper bound c13.1:
[TABLE]
Let ν be an RPF such that ν(Ssde)=1.
Let ν[m] be the m-reduced Campbell measure of ν for m∈N and ν[0]=ν.
Let l:Ss,i→(Rd)N be a label.
Let {Qs} be a family of probability measures on (Ω,F,{Ft}) such that (X,B) defined on (Ω,F,Qs,{Ft})
is a weak solution of (13.7) starting at s=l(s) for ν-a.s. s.
We assume {Qs} is a measurable family in the following sense.
{MF}Qs(A) is
B(S)ν-measurable in s for each A∈F.
Assume {MF} and let Qν=∫SQsdν.
Then (X,B) under Qν is a solution of (13.7) with the initial distribution ν∘l−1.
{BX}σ[Bs;s≤t]⊂σ[Xs;s≤t] for all t under Qν.
Let Qx,s[m] be the distribution of X[m]=(Xm,Xm∗) under
Qu(x)+s. We define
Qν[m]=∫(Rd)m×SQx,s[m]dν[m] for m∈N and Qν[0]=Qν∘X−1.
{Sν}
For each m∈{0}∪N, X[m] under (Ω,F,{Qu(x)+s},{Ft})
gives a symmetric, Markovian semi-group Tt[m] on L2(ν[m]) defined by
[TABLE]
Furthermore, ν[m] is an invariant measure of Tt[m].
We label s=∑iδsi in such a way that ∣si∣≤∣si+1∣ for all i.
Let aq={aq(r)}r∈N be the increasing sequences in (13.3).
For Q∈N∪{∞}, define
[TABLE]
Let θ∈C∞(R) such that 0≤θ(t)≤1 for all t∈R,
θ(t)=0 for t≤ϵ, and
θ(t)=1 for t≥1−ϵ for a sufficiently small ϵ>0.
Assume that 0≤θ′(t)≤2 for all t∈R.
Let
[TABLE]
The limit in (13.15) exists since θ∘dqQ≥0 is non-decreasing in Q.
Assume {UB},{MF},{BX},{Sν}, and (13.16).
Let X=upath(X) be the unlabeled dynamics under Qν.
Then
[TABLE]
Here κq=inf{t>0;Xt∈/K[aq]}.
Proof.
The proof follows from [18, Prop. 6.2], where an additional assumption {D} was imposed.
This assumption was not used in [18, Prop. 6.3]. (It was used in [18, Prop. 6.2]).
∎
Bibliography46
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] Armstrong, S., Serfaty, S., Local laws and rigidity for Coulomb gases at any temperature, Ann. Probab. 49 , 46–121 (2022)
2[2] Assiotis, T., Mirsajjadi, Z.S., ISDE with logarithmic interaction and characteristic polynomials, preprint
3[3] Bouleau, N., Hirsch, F., Dirichlet Forms and Analysis on Wiener Space , Walter de Gruyter, Berlin (1991)
4[4] Bufetov, A.I., Kawamoto, Y., The Intertwining Property for Laguerre Processes with a Fixed Parameter. J. Stat. Phys. 192 , 58 (2025)
5[5] Chen, Z.-Q., Fukushima, M., Symmetric Markov Processes, Time Change, and Boundary Theory , Princeton University Press, Princeton, NJ (2012)
6[6] Esaki, S., Tanemura, H., Stochastic differential equations for infinite particle systems of jump type with long-range interactions, J. Math. Soc. Japan 76 , 283–336 (2024)
7[7] Fritz, J., Gradient dynamics of infinite point systems, Ann. Probab. 15 , 478–514 (1987)
8[8] Fukushima, M., Oshima, Y., Takeda, M., Dirichlet Forms and Symmetric Markov Processes , 2nd ed., Walter de Gruyter, Berlin (2011)