We study equations from the area of peridynamics, which is a nonlocal extension of elasticity. The governing equations form a system of nonlocal wave equations. We take a novel approach by applying operator theory methods in a systematic way. On the unbounded domain ℝ^{n}, we present three main results. As main result 1, we find that the governing operator is a bounded function of the governing operator of classical elasticity. As main result 2, a consequence of main result 1, we prove that the peridynamic solutions strongly converge to the classical solutions by utilizing, for the first time, strong resolvent convergence. In addition, main result 1 allows us to incorporate local boundary conditions, in particular, into peridynamics. This avenue of research is developed in companion papers, providing a remedy for boundary effects. As main result 3, employing spherical Bessel functions, we give a new practical series representation of the solution which allows straightforward numerical treatment with symbolic computation.

## I. INTRODUCTION AND MOTIVATION

Classical elasticity has been successful in characterizing and measuring the resistance of materials to crack growth. On the other hand, peridynamics (PD), a nonlocal extension of continuum mechanics developed by Silling,^{55} is capable of quantitatively predicting the dynamics of propagating cracks, including bifurcation. Its effectiveness has been established in sophisticated applications such as Kalthoff-Winkler experiments of the fracture of a steel plate with notches,^{31,58} fracture and failure of composites, nanofiber networks, and polycrystal fracture.^{34,44,60,59} Further applications are in the context of multiscale modeling, where PD has been shown to be an upscaling of molecular dynamics^{52,54} and has been demonstrated as a viable multiscale material model for length scales ranging from molecular dynamics to classical elasticity.^{10} Also see other related engineering applications,^{13,33,35,46,45} the review, and news articles^{14,16,21,57} for a comprehensive discussion, and the book.^{38}

We study a class of nonlocal wave equations. The driving application is PD whose equation of motion corresponds exactly to the nonlocal wave equation under consideration. The same operator is also employed in nonlocal diffusion.^{9,14,51} Similar classes of operators are used in numerous applications such as population models, particle systems, phase transition, and coagulation. Nonlocal operators have also been used in image processing.^{27,36} In addition, we witness a strong interest for PD, its applications, and related nonlocal problems addressing, for instance, conditioning analysis, domain decomposition, and variational theory,^{5–7} discretization,^{1,7,25,62} numerical methods,^{15,17,50} nonlinear PD,^{26,37} nonlinearity in nonlocal wave equations,^{19,20} well-posedness in various forms,^{5,6,8,18,22–24,39,40,67} and other aspects.^{12,28,41,53}

It is part of the folklore in physics that the point particle model, which is the root for *locality* in physics, is the cause of unphysical singular behavior in the description of the underlying phenomena. This fact is a strong indication that, in the long run, the development of nonlocal theories is necessary for description of natural phenomena. Similar set of operator theory tools used for local problems can be transferred to study nonlocal problems because operator theory does not discern the locality or nonlocality of the governing operator. This article adds valuable tools to the arsenal of methods to analyze nonlocal problems.

The equation of interest falls into the class of abstract evolution equations, more precisely, abstract linear wave equations. Methods from operator theory are ideal to treat such equations. One can directly gain access to powerful tools such as functional calculus for bounded self-adjoint operators, spectral theorems (of densely defined, linear, and self-adjoint operators in Hilbert spaces), and strong resolvent convergence of operators. Well-posedness of the initial value problem, conservation of energy, stability of solutions (merely determined by the spectrum) are all immediately available; see Theorem 1 and Corollary 2. Furthermore, representation of solutions can easily be constructed through functional calculus for bounded operators; see our main result 3.

There are many studies^{6–8,18,24,39,40,56,67} that aim to establish a connection between the nonlocal operator and the classical one. We denote this nonlocal-to-local connection by *NL*2*L*. PD formulation utilizes a distance parameter *δ*, called the horizon, in order to limit interactions to nearby points. The source of nonlocality is the horizon and it is embedded in the support of the micromodulus function. Studies turn to the limit of vanishing nonlocality, i.e., *δ* → 0 when they aim to establish *NL*2*L*. The exact expression of *NL*2*L* is lost due to taking limit. We take general micromodulus functions in *L*^{1}(ℝ^{n}), not necessarily parametrized by *δ*. Without resorting to a *δ* → 0 argument, we identify the exact expression of the *NL*2*L* connection, what we call as main result 1.

Identifying the *NL*2*L* connection leads to a fruitful research direction because it helped us to extend the construction for unbounded domains to the bounded domain case by using the same function of the classical operator used in PD. That way, we kept a close proximity to PD in the bounded domain case. More importantly, *NL*2*L* connection turned out to be the notable result that *the governing operator is a bounded function of the classical (local) operator*.^{69} This has far reaching consequences. It enables the incorporation of local boundary conditions into nonlocal theories; see the companion papers.^{2,3} It takes a lot of effort to mitigate the boundary effects. See the comprehensive discussion [Ref. 38, Chaps. 4, 5, 7, and 12] and the numerical study.^{32} Incorporation of local boundary conditions provides a remedy for boundary effects seen in PD.

The rest of the article is structured as follows. We provide mathematical introduction in Section I. In Section II, we set the operator theory framework to treat the nonlocal wave equation. We prove basic properties of the solutions such as well-posedness of the initial value problem and conservation of energy. We provide a representation of the solutions in terms of bounded functions of the governing operator. We study the stability of solutions. We also consider the class of inhomogeneous wave equations and prove well-posedness of the corresponding initial value problem as well as a representation of the solutions in terms of bounded functions of the governing operator.

