Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Matrix representations of multidimensional integral and ergodic operators

Matrix representations of multidimensional integral and ergodic operators We provide a representation of the C -algebra generated by multidimensional integral op- erators with piecewise constant kernels and discrete ergodic operators. This representation allows us to nd the spectrum and to construct the explicit functional calculus on this algebra. The method can be useful in various applications, since many discrete approxima- tions of integral and di erential operators belong to this algebra. Some examples are also presented: 1) we construct an explicit functional calculus for extended Fredholm integral operators with piecewise constant kernels, 2) we nd a wave function and spectral estimates for 3D discrete Schr odinger equation with planar, guided, local potential defects, and point sources. The accuracy of approximation of continuous multi-kernel integral operators by the operators with piecewise constant kernels is also discussed. Keywords: integral equations, functional calculus, Schr odnger operator with defects, operator algebras 1. Introduction N 2 2 M Let R be some domain. Consider the Hilbert space L := L ( ! C ; ) of vector-valued functions acting on ( is the Lebesgue measure). The goal of the paper is the study of a discrete analogue of the following algebra of operators acting on L : Z Z A = Alg A; dx ; :::; dx ;T : (1) c 1 N In other words, A is generated by: 1) multiplication operators A (the dot  denotes the M 2 place of the operator argument u = (u ) 2 L ), where A = A(x ; :::; x ) are bounded m 1 N m=1 measurable M  M matrix-valued functions with complex entries de ned on ; 2) integral R R operators dx = dx , where I is the one-dimensional set (line) de ned by n n \I I = I (x ; :::; x ; x ; :::; x ) = f(x ; :::; x ; y ; x ; :::; x ) : y 2 Rg; (2) 1 n1 n+1 N 1 n1 n n+1 N n Preprint submitted to Elsevier May 22, 2020 arXiv:1805.09100v2 [math-ph] 21 May 2020 R so, the operator u(x ; :::; x )dx integrates u in one coordinate x producing the function 1 N n n that is constant along this coordinate; 3) operators T u(x) = (u (T (x))) ; x 2 ; where T : (3) m m m m=1 are measurable mappings. In particular, the algebra A contains various integral operators of Fredholm type, ergodic operators (which are based on T ), and their combinations. All such operators have di erent applications in mathematical physics, e.g., they describe a propagation of waves and other phenomena in complex structures with defects [1, 2, 3, 4, 5, 6], di usion [7, 8], thermodynamic processes [9, 10], random Schr odinger operators and various operators on discrete graphs [11, 12, 13], electromagnetic scattering [14]. Some general aspects of the connection between integral and ergodic operators are discussed in [15, 16, 17]. Integral operators and some of their nite-dimensional approximations are discussed in [18, 19, 20]. There are also useful purely algebraic approaches for various integro-di erential modules over C (R) developed in, e.g., [21], [22], [23]. The main problems for operators from A are to nd the spectrum, to nd the inverse operators, square roots, or, more generally, to construct the functional calculus on this algebra. The diculty is that A is very complex. We try to nd some discrete analogue A of A for which the functional calculus can be constructed explicitly. One of the most important requirements to A is to be nite dimensional. Because in this case A can be expressed in terms of matrix algebras for which the functional calculus is well known. If A is nite dimensional then, due to the Stone-Weierstrass theorem, all the matrix-valued functions A should be piecewise constant, otherwise the subalgebra generated by A has an in nite dimension. This tells us how the operator algebra A should be arranged. It is natural to suppose that is a union of a nite number of shifted copies of a cube H = [0; h) (h > 0) = a + H; (4) i i i i=1 N 2 2 where a 2 R are some vertices such that are disjoint. Recall that L := L ( i i M 2 M C ; ) = L ( ) is the Hilbert space of vector-valued functions acting on , square- MM integrable with respect to the Lebesgue measure . For any A 2 C ,  f1; :::; Ng and 2 2 i; j 2 f1; :::; Sg introduce the following elementary operators E [A] : L ! L de ned by ij j j h A u(x a + a )dx ; x 2 i j i E [A]u(x) = i (5) ij 0; x 2 where dx = dx and the number of elements in is denoted by j j. Note that can n2 be the empty set ;, in this case there is no in (5). Also note that E [A]u(x) is constant ij along x in for n 2 . In fact, the operator E [A] translates the values of u(x) from the n i ij cube to the cube , then, it takes the average in along x , n 2 and multiplies the j i i n average by A, and, nally, it puts zero values inside other cubes , r 6= i. The examples of 2 (a + h, b + h) u u u j j j u u u i i i (a, b) { 1} { 2} E [1] hE [1] hE [1] ij ij ij 0 0 0 R R a+h b+h u dx u dx j 1 j 2 a b Figure 1: Schematic view of three di erent operators E in 2D case. While u(x) and Eu(x) are de ned for all x 2 = [ , the action of E is shown for x 2 , for simplicity, because Eu(x) = 0, x 62 . Here, r i j i u means uj , r = i; j. action of di erent operators E is demonstrated in Fig. 1. The operators E [A] provide the interaction between the components of u in the var- ij ious micro-domains . Roughly speaking, if h > 0 is suciently small then can be approximately represented by (4) and the discrete analogue of A can be chosen as A = Alg(E [A]): (6) ij MM In other words, A is generated by E [A] for all A 2 C ,  f1; :::; Ng and i; j 2 ij f1; :::; Sg. Let us discuss why A is a discrete analogue of A de ned by (1). The discrete analogues of the multiplication operators A are the multiplication operators with piecewise constant functions A = A(x), i.e. A(x) = A = const; x 2 ; i = 1; :::; S: (7) i i The corresponding operators are expressed in terms of E as follows: A = E [A ]: (8) ii i=1 The integral operators, see (1) and (2), are also expressed in terms of E : XX fng dx = hE [I]; (9) ij i=1 j2 3 where I is the identity matrix, = fj : a (m) = a (m); 8m 6= ng (10) i j i and a(m) is the m-th entry of the vector a. The discrete analogues of change-of-variables operators T : (see (3)) should be based on discrete mappings p : f1; :::; Sg ! m m f1; :::; Sg, since = [ . Taking arbitrary p we construct i m i=1 T x = x a + a ; x 2 ; i = 1; :::; S: (11) m i i p (i) Then the discrete analogues of T (3) are expressed in terms of E as follows: M S X X ; M MM T = E [I ]; where I = (  ) 2 C (12) m m im jm i;j=1 i;p (i) m=1 i=1 and  is the Kronecker delta. The above arguments show that A can be considered as the discrete approximation of A . As it is shown in Theorem 1.3 and Example 2 below, if h ! 0 then the approximation becomes better and better. Moreover, A is a closed C -subalgebra of A and, hence, deserves its own study. The properties of A are enough for most practical applications. Denoting by the Hermitian conjugation, one can easily check the fundamental relations E [A] = E [A ]; E [A] +E [B] = E [A + B]; E [A]E [B] =  E [AB]: (13) jk ij ji ij ij ij ij kl il Hence, E [A] are basis elements and any operator A 2 A has the form ij X X A = E [A ]; (14) ij ij i;j=1 f1;:::;Ng where A are some M  M matrices. In practice, form (14) is available after taking an ij approximation of the initial operator A 2 A . The question is how to nd explicitly the c c spectrum of A, inverse A , square root A, etc. As mentioned above, if we provide a representation of A in terms of simple matrix algebras, then the answers on all these nn questions become explicit. We denote the simple matrix algebras as C , n > 1. Introduce the following matrices S MSMS A = (A ) 2 C ; B = A (15) ij i;j=1 and the following mapping MSMS 2 : A ! (C ) ; (A) = (B ) ; (16) f1;:::;Ng where A is given by (14). The next theorem is our main result. 4 N MSMS 2 Theorem 1.1. The mapping  is the C -isomorphism between C -algebras A and (C ) . The inverse mapping has the form X X ((B ) ) = E [A ]; (17) f1;:::;Ng ij ij i;j=1 f1;:::;Ng MM S where A 2 C are the blocks of the matrix A = (A ) given by ij ij i;j=1 j n j A = (1) B : (18) Note that while the most of operators from A are in nite-dimensional and even non- compact, the algebra A has a nite dimension. We immediately obtain the following Corollary 1.2. i) The operator A is invertible if and only if all the matrices B are invert- ible. In this case, A can be computed explicitly 1 1 1 A =  ((B ) ): (19) f1;:::;Ng ii) Generalizing i) we can take rational functions f and write f (A) =  ((f (B )) ): (20) f1;:::;Ng The extension to algebraic and transcendent functions f is also obvious. iii) The trace and the determinant of A can be de ned as Y X det(A) = det B ; tr(A) = trB ; (21) f1;:::;Ng f1;:::;Ng they satisfy the usual properties det(AB) = det(A) det(B); tr(aA + bB) = atr(A) + btr(B); (22) A tr(A) det(e ) = e ; sp(A) = f : det(A I ) = 0g; (23) where I =  ((I) ) is the identity operator, I is MS  MS identity matrix, and f1;:::;Ng a; b 2 C, A;B 2 A . iv) Since  is the C -isomorphism, the operator norm of A 2 A can be computed explic- itly kAk 2 2 = maxf :  are s values of B g; (24) L !L where, recall that (A) = (B ) . f1;:::;Ng 5 Remark. The algebra Z Z 1 1 B = Alg A; dx ; :::; dx 1 n 0 0 with piecewise constant functions A acting on = [0; 1) is considered in [24]. The di er- ence between A and B is the presence of more general class of sets and the addition of change-of-variables operators T . In the case = [0; 1) with the uniform partition of each of the segments [0; 1) on p intervals (i.e. S = p ), B becomes a C -subalgebra of A : Nn n n N N N Mp Mp p Mp Mp 2 ( ) B (C )  (C ) A ; = = n=0 where are binomial coecients. It is interesting to note that the generalization of the structure of A simpli es the proof of main results. Let us discuss a norm of approximation of continuous operators from A by discrete operators from A . For simplicity, consider the case = [0; 1) and M = 1. The general- 2 2 ization to M > 1 is similar. Consider a multi-kernel operator A : L ( ) ! L ( ) of the form, common in applications, A u(k) = A (k; x )u(k + x )dx ; k 2 ; u 2 L ( ): (25) j j [0;1) f1;:::;Ng Here, we also use the notation x ; n 2 ; x = (xe ); xe = n n 0; n 62 : Note that x = x + x , where = f1; :::; Ngn is the complement to the set . Consider the uniform partition of onto p identical cubes = a + [0; 1=p) ; (26) i i i i=1 where N j1 a = (b ) ; i 1 = p b ; b 2 f0; :::; p 1g (27) i j j j j=1 j=1 is the representation of i 1 in the base p numeral system. Let us take the approximation of A by N N p p X XX A = E [A ]; (28) ij ij j j i=1 j=1 f1;:::;Ng where " " A (a + ""; (a + "") ); (a ) = (a ) ; i j i j A = and """ = (1=(2p)) : (29) ij r=1 0; otherwise; The next proposition shows us that A ! A in the operator norm. 6 Theorem 1.3. Suppose that A (k; x ) in (25) are real functions with bounded rst deriva- tives. Then 2 2 1 kA Ak 6 krA k ! 0 for p ! 1: (30) c L L !L 2p f1;:::;Ng Along with Theorem 1.1 and Corollary 1.2, we can use Theorem 1.3 to determine the spec- trum and the inverse operator. Corollary 1.4. Under the assumptions of Theorem 1.3, we assume also that A is self- adjoint, i.e. A (k; x ) = A (k + x ; k ). Then A is self-adjoint and sp(A )  B (sp(A)); sp(A)  B (sp(A )); (31) c   c where B means -neighbourhood and we can set  = krA k 1 . If 0 62 sp(A ) L c f1;:::;Ng 2p 1 1 2 2 then A is invertible and there is p such that 0 62 sp(A) and kA k 6  , and the L !L following estimate is true 1 2 kA k 2 2 1 1 L !L kA A k 2 2 6 : (32) L !L 2 2 1 kA k L !L Moreover, RHS of (32) tends to 0 for p ! 1. The rest of the paper is organized as follows. Section 2 contains two examples: 1) new formulas for the functions of 1D Fredholm integral operators with step kernels; 2) the application of the method for obtaining a solution (with arbitrary precision) and spectral estimates for 3D discrete Schr odinger equation with planar, guided, and local potential defects. A short proof of the main result based on the explicit representation of a semigroup algebra of subsets is given in Section 3. We conclude in Section 4. 2. Examples Example 1. Consider the case N; M = 1, and the classical Fredholm integral operators (see [18]) 2 2 A : L ([0; 1)) ! L ([0; 1)); Au(x) = B(x; y)u(y)dy (33) with S-step (piecewise constant) kernels i1 i 1; x 2 [ ; ); S S B(x; y) = S B  (x) (y); B 2 C;  (x) = (34) ij i j ij i 0; otherwise: i;j=1 SS Such operators form an algebra isomorphic to C (see, e.g., [19]). But, this algebra does not contain the identity operator (Iu = u). Let us supplement it by adding new operators Au(x) = A(x)u(x) + B(x; y)u(y)dy; where A(x) = A  (x): (35) i i i=1 7 In other words > > A(x) =  (x)A(x); B(x; y) = S (x)B(y); where (36) S S = ( ) ; A = diag(A ); B = (B ) : (37) i i ij i=1 i;j=1 Taking a = (i 1)=S and using notations (4)-(5), we obtain S S X X f1g A = E [A ] + E [B ]: (38) i ij ii ij i=1 i;j=1 Applying the results of Theorem 1.1 and Corollary 1.2 along with (35)-(37) we obtain that 1 > 1 > 1 1 A u(x) =  (x)A (x)u(x) + S  (x)((A + B) A )(y)u(y)dy; (39) if and only if A and A + B are invertible (otherwise A is non-invertible). The spectrum of A consists of all the eigenvalues of A and A + B. In general, > > f (A)u(x) =  (x)f (A)(x)u(x) + S  (x)(f (A + B) f (A))(y)u(y)dy (40) for various functions f (e.g. rational/algebraic/transcendent:exp; ; :::). We have reduced explicitly the functional calculus on integral operators to the functional calculus on matrices. Let us consider a concrete example with some oscillating non-symmetric integral kernel, say A u(x) = A(x)u(x) + B(x; y)u(y)dy; (41) where A(x) = x sin(x) 1; B(x; y) = 4 sin(6xy) + exp(x) y + 2: (42) Taking the step approximation of A and B, see (34)-(37), with suciently large S, and using Corollary 1.2, we can estimate the eigenvalue  closest to 0. It is   0:1 6= 0. Hence, 0 0 the operator A has an inverse. Note that we can not use a Neumann series for A , since the kernel B(x; y) is not small relative to the rst term A(x). Let us approximately compute 1 3 A and, say, A , where the real branch of cubic root is assumed. To do this, we apply (40) to the step approximations of A and B. Thus, Z Z 1 1 u(x) 1 1 3 A u(x) = + C (x; y)u(y)dy; A u(x) = (A(x)) u(x) + D(x; y)u(y)dy: (43) A(x) 0 0 Computations show a good convergence of step kernels obtained by (40). For example, the relative deviations computed for S = 60, 300, and 1500 are suciently small kC C k 1 kC C k 1 60 300 L 300 1500 L = 0:1; = 0:02; kC k 1 kC k 1 300 L 300 L 8 1 1 kD D k kD D k 60 300 L 300 1500 L = 0:07; = 0:01; 1 1 kD k kD k 300 L 300 L while kC k 1  kC k 1  kC k 1  30; kD k 1  kD k 1  kD k 1  5: 60 L 300 L 1500 L 60 L 300 L 1500 L The corresponding integral kernels are plotted in Fig. 2. Figure 2: The original integral kernel B(x; y), see (42), and the step approximations of integral kernels C (x; y), D(x; y), see (43), computed by (40) for S = 60 and S = 300. Example 2. Actually, the method based on the expansion in E works well for more nm complex mixed multidimensional integral operators. Let us consider a discrete normalized 3D time-dependent Schr odinger operator with three planar potentials V , V , V with their 1 2 3 (i+")t intersections (see Fig. 3) and with attenuated harmonic source term e ( 2 R is the frequency, " > 0 is the attenuation factor, and the amplitude is 1) located at the origin. Then the corresponding equation on the wave function is (i+")t 3 3 i  V =  e ; n = (n ) 2 Z ; (44) discr n n n n0 i i=1 @t 9 3 Figure 3: The lattice Z with planar potential defects V , line potential defects V + V , and point defect i i j V + V + V . 1 2 3 where V ; n = 0; n 6= 0; i i j j6=i V + V ; n = n = 0; n 6= 0 (k 62 fi; jg) i j i j k V = ;  = ( ) (45) n discr n n n e V + V + V ; n = 0; 1 2 3 > n e 0; otherwise and ne denotes points adjacent to n. Our goal is to solve the Schr odinger equation and to (i+")t ^ ^ nd . Assuming (t) = e with time-independent and taking Fourier series n n n n 2in k 3 3 ^ ^ (k) = e , k = (k ) 2 [0; 1) we can rewrite (44) in the integral form n i n2Z i=1 A(k) (k) Z Z Z 1 1 1 ^ ^ ^ V (x ; k ; k )dx V (k ; x ; k )dx V (k ; k ; x )dx = 1; (46) 1 1 2 3 1 2 1 2 3 2 3 1 2 3 3 0 0 0 where A(k) = i" +  4 sin k :This 3D system contains non-parallel potential i=1 defects and, hence, does not admit explicit procedure for nding (see, e.g., [25] for 2D case). There are various methods to obtain the approximated solution of (46). We will use the method based on the expansion of the operator (46) in the basis E to obtain the nm 10 ^ approximation of with an arbitrary precision. Let p 2 N, h = 1=p. Then, following (4) we have 3 3 3 [0; 1) = fa + [0; h) g; N = f0; ::; p 1g ; a = hn: (47) n n n2N It is convenient to use the indices n; m 2 N in (5) instead of numbers i; j. Introducing A = A(a + 1h=2) (where 1 = (1; 1; 1)) we can write the approximation of the equation n n ^ ee A = 1 (A denotes LHS of (46)) as A = 1, where X X X X ; f1g f2g f3g A = E [A ] E [hV ] E [hV ] E [hV ]: (48) n 1 2 3 nn nm nm nm n2N n;m2N n;m2N n;m2N n =m n =m n =m 2 2 1 1 1 1 n =m n =m n =m 3 3 3 3 2 2 3 3 p p Then the matrices A 2 C , see (14), (15), are de ned by A = diag(A ); A = hV (  ) ; A = hV (  ) ; ; n f1g 1 n m n m n;m2N f2g 2 n m n m n;m2N 2 2 3 3 1 1 3 3 A = hV (  ) ; A = A = A = A = 0: (49) f3g 3 n m n m n;m2N f1;2g f1;3g f2;3g f1;2;3g 1 1 2 2 The matrices B , see (15), are B = A ; B = A + A (i = 1; 2; 3); ; ; fig ; fig B = A + A + A ; B = A + A + A ; f1;2g ; f1g f2g f1;3g ; f1g f3g B = A + A + A ; ; B = A + A + A + A : (50) f2;3g ; f2g f3g f1;2;3g ; f1g f2g f3g To compute A we follow (19) and (17), (18). De ne the matrices C by 1 1 1 1 1 1 1 C = B ; C = B B (i = 1; 2; 3); C = B B B + B ; ; fig f1;2g ; fig ; f1;2g f1g f2g ; 1 1 1 1 1 1 1 1 C = B B B + B ; C = B B B + B ; f1;3g f2;3g f1;3g f1g f3g ; f2;3g f2g f3g ; 1 1 1 1 1 1 1 1 C = B B B B + B + B + B B : (51) f1;2;3g f1;2;3g f1;2g f1;3g f2;3g f1g f2g f3g ; e e e All the matrices B are invertible since A exists (because A + i" is self-adjoint). Then has the form Z Z 1 1 1 > e e = A 1 =  C  + pC (x ; k ; k )dx + pC (k ; x ; k )dx + ; f1g 1 2 3 1 f2g 1 2 3 2 0 0 Z Z Z 1 1 1 pC (k ; k ; x )dx + p C (x ; x ; k )dx dx + f3g 1 2 3 3 f1;2g 1 2 3 1 2 0 0 0 Z Z Z Z 1 1 1 1 2 2 p C (x ; k ; x )dx dx + p C (k ; x ; x )dx dx + f1;3g 1 2 3 1 3 f2;3g 1 2 3 2 3 0 0 0 0 Z Z Z 1 1 1 p C (x ; x ; x )dx dx dx ; (52) f1;2;3g 1 2 3 1 2 3 0 0 0 11 Q where the vector-valued function  = (  (k )) (see also (34) for  ). Due to the n i n2N n i=1 i simple form of , all the integrals in (52) can be computed explicitly: for example, Z Z 1 1 p (x ; k ; x )dx dx = ( (k )) : 1 2 3 1 3 n 2 n2N 0 0 e e So, in (52) is given explicitly. The norm of di erence between A (46) and A (48) depends on how much A approximates A, at least kA Ak 2 2 6 24h, see also Theorem 1.3. n L !L 1 1 1 1 The norm of the inverse operators kA k 2 2 6 " and kA k 2 2 6 " since A + i" L !L L !L 1 1 2 e e and A + i" are self-adjoint. Hence kA A k 2 2 6 24h=" . In particular, L !L 2 2 2 2 ^ e ^ e k k 6 24h=" ; j j 6 (24h=" ) ; (53) L n n n2Z e e where are Fourier coecients of . On the side, we obtain the following spectral estimates dist(sp(A); sp(B )) 6 24h: (54) f1;2;3g Estimates (53), (54) become better and better for h ! 0. Depending on the number j j, the components sp(B )n ([ sp(B )) approximate the volume, planar, guided, and local j j<j j isolated spectral components of A which correspond to bulk, planar, guided, and local wave functions. These waves propagate along the potential defects of the corresponding dimension 3j j and exponentially attenuate in the perpendicular directions to the defect, see, e.g., the corresponding discussion in [25]. 3. Proof of Theorems 1.1 and 1.3 At rst, let us consider the semigroup of subsets G = fe :  f1; :::; Ngg; e e = e (55) and the corresponding C -algebra M = A e : A 2 C : (56) f1;:::;Ng The identity element in this algebra is 1 = e , where ; is the empty set. All the basis elements e = e are self-adjoint. De ne the mapping : M ! C ; (57) X X A e = (B ) ; B = A : (58) 1 f1;:::;Ng f1;:::;Ng Lemma 3.1. The mapping  is the C -isomorphism. The inverse mapping is de ned by X X 1 j n j ((B ) ) = A e ; A = (1) B : (59) f1;:::;Ng f1;:::;Ng Proof. Consider the following basis in M f = e (1 e ): (60) fng n62 Direct calculations give us f = f ; f f =  f : (61) This means that ff g can be considered as an orthogonal (in the algebraic sense) f1;:::;Ng 2 N basis in C because the number of elements #f :  f1; :::; Ngg = 2 . Hence, M is isomorphic to C with the isomorphism e $ f . Using X X Y X X Y A e = A e (1 e + e ) = A e (1 e ) (62) fng fng fng n62  n62 f1;:::;Ng f1;:::;Ng f1;:::;Ng X X X X X = A f = ( A )f = B f (63) f1;:::;Ng f1;:::;Ng f1;:::;Ng we obtain that the isomorphism e $ f coincides with  (57)-(58). Similarly, identities X X Y X X j n j B f = B e (1 e ) = B (1) e (64) fng n62 f1;:::;Ng f1;:::;Ng f1;:::;Ng X X X j n j = ( (1) B )e = A e (65) f1;:::;Ng f1;:::;Ng give us the form of the inverse mapping  (59). MM By Lemma 3.1, the C -algebra M = M of MM matrices with entries belonging MM 2 to M is isomorphic to (C ) . The corresponding isomorphism is X X MM 2 : M ! (C ) ;  ( A e ) = (B ) ; B = A ; (66) M M M f1;:::;Ng f1;:::;Ng X X 1 j n j ((B ) ) = A e ; A = (1) B : (67) f1;:::;Ng f1;:::;Ng SS MSMS MM The same result holds for the C -algebra M M M . For any A 2 C = = MS and  f1; :::; Ng, the basis elements F [A] = (  e A) 2 M satisfy the same ~ ~ MS ij ii jj ~ ~ i;j=1 equations as the basis elements E [A] 2 A (see (13)) ij F [A] = F [A ]; F [A] +F [B] = F [A + B]; F [A]F [B] =  F [AB]: (68) jk ij ji ij ij ij ij kl il This means that A = M with the natural isomorphism E [A] $ F [A]. Then, com- MS ij ij paring (14)-(16) and (66)-(67) (with MS instead of M ) we deduce that  in (16) is the isomorphism satisfying (17)-(18). Now, let us prove Theoem 1.3. It is seen that X X (A A)u(k) = (A (k; x ) A (k; (xe ) ))u(k + x )dx ; c i j j [0;1) \ f1;:::;Ng i=1 where k is the center of the cube for which k 2 , and xe is the center of . Then, the j j i i mean value estimates maxjA (k; x ) A (k; (xe ) )j 6 krA k i L x2 i 2p nish the proof of (30). Estimates (31) are classical spectral estimates for self-adjoint op- erator A and self-adjoint perturbed operator A = A + A A in terms of the norm of c c c perturbation  = kAA k 2 2 , see, e.g., [26]. If A is invertible then A is invertible for c L !L c 1 1 all suciently large p and kA k 2 2 ! kA k 2 2 for p ! 1. Moreover, using the L !L L !L Neumann series for the inverse of perturbed operator, we deduce that 1 2 kA k 2 2 1 1 1 n 1 L !L 2 2 2 2 kA A k = k (A (A A)) A k 6 : L !L L !L 1 kA k 2 2 L !L n=1 The last term tends to 0 since  tends to 0 for p ! 1. 4. Conclusion We have shown that the analysis of mixed multidimensional integral and some type of ergodic operators can be explicitly reduced to the analysis of special matrices. This allows us to compute functions of such operators and their spectra explicitly with an arbitrary precision. Funding statement This paper is a contribution to the project M3 of the Collaborative Research Centre TRR 181 "Energy Transfer in Atmosphere and Ocean" funded by the Deutsche Forschungs- gemeinschaft (DFG, German Research Foundation) - Projektnummer 274762653. This work is also supported by the RFBR (RFFI) grant No. 19-01-00094. References [1] A. B. Movchan, and L. I. Slepyan, Band gap green's functions and localized oscillations, Proc. R. Soc. A 463 (2007) 2709{2727. [2] D. J. Colquitt, M. J. Nieves, I. S. Jones, A. B. Movchan, and N. V. Movchan, Waves in lattices with imperfect junctions and localised defect modes, Proc. R. Soc. A 469 (2150) (2013) 20120579. 14 [3] M. Makwana, and R. V. Craster, Localised point defect states in asymptotic models of discrete lattices, Q. J. Mechanics Appl. Math. 66 (2013) 289{316. [4] A. A. Kutsenko, Algebra of multidimensional periodic operators with defects, J. Math. Anal. Appl. 428 (2015) 221{230. [5] A. A. Kutsenko, Recovery of defects from the information at detectors, Inverse Problems 32 (2016) [6] V. Caudrelier, Multisymplectic approach to integrable defects in the sine-gordon model, J Phys A Math Gen 48 (19) (2015) 195203. [7] M. F. Norman, Ergodicity of di usion and temporal uniformity of di usion approximation, J. Appl. Prob. 14 (1977) 399{404. [8] J. Schenker, Di usion in the mean for an ergodic Schrodinger equation perturbed by a uctuating potential, Comm. Math. Phys. 339 (2015) 859{901. [9] I. Prigogine, A. P. Grecos, and Cl. George , Kinetic theory and ergodic properties, Proc. Nati. Acad. Sci. USA 14 (1976) 1802{1805. [10] E. Kreyszig, Advanced Engineering Mathematics, Wiley, 2011. [11] V. Kirsh, and L. A. Pastur , Ergodic theory and discrete one-dimensional random Schrodinger operators: uniform existence of the Lyapunov exponent , Contemp. Math. 327 (2003) 1{23. [12] D. Lenz, F. Schwarzenberger, and I. Veselic, A Banach space-valued ergodic theorem and the uniform approximation of the integrated density of states , Geom. Dedicata 150 (2011) 1{34. [13] E. Korotyaev, N. Saburova, E ective masses for laplacians on periodic graphs, J. Math. Anal. Appl. 436 (2016) 104{130. [14] F. Vico, M. Ferrando, L. Greengard, Z. Gimbutas, The decoupled potential integral equation for time- harmonic electromagnetic scattering, Comm. Pure Appl. Math. 69 (2016) 771{812. [15] T. Eisner, B. Farkas, M. Haase, and R. Nagel, Operator theoretic aspects of ergodic theory, Springer, [16] V. Kirsh, and L. A. Pastur , Analogues of Szego's theorem for ergodic operators , Mat. Sb. 206 (2015) 93{119. [17] C. Houdayer, and Y. Isono, Bi-exact groups, strongly ergodic actions and group measure space type III factors with no central sequence , Comm. Math. Phys. 348 (2016) 991{1015. [18] E. I. Fredholm, Sur une classe d'equations fonctionnelles, Acta Math. 27 (1903) 365{390. [19] M. S. Gockenbach, Finite-Dimensional Linear Algebra, Discrete mathematics and its applications., CRC Press, 2010. [20] S. Albeverio, E. I. Gordon, A. Yu. Hrennikov, Finite-dimensional approximations of operators in the hilbert spaces of functions on locally compact abelian groups, Acta Appl. Math. 64 (2000) 33{73. [21] M. Rosenkranz, A new symbolic method for solving linear two-point boundary value problems on the level of operators, Journal of Symbolic Computation 39 (2005) 171{199. [22] V. V. Bavula, The algebra of integro-di erential operators on an ane line and its modules, J. Pure Appl. Algebra 217 (2013) 495{529. [23] L. Guo, G. Regensburger, M. Rosenkranz, On integro-di erential algebras, J. Pure Appl. Algebra 218 (2014) 456{473. [24] A. Kutsenko, Mixed multidimensional integral operators with piecewise constant kernels and their representations, Linear and Multilinear Algebra 0 (2017) 0. doi:10.1080/03081087.2017.1415294. [25] A. A. Kutsenko, Algebra of 2d periodic operators with local and perpendicular defects, J. Math. Anal. Appl. 442 (2016) 796{803. [26] T. Kato, Perturbation theory for linear operators, Springer-Verlag Berlin Heidelberg, 1995. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Mathematical Physics arXiv (Cornell University)

Matrix representations of multidimensional integral and ergodic operators

Kutsenko, Anton A.
Mathematical Physics , Volume 2018 (1805) – May 23, 2018
Download PDF
15 pages

Loading next page...
 
/lp/arxiv-cornell-university/matrix-representations-of-multidimensional-integral-and-ergodic-Tjb65g2q7W

References

References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.

ISSN
0096-3003
eISSN
ARCH-3337
DOI
10.1016/j.amc.2019.124818
Publisher site
See Article on Publisher Site

Abstract

We provide a representation of the C -algebra generated by multidimensional integral op- erators with piecewise constant kernels and discrete ergodic operators. This representation allows us to nd the spectrum and to construct the explicit functional calculus on this algebra. The method can be useful in various applications, since many discrete approxima- tions of integral and di erential operators belong to this algebra. Some examples are also presented: 1) we construct an explicit functional calculus for extended Fredholm integral operators with piecewise constant kernels, 2) we nd a wave function and spectral estimates for 3D discrete Schr odinger equation with planar, guided, local potential defects, and point sources. The accuracy of approximation of continuous multi-kernel integral operators by the operators with piecewise constant kernels is also discussed. Keywords: integral equations, functional calculus, Schr odnger operator with defects, operator algebras 1. Introduction N 2 2 M Let R be some domain. Consider the Hilbert space L := L ( ! C ; ) of vector-valued functions acting on ( is the Lebesgue measure). The goal of the paper is the study of a discrete analogue of the following algebra of operators acting on L : Z Z A = Alg A; dx ; :::; dx ;T : (1) c 1 N In other words, A is generated by: 1) multiplication operators A (the dot  denotes the M 2 place of the operator argument u = (u ) 2 L ), where A = A(x ; :::; x ) are bounded m 1 N m=1 measurable M  M matrix-valued functions with complex entries de ned on ; 2) integral R R operators dx = dx , where I is the one-dimensional set (line) de ned by n n \I I = I (x ; :::; x ; x ; :::; x ) = f(x ; :::; x ; y ; x ; :::; x ) : y 2 Rg; (2) 1 n1 n+1 N 1 n1 n n+1 N n Preprint submitted to Elsevier May 22, 2020 arXiv:1805.09100v2 [math-ph] 21 May 2020 R so, the operator u(x ; :::; x )dx integrates u in one coordinate x producing the function 1 N n n that is constant along this coordinate; 3) operators T u(x) = (u (T (x))) ; x 2 ; where T : (3) m m m m=1 are measurable mappings. In particular, the algebra A contains various integral operators of Fredholm type, ergodic operators (which are based on T ), and their combinations. All such operators have di erent applications in mathematical physics, e.g., they describe a propagation of waves and other phenomena in complex structures with defects [1, 2, 3, 4, 5, 6], di usion [7, 8], thermodynamic processes [9, 10], random Schr odinger operators and various operators on discrete graphs [11, 12, 13], electromagnetic scattering [14]. Some general aspects of the connection between integral and ergodic operators are discussed in [15, 16, 17]. Integral operators and some of their nite-dimensional approximations are discussed in [18, 19, 20]. There are also useful purely algebraic approaches for various integro-di erential modules over C (R) developed in, e.g., [21], [22], [23]. The main problems for operators from A are to nd the spectrum, to nd the inverse operators, square roots, or, more generally, to construct the functional calculus on this algebra. The diculty is that A is very complex. We try to nd some discrete analogue A of A for which the functional calculus can be constructed explicitly. One of the most important requirements to A is to be nite dimensional. Because in this case A can be expressed in terms of matrix algebras for which the functional calculus is well known. If A is nite dimensional then, due to the Stone-Weierstrass theorem, all the matrix-valued functions A should be piecewise constant, otherwise the subalgebra generated by A has an in nite dimension. This tells us how the operator algebra A should be arranged. It is natural to suppose that is a union of a nite number of shifted copies of a cube H = [0; h) (h > 0) = a + H; (4) i i i i=1 N 2 2 where a 2 R are some vertices such that are disjoint. Recall that L := L ( i i M 2 M C ; ) = L ( ) is the Hilbert space of vector-valued functions acting on , square- MM integrable with respect to the Lebesgue measure . For any A 2 C ,  f1; :::; Ng and 2 2 i; j 2 f1; :::; Sg introduce the following elementary operators E [A] : L ! L de ned by ij j j h A u(x a + a )dx ; x 2 i j i E [A]u(x) = i (5) ij 0; x 2 where dx = dx and the number of elements in is denoted by j j. Note that can n2 be the empty set ;, in this case there is no in (5). Also note that E [A]u(x) is constant ij along x in for n 2 . In fact, the operator E [A] translates the values of u(x) from the n i ij cube to the cube , then, it takes the average in along x , n 2 and multiplies the j i i n average by A, and, nally, it puts zero values inside other cubes , r 6= i. The examples of 2 (a + h, b + h) u u u j j j u u u i i i (a, b) { 1} { 2} E [1] hE [1] hE [1] ij ij ij 0 0 0 R R a+h b+h u dx u dx j 1 j 2 a b Figure 1: Schematic view of three di erent operators E in 2D case. While u(x) and Eu(x) are de ned for all x 2 = [ , the action of E is shown for x 2 , for simplicity, because Eu(x) = 0, x 62 . Here, r i j i u means uj , r = i; j. action of di erent operators E is demonstrated in Fig. 1. The operators E [A] provide the interaction between the components of u in the var- ij ious micro-domains . Roughly speaking, if h > 0 is suciently small then can be approximately represented by (4) and the discrete analogue of A can be chosen as A = Alg(E [A]): (6) ij MM In other words, A is generated by E [A] for all A 2 C ,  f1; :::; Ng and i; j 2 ij f1; :::; Sg. Let us discuss why A is a discrete analogue of A de ned by (1). The discrete analogues of the multiplication operators A are the multiplication operators with piecewise constant functions A = A(x), i.e. A(x) = A = const; x 2 ; i = 1; :::; S: (7) i i The corresponding operators are expressed in terms of E as follows: A = E [A ]: (8) ii i=1 The integral operators, see (1) and (2), are also expressed in terms of E : XX fng dx = hE [I]; (9) ij i=1 j2 3 where I is the identity matrix, = fj : a (m) = a (m); 8m 6= ng (10) i j i and a(m) is the m-th entry of the vector a. The discrete analogues of change-of-variables operators T : (see (3)) should be based on discrete mappings p : f1; :::; Sg ! m m f1; :::; Sg, since = [ . Taking arbitrary p we construct i m i=1 T x = x a + a ; x 2 ; i = 1; :::; S: (11) m i i p (i) Then the discrete analogues of T (3) are expressed in terms of E as follows: M S X X ; M MM T = E [I ]; where I = (  ) 2 C (12) m m im jm i;j=1 i;p (i) m=1 i=1 and  is the Kronecker delta. The above arguments show that A can be considered as the discrete approximation of A . As it is shown in Theorem 1.3 and Example 2 below, if h ! 0 then the approximation becomes better and better. Moreover, A is a closed C -subalgebra of A and, hence, deserves its own study. The properties of A are enough for most practical applications. Denoting by the Hermitian conjugation, one can easily check the fundamental relations E [A] = E [A ]; E [A] +E [B] = E [A + B]; E [A]E [B] =  E [AB]: (13) jk ij ji ij ij ij ij kl il Hence, E [A] are basis elements and any operator A 2 A has the form ij X X A = E [A ]; (14) ij ij i;j=1 f1;:::;Ng where A are some M  M matrices. In practice, form (14) is available after taking an ij approximation of the initial operator A 2 A . The question is how to nd explicitly the c c spectrum of A, inverse A , square root A, etc. As mentioned above, if we provide a representation of A in terms of simple matrix algebras, then the answers on all these nn questions become explicit. We denote the simple matrix algebras as C , n > 1. Introduce the following matrices S MSMS A = (A ) 2 C ; B = A (15) ij i;j=1 and the following mapping MSMS 2 : A ! (C ) ; (A) = (B ) ; (16) f1;:::;Ng where A is given by (14). The next theorem is our main result. 4 N MSMS 2 Theorem 1.1. The mapping  is the C -isomorphism between C -algebras A and (C ) . The inverse mapping has the form X X ((B ) ) = E [A ]; (17) f1;:::;Ng ij ij i;j=1 f1;:::;Ng MM S where A 2 C are the blocks of the matrix A = (A ) given by ij ij i;j=1 j n j A = (1) B : (18) Note that while the most of operators from A are in nite-dimensional and even non- compact, the algebra A has a nite dimension. We immediately obtain the following Corollary 1.2. i) The operator A is invertible if and only if all the matrices B are invert- ible. In this case, A can be computed explicitly 1 1 1 A =  ((B ) ): (19) f1;:::;Ng ii) Generalizing i) we can take rational functions f and write f (A) =  ((f (B )) ): (20) f1;:::;Ng The extension to algebraic and transcendent functions f is also obvious. iii) The trace and the determinant of A can be de ned as Y X det(A) = det B ; tr(A) = trB ; (21) f1;:::;Ng f1;:::;Ng they satisfy the usual properties det(AB) = det(A) det(B); tr(aA + bB) = atr(A) + btr(B); (22) A tr(A) det(e ) = e ; sp(A) = f : det(A I ) = 0g; (23) where I =  ((I) ) is the identity operator, I is MS  MS identity matrix, and f1;:::;Ng a; b 2 C, A;B 2 A . iv) Since  is the C -isomorphism, the operator norm of A 2 A can be computed explic- itly kAk 2 2 = maxf :  are s values of B g; (24) L !L where, recall that (A) = (B ) . f1;:::;Ng 5 Remark. The algebra Z Z 1 1 B = Alg A; dx ; :::; dx 1 n 0 0 with piecewise constant functions A acting on = [0; 1) is considered in [24]. The di er- ence between A and B is the presence of more general class of sets and the addition of change-of-variables operators T . In the case = [0; 1) with the uniform partition of each of the segments [0; 1) on p intervals (i.e. S = p ), B becomes a C -subalgebra of A : Nn n n N N N Mp Mp p Mp Mp 2 ( ) B (C )  (C ) A ; = = n=0 where are binomial coecients. It is interesting to note that the generalization of the structure of A simpli es the proof of main results. Let us discuss a norm of approximation of continuous operators from A by discrete operators from A . For simplicity, consider the case = [0; 1) and M = 1. The general- 2 2 ization to M > 1 is similar. Consider a multi-kernel operator A : L ( ) ! L ( ) of the form, common in applications, A u(k) = A (k; x )u(k + x )dx ; k 2 ; u 2 L ( ): (25) j j [0;1) f1;:::;Ng Here, we also use the notation x ; n 2 ; x = (xe ); xe = n n 0; n 62 : Note that x = x + x , where = f1; :::; Ngn is the complement to the set . Consider the uniform partition of onto p identical cubes = a + [0; 1=p) ; (26) i i i i=1 where N j1 a = (b ) ; i 1 = p b ; b 2 f0; :::; p 1g (27) i j j j j=1 j=1 is the representation of i 1 in the base p numeral system. Let us take the approximation of A by N N p p X XX A = E [A ]; (28) ij ij j j i=1 j=1 f1;:::;Ng where " " A (a + ""; (a + "") ); (a ) = (a ) ; i j i j A = and """ = (1=(2p)) : (29) ij r=1 0; otherwise; The next proposition shows us that A ! A in the operator norm. 6 Theorem 1.3. Suppose that A (k; x ) in (25) are real functions with bounded rst deriva- tives. Then 2 2 1 kA Ak 6 krA k ! 0 for p ! 1: (30) c L L !L 2p f1;:::;Ng Along with Theorem 1.1 and Corollary 1.2, we can use Theorem 1.3 to determine the spec- trum and the inverse operator. Corollary 1.4. Under the assumptions of Theorem 1.3, we assume also that A is self- adjoint, i.e. A (k; x ) = A (k + x ; k ). Then A is self-adjoint and sp(A )  B (sp(A)); sp(A)  B (sp(A )); (31) c   c where B means -neighbourhood and we can set  = krA k 1 . If 0 62 sp(A ) L c f1;:::;Ng 2p 1 1 2 2 then A is invertible and there is p such that 0 62 sp(A) and kA k 6  , and the L !L following estimate is true 1 2 kA k 2 2 1 1 L !L kA A k 2 2 6 : (32) L !L 2 2 1 kA k L !L Moreover, RHS of (32) tends to 0 for p ! 1. The rest of the paper is organized as follows. Section 2 contains two examples: 1) new formulas for the functions of 1D Fredholm integral operators with step kernels; 2) the application of the method for obtaining a solution (with arbitrary precision) and spectral estimates for 3D discrete Schr odinger equation with planar, guided, and local potential defects. A short proof of the main result based on the explicit representation of a semigroup algebra of subsets is given in Section 3. We conclude in Section 4. 2. Examples Example 1. Consider the case N; M = 1, and the classical Fredholm integral operators (see [18]) 2 2 A : L ([0; 1)) ! L ([0; 1)); Au(x) = B(x; y)u(y)dy (33) with S-step (piecewise constant) kernels i1 i 1; x 2 [ ; ); S S B(x; y) = S B  (x) (y); B 2 C;  (x) = (34) ij i j ij i 0; otherwise: i;j=1 SS Such operators form an algebra isomorphic to C (see, e.g., [19]). But, this algebra does not contain the identity operator (Iu = u). Let us supplement it by adding new operators Au(x) = A(x)u(x) + B(x; y)u(y)dy; where A(x) = A  (x): (35) i i i=1 7 In other words > > A(x) =  (x)A(x); B(x; y) = S (x)B(y); where (36) S S = ( ) ; A = diag(A ); B = (B ) : (37) i i ij i=1 i;j=1 Taking a = (i 1)=S and using notations (4)-(5), we obtain S S X X f1g A = E [A ] + E [B ]: (38) i ij ii ij i=1 i;j=1 Applying the results of Theorem 1.1 and Corollary 1.2 along with (35)-(37) we obtain that 1 > 1 > 1 1 A u(x) =  (x)A (x)u(x) + S  (x)((A + B) A )(y)u(y)dy; (39) if and only if A and A + B are invertible (otherwise A is non-invertible). The spectrum of A consists of all the eigenvalues of A and A + B. In general, > > f (A)u(x) =  (x)f (A)(x)u(x) + S  (x)(f (A + B) f (A))(y)u(y)dy (40) for various functions f (e.g. rational/algebraic/transcendent:exp; ; :::). We have reduced explicitly the functional calculus on integral operators to the functional calculus on matrices. Let us consider a concrete example with some oscillating non-symmetric integral kernel, say A u(x) = A(x)u(x) + B(x; y)u(y)dy; (41) where A(x) = x sin(x) 1; B(x; y) = 4 sin(6xy) + exp(x) y + 2: (42) Taking the step approximation of A and B, see (34)-(37), with suciently large S, and using Corollary 1.2, we can estimate the eigenvalue  closest to 0. It is   0:1 6= 0. Hence, 0 0 the operator A has an inverse. Note that we can not use a Neumann series for A , since the kernel B(x; y) is not small relative to the rst term A(x). Let us approximately compute 1 3 A and, say, A , where the real branch of cubic root is assumed. To do this, we apply (40) to the step approximations of A and B. Thus, Z Z 1 1 u(x) 1 1 3 A u(x) = + C (x; y)u(y)dy; A u(x) = (A(x)) u(x) + D(x; y)u(y)dy: (43) A(x) 0 0 Computations show a good convergence of step kernels obtained by (40). For example, the relative deviations computed for S = 60, 300, and 1500 are suciently small kC C k 1 kC C k 1 60 300 L 300 1500 L = 0:1; = 0:02; kC k 1 kC k 1 300 L 300 L 8 1 1 kD D k kD D k 60 300 L 300 1500 L = 0:07; = 0:01; 1 1 kD k kD k 300 L 300 L while kC k 1  kC k 1  kC k 1  30; kD k 1  kD k 1  kD k 1  5: 60 L 300 L 1500 L 60 L 300 L 1500 L The corresponding integral kernels are plotted in Fig. 2. Figure 2: The original integral kernel B(x; y), see (42), and the step approximations of integral kernels C (x; y), D(x; y), see (43), computed by (40) for S = 60 and S = 300. Example 2. Actually, the method based on the expansion in E works well for more nm complex mixed multidimensional integral operators. Let us consider a discrete normalized 3D time-dependent Schr odinger operator with three planar potentials V , V , V with their 1 2 3 (i+")t intersections (see Fig. 3) and with attenuated harmonic source term e ( 2 R is the frequency, " > 0 is the attenuation factor, and the amplitude is 1) located at the origin. Then the corresponding equation on the wave function is (i+")t 3 3 i  V =  e ; n = (n ) 2 Z ; (44) discr n n n n0 i i=1 @t 9 3 Figure 3: The lattice Z with planar potential defects V , line potential defects V + V , and point defect i i j V + V + V . 1 2 3 where V ; n = 0; n 6= 0; i i j j6=i V + V ; n = n = 0; n 6= 0 (k 62 fi; jg) i j i j k V = ;  = ( ) (45) n discr n n n e V + V + V ; n = 0; 1 2 3 > n e 0; otherwise and ne denotes points adjacent to n. Our goal is to solve the Schr odinger equation and to (i+")t ^ ^ nd . Assuming (t) = e with time-independent and taking Fourier series n n n n 2in k 3 3 ^ ^ (k) = e , k = (k ) 2 [0; 1) we can rewrite (44) in the integral form n i n2Z i=1 A(k) (k) Z Z Z 1 1 1 ^ ^ ^ V (x ; k ; k )dx V (k ; x ; k )dx V (k ; k ; x )dx = 1; (46) 1 1 2 3 1 2 1 2 3 2 3 1 2 3 3 0 0 0 where A(k) = i" +  4 sin k :This 3D system contains non-parallel potential i=1 defects and, hence, does not admit explicit procedure for nding (see, e.g., [25] for 2D case). There are various methods to obtain the approximated solution of (46). We will use the method based on the expansion of the operator (46) in the basis E to obtain the nm 10 ^ approximation of with an arbitrary precision. Let p 2 N, h = 1=p. Then, following (4) we have 3 3 3 [0; 1) = fa + [0; h) g; N = f0; ::; p 1g ; a = hn: (47) n n n2N It is convenient to use the indices n; m 2 N in (5) instead of numbers i; j. Introducing A = A(a + 1h=2) (where 1 = (1; 1; 1)) we can write the approximation of the equation n n ^ ee A = 1 (A denotes LHS of (46)) as A = 1, where X X X X ; f1g f2g f3g A = E [A ] E [hV ] E [hV ] E [hV ]: (48) n 1 2 3 nn nm nm nm n2N n;m2N n;m2N n;m2N n =m n =m n =m 2 2 1 1 1 1 n =m n =m n =m 3 3 3 3 2 2 3 3 p p Then the matrices A 2 C , see (14), (15), are de ned by A = diag(A ); A = hV (  ) ; A = hV (  ) ; ; n f1g 1 n m n m n;m2N f2g 2 n m n m n;m2N 2 2 3 3 1 1 3 3 A = hV (  ) ; A = A = A = A = 0: (49) f3g 3 n m n m n;m2N f1;2g f1;3g f2;3g f1;2;3g 1 1 2 2 The matrices B , see (15), are B = A ; B = A + A (i = 1; 2; 3); ; ; fig ; fig B = A + A + A ; B = A + A + A ; f1;2g ; f1g f2g f1;3g ; f1g f3g B = A + A + A ; ; B = A + A + A + A : (50) f2;3g ; f2g f3g f1;2;3g ; f1g f2g f3g To compute A we follow (19) and (17), (18). De ne the matrices C by 1 1 1 1 1 1 1 C = B ; C = B B (i = 1; 2; 3); C = B B B + B ; ; fig f1;2g ; fig ; f1;2g f1g f2g ; 1 1 1 1 1 1 1 1 C = B B B + B ; C = B B B + B ; f1;3g f2;3g f1;3g f1g f3g ; f2;3g f2g f3g ; 1 1 1 1 1 1 1 1 C = B B B B + B + B + B B : (51) f1;2;3g f1;2;3g f1;2g f1;3g f2;3g f1g f2g f3g ; e e e All the matrices B are invertible since A exists (because A + i" is self-adjoint). Then has the form Z Z 1 1 1 > e e = A 1 =  C  + pC (x ; k ; k )dx + pC (k ; x ; k )dx + ; f1g 1 2 3 1 f2g 1 2 3 2 0 0 Z Z Z 1 1 1 pC (k ; k ; x )dx + p C (x ; x ; k )dx dx + f3g 1 2 3 3 f1;2g 1 2 3 1 2 0 0 0 Z Z Z Z 1 1 1 1 2 2 p C (x ; k ; x )dx dx + p C (k ; x ; x )dx dx + f1;3g 1 2 3 1 3 f2;3g 1 2 3 2 3 0 0 0 0 Z Z Z 1 1 1 p C (x ; x ; x )dx dx dx ; (52) f1;2;3g 1 2 3 1 2 3 0 0 0 11 Q where the vector-valued function  = (  (k )) (see also (34) for  ). Due to the n i n2N n i=1 i simple form of , all the integrals in (52) can be computed explicitly: for example, Z Z 1 1 p (x ; k ; x )dx dx = ( (k )) : 1 2 3 1 3 n 2 n2N 0 0 e e So, in (52) is given explicitly. The norm of di erence between A (46) and A (48) depends on how much A approximates A, at least kA Ak 2 2 6 24h, see also Theorem 1.3. n L !L 1 1 1 1 The norm of the inverse operators kA k 2 2 6 " and kA k 2 2 6 " since A + i" L !L L !L 1 1 2 e e and A + i" are self-adjoint. Hence kA A k 2 2 6 24h=" . In particular, L !L 2 2 2 2 ^ e ^ e k k 6 24h=" ; j j 6 (24h=" ) ; (53) L n n n2Z e e where are Fourier coecients of . On the side, we obtain the following spectral estimates dist(sp(A); sp(B )) 6 24h: (54) f1;2;3g Estimates (53), (54) become better and better for h ! 0. Depending on the number j j, the components sp(B )n ([ sp(B )) approximate the volume, planar, guided, and local j j<j j isolated spectral components of A which correspond to bulk, planar, guided, and local wave functions. These waves propagate along the potential defects of the corresponding dimension 3j j and exponentially attenuate in the perpendicular directions to the defect, see, e.g., the corresponding discussion in [25]. 3. Proof of Theorems 1.1 and 1.3 At rst, let us consider the semigroup of subsets G = fe :  f1; :::; Ngg; e e = e (55) and the corresponding C -algebra M = A e : A 2 C : (56) f1;:::;Ng The identity element in this algebra is 1 = e , where ; is the empty set. All the basis elements e = e are self-adjoint. De ne the mapping : M ! C ; (57) X X A e = (B ) ; B = A : (58) 1 f1;:::;Ng f1;:::;Ng Lemma 3.1. The mapping  is the C -isomorphism. The inverse mapping is de ned by X X 1 j n j ((B ) ) = A e ; A = (1) B : (59) f1;:::;Ng f1;:::;Ng Proof. Consider the following basis in M f = e (1 e ): (60) fng n62 Direct calculations give us f = f ; f f =  f : (61) This means that ff g can be considered as an orthogonal (in the algebraic sense) f1;:::;Ng 2 N basis in C because the number of elements #f :  f1; :::; Ngg = 2 . Hence, M is isomorphic to C with the isomorphism e $ f . Using X X Y X X Y A e = A e (1 e + e ) = A e (1 e ) (62) fng fng fng n62  n62 f1;:::;Ng f1;:::;Ng f1;:::;Ng X X X X X = A f = ( A )f = B f (63) f1;:::;Ng f1;:::;Ng f1;:::;Ng we obtain that the isomorphism e $ f coincides with  (57)-(58). Similarly, identities X X Y X X j n j B f = B e (1 e ) = B (1) e (64) fng n62 f1;:::;Ng f1;:::;Ng f1;:::;Ng X X X j n j = ( (1) B )e = A e (65) f1;:::;Ng f1;:::;Ng give us the form of the inverse mapping  (59). MM By Lemma 3.1, the C -algebra M = M of MM matrices with entries belonging MM 2 to M is isomorphic to (C ) . The corresponding isomorphism is X X MM 2 : M ! (C ) ;  ( A e ) = (B ) ; B = A ; (66) M M M f1;:::;Ng f1;:::;Ng X X 1 j n j ((B ) ) = A e ; A = (1) B : (67) f1;:::;Ng f1;:::;Ng SS MSMS MM The same result holds for the C -algebra M M M . For any A 2 C = = MS and  f1; :::; Ng, the basis elements F [A] = (  e A) 2 M satisfy the same ~ ~ MS ij ii jj ~ ~ i;j=1 equations as the basis elements E [A] 2 A (see (13)) ij F [A] = F [A ]; F [A] +F [B] = F [A + B]; F [A]F [B] =  F [AB]: (68) jk ij ji ij ij ij ij kl il This means that A = M with the natural isomorphism E [A] $ F [A]. Then, com- MS ij ij paring (14)-(16) and (66)-(67) (with MS instead of M ) we deduce that  in (16) is the isomorphism satisfying (17)-(18). Now, let us prove Theoem 1.3. It is seen that X X (A A)u(k) = (A (k; x ) A (k; (xe ) ))u(k + x )dx ; c i j j [0;1) \ f1;:::;Ng i=1 where k is the center of the cube for which k 2 , and xe is the center of . Then, the j j i i mean value estimates maxjA (k; x ) A (k; (xe ) )j 6 krA k i L x2 i 2p nish the proof of (30). Estimates (31) are classical spectral estimates for self-adjoint op- erator A and self-adjoint perturbed operator A = A + A A in terms of the norm of c c c perturbation  = kAA k 2 2 , see, e.g., [26]. If A is invertible then A is invertible for c L !L c 1 1 all suciently large p and kA k 2 2 ! kA k 2 2 for p ! 1. Moreover, using the L !L L !L Neumann series for the inverse of perturbed operator, we deduce that 1 2 kA k 2 2 1 1 1 n 1 L !L 2 2 2 2 kA A k = k (A (A A)) A k 6 : L !L L !L 1 kA k 2 2 L !L n=1 The last term tends to 0 since  tends to 0 for p ! 1. 4. Conclusion We have shown that the analysis of mixed multidimensional integral and some type of ergodic operators can be explicitly reduced to the analysis of special matrices. This allows us to compute functions of such operators and their spectra explicitly with an arbitrary precision. Funding statement This paper is a contribution to the project M3 of the Collaborative Research Centre TRR 181 "Energy Transfer in Atmosphere and Ocean" funded by the Deutsche Forschungs- gemeinschaft (DFG, German Research Foundation) - Projektnummer 274762653. This work is also supported by the RFBR (RFFI) grant No. 19-01-00094. References [1] A. B. Movchan, and L. I. Slepyan, Band gap green's functions and localized oscillations, Proc. R. Soc. A 463 (2007) 2709{2727. [2] D. J. Colquitt, M. J. Nieves, I. S. Jones, A. B. Movchan, and N. V. Movchan, Waves in lattices with imperfect junctions and localised defect modes, Proc. R. Soc. A 469 (2150) (2013) 20120579. 14 [3] M. Makwana, and R. V. Craster, Localised point defect states in asymptotic models of discrete lattices, Q. J. Mechanics Appl. Math. 66 (2013) 289{316. [4] A. A. Kutsenko, Algebra of multidimensional periodic operators with defects, J. Math. Anal. Appl. 428 (2015) 221{230. [5] A. A. Kutsenko, Recovery of defects from the information at detectors, Inverse Problems 32 (2016) [6] V. Caudrelier, Multisymplectic approach to integrable defects in the sine-gordon model, J Phys A Math Gen 48 (19) (2015) 195203. [7] M. F. Norman, Ergodicity of di usion and temporal uniformity of di usion approximation, J. Appl. Prob. 14 (1977) 399{404. [8] J. Schenker, Di usion in the mean for an ergodic Schrodinger equation perturbed by a uctuating potential, Comm. Math. Phys. 339 (2015) 859{901. [9] I. Prigogine, A. P. Grecos, and Cl. George , Kinetic theory and ergodic properties, Proc. Nati. Acad. Sci. USA 14 (1976) 1802{1805. [10] E. Kreyszig, Advanced Engineering Mathematics, Wiley, 2011. [11] V. Kirsh, and L. A. Pastur , Ergodic theory and discrete one-dimensional random Schrodinger operators: uniform existence of the Lyapunov exponent , Contemp. Math. 327 (2003) 1{23. [12] D. Lenz, F. Schwarzenberger, and I. Veselic, A Banach space-valued ergodic theorem and the uniform approximation of the integrated density of states , Geom. Dedicata 150 (2011) 1{34. [13] E. Korotyaev, N. Saburova, E ective masses for laplacians on periodic graphs, J. Math. Anal. Appl. 436 (2016) 104{130. [14] F. Vico, M. Ferrando, L. Greengard, Z. Gimbutas, The decoupled potential integral equation for time- harmonic electromagnetic scattering, Comm. Pure Appl. Math. 69 (2016) 771{812. [15] T. Eisner, B. Farkas, M. Haase, and R. Nagel, Operator theoretic aspects of ergodic theory, Springer, [16] V. Kirsh, and L. A. Pastur , Analogues of Szego's theorem for ergodic operators , Mat. Sb. 206 (2015) 93{119. [17] C. Houdayer, and Y. Isono, Bi-exact groups, strongly ergodic actions and group measure space type III factors with no central sequence , Comm. Math. Phys. 348 (2016) 991{1015. [18] E. I. Fredholm, Sur une classe d'equations fonctionnelles, Acta Math. 27 (1903) 365{390. [19] M. S. Gockenbach, Finite-Dimensional Linear Algebra, Discrete mathematics and its applications., CRC Press, 2010. [20] S. Albeverio, E. I. Gordon, A. Yu. Hrennikov, Finite-dimensional approximations of operators in the hilbert spaces of functions on locally compact abelian groups, Acta Appl. Math. 64 (2000) 33{73. [21] M. Rosenkranz, A new symbolic method for solving linear two-point boundary value problems on the level of operators, Journal of Symbolic Computation 39 (2005) 171{199. [22] V. V. Bavula, The algebra of integro-di erential operators on an ane line and its modules, J. Pure Appl. Algebra 217 (2013) 495{529. [23] L. Guo, G. Regensburger, M. Rosenkranz, On integro-di erential algebras, J. Pure Appl. Algebra 218 (2014) 456{473. [24] A. Kutsenko, Mixed multidimensional integral operators with piecewise constant kernels and their representations, Linear and Multilinear Algebra 0 (2017) 0. doi:10.1080/03081087.2017.1415294. [25] A. A. Kutsenko, Algebra of 2d periodic operators with local and perpendicular defects, J. Math. Anal. Appl. 442 (2016) 796{803. [26] T. Kato, Perturbation theory for linear operators, Springer-Verlag Berlin Heidelberg, 1995.

Journal

Mathematical Physics – arXiv (Cornell University)

Published: May 23, 2018

References

  • Band gap Green’s functions and localized oscillations
    Movchan
  • Waves in lattices with imperfect junctions and localised defect modes
    Colquitt
  • Localised point defect states in asymptotic models of discrete lattices
    Makwana
  • Algebra of multidimensional periodic operators with defects
    Kutsenko
  • Recovery of defects from the information at detectors
    Kutsenko
  • Multisymplectic approach to integrable defects in the Sine–Gordon model
    Caudrelier
  • Ergodicity of diffusion and temporal uniformity of diffusion approximation
    Norman
  • Diffusion in the mean for an ergodic schrodinger equation perturbed by a fluctuating potential
    Schenker
  • Kinetic theory and ergodic properties
    Prigogine