In Section III, in the vector-valued case, we note that the governing operator becomes an operator matrix. The generality of operator theory allows a simple extension of the results established for the scalar-valued functions to the vector-valued ones. We prove the boundedness of the entries of the governing operator matrix. The proof is natural due to operator theory again, because it relies on a well-known criterion for integral operators. We present a “diagonalization” of the matrix entries. This is accomplished by employing the unitary Fourier transform and connecting the entries to maximal multiplication operators.

In Section IV, we collect the statements of the three main results:

**Main Result 1:**Nonlocal operator is a bounded function of the classical operator.**Main Result 2:**Strong convergence of nonlocal solutions to classical ones through strong resolvent convergence.**Main Result 3:**Representation of the solution in terms of spherical Bessel functions.

In Section V, we construct two examples to demonstrate the usefulness of main result 3. We give explicit representation of the solutions in terms of Bessel functions. Since the explicit representation is available, we easily compute the solutions by symbolic computation, depict the resulting solutions of PD wave equation, and compare to the classical solutions.

In Section VI, we collect the proofs of the three main results together with related supporting material. In particular, in Section VI A, utilizing Fourier transforms, we turn the nonlocal governing operator into maximal multiplication operator. This process can be viewed as a form of “diagonalization.” The spectra of maximal multiplication operators are well understood. In addition, the functional calculus associated to maximal multiplication operators is known. Spectra and functional calculus allow the construction of the functional calculi for the governing operator, which in turn is used to prove that the governing operators with spherically symmetric micromoduli are functions of the Laplace operator. In Section VI B, we prove strong convergence of nonlocal solutions to classical ones through the notion of strong resolvent convergence of functions of the classical operator. Different types of resolvent convergences, i.e., norm, strong, and weak, are available; see Ref. 65 for comparison. We utilize the notion of strong resolvent convergence in order to obtain strong convergence of solution operators. This allows us to prove main result 3. We give examples of sequences of micromoduli that are instances of this result. In Section VI C, we consider the calculation of the solution of the wave equation. Since the governing operator is bounded, holomorphic functions of that operator can be represented in the form of power series in the operator. Then, we give a representation of holomorphic functions, present in the solution of the initial value problem, utilizing the fact that the governing operator is a sum of two commuting operators. We establish main result 3 by discovering that the corresponding power series can be given in terms of a series of Bessel functions. In addition, we provide an error estimate for the series representation. We use this estimate in the plots of solutions of the PD wave equation; see Figures 1 and 2. In Section VI D, we give the explicit representation of the solutions of inhomogeneous wave equations and prove well-posedness of the corresponding initial value problem. We conclude in Section VII.

Next, we present the formal system of linear PD wave equations in *n*-space dimensions [Ref. 55, Eq. (54)], *n* ∈ ℕ^{∗},

Here, “⋅” indicates matrix multiplication, or equivalently by the system

where *x* ∈ ℝ^{n}, *t* ∈ ℝ, *C* : ℝ^{n} → *M*(*n* × *n*, ℝ) is the micromodulus tensor, assumed to be even and assuming values inside the subspace of symmetric matrices, *ρ* > 0 is the mass density, *b* : ℝ^{n} × ℝ → ℝ^{n} is the prescribed body force density, and *u* : ℝ^{n} × ℝ → ℝ^{n} is the displacement field.

For comparison, e.g., the corresponding wave equation in classical elasticity in 1-space dimension is given by

where *E* > 0 is the so called “Young’s modulus,” and describing compression waves in a rod.

If *j*, *k* ∈ {1, …, *n*} and *C _{jk}* ∈

*L*

^{1}(ℝ

^{n}), we can rewrite (1.1) as

## II. OPERATOR-THEORETIC TREATMENT OF SYSTEMS OF WAVE EQUATIONS

Analogous to the majority of evolution equations from classical and quantum physics, (1.3) can be treated with methods from operator theory, see, e.g., Refs. 11 and 47 for substantiation of this claim and Ref. 29 for applications of operator theory in engineering. More specifically, this system falls into the class of abstract linear wave equations from Theorem 1. For the proof of this theorem see, e.g., Ref. 11 [Theorem 2.2.1 and Corollary 2.2.2]. Special cases of this theorem are proved in Refs. 30 and 42 and Ref. 48 [Vol. II]. Statements and proofs make use of the spectral theorems of (densely defined, linear and) self-adjoint operators in Hilbert spaces, including the concept of functions of such operators, see, e.g., Ref. 48 [Vol. I], or standard books on functional analysis, such as Refs. 49 and 66. These methods are also used throughout the paper.

This section provides the basic properties of the solutions of abstract wave equations of the form (2.1). Theorem 1 gives the well-posedness of the initial value problem for a class of abstract wave equations, conservation of energy, and a representation of the solutions in terms of bounded functions of the governing operator. Theorem 3 provides special solutions of the associated class of inhomogeneous wave equations. Together with Theorem 1, these solutions provide the well-posedness of the initial value problem of the latter equations as well as a representation of the solutions in terms of bounded functions of the governing operator.

*Let* $X,\u3008|\u3009$ *be some non-trivial complex Hilbert space. Furthermore, let* *A* : *D*(*A*) → *X* *be some densely defined, linear, semibounded self-adjoint operator in* *X* *with spectrum* *σ*(*A*)*. Finally, let* *ξ*, *η* ∈ *D*(*A*)*.*

*Then there is a unique twice continuously differentiable map**u*: ℝ →*X**assuming values in**D*(*A*)*and satisfying*$u\u2033(t)=\u2212A\u2009u(t)$*for all**t*∈ ℝ*as well as*$u(0)=\xi ,\u2009\u2009u\u2032(0)=\eta .$*For this**u, the corresponding energy function**E*: ℝ → ℝ_{u}*, defined by*$Eu(t)\u225412\u2009\u2009\u3008u\u2032(t)|u\u2032(t)\u3009+\u3008u(t)|Au(t)\u3009\u2009$*for all**t*∈ ℝ*, is constant.**Moreover, this**u**is represented by the following solution operators for all**t*∈ ℝ*:*^{70}$u(t)=costA\u2009\xi +sintAA\u2009\eta .$

See the proofs of Ref. 11 [Theorem 2.2.1 and Corollary 2.2.2]. □

Moreover, if *A* is positive, the solutions of (2.1) are stable, i.e., there are no solutions that are growing exponentially in the norm.

*If*

*A*

*is positive, then for every*

*t*∈ ℝ

*, we have*

The proof is obvious. □

Duhamel’s principle leads to a solution of (2.1) for vanishing data, the proof of the well-posedness and a representation of the solutions of the initial value problem of the inhomogeneous equation,

For simplicity, the corresponding subsequent Theorem 3 assumes that *A* is in addition positive, which is the most relevant case for applications because otherwise there are exponentially growing solutions, indicating that the system is unstable. The same statement is true if *σ*(*A*) is only bounded from below. On the other hand, Theorem 3 can also be obtained by application of the corresponding well-known more general theorem for strongly continuous semigroups; see, e.g., Ref. 11 [Theorem 4.6.2]. We give a direct proof of Theorem 3 in Section VI D, which does not rely on methods from the theory of strongly continuous semigroups. For the definition of weak integration; see, e.g., Ref. 11 [Sec. 3.2].

*Let*$X,\u3008|\u3009$

*,*

*A*:

*D*(

*A*) →

*X,*

*σ*(

*A*)

*be as in Theorem 1 and, in addition,*

*A*

*be positive. Let*

*b*: ℝ →

*X*

*be a continuous map, assuming values in*

*D*(

*A*

^{2})

*such that*

*Ab,*

*A*

^{2}

*b*

*are continuous. Then,*

*v*: ℝ →

*X, for every*

*t*∈ ℝ

*defined by*

*where*∫

*denotes weak integration in*

*X,*

*is twice continuously differentiable, assumes values in*

*D*(

*A*)

*, is such that*

See Section VI D. □

## III. THE GOVERNING OPERATOR AND PROPERTIES

The standard data space for the classical wave equation is a *L*^{2}-space with constant weight, on a non-empty open subset of ℝ^{n}, *n* ∈ ℕ^{∗}, for instance, $LC2(R)$ for a bar of infinite extension in 1-space dimension. It turns out that the classical data spaces are suitable also as data spaces for peridynamics, for instance, again $LC2(R)$ for a bar of infinite extension in 1-space dimension, composed of a “linear peridynamic material.” This simplifies the discussion of the convergence of peridynamic solutions to classical solutions.

In the following, we represent (1.3) in the form of (2.1), where the governing operator *A* is an “operator matrix,” consisting of sums of multiples of the identity and convolution operators, as indicated in (1.3). These matrix entries will turn out to be pairwise commuting. The following remark provides some known relevant information on operator matrices of bounded operators. On the other hand, we avoid explicit matrix notation.

*If* 𝕂 ∈ {ℝ, ℂ}*,* (*X*, 〈|〉) *a non-trivial* 𝕂*-Hilbert space,* (*A _{jk}*)

_{j,k∈{1,…,n}}

*a family of elements of*

*L*(

*X*,

*X*)

*.*

*Then by*$A(\xi 1,\u2026,\xi n)\u2254\u2009\u2211k=1nA1k\xi k,\u2026,\u2211k=1nAnk\xi k$*for every*(*ξ*_{1}, …,*ξ*) ∈_{n}*X*^{n}, there is defined a bounded linear operator with adjoint*A*^{∗}*given by*$A\u2217(\xi 1,\u2026,\xi n)=\u2009\u2211k=1nAk1\u2217\xi k,\u2026,\u2211k=1nAkn\u2217\xi k$*for every*(*ξ*_{1}, …,*ξ*) ∈_{n}*X*^{n}.*If the members of*(*A*)_{jk}_{j,k∈{1,…,n}}*are pairwise commuting, then**A**is bijective if and only if*det(*A*)*is bijective, where*$det(A)\u2254\u2211\sigma \u2208Snsign(\sigma )\u2009A1\sigma (1)\cdots An\sigma (n),$*S*_{n}*denotes the set of permutations of*{1, …,*n*}*,*$sign(\sigma )\u2254\u220fi,j=1,i<jnsign(\sigma (j)\u2212\sigma (i))$*for all**σ*∈*S*_{n}*and sign denotes the signum function.*

The basic properties of the entries of the operator matrix are given in the following lemma. In fact, these operators turn out to be bounded linear operators on $LC2(Rn)$. Hence, the boundedness and self-adjointness of *A* follows from those of *A _{C}*. The boundedness of

*A*has appeared in various forms

^{5,6,8,18,23,24,67}and sometimes for special class of kernel functions. We give a result using kernel functions that are in

*L*

^{1}(ℝ

^{n}) by utilizing a well-known criterion for integral operators; see, e.g., corollary to Ref. 64 [Theorem 6.24].

*Let*

*n*∈ ℕ

^{∗}

*,*

*ρ*> 0

*and*

*C*∈

*L*

^{1}(ℝ

^{n})

*be even. Then,*

*for every*$f\u2208LC2(Rn)$

*, where*∗

*denotes the convolution product, there is defined a self-adjoint bounded linear operator on*$LC2(Rn)$

*with operator norm*$AC$

*satisfying*

By a consequence of a well-known criterion for integral operators on *L*^{2}-spaces, see, e.g., Corollary to Ref. 64 [Theorem 6.24]. □

For other governing operators related to *A _{C}* on bounded domains, see Ref. 4.

## IV. MAIN RESULTS

The following are the three main results in this paper:

**Main Result 1:**Nonlocal operator is a bounded function of the classical operator.**Main Result 2:**Strong convergence of nonlocal solutions to classical ones through strong resolvent convergence.**Main Result 3:**Representation of the solution in terms of spherical Bessel functions.

Our first main result establishes the connection between the nonlocal operator and the classical one. We identify the exact expression of this connection. Namely, the governing operator *A _{C}* in (3.1) of the peridynamic wave equation is a bounded function of the classical governing operator in (1.2). This observation enables the comparison of peridynamic solutions to those of classical elasticity. In addition, it enables the generalization of peridynamic-type operators from ℝ

^{n}to bounded domains as functions of the corresponding classical operator. This is the subject of our companion papers.

^{2,3}

*Let*

*n*∈ ℕ

^{∗}

*,*$Ln$

*be the closure of the positive symmetric, essentially self-adjoint operator in*$LC2(Rn)$

*, given by*

*where*

*ρ*> 0

*and*

*E*> 0

*. In addition, for*

*n*> 1

*, let*

*C*

*be spherically symmetric, i.e., such that*

*for every*

*R*∈

*SO*(

*n*)

*, where*

*SO*(

*n*)

*denotes the map of group of special orthogonal transformations on*ℝ

^{n}

*. Then*

*where*

*ι*: [0, ∞) → ℝ

^{n}

*is defined by*

*for every*

*s*⩾ 0

*and*

*e*

_{1}, …,

*e*

_{n}*denotes the canonical basis of*ℝ

^{n}

*.*

Employing the notion of strong resolvent convergence, our second main result gives the strong convergence of solutions of the governing equation to that of the classical equation. There are results that establish the convergence of the peridynamic operator to the classical operator with vanishing nonlocality; see Refs. 6, 7, 24, and 56. More important for applications is the corresponding convergence of solutions. This question has also been pointed out in Ref. 24 [p. 862] for the convergence of PD solutions to that of Navier equation. The studies^{8,18,67} considered convergence of solutions. In particular, Refs. 39 and 40 put an effort to establish strong convergence of solutions. One can find a survey of existing convergence results in Ref. 21 [Chap. 4]. The tool that we develop for this purpose is the notion of strong resolvent convergence^{64,65} used in Theorem 7. To the best of our knowledge, our study provides the first exploitation of strong resolvent convergence of the functions of the classical operator to obtain strong convergence of solutions.

*Let*(

*X*, 〈|〉)

*be a non-trivial complex Hilbert space and*

*A*:

*D*(

*A*) →

*X*

*a densely defined, linear, and self-adjoint operator with spectrum*

*σ*(

*A*)

*. Furthermore, let*

*f*

_{1},

*f*

_{2}, …

*be a sequence of real-valued functions in*$UCs(\sigma (A))$

^{71}

*that is everywhere on*

*σ*(

*A*)

*pointwise convergent to*

*id*

_{σ(A)}

*, and for which there is*

*M*> 0

*such that*

*for all*

*ν*∈ ℝ

*. Then for every*

*g*∈

*BC*(ℝ, ℂ)

*, where*

*BC*(ℝ

^{n}, ℂ)

*is the space of complex-valued bounded continuous functions on*ℝ

^{n}

*,*

*where for every*

*ν*∈ ℕ

^{∗}

*,*

*σ*(

*f*

_{ν}(

*A*))

*denotes the spectrum of*

*f*

_{ν}(

*A*)

*.*

*The mechanism to obtain strong convergence of nonlocal solutions to classical ones is as follows. We first choose a sequence of micromoduli*

*C*

_{ν}

*so that*

*is pointwise convergent to*$id\sigma (L1)$

*and that satisfies condition (6.1), where*$L1$

*is the classical governing operator in 1 dimension defined in Theorem 6. We choose*$A=L1$

*and*

*f*

_{Cν}

*in (4.1). Then, a well-known result [Ref.*

*64*

*, Theorem 9.16] indicates that Lemma 15 implies strong resolvent convergence of*$fC\nu (L1)$

*to*$L1$

*. Then, Theorem 7 implies that for every*

*g*∈

*BC*(ℝ, ℂ)

*, we obtain the following strong convergence:*

*We denote the nonlocal operator corresponding to micromoduli*

*C*

_{ν}

*by*$AC\nu \u2254fC\nu (L1)$

*. Then, we apply (4.2) to solution operators*

*g*

*in (2.2), namely,*$cos(tAC\nu \u2009)$

*and*$sin(tAC\nu )/AC\nu $

*, both of which are in*

*BC*(ℝ, ℂ)

*. Consequently, we obtain strong convergence of nonlocal solution operators*$[g|\sigma (fC\nu (L1)](fC\nu (L1))$

*in (3.1) to the classical solution operators*$[g|\sigma (L1)](L1).$

*As an application, we establish strong convergence of solutions generated by standard examples of micromoduli* *C*_{ν} *in the literature, given in (6.3) and (6.4); see Lemmata 16 and 18.*

Our third main result is the discovery that the solution can be represented in terms of a series of Bessel functions. This property can be exploited in obtaining numerical representation of the solutions. We utilized this for the solutions depicted in Section V.

*Let*

*n*∈ ℕ

^{∗}

*,*

*ρ*> 0

*,*

*C*∈

*L*

^{1}(ℝ

^{n})

*be even such that*

*Then, for*

*A*

_{C}*in (3.1) and*

*t*∈ ℝ

*, the solutions have the following representation:*

*for every*$f\u2208LC2(Rn)$

*, where the spherical Bessel functions*

*j*

_{0},

*j*

_{1}, …

*are defined as in Ref.*

*43*

*,*

*and the members of the sums are defined for*

*t*= 0

*by continuous extension.*

## V. EXAMPLES AND ILLUSTRATIONS OF MAIN RESULT 3

As an application of main result 3, we construct examples of the nonlocal wave equation to illustrate the separation of waves and propagation of discontinuities, Examples 10 and 11, respectively. The depicted solutions in Figures 1 and 2 were computed symbolically using spherical Bessel function representation of the solutions. The infinite series in (5.1) and (5.2), respectively, are truncated after adding 46 terms. Visually, we do not observe a difference in the solutions if more terms are added.

For constructing the case that characterizes separation of waves, we choose an example that involves a micromodulus and input function both of which are normal distributions with mean value zero and standard deviation *σ* and *σ _{d}*, respectively. Similar micromodulus functions have also been used in Refs. 41, 63, and 68.

*For*

*ρ*,

*σ*,

*σ*,

_{d}*a*> 0

*, we define*

*C*

_{σ}∈

*L*

^{1}(ℝ)

*and*$f\u2208LC2(R)$

*by*

*Then for*

*k*∈ ℕ

^{∗}

*and hence*

*Since*

*we conclude from Theorem 9 that for*

*A*

_{C}*in (3.1) and*

*t*∈ ℝ

*,*

We depict and compare the solutions of the classical and nonlocal wave equations in Figures 1 and 2. In the classical case, as expected, we observe the propagation of waves along characteristics; see Figures 1(a) and 1(c) for vanishing initial velocity and displacement, respectively. In the nonlocal case, we observe repeated separation of waves and an oscillation at the center of the initial pulse; see Figures 1(b) and 1(d) for vanishing initial velocity and displacement, respectively. Solution wave patterns were also reported earlier^{23,25} for the same micromodulus function in Example 11.

Now, we study the propagation of discontinuity in the data for classical and nonlocal wave equations in the following example.

*As in the previous example, for*

*ρ*,

*σ*,

*a*,

*b*, ε > 0

*, we define*

*C*

_{σ}∈

*L*

^{1}(ℝ)

*by*

*Then*

*and for*

*k*∈ ℕ

^{∗}

*,*

*and hence*

*Furthermore, we define*$f\u2208LC2(R)$

*by*

*Then for*

*x*∈ ℝ

*,*

*where*erfc

*denotes the error function defined according to DLMF.*

^{43}

*We note for*

*x*∈ ℝ

*that*

*Since*

*we conclude from Theorem 9 that for*

*t*∈ ℝ

*,*

*for*

*A*

_{C}*in (3.1), where for every*

*x*∈ ℝ

*and*

*k*∈ ℕ

^{∗}

*,*

In the classical wave equation, as expected, discontinuities propagate along the characteristics; see Figures 2(a) and 2(c) for vanishing initial velocity and displacement, respectively. On the other hand, in the nonlocal case, the discontinuity remains in the same place for all time; see Figures 2(b) and 2(d) for vanishing initial velocity and displacement, respectively. This confirms the results given in Ref. 63.

## VI. PROOFS AND RELATED RESULTS

### A. Proof of main result 1 (nonlocal operator is a bounded function of the classical operator)

For the study of the spectral properties of the matrix entries, needed for the application of the results from Section II, we use Fourier transformations. This step parallels the common procedure for constant coefficient differential operators on ℝ^{n}, *n* ∈ ℕ^{∗}. With the help of the unitary Fourier transform *F*_{2}, Theorem 13 represents the matrix entries as maximal multiplication operators. This process can be viewed as a form of “diagonalization” of the entries. Also, since bounded maximal multiplication operators commute, the entries commute pairwise. The spectra of maximal multiplication operators are well understood, leading to Corollary 14. Also, the functional calculus which is associated to maximal multiplication operators is known and allows the construction of the functional calculi of the entries. The latter is used in the proof of Theorem 6 which proves that matrix entries corresponding to spherically symmetric micromoduli are functions of the Laplace operator.

*In the following, for*

*n*∈ ℕ

^{∗}

*,*

*F*

_{2}

*denotes the unitary Fourier transformation on*$LC2(Rn)$

*which, for every rapidly decreasing test function*

*f*∈ 𝒮

_{ℂ}(ℝ)

*, is defined by*

*Also, we denote by*

*F*

_{1}

*the map from*$LC1(Rn)$

*to*

*C*

_{∞}(ℝ

^{n}, ℂ)

*, the space of continuous functions on*ℝ

^{n}

*vanishing at infinity, which for every*$f\u2208LC1(Rn)$

*, is defined by*

*Let*

*denote the maximal multiplication operator by the bounded continuous function*

*on*$LC2(Rn)$

*. Then*

*K*≔

*C*∘(

*p*

_{1}−

*p*

_{2}), with projections

*p*

_{1},

*p*

_{2}: ℝ

^{2n}→ ℝ

^{n}, where

*x*

_{1}, …,

*x*,

_{n}*y*

_{1}, …,

*y*) ∈ ℝ

_{n}^{2n}. Since

*C*is even,

*K*is symmetric. Furthermore, for every

*x*∈ ℝ

^{n}and

*y*∈ ℝ

^{n},

*K*induces a self-adjoint bounded linear integral operator Int(

*K*) on $LC2(Rn)$ defined by

*F*

_{2}∘Int(

*K*) and

*T*

_{F1C}∘

*F*

_{2}coincide on a dense subspace of $LC2(Rn)$, they coincide on the whole of $LC2(Rn)$. Consequently, the result follows from

We finally arrive at the proof of main result 1.

*F*

_{1}

*C*)(0) −

*F*

_{1}

*C*∈

*BC*(ℝ

^{n}, ℝ) and that (

*F*

_{1}

*C*)(0) −

*F*

_{1}

*C*is even, since for every

*k*∈ ℝ

^{n}

*n*> 1, we note that

*R*∈

*SO*(

*n*) and

*k*∈ ℝ

^{n}and hence that

*k*∈ ℝ

^{n}. In particular,

We give the spectrum and point spectrum of *A _{C}*. The spectrum is an essential ingredient in constructing functional calculus for

*A*; see Lemma 20.

_{C}*A*

_{C}*where the overline denotes the closure in*ℝ

*. Finally, for every*

*λ*∈

*σ*(

*A*)

_{C}*,*

*A*−

_{C}*λ*

*is not surjective.*

The result is a consequence of well-known properties of multiplication operators. □

### B. Proof of main result 2 (strong convergence of nonlocal solutions to classical ones through strong resolvent convergence)

Lemma 15 gives conditions for the convergence of bounded functions of a self-adjoint operator to converge to that operator, which implies strong resolvent convergence and also the strong convergence of the same bounded continuous function of each member of the sequence against that bounded continuous function of the self-adjoint operator [Ref. 64, Theorem 9.16]; see Theorem 7.

*Let*(

*X*, 〈|〉)

*be a non-trivial complex Hilbert space and*

*A*:

*D*(

*A*) →

*X*

*a densely defined, linear, and self-adjoint operator with spectrum*

*σ*(

*A*)

*. Furthermore, let*

*f*

_{1},

*f*

_{2}, …

*be a sequence in*$UCs(\sigma (A))$

*that is everywhere on*

*σ*(

*A*)

*pointwise convergent to*

*id*

_{σ(A)}

*, and for which there is*

*M*> 0

*such that*

*for all*

*ν*∈ ℝ

*. Then*

*ξ*∈

*D*(

*A*) and

*ψ*

_{ξ}the corresponding spectral measure. According to the spectral theorem for densely defined, self-adjoint linear operators in Hilbert spaces, $idR2$ is

*ψ*

_{ξ}-summable and

*μ*,

*ν*∈ ℕ

^{∗}. As a consequence of the pointwise convergence of

*f*

_{1},

*f*

_{2}, … on

*σ*(

*A*) to id

_{σ(A)}, (6.1) and Lebesgue’s dominated convergence theorem, it follows that

*f*

_{1}(

*A*)

*ξ*,

*f*

_{2}(

*A*)

*ξ*, … is a Cauchy sequence in

*X*. Since $(X,\u2009\u2009)$ is in particular complete, the latter implies that

*f*

_{1}(

*A*)

*ξ*,

*f*

_{2}(

*A*)

*ξ*, … is convergent in $(X,\u2009\u2009)$. Furthermore,

*f*

_{1},

*f*

_{2}, … on

*σ*(

*A*) to id

_{σ(A)}, (6.1), Lebesgue’s dominated convergence theorem and the spectral theorem for densely defined, self-adjoint linear operators in Hilbert spaces has been applied. From the polarization identity for 〈|〉, it follows that

*ξ*,

*η*∈

*D*(

*A*). Since

*D*(

*A*) is dense in

*X*, the latter implies that

*ξ*∈

*X*,

*η*∈

*D*(

*A*) and hence for every

*η*∈

*D*(

*A*) that

Now, main result 2, Theorem 7, is a consequence of Lemma 15, and, for example, Ref. 48 [Vol. I, Theorems 8.20 and 8.25].

We apply our main result 2 to standard examples treated in the literature. Hence, we easily obtain strong resolvent convergence for these cases, i.e., Lemmas 16 and 18. They provide sequences of micromoduli which satisfy the conditions of Lemma 15. Lemma 16 has also been treated in Refs. 41 and 61 and Lemma 18 has been treated in Ref. 41. Remark 8 applies Theorem 7 to the sequences of micromoduli from Lemmas 16 and 18. As a consequence, for fixed data and *t* ∈ ℝ, the solutions of the initial value problem at time *t* corresponding to the members of each sequence of micromoduli converge in $LC2(R)$ to the corresponding classical solution at time *t*.

*For every*

*ν*∈ ℕ

^{∗}

*, we define*

*C*

_{ν}∈

*L*

^{1}(ℝ)

*by*

*Then, the result (6.2) preparing for strong convergence of solutions is satisfied*

*By applying Theorem 7 as in Remark 8, we obtain that the solutions of the nonlocal operator* *A*_{Cν} *strongly converge to that of the local operator for the micromoduli* *C*_{ν} *in (6.3).*

*ν*∈ ℕ

^{∗},

*ν*∈ ℕ

^{∗},

*λ*⩾ 0,

*k*> 0

*λ*⩾ 0 that

*For every*

*ν*∈ ℕ

^{∗}

*, we define*

*C*

_{ν}∈

*L*

^{1}(ℝ)

*by*

*Then, the result (6.2) preparing for strong convergence of solutions is satisfied,*

*By applying Theorem 7 as in Remark 8, we obtain that the solutions of the nonlocal operator* *A*_{Cν} *strongly converge to that of the local operator for the micromoduli* *C*_{ν} *in (6.4).*

*ν*∈ ℕ

^{∗},

*λ*⩾ 0,

*k*⩾ 0

*k*⩾ 0, with the help of Lebesgue’s dominated convergence theorem, that

*λ*⩾ 0 that

### C. Proof of main result 3 (representation of the solution in terms of spherical Bessel functions)

In this section, the goal is to prove that the solutions of the nonlocal wave equation given in (2.2) can be expressed in terms of spherical Bessel functions. This result has a practical implication because it allows symbolic computation of solutions. Indeed, using symbolic computation, we generate solutions of the nonlocal wave equation and illustrate how separation of waves and propagation of discontinuities occur; see Figures 1 and 2. In addition, one can use the explicit representations of solutions for benchmarking numerical results which helps in the development of numerical methods for computing approximate solutions. Such benchmarking could be used, for instance, for verifying and validating numerical solutions.

Main result 1 indicates that the governing peridynamic operator is bounded. Hence, functions of that operator in (2.2) can be represented in the form of power series in the governing operator due to functional calculus of self-adjoint and bounded operators. In Lemma 20, we prove that holomorphic functional calculus is available for power series representation. In Lemma 21, we apply this representation to the functions present in the solution of the initial value problem of the homogeneous wave equation. Lemmas 20 and 21 can be viewed as straightforward applications of the spectral theorems for densely defined, self-adjoint linear operators in Hilbert spaces.

*Let*(

*X*, 〈|〉)

*be a non-trivial complex Hilbert space,*

*A*∈

*L*(

*X*,

*X*)

*self-adjoint, and*

*σ*(

*A*) ⊂ ℝ

*the (non-empty, compact) spectrum of*

*A. Furthermore, let*$R>A$

*and*

*f*:

*U*(0) → ℂ

_{R}*be holomorphic. Then, the sequence*

*is absolutely summable in*

*L*(

*X*,

*X*)

*and*

*σ*(

*A*) that

*z*∈

*U*(0) as well as, since $\sigma (A)\u2282BA(0)\u2282UR(0),$ that the sequence of continuous functions

_{R}*f*|

_{σ(A)}. In particular, since $A<R$, this implies that the sequence

*L*(

*X*,

*X*), and it follows from the spectral theorem for bounded self-adjoint operators in Hilbert spaces that

Now, we give the exact representation of the terms present in (2.2) in terms of power series in *A*.

*Let*(

*X*, 〈|〉)

*be a non-trivial complex Hilbert space,*$ij$

*the complex square-root function, with domain*ℂ∖((−∞, 0] × {0})

*.*

*A*∈

*L*(

*X*,

*X*)

*self-adjoint and*

*σ*(

*A*) ⊂ ℝ

*the (non-empty, compact) spectrum of*

*A. For every*

*t*∈ ℝ

*, the sequences*

*are absolutely summable in*

*L*(

*X*,

*X*)

*and*

*t*∈ ℝ,

*z*∈ ℂ∖((−∞, 0] × {0}) and

*z*∈ ℂ∖[0, ∞) × {0}), respectively. Furthermore,

*z*∈ ℂ∖((−∞, 0] × {0}) and

*z*∈ ℂ∖[0, ∞) × {0}), respectively. Then, it follows from Lemma 20 that the sequences

*L*(

*X*,

*X*) and that

The goal here was to show that (2.2) can be expressed by spherical Bessel functions. In order to arrive at a Bessel function representation, we first connect (2.2) to generalized hypergeometric functions. In Theorem 22, we give a general result which involves two arbitrary bounded commuting operators *A* and *B*. In Theorem 24, we will apply this result to special commuting operators, multiple of the identity operator, and a convolution, i.e., *c* − *C*, because the governing operator is of this form. In addition, we expect that the expansion below can be used for large time asymptotic of the solutions of the nonlocal wave equation, a research direction beyond the scope of this paper.

*Let*(

*X*, 〈|〉)

*be a non-trivial complex Hilbert space,*$ij$

*the complex square-root function, with domain*ℂ∖((−∞, 0] × {0})

*.*

*A*,

*B*∈

*L*(

*X*,

*X*)

*self-adjoint such that*[

*A*,

*B*] = 0

*and*

*σ*(

*A*),

*σ*(

*A*+

*B*) ⊂ ℝ

*the (non-empty, compact) spectra of*

*A*

*and*

*A*+

*B, respectively. Then*

*where*

_{0}

*F*

_{1}

*denotes the generalized hypergeometric function, defined as in Ref.*

*43*

*.*

*t*∈ ℝ that the family

*L*(

*X*,

*X*), since for (

*k*,

*l*) ∈ ℕ

^{2}

*S*⊂ ℕ

^{2},

*L*(

*X*,

*X*), since for (

*k*,

*l*) ∈ ℕ

^{2},

*S*⊂ ℕ

^{2}. Hence, we conclude the following:

*k*,

*l*∈ ℕ,

*k*. First, we note that

*k*∈ ℕ. Then

*k*+ 1. The equality (6.5) implies for every

*k*,

*l*∈ ℕ that

_{0}

*F*

_{1}, for every

*l*∈ ℕ,

*z*∈ ℂ,

*k*,

*l*∈ ℕ,

We give a connection between generalized hypergeometric and spherical Bessel functions which will lead to main result 3.

*For every*

*k*∈ ℕ

*and*

*x*> 0

*,*

*where the spherical Bessel functions*

*j*

_{0},

*j*

_{1}, …

*are defined as in Ref.*

*43*

*and*

*ν*∈ (0, ∞),

*k*∈ ℕ, and

*x*> 0,

*k*∈ ℕ,

*x*> 0

*k*∈ ℕ

^{∗},

*x*> 0

*x*> 0

*k*= 0, if we define

Eventually, we have a representation involving two commuting operators, with one of the operators being a multiple of the identity and a general *C* which is not necessarily a convolution operator.

*Let*(

*X*, 〈|〉)

*be a non-trivial complex Hilbert space,*$ij$

*the complex square-root function, with domain*ℂ∖((−∞, 0] × {0})

*.*

*c*> 0

*,*

*C*∈

*L*(

*X*,

*X*)

*self-adjoint and*

*σ*(

*c*−

*C*) ⊂ ℝ

*the (non-empty, compact) spectrum of*

*c*−

*C. Then for every*

*t*∈ ℝ

*,*

*where the spherical Bessel functions*

*j*

_{0},

*j*

_{1}, …

*are defined as in Ref.*

*43*

*and*

*and the members of the sums are defined for*

*t*= 0

*by continuous extension.*

Direct consequence of Theorem 22 and Lemma 23. □

Now, main result 3, Theorem 9, is a consequence of Theorem 24. The representations given in Theorems 22, 24, and 9 were used to express the solutions in Examples 10 and 11.

As a side result, we provide an error estimate for the representation in Theorem 24. For the solution plots, Figures 1 and 2, the infinite series is truncated after adding 46 terms, i.e., *N* = 45. The below error estimate has been used to monitor the error. For *N* = 45 and *t* = 20, the error is in the order of 10^{−13}.

*Let*$(X,\u3008|\u3009),ij,c,C,\sigma (c\u2212C),j\u22121,j0,j1,\u2026$

*as in Theorem 24 and*

*N*∈ ℕ

*. Then for every*

*t*∈ ℝ

*,*

*t*∈ ℝ,

### D. Proof of Theorem 3 (solutions of inhomogeneous wave equations)

We give a proof of Theorem 3.

*λ*> 0 that

*t*,

*τ*∈ ℝ, where

*α*,

*β*: ℝ →

*L*(

*X*,

*X*) are defined by

*t*∈ ℝ. In the following, for

*ξ*∈

*D*(

*A*), we are going to use that the maps

*a*,

*b*are strongly continuous and that

*t*∈ ℝ. Also for every

*k*∈

*U*

_{ℂ}(

*σ*(

*A*))

^{s},

*k*(

*A*)

*D*(

*A*) ⊂

*D*(

*A*), and for

*ξ*∈

*D*(

*A*),

*a*,

*b*induce strongly continuous maps from ℝ to

*X*, which we indicate with the same symbols, and where $XA\u2254(D(A),\u2009\u2009A)$. In addition, we note that the inclusion

_{A}*ι*of

*X*into

_{A}*X*is continuous. In the next step, we observe for a strongly continuous

*c*: ℝ →

*L*(

*X*,

_{A}*X*) and a continuous

_{A}*g*: ℝ →

*X*that

_{A}*X*,

_{A}*t*↦

*c*(

*t*)

*g*(

*t*)) is continuous as well as that

^{A}denotes weak integration in

*X*, is differentiable with derivative

_{A}*t*∈ ℝ that

*c*: ℝ →

*L*(

*X*,

*X*),

*g*: ℝ →

*X*such that Ran(

*g*) ⊂

*D*(

*A*),

*t*∈ ℝ, and

*h*∈ ℝ

^{∗}that

*v*gives for

*t*∈ ℝ,

*X*, and that

## VII. CONCLUSION

Our result that the governing operator is a bounded function of the classical local operator for scalar-valued functions should be generalizable to vector-valued case. Our notable result that the governing operator *A _{C}* of the peridynamic wave equation is a bounded function of the classical governing operator has far reaching consequences. It enables the comparison of peridynamic solutions to those of classical elasticity. The remarkable implication is that it opens the possibly of defining peridynamic-type operators on bounded domains as functions of the corresponding classical operator. Since the classical operator is defined through

*local*boundary conditions, the functions inherit this knowledge. This observation opens a gateway to incorporate local boundary conditions into nonlocal theories, which has vital implications for numerical treatment of nonlocal problems. This is the subject of our companion papers.

^{2,3}

We expect that the expansions in Theorems 22 and 24 can be used for obtaining the large time asymptotic of solutions of the nonlocal wave equation. In the classical case, as expected, we observe the propagation of waves along characteristics. In the nonlocal case, we observe oscillatory recurrent wave separation. We think that this phenomenon is worth investigating. On the other hand, we observe that discontinuity remains stationary in the nonlocal case, whereas, it is well-known that discontinuities propagate along characteristics. We hold that this fundamental difference is one of the most distinguishing features of PD. In conclusion, we believe that we added valuable tools to the arsenal of methods to analyze nonlocal problems.

## Acknowledgments

Burak Aksoylu was supported in part by National Science Foundation Grant No. DMS 1016190, European Commission Marie Curie Career Integration Grant No. 293978, and Scientific and Technological Research Council of Turkey (TÜBİTAK) Grant Nos. TBAG 112T240 and MAG 112M891. Research visit of Horst R. Beyer was supported in part by No. TÜBİTAK 2221 Fellowship for Visiting Scientist Program. Sabbatical visit of Fatih Celiker was supported in part by No. TÜBİTAK 2221 Fellowship for Scientist on Sabbatical Leave Program. Fatih Celiker was supported in part by National Science Foundation Grant No. DMS 1115280.

## REFERENCES

Throughout the paper, we refer to the Laplace operator as the local or the classical operator interchangeably.

Solution operators, $cos(tA\u2009)$ and $sin(tA)/A,$ denote operators corresponding to unique extensions to entire holomorphic functions. A pedantic notation would be $cost\xaf\u2009|\sigma (A)(A)$ and $\u2009sint\xaf\u2009|\sigma (A)(A)$, respectively. For brevity, we prefer to use the former.

$UCs(\sigma (A),C)$ denotes the space of bounded complex-valued functions on *σ*(*A*) that are strongly measurable in the sense that they are everywhere on *σ*(*A*) limit of a sequence of step functions.