Rational approximations to the zeta function

Size: px
Start display at page:

Download "Rational approximations to the zeta function"

Transcription

1 Rational appoximations to the zeta function Keith Ball June 29, 27 Abstact This aticle descibes a sequence of ational functions which conveges locally unifomly to ζ The numeatos (and denominatos) of these ational functions can be expessed as chaacteistic polynomials of matices that ae on the face of it vey simple As a consequence, the Riemann hypothesis can be estated as what looks like a athe conventional spectal poblem but which is elated to the one found by Connes in his analysis of the zeta function Howeve the point hee is that the ational appoximations look to be susceptible of quantitative estimation Intoduction The Riemann hypothesis [R] states that the Riemann zeta function has all its nontivial zeos on the citical line {s : Rs = /2} (the tivial ones being at 2, 4, 6 and so on) The eal content of the conjectue is that ζ has no zeoes to the ight of the citical line The customay fomulation is equivalent by vitue of the functional equation poved by Riemann In ode to show that a holomophic function has no zeoes in (say) a half-plane it suffices to expess the function as a locally unifom limit of holomophic functions with no zeoes thee This aticle descibes a sequence of ational functions which begins (s ), s + 2(s ), 4s 2 + s + 9 6(s + 3)(s ), (s + 2)(3s 2 + s + ) 4(s 2 + 6s + )(s ), (s + 2)(72s s s + 25) 3(3s s 2 + 6s + 5)(s ), ()

2 We define the sequence as follows Fo each intege m we define and the coefficients (a m,j ) by ( p m (t) = ( t) t ) ( t ) 2 m p m (t) = ( ) j a m,j t j These coefficients ae escaled Stiling numbes of the fist kind but it is moe convenient fo us to use the diffeent indexation and scaling hee We then set F m (s) = a m,j B j s + j whee the B j ae the usual Benoulli numbes given by the geneating function and z e z = B j j! zj G m (s) = ( ) j a m,j j= s + j The ational functions in question ae the atios Fo example and F m (s) (s )G m (s) F 3 (s) = s 2s + 6(s + ) = 3s2 + s + 2(s )s(s + ) G 3 (s) = s 6s + s + 6(s + 2) = s 2 + 6s + 3(s )s(s + )(s + 2) It is immediate that the m th atio intepolates ζ at the points,, 2,, m and has a simple pole with esidue at s = We shall show that the sequence conveges locally unifomly to ζ, at least to the ight of the line {s : Rs = }, (with the obvious convention at s = ) The numeatos (and denominatos) of these ational functions can be expessed as chaacteistic polynomials of matices that 2

3 ae on the face of it vey simple One way to state this is that the numeato of the m th function is the deteminant of 2 3 m m m + ( s) m m+ m m 2 As a consequence, the Riemann hypothesis can be estated as what looks like a athe conventional spectal poblem: to show that the specta of cetain matices stay to the left of the citical line, at least asymptotically as m Needless to say, I spent some time thinking about this poblem without success but am cetain that I have not exhausted the possible lines of attack Even if the simplicity of these matices is indeed an illusion, as one would expect, thee ae concete easons to think that these ational appoximations to zeta might be useful fo estimating the size of zeta, in the sense of the Lindelöf hypothesis: see fo example the book of Patteson [P] Polyá suggested that the Riemann hypothesis should be poved by expessing the zeoes of zeta (otated onto the eal line) as eigenvalues of a self-adjoint opeato (The statement is often cedited to Hilbet: it was Tey Tao who pointed out the mistake to me) A numbe of candidates fo such opeatos have been poposed, coming mainly fom quantum theoy The best known of these was found by Connes [C] Thee is a connection between Connes infinite-dimensional opeato and the finite-dimensional ones descibed hee, which will be explained biefly in Section 5 This became appaent to me fom the vey eadable aticle of Lachaud [L] My hope is that the finite-dimensional opeatos and the ational functions they coespond to, ae susceptible of quantitative estimation that would not make sense fo the infinite-dimensional opeato It is well known that the zeoes of ζ should be modelled by the eigenvalues of cetain andom matices This oiginated in the wok of Dyson and Montgomey [M] and was expeimentally confimed by emakable calculations of Odlyzko [O] Katz and Sanak extended the model to othe L-functions [KaSa] In the past two decades a huge amount of wok has been done on this connection in paticula by 3

4 Keating, Snaith and thei collaboatos [KeSn] While stictly speaking this is only indiectly elated to the esults in this aticle the andom model is clealy a cucial inspiation In ode to pove the convegence of F m(s) (s )G m(s) F m (s) h s m Γ(s)ζ(s) to ζ(s) we shall show that and (s )G m (s) h s m Γ(s) whee h m is the patial sum m j= /j of the hamonic seies We can get a sense of why this is, quite easily Using the following vaiant of Konecke s fomula B j = ( ) j j k= it is easy to show that fo Rs > F m (s) = ( m k= k + k + k = k = If x is close to zeo then p m (x) is appoximately ( ) k ( ) ( + ) j (2) ( ) ) k ( ) p m (( + )x) x s 2 dx (3) m e x/i = e hmx i= So the sum ove in equation (3) is appoximately k =( ) ( k ) e hm(+)x = e hmx ( e hmx ) k Thus fo small values of x the integand in equation (3) is appoximately ( m k= ) k + e hmx ( e hmx ) k x s 2 If the appoximation wee good fo all x between and then F m (s) would be close to k= k + e hmx ( e hmx ) k x s 2 dx = hm s hm k= k + e y ( e y ) k y s 2 dy 4

5 The last integal plainly conveges to k= k + e y ( e y ) k y s 2 y dy = e y e y y s 2 dy as m povided Rs > and the latte is easily seen to be Γ(s)ζ(s) Ou fist aim will be to show that indeed h s m F m (s) Γ(s)ζ(s) locally unifomly fo Rs > (not just Rs > ) as m This looks like a tall ode Cossing the pole at s = is not the poblem The difficulty is that unless x is vey close to, the expessions m,k (x) = k = ( ) k ( ) p m (( + )x) involve values of p m at points well outside the inteval [, ], whee p m looks nothing like a negative exponential Indeed m,k (x) is a divided diffeence of p m and consequently equal to ( x) k p (k) m (u) fo some u between x and (k + )x that is not easily specified Since p m oscillates epeatedly on the inteval [, m + ] it would seem that m,k (x) could be vey lage in size and of moe o less andom sign So the following lemma comes as something of a shock Lemma If m is a non-negative intege and p m (x) = ( x)( x/2) ( x/m) then fo each intege k and each x [, ] ( ) k k m,k (x) = ( ) p m (( + )x) = It is tivial to check that m,k (x) = k= fo all x, so the lemma shows that fo each m the m,k fom a patition of unity on [, ] and thus automatically contols the sizes of the m,k as well as thei signs Once the lemma is established the convegence poof is faily staightfowad: this will be the content of Section below The obvious way to pove the Konecke fomula (2) mentioned above is to use the expansion y e y = k= 5 ( e y ) k k +

6 that aleady appeaed in the integal fomula fo Γ(s)ζ(s) So it might be logically moe easonable to define the F m by using the fomula ( m ) F m (s) = k + m,k(x) x s 2 dx k= and simply avoid mention of the Benoulli numbes Howeve it seemed a little odd to define a ational function with known poles and esidues as the analytic continuation of an integal The convegence poof just alluded to elies on the fact that the F m ae defined as sums which we can pass though integal signs The point of the second section of the aticle will be to povide a bidge between the definition of the F m and thei epesentation as chaacteistic polynomials: in othe wods to epesent the F m as something moe like a poduct than a sum The main fomula in Section 2 is the following ecuence fo the F m : fo each m Thus and so on (s + m )F m (s) = m m + + (m + ) F m j (s) j(j + ) (s )F (s) = sf (s) = 2 + F (s) j= (s + )F 2 (s) = F (s) F (s) If we teat the fist m + of these elations as a linea system fo the values F (s), F (s),, F m (s) we can expess the fact that F m (s) = by the vanishing of a cetain deteminant In Section 3 shall show that this deteminant can be witten as det(l m + ( s)u m ) whee as mentioned ealie L m is the (m + ) (m + ) Toeplitz matix m+ m m 2 6

7 and U m is the matix 2 3 m 2 3 m 3 m m If we set s = z/(z ) the deteminant becomes det(zl m (L m +U m )) (apat fom a facto of ( z) m+ ) The Riemann hypothesis would follow if this deteminant vanishes only at points with z o equivalently that the matix I m+ + L m U m has spectal adius at most, whee I m+ is the (m+) (m+) identity matix Fo small values of m this is tue In the fist vesion of this aticle I stated that thee ae good easons to believe that the zeoes of the F m do leak acoss the citical line (and then come back again) Pace Nielsen quickly confimed that when m = 643 thee is a zeo to the ight of the citical line He also infomed me that when m = 27 thee is a zeo with eal pat lage than (which I had not expected) The Riemann hypothesis is equivalent to the statement that the spectal adius of I m+ + L m U m is at most + o() as m In Section 4 of the aticle I shall explain why I think that the appoximations F m (s) might be useful to estimate the size of ζ The main point is that wheeas appoximations to zeta that ae sums of powes oscillate wildly all the way up the citical line, a polynomial of degee m cannot oscillate too often Wheneve one has a new sequence of appoximations to ζ it is natual to ask whethe they can help to pove Diophantine popeties (iationality o tanscendence) of values of the zeta function and most especially Eule s constant ( γ = lim ζ(s) ) s s Since the appoximations descibed hee ae ational functions (with intege coefficients) they do povide ational appoximations to Eule s constant but fo this paticula value of zeta the appoximations ae not new 7

8 Acknowledgements I am extemely gateful to David Peiss fo his advice duing this wok My thanks also to Tey Tao and Pete Sanak who made vey helpful suggestions concening the pesentation The key lemma and convegence In the intoduction we defined, fo each m, p m (x) = ( x)( x/2) ( x/m), and fo each k m,k (x) = k = ( ) k ( ) p m (( + )x) Note that the sum makes sense and is zeo if k < o k > m We also intoduced the function F m (s) as a ational function, F m (s) = We also mentioned the Konecke fomula B j = ( ) j j k= k + a m,j B j s + j k = ( ) k ( ) ( + ) j It is a consequence of standad popeties of the binomial coefficients that fo all j between and m, the sum on the ight is unchanged if the uppe limit is inceased fom j to m This implies that F m (s) = = Fo Rs > the sum ove j can be witten as ( ) ( ) j a m,j k k ( ) ( + ) j j= s + j k= k + = ( ) k k ( ) m ( ) j a mj( + ) j k= k + = j= s + j ( ) j a mj ( + ) j x s+j 2 dx = p m (( + )x)x s 2 dx j= and so fo Rs > we have F m (s) = k= m,k (x) k + xs 2 dx 8

9 Fom now on we use f m (x) to denote m k= m,k (x) k+ The aim of this section is to pove the following theoem Theoem 2 hm s (s )F m (s) (s )Γ(s)ζ(s) locally unifomly fo Rs > (with the obvious convention at s = ) It is easy to check that fo Rs > Γ(s)ζ(s) = y e y e y y s 2 dy If we set u = e y fo y >, then since < u <, Theefoe y log( u) = e y u Γ(s)ζ(s) = k= = k= u k k + = k= ( e y ) k k + k + e y ( e y ) k y s 2 dy Most of the effot in poving Theoem 2 will go into showing that the tuncated functions f m (x/h m ) [,hm] convege to x e x = k + e x ( e x ) k k= on [, ) with the convegence dominated by a negative exponential function It is clea that fo each fixed x, p m (x/h m ) e x and hence that fo each fixed k and x m,k (x/h m ) e x ( e x ) k (4) as m We need to establish two types of dominance: one to confim that the sum f m (x) = k= m,k (x) k + conveges pointwise in x to k= k+ e x ( e x ) k and one to check that this convegence is dominated on [, ) Fo almost evey estimate we make it is essential to have the key lemma stated in the intoduction: 9

10 Lemma If m is a non-negative intege and p m (x) = ( x)( x/2) ( x/m) then fo each intege k and each x [, ] m,k (x) = k = ( ) k ( ) p m (( + )x) We also need a simple popety of the divided diffeences that depends only upon the fact that p m is a polynomial of degee at most m Lemma 3 If j is a non-negative intege then fo evey x, (k + )(k + 2) (k + j) m,k (x) = j!p m ( jx) k= In paticula fo j = m,k (x) = p m () = k= Poof We shall confim that fo any polynomial q of degee at most m ) k (k + )(k + 2) (k + j) k= =( ) ( k q(( + )x) = j!q( jx) and in checking this we may assume that x = So ou aim is to veify that fo each such q ) k (k + )(k + 2) (k + j) k= =( ) ( k q( + ) = j!q( j) It suffices to check this fo each polynomial of the fom q n : t (t )(t 2) (t n) with n m The intenal sum vanishes if q has degee less than k and hence it vanishes fo q n if k > n It also vanishes if k < n because of the fom of q n The only emaining case is k = n and in that case the intenal sum has value ( ) n n! So the double sum is (n + )(n + 2) (n + j)( ) n n! = ( ) n (n + j)! = j!q n ( j)

11 The poof of Lemma involves the intoduction of an additional paamete as follows Fo each v define and P m (v, x) = m!p m (x v) = (v + x)(v + 2 x) (v + m x) m,k (v, x) = k = ( ) k ( ) P m (v, ( + )x) Obseve that P m (, x) = m!p m (x) and m,k (, x) = m! m,k (x) So Lemma follows fom: Lemma 4 If m is a non-negative intege, k is an intege, v and x m,k (v, x) Poof We use induction on m When m =, m,k (v, x) is zeo unless k = in which case it is We claim that fo m > m,k (v, x) = (v + x) m,k (v +, x) + kx m,k (v + x, x) Once this is established the inductive step is clea because we can assume that k and fo the given ange of v and x, the numbe v + x is also at least Now fo any v and x P m (v, x) = (v + x)p m (v +, x) and so m,k (v, x) = k = ( ) k ( ) (v + ( + )x)p m (v +, ( + )x) ( ) k k = ( ) (v + x)p m (v +, ( + )x) = ( ) k k ( ) xp m (v +, ( + )x) = ( ) = (v + x)) k k m,k (v +, x) k ( ) xp m (v +, ( + )x) = (v + x)) m,k (v +, x) + kx = k = ( ) k ( ) P m (v +, ( + 2)x) = (v + x) m,k (v +, x) + kx m,k (v + x, x)

12 whee the last step follows fom the fact that fo all m, v and x, P m (v +, ( + 2)x) = P m (v + x, ( + )x) By combining Lemmas and 3 we can immediately make some estimates fo the m,k (x) that will give us pat of the dominance we need to get convegence Lemma 5 Fo each m, each k and each x [, h m ] and fo each m m,k (x/h m ) f m (x/h m ) p m (x/h m ) = k= (k + ) ex k + m,k(x/h m ) e x Poof Fo the fist one we apply Lemma 3 with j = and use the positivity of the m,k to deduce that fo each k m,k (x) k + p m( x) = k + ( + x)( + x/2) ( + x/m) (k + ) ehmx Fo the second one we obseve that k= k + m,k(x) k m,k (x) = p m ( x) p m () e hmx k= As aleady emaked it is clea that fo each fixed x, p m (x/h m ) e x and hence that fo each fixed k and x as m Fom Lemma 5 we have m,k (x/h m ) e x ( e x ) k (5) m,k (x/h m ) (k + ) ex 2

13 so the convegence in (5) is dominated (on the space of non-negative integes with counting measue) by a sequence summable against (/(k + )) Hence fo each x > f m (x/h m ) = k= k + m,k(x/h m ) x e x We have that hm s hm F m (s) = f m (x/h m )x s 2 dx In ode to use dominated convegence on [, ) we need an estimate fo f m which we pove by intoducing anothe exta paamete Fo each m and p define f m (p, x) = and obseve that f m (x) = f m (, x) Lemma 6 Fo each m and p k= k + + p m,k(x) f m (p, x) = f m (p, x) + px m f m (p, x) (p + )x m f m (p +, x) Poof As long as p > we have f m (p, x) = = k= m,k (x) ( + u) K m,x (u) ( + u) p+ du k+2+p du whee K m,x (u) = k= ( + u) k+ m,k(x) This function is holomophic on the plane apat fom its pole at Its deivatives at ae successively m,k (x) k= (k + ) m,k (x) k= (k + )(k + 2) m,k (x) k= 3

14 and so on and theefoe by Lemma 3 its powe seies expansion at is Theefoe, fo u < ( ) j p m ( jx)u j j= K m,x (u) = ( ) j p m ( jx)u j j= = ( ( ) j p m ( jx) + jx ) u j j= m = K m,x (u) + x ( ) j p m ( jx) ju j m j= = K m,x (u) + x m uk m,x(u) The outemost identity continues analytically so it holds fo all u > Theefoe f m (p, x) = K m,x (u) ( + u) = f m (p, x) + x m = f m (p, x) x m = f m (p, x) + x m du p+ uk m,x(u) ( + u) K m,x (u) d u du ( + u) ( K m,x (u) p+ du p+ du p ( + u) p + p+ ( + u) p+2 = f m (p, x) + p m x f m (p, x) p + m x f m (p +, x) ) du Now we can estimate f m (x) as follows Using the key lemma we have an inequality f m (, x) = k= k + 2 m,k(x) 2 k= k + m,k(x) = 2 f m(x) povided x Then by Lemma (6) f m (x) = f m (, x) = f m (, x) x m f ( m (, x) f m (x) x ) e x/(2m) f m (x) 2m So by induction f m (x) e hmx/2 4

15 and we get the negative exponential dominance f m (x/h m ) e x/2 on the ange of integation [, h m ] This suffices to guaantee that fo Rs > hm s hm F m (s) = f m (x/h m )x s 2 dx Γ(s)ζ(s) We wish to coss the pole and so we need to modify the integand Fo Rs > we have ( Γ(s) ζ(s) ) s = = x e x xs 2 dx ( x e x e x ) x s 2 dx e x x s 2 dx The last integand behaves like x nea so the integal conveges locally unifomly fo Rs > and epesents the holomophic function Γ(s) ( ) ζ(s) s on this lage egion As in the intoduction we set and obseve that fo Rs > G m (s) = ( ) j a m,j j= s + j So on this half-plane G m (s) = p m (x)x s 2 dx h s hm m G m (s) = p m (x/h m )x s 2 dx Now p m (x/h m ) e x fo each fixed x and p m (x/h m ) e x as long as x h m so fo Rs > hm p m (x/h m )x s 2 dx Γ(s) s Theefoe, still only fo Rs >, we have h s hm ( m (F m (s) G m (s)) = (f m (x/h m ) p m (x/h m )) x s 2 dx Γ(s) ζ(s) ) s The integand is dominated as x by e x/2 but also, as x by e x owing to Lemma 5 Moeove, F m and G m both have esidue at s = so the diffeence 5

16 is holomophic fo Rs > So the integal epesents h s m (F m (s) G m (s)) on the lage egion and also conveges to Γ(s) ( ) ζ(s) s on this egion To complete the poof of Theoem 2 it suffices to show that (s )h s m G m (s) Γ(s) locally unifomly fo Rs > We have that fo Rs >, (s )h s hm m G m (s) = (s ) p m (x/h m )x s 2 hm dx = /h m p m(x/h m )x s dx because p m () = The latte integal conveges as long as Rs > so it epesents (s )hm s G m (s) on the lage egion It suffices to show that /h m p m(x/h m ) [,hm] e x fo x > and that the convegence is dominated by a negative exponential Obseve that p m is deceasing on [, ] so p m(x) is positive fo x < On the othe hand p m m(x) = p m (x) j= j x ( x/j) + p m (x) j=2 j=2 j e hmx ( + h m ) Thus fo x h m giving the equied dominance Also /h m p m(x/h m ) = h m /h m p m(x/h m ) e x (/h m + ) 2e x m ( x/(jh m )) + p m (x/h m ) j=2 h m j=2 j x/h m The fist tem is at most h m e (hm )x/hm and tends to while the second tem behaves like e x (h m )/h m and tends to e x This establishes Theoem 2 Finally we have that the atios F m (s) (s )G m (s) convege locally unifomly to ζ(s) fo Rs > My guess is that they do so on the entie complex plane 6

17 2 The bidge fom sum to deteminant The pupose of this section is to establish the ecuence (s + m )F m (s) = m m + + (m + ) F m j (s) j(j + ) fo m which will enable us to expess the numeato of F m as a deteminant A simila ecuence holds fo the functions G m : if m (s + m )G m (s) = (m + ) j= j= G m j (s) j(j + ) A small modification of the poof below actually yields this as well We begin with a simple emak Lemma 7 Fo each non-negative intege m we have lim sf m(s) = s m + Poof lim sf m(s) = f m () s All the m,k vanish at x = apat fom m,m since they involve only values of p m at the integes, 2,, m By Lemma 3 the m,k add up to so m,m () = So f m () = k= m,k () k + = m + Now fo the main lemma of the section Lemma 8 Fo each non-negative intege m (s + m )F m (s) = m m + + (m + ) F m j (s) j(j + ) j= Poof The two sides have the same limits at infinity so it suffices to check that they have the same esidues at each of the points,,,, 2 m At s = the esidue on the left is (m )a m, B while the esidue on the ight is (m + ) j= 7 a m j, B j(j + )

18 It thus suffices to check that fo each m, (m )a m, = (m + ) j= a m j, j(j + ) (The case = m is also obvious) Multiplying by x and summing ove m it is enough to check that m p m (x) xp m(x) = (m + ) j= p m j (x) j(j + ) Both sides ae polynomials of degee m so we need only check the values at x =, 2,, m Let k be one of these integes p m (k) = and it is easy to check that On the othe hand (m + ) j= p m j (k) j(j + ) k p m(k) = ( ) k ( m k ) = (m + ) m j=m k+ j(j + ) m j because p m j vanishes at k if m j k The latte expession is (m + ) j=m k+ ( ) m j j(j + ) = (m + ) = (m + ) ( ) k m j k = = ( ) k (m + ) i= ( k i k ( ) ( ) k = (m + ) = (m )(m + ) ( ) k ( ) (x m x m ) dx ( x m x) k dx x m k ( x) k dx and the beta integal gives the appopiate ecipocal of the binomial coefficient The ecuence elation given by Lemma 8 descibes a dynamical system fo the sequence F m (s) Numeically this system appeas to evolve vey slowly and indeed the convegence of the sequence is vey slow This makes the appoximations useless fo effective calculation of the value ζ(s) but suggests that it might be possible to tack the dynamical system: it is almost a continuous-time system 8 )

19 3 The spectum Fom the pevious section we have that fo each m (s + m )F m (s) = m m + + (m + ) F m j (s) j(j + ) (6) j= The fist m + of these elations give us the linea system s 2 2 s 3 3 s m+ m(m+) m+ 6 m+ 2 s + m F (s) F (s) F m (s) = 2 m+ So F m (s) can be witten as the atio of two deteminants The denominato is the deteminant of the (m+) (m+) matix on the left, namely (s )s(s+) (s+ m ) The numeato is the deteminant of the matix obtained by eplacing the last column of the oiginal matix with the vecto (, /2,, /(m + )) It will be moe convenient to move this vecto to the fist column of the matix thus intoducing a facto of ( ) m into the deteminant, and to change the sign of all the othe columns, thus emoving the facto again Then the numeato can be witten det s 2 3 m m+ 2 s s 6 2 m m (m )m 6 m+ m+ m(m+) 6 m s m Regad the columns as labelled,,, m We leave the zeo and columns unchanged as (, /2, /3,, /(m + )) and ( s,, /2,, /m) espectively We add the column to the 2 column to get (( s), ( s), 2,, 2/3,, 2/(m )) m+ 2 9

20 We add the (new) 2 column to the 3 column and we get (( s), ( s), ( s), 3, 3/2, 3/3,, 3/(m 2)) Continue in this way and afte all the additions divide the 2 column by 2, the 3 column by 3 and so on to get F m (s) = m! det(l m + ( s)u m ) (s )s (s + m ) (7) whee and L m = m m 2 m+ 2 3 m 2 3 m 3 m U m = (9) m Notice that F m does not have poles at 2, 4, because the coesponding Benoulli numbes vanish So the deteminant in the numeato picks up the tivial zeos of zeta at these numbes: indeed the facto s + 2 appeas as soon as m 3, the facto s + 4 as soon as m 5 and so on The zeoes of F m ae thus elated in a simple way to the spectum of L m U m This fomulation fo F m is pehaps the most elegant one in tems of a deteminant but it is inteesting to expess F m in a fom in which the poblem looks moe like a conventional spectal poblem To begin with we multiply column j by j fo each j but leave the zeo column unchanged This includes the m! facto appeaing in (7) into the deteminant We now subtact the m ow fom the last, the m 2 ow fom the m and so on (8) 2

21 to poduce the matix m(m+) (m )m 2 m (m 2)(m ) 2 +( s) We now add all the ows below the top one, to the top one, so that the second matix now has a zeo top ow The fist matix now has top ow ( m +, m, 2 m,, m ) Since the vaiable s now appeas on the diagonal in all places except the fist, ou aim is to educe the dimension by one so as to ceate a chaacteistic polynomial pope We add multiples (m + )/2, (m + )/6 and so on of the top ow to the successive ows below Since this eliminates the fist column below the fist ow, the deteminant is now the top left enty /(m + ) multiplied by the deteminant of the emaining squae So we get whee T m and R m ae as follows T m = det(t m + R m + (s )I m ) m (m )m 2 m m (m 2)(m ) 2 which is lowe tiangula with enty j in the jj diagonal place and enty j/((i j)(i j + )) in the ij place, if i > j R m is the ank one matix given by ( ) (m + )j R m = i(i + )(m j + ) 2 ij

22 We have F m (s) = m + det(t m + R m + (s )I m ) mj= (s + j ) We ae thus inteested in the spectum of the matix T m + R m whee T m is a cetain lowe tiangula matix and R m has ank It will be seen that this fomulation is actually somewhat close to the ecuence elation (6) If we set A m = T m + R m and B m = T m + R m I m then the deteminant in question is det(sa m (s )B m ) The complex numbes to the ight of the citical line ae those fo which s > s so a natual way to tackle the spectal poblem would be to ty to find a nom on C m with the popety that fo evey z C m A m z B m z The obvious choice would be an Hilbetian nom So we look fo a positive definite matix H fo which A mha m BmHB m is also positive definite o altenatively one fo which A m HA m B m HBm is positive definite The fom of the matix T m makes this vey tempting If we ignoe the ank one matix R m then we can cetainly compute the spectum of T m since it is lowe tiangula Howeve thee is a natual choice of nom which shows that the spectum is to the left of the citical line and theefoe povides a moe obust agument that one could ty to petub If H is the diagonal matix with enties, /2, /3,, /m on the diagonal then T m HTm (T m I m )H(T m I m ) = T m H + HTm H 22

23 is the matix (m )m (m )m (m 2)(m ) (m 2)(m ) 2 2 m This is obviously positive definite because it has negative off diagonal enties and ow sums that ae positive because of the familia telescoping sum n= n(n + ) = m + As emaked ealie thee ae good easons to think that the zeoes of the F m do leak though the citical line so that the best one can hope fo is to find matices H m with positive definite and ɛ m ( + ɛ m )A mh m A m B mh m B m Once one is in possession of the matices T m and R m one could confim that they yield the F m in a diect way Diagonalise T m as P DP and check that when you apply P and P to the components of R m you ecove the Benoulli and Stiling numbes Such an agument would be a bit hash on the eade since thee would be little motivation fo intoducing these paticula matices Moe impotantly the dynamical system descibed by the ecuence elation (6) is of inteest in itself 4 Estimating the size of ζ Numeical evidence indicates that the function f m (x/h m ) diffes fom x/(e x ) by only about h m /m at any point of [, h m ] and so we expect the atio F m (s) (s )G m (s) to povide a good appoximation to ζ at s = /2 + it as long as Γ(s) is as lage as h m /m This happens if t is at most a bit less than 2 log m In fact, numeical π evidence and ough calculations indicate that the atio is not too fa fom ζ fo t all the way up to log m At the same time thee ae good easons to think that F m (s) 23

24 does not oscillate significantly fo t lage than log m So we have the tantalising possibility that the two egions ovelap: the t < log m egion whee F m tells us about ζ and the t > log m egion whee F m is smooth enough to be estimated This discussion suggests that one should look at the asymptotic expansion fo F m (s) which stats off F m (s) = (m + )(s ) + c m (m + )(s )s + whee the coefficient c m gows logaithmically with m, the next coefficient like (log m) 2 and so on Howeve my feeling is that the moe pomising appoach is the usual one: to look at an integal (say) F m (s) G m (s) = (f m (x) p m (x))x s 2 dx and move the contou into the egion whee x s 2 is vey small if s = /2 + it with t lage Fo the genuine ζ integal x e x e x x s 2 dx this appoach is hopeless because the contou is foced up against the imaginay axis and hence picks up the poles of z z e z Being a polynomial, z f m (z) has no poles so the issue does not aise The poblem is to estimate f m off the eal line 5 The connection with Connes opeato The Toeplitz matix L m given in equation (8) can be thought of as acting on polynomials a + a x + a 2 x a m x m athe than sequences (a,, a m ) It does so by multiplication by the patial sum x j j + 24

25 of the seies fo log( x) (followed by tuncation back to a polynomial of degee m) x In this context the uppe tiangula matix U m in (9) maps the constant function to and fo each k the monomial x k to the sum + x + + x k k = xk k( x) Thus fo any polynomial q of degee m the image is q(t) q() dt x x t The opeato of Connes is built fom a multiplication opeato and an integal opeato much like these, acting on an infinite-dimensional function space Refeences [C] [KaSa] [KeSn] [L] [M] [O] [P] A Connes, Tace fomula in noncommutative geomety and the zeos of the Riemann zeta function, Sel Math, New Se 5 (999) N M Katz and P Sanak, Random matices, Fobenius eigenvalues, and monodomy, Colloquium publications, bf 45, Ameican Mathematical Society, Povidence, RI, (999) JP Keating and NC Snaith, Random matix theoy and ζ(/2 + it), Comm Math Phys, 24 (2) G Lachaud, Spectal analysis and the Riemann hypothesis, J Comp Appl Math, 6 (23) HL Montgomey, The pai coelation of the zeta function, Poc Symp Pue Math, 24 (973) A M Odlyzko, The 22 -nd zeo of the Riemann zeta function, Dynamical, Spectal, and Aithmetic Zeta Functions (M van Fankenhuysen and ML Lapidus, eds), Contempoay Math, Ame Math Soc, Povidence, RI, (2) S J Patteson, An intoduction to the theoy of the Riemann Zeta- Function, Cambidge studies in advanced mathematics 4, CUP (988) 25

26 [R] B Riemann, Übe die Anzahl de Pimzahlen unte eine gegebenen Gösse, Monatsb de Beline Akad, (858) 26

Numerical approximation to ζ(2n+1)

Numerical approximation to ζ(2n+1) Illinois Wesleyan Univesity Fom the SelectedWoks of Tian-Xiao He 6 Numeical appoximation to ζ(n+1) Tian-Xiao He, Illinois Wesleyan Univesity Michael J. Dancs Available at: https://woks.bepess.com/tian_xiao_he/6/

More information

ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},

ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0}, ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION E. J. IONASCU and A. A. STANCU Abstact. We ae inteested in constucting concete independent events in puely atomic pobability

More information

arxiv: v1 [math.co] 4 May 2017

arxiv: v1 [math.co] 4 May 2017 On The Numbe Of Unlabeled Bipatite Gaphs Abdullah Atmaca and A Yavuz Ouç axiv:7050800v [mathco] 4 May 207 Abstact This pape solves a poblem that was stated by M A Haison in 973 [] This poblem, that has

More information

Journal of Number Theory

Journal of Number Theory Jounal of umbe Theoy 3 2 2259 227 Contents lists available at ScienceDiect Jounal of umbe Theoy www.elsevie.com/locate/jnt Sums of poducts of hypegeometic Benoulli numbes Ken Kamano Depatment of Geneal

More information

Math 301: The Erdős-Stone-Simonovitz Theorem and Extremal Numbers for Bipartite Graphs

Math 301: The Erdős-Stone-Simonovitz Theorem and Extremal Numbers for Bipartite Graphs Math 30: The Edős-Stone-Simonovitz Theoem and Extemal Numbes fo Bipatite Gaphs May Radcliffe The Edős-Stone-Simonovitz Theoem Recall, in class we poved Tuán s Gaph Theoem, namely Theoem Tuán s Theoem Let

More information

Solution to HW 3, Ma 1a Fall 2016

Solution to HW 3, Ma 1a Fall 2016 Solution to HW 3, Ma a Fall 206 Section 2. Execise 2: Let C be a subset of the eal numbes consisting of those eal numbes x having the popety that evey digit in the decimal expansion of x is, 3, 5, o 7.

More information

A generalization of the Bernstein polynomials

A generalization of the Bernstein polynomials A genealization of the Benstein polynomials Halil Ouç and Geoge M Phillips Mathematical Institute, Univesity of St Andews, Noth Haugh, St Andews, Fife KY16 9SS, Scotland Dedicated to Philip J Davis This

More information

Method for Approximating Irrational Numbers

Method for Approximating Irrational Numbers Method fo Appoximating Iational Numbes Eic Reichwein Depatment of Physics Univesity of Califonia, Santa Cuz June 6, 0 Abstact I will put foth an algoithm fo poducing inceasingly accuate ational appoximations

More information

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr. POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and

More information

3.1 Random variables

3.1 Random variables 3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated

More information

Math 124B February 02, 2012

Math 124B February 02, 2012 Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial

More information

On the integration of the equations of hydrodynamics

On the integration of the equations of hydrodynamics Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious

More information

Σk=1. g r 3/2 z. 2 3-z. g 3 ( 3/2 ) g r 2. = 1 r = 0. () z = ( a ) + Σ. c n () a = ( a) 3-z -a. 3-z. z - + Σ. z 3, 5, 7, z ! = !

Σk=1. g r 3/2 z. 2 3-z. g 3 ( 3/2 ) g r 2. = 1 r = 0. () z = ( a ) + Σ. c n () a = ( a) 3-z -a. 3-z. z - + Σ. z 3, 5, 7, z ! = ! 09 Maclauin Seies of Completed Riemann Zeta 9. Maclauin Seies of Lemma 9.. ( Maclauin seies of gamma function ) When is the gamma function, n is the polygamma function and B n,kf, f, ae Bell polynomials,

More information

6 PROBABILITY GENERATING FUNCTIONS

6 PROBABILITY GENERATING FUNCTIONS 6 PROBABILITY GENERATING FUNCTIONS Cetain deivations pesented in this couse have been somewhat heavy on algeba. Fo example, detemining the expectation of the Binomial distibution (page 5.1 tuned out to

More information

Stanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012

Stanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012 Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,

More information

k. s k=1 Part of the significance of the Riemann zeta-function stems from Theorem 9.2. If s > 1 then 1 p s

k. s k=1 Part of the significance of the Riemann zeta-function stems from Theorem 9.2. If s > 1 then 1 p s 9 Pimes in aithmetic ogession Definition 9 The Riemann zeta-function ζs) is the function which assigns to a eal numbe s > the convegent seies k s k Pat of the significance of the Riemann zeta-function

More information

On a quantity that is analogous to potential and a theorem that relates to it

On a quantity that is analogous to potential and a theorem that relates to it Su une quantité analogue au potential et su un théoème y elatif C R Acad Sci 7 (87) 34-39 On a quantity that is analogous to potential and a theoem that elates to it By R CLAUSIUS Tanslated by D H Delphenich

More information

MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form

MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE ANDRAS VASY We conside second ode constant coefficient scala linea PDEs on R n. These have the fom Lu = f L = a ij xi xj + b i xi + c i whee a ij b i and

More information

A Bijective Approach to the Permutational Power of a Priority Queue

A Bijective Approach to the Permutational Power of a Priority Queue A Bijective Appoach to the Pemutational Powe of a Pioity Queue Ia M. Gessel Kuang-Yeh Wang Depatment of Mathematics Bandeis Univesity Waltham, MA 02254-9110 Abstact A pioity queue tansfoms an input pemutation

More information

arxiv: v1 [math.co] 6 Mar 2008

arxiv: v1 [math.co] 6 Mar 2008 An uppe bound fo the numbe of pefect matchings in gaphs Shmuel Fiedland axiv:0803.0864v [math.co] 6 Ma 2008 Depatment of Mathematics, Statistics, and Compute Science, Univesity of Illinois at Chicago Chicago,

More information

Exceptional regular singular points of second-order ODEs. 1. Solving second-order ODEs

Exceptional regular singular points of second-order ODEs. 1. Solving second-order ODEs (May 14, 2011 Exceptional egula singula points of second-ode ODEs Paul Gaett gaett@math.umn.edu http://www.math.umn.edu/ gaett/ 1. Solving second-ode ODEs 2. Examples 3. Convegence Fobenius method fo solving

More information

An Exact Solution of Navier Stokes Equation

An Exact Solution of Navier Stokes Equation An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in

More information

Journal of Inequalities in Pure and Applied Mathematics

Journal of Inequalities in Pure and Applied Mathematics Jounal of Inequalities in Pue and Applied Mathematics COEFFICIENT INEQUALITY FOR A FUNCTION WHOSE DERIVATIVE HAS A POSITIVE REAL PART S. ABRAMOVICH, M. KLARIČIĆ BAKULA AND S. BANIĆ Depatment of Mathematics

More information

6 Matrix Concentration Bounds

6 Matrix Concentration Bounds 6 Matix Concentation Bounds Concentation bounds ae inequalities that bound pobabilities of deviations by a andom vaiable fom some value, often its mean. Infomally, they show the pobability that a andom

More information

Math 2263 Solutions for Spring 2003 Final Exam

Math 2263 Solutions for Spring 2003 Final Exam Math 6 Solutions fo Sping Final Exam ) A staightfowad appoach to finding the tangent plane to a suface at a point ( x, y, z ) would be to expess the cuve as an explicit function z = f ( x, y ), calculate

More information

EFFECTS OF FRINGING FIELDS ON SINGLE PARTICLE DYNAMICS. M. Bassetti and C. Biscari INFN-LNF, CP 13, Frascati (RM), Italy

EFFECTS OF FRINGING FIELDS ON SINGLE PARTICLE DYNAMICS. M. Bassetti and C. Biscari INFN-LNF, CP 13, Frascati (RM), Italy Fascati Physics Seies Vol. X (998), pp. 47-54 4 th Advanced ICFA Beam Dynamics Wokshop, Fascati, Oct. -5, 997 EFFECTS OF FRININ FIELDS ON SINLE PARTICLE DYNAMICS M. Bassetti and C. Biscai INFN-LNF, CP

More information

F-IF Logistic Growth Model, Abstract Version

F-IF Logistic Growth Model, Abstract Version F-IF Logistic Gowth Model, Abstact Vesion Alignments to Content Standads: F-IFB4 Task An impotant example of a model often used in biology o ecology to model population gowth is called the logistic gowth

More information

Compactly Supported Radial Basis Functions

Compactly Supported Radial Basis Functions Chapte 4 Compactly Suppoted Radial Basis Functions As we saw ealie, compactly suppoted functions Φ that ae tuly stictly conditionally positive definite of ode m > do not exist The compact suppot automatically

More information

Berkeley Math Circle AIME Preparation March 5, 2013

Berkeley Math Circle AIME Preparation March 5, 2013 Algeba Toolkit Rules of Thumb. Make sue that you can pove all fomulas you use. This is even bette than memoizing the fomulas. Although it is best to memoize, as well. Stive fo elegant, economical methods.

More information

A Multivariate Normal Law for Turing s Formulae

A Multivariate Normal Law for Turing s Formulae A Multivaiate Nomal Law fo Tuing s Fomulae Zhiyi Zhang Depatment of Mathematics and Statistics Univesity of Noth Caolina at Chalotte Chalotte, NC 28223 Abstact This pape establishes a sufficient condition

More information

A Relativistic Electron in a Coulomb Potential

A Relativistic Electron in a Coulomb Potential A Relativistic Electon in a Coulomb Potential Alfed Whitehead Physics 518, Fall 009 The Poblem Solve the Diac Equation fo an electon in a Coulomb potential. Identify the conseved quantum numbes. Specify

More information

Surveillance Points in High Dimensional Spaces

Surveillance Points in High Dimensional Spaces Société de Calcul Mathématique SA Tools fo decision help since 995 Suveillance Points in High Dimensional Spaces by Benad Beauzamy Januay 06 Abstact Let us conside any compute softwae, elying upon a lage

More information

Lecture 8 - Gauss s Law

Lecture 8 - Gauss s Law Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.

More information

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two

More information

SOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF 2 2 OPERATOR MATRICES

SOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF 2 2 OPERATOR MATRICES italian jounal of pue and applied mathematics n. 35 015 (433 44) 433 SOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF OPERATOR MATRICES Watheq Bani-Domi Depatment of Mathematics

More information

2 x 8 2 x 2 SKILLS Determine whether the given value is a solution of the. equation. (a) x 2 (b) x 4. (a) x 2 (b) x 4 (a) x 4 (b) x 8

2 x 8 2 x 2 SKILLS Determine whether the given value is a solution of the. equation. (a) x 2 (b) x 4. (a) x 2 (b) x 4 (a) x 4 (b) x 8 5 CHAPTER Fundamentals When solving equations that involve absolute values, we usually take cases. EXAMPLE An Absolute Value Equation Solve the equation 0 x 5 0 3. SOLUTION By the definition of absolute

More information

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects Pulse Neuton Neuton (PNN) tool logging fo poosity Some theoetical aspects Intoduction Pehaps the most citicism of Pulse Neuton Neuon (PNN) logging methods has been chage that PNN is to sensitive to the

More information

KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS

KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS Jounal of Applied Analysis Vol. 14, No. 1 2008), pp. 43 52 KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS L. KOCZAN and P. ZAPRAWA Received Mach 12, 2007 and, in evised fom,

More information

arxiv: v1 [math.na] 8 Feb 2013

arxiv: v1 [math.na] 8 Feb 2013 A mixed method fo Diichlet poblems with adial basis functions axiv:1302.2079v1 [math.na] 8 Feb 2013 Nobet Heue Thanh Tan Abstact We pesent a simple discetization by adial basis functions fo the Poisson

More information

12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,

12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31, th WSEAS Int. Conf. on APPLIED MATHEMATICS, Caio, Egypt, Decembe 9-3, 7 5 Magnetostatic Field calculations associated with thick Solenoids in the Pesence of Ion using a Powe Seies expansion and the Complete

More information

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in

More information

Brief summary of functional analysis APPM 5440 Fall 2014 Applied Analysis

Brief summary of functional analysis APPM 5440 Fall 2014 Applied Analysis Bief summay of functional analysis APPM 5440 Fall 014 Applied Analysis Stephen Becke, stephen.becke@coloado.edu Standad theoems. When necessay, I used Royden s and Keyzsig s books as a efeence. Vesion

More information

A STABILITY RESULT FOR p-harmonic SYSTEMS WITH DISCONTINUOUS COEFFICIENTS. Bianca Stroffolini. 0. Introduction

A STABILITY RESULT FOR p-harmonic SYSTEMS WITH DISCONTINUOUS COEFFICIENTS. Bianca Stroffolini. 0. Introduction Electonic Jounal of Diffeential Equations, Vol. 2001(2001), No. 02, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.swt.edu o http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) A STABILITY RESULT

More information

Multiple Criteria Secretary Problem: A New Approach

Multiple Criteria Secretary Problem: A New Approach J. Stat. Appl. Po. 3, o., 9-38 (04 9 Jounal of Statistics Applications & Pobability An Intenational Jounal http://dx.doi.og/0.785/jsap/0303 Multiple Citeia Secetay Poblem: A ew Appoach Alaka Padhye, and

More information

EM Boundary Value Problems

EM Boundary Value Problems EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do

More information

On the ratio of maximum and minimum degree in maximal intersecting families

On the ratio of maximum and minimum degree in maximal intersecting families On the atio of maximum and minimum degee in maximal intesecting families Zoltán Lóánt Nagy Lale Özkahya Balázs Patkós Máté Vize Septembe 5, 011 Abstact To study how balanced o unbalanced a maximal intesecting

More information

Goodness-of-fit for composite hypotheses.

Goodness-of-fit for composite hypotheses. Section 11 Goodness-of-fit fo composite hypotheses. Example. Let us conside a Matlab example. Let us geneate 50 obsevations fom N(1, 2): X=nomnd(1,2,50,1); Then, unning a chi-squaed goodness-of-fit test

More information

OSCILLATIONS AND GRAVITATION

OSCILLATIONS AND GRAVITATION 1. SIMPLE HARMONIC MOTION Simple hamonic motion is any motion that is equivalent to a single component of unifom cicula motion. In this situation the velocity is always geatest in the middle of the motion,

More information

Encapsulation theory: the transformation equations of absolute information hiding.

Encapsulation theory: the transformation equations of absolute information hiding. 1 Encapsulation theoy: the tansfomation equations of absolute infomation hiding. Edmund Kiwan * www.edmundkiwan.com Abstact This pape descibes how the potential coupling of a set vaies as the set is tansfomed,

More information

-Δ u = λ u. u(x,y) = u 1. (x) u 2. (y) u(r,θ) = R(r) Θ(θ) Δu = 2 u + 2 u. r = x 2 + y 2. tan(θ) = y/x. r cos(θ) = cos(θ) r.

-Δ u = λ u. u(x,y) = u 1. (x) u 2. (y) u(r,θ) = R(r) Θ(θ) Δu = 2 u + 2 u. r = x 2 + y 2. tan(θ) = y/x. r cos(θ) = cos(θ) r. The Laplace opeato in pola coodinates We now conside the Laplace opeato with Diichlet bounday conditions on a cicula egion Ω {(x,y) x + y A }. Ou goal is to compute eigenvalues and eigenfunctions of the

More information

Quantum Mechanics II

Quantum Mechanics II Quantum Mechanics II Pof. Bois Altshule Apil 25, 2 Lectue 25 We have been dicussing the analytic popeties of the S-matix element. Remembe the adial wave function was u kl () = R kl () e ik iπl/2 S l (k)e

More information

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx. 9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can

More information

I. CONSTRUCTION OF THE GREEN S FUNCTION

I. CONSTRUCTION OF THE GREEN S FUNCTION I. CONSTRUCTION OF THE GREEN S FUNCTION The Helmohltz equation in 4 dimensions is 4 + k G 4 x, x = δ 4 x x. In this equation, G is the Geen s function and 4 efes to the dimensionality. In the vey end,

More information

Nuclear and Particle Physics - Lecture 20 The shell model

Nuclear and Particle Physics - Lecture 20 The shell model 1 Intoduction Nuclea and Paticle Physics - Lectue 0 The shell model It is appaent that the semi-empiical mass fomula does a good job of descibing tends but not the non-smooth behaviou of the binding enegy.

More information

Chapter 3: Theory of Modular Arithmetic 38

Chapter 3: Theory of Modular Arithmetic 38 Chapte 3: Theoy of Modula Aithmetic 38 Section D Chinese Remainde Theoem By the end of this section you will be able to pove the Chinese Remainde Theoem apply this theoem to solve simultaneous linea conguences

More information

Chapter 5 Linear Equations: Basic Theory and Practice

Chapter 5 Linear Equations: Basic Theory and Practice Chapte 5 inea Equations: Basic Theoy and actice In this chapte and the next, we ae inteested in the linea algebaic equation AX = b, (5-1) whee A is an m n matix, X is an n 1 vecto to be solved fo, and

More information

Multiple Experts with Binary Features

Multiple Experts with Binary Features Multiple Expets with Binay Featues Ye Jin & Lingen Zhang Decembe 9, 2010 1 Intoduction Ou intuition fo the poect comes fom the pape Supevised Leaning fom Multiple Expets: Whom to tust when eveyone lies

More information

arxiv: v1 [math.co] 1 Apr 2011

arxiv: v1 [math.co] 1 Apr 2011 Weight enumeation of codes fom finite spaces Relinde Juius Octobe 23, 2018 axiv:1104.0172v1 [math.co] 1 Ap 2011 Abstact We study the genealized and extended weight enumeato of the - ay Simplex code and

More information

Bounds for Codimensions of Fitting Ideals

Bounds for Codimensions of Fitting Ideals Ž. JOUNAL OF ALGEBA 194, 378 382 1997 ATICLE NO. JA966999 Bounds fo Coensions of Fitting Ideals Michał Kwiecinski* Uniwesytet Jagiellonski, Instytut Matematyki, ul. eymonta 4, 30-059, Kakow, Poland Communicated

More information

Geometry of the homogeneous and isotropic spaces

Geometry of the homogeneous and isotropic spaces Geomety of the homogeneous and isotopic spaces H. Sonoda Septembe 2000; last evised Octobe 2009 Abstact We summaize the aspects of the geomety of the homogeneous and isotopic spaces which ae most elevant

More information

Appendix A. Appendices. A.1 ɛ ijk and cross products. Vector Operations: δ ij and ɛ ijk

Appendix A. Appendices. A.1 ɛ ijk and cross products. Vector Operations: δ ij and ɛ ijk Appendix A Appendices A1 ɛ and coss poducts A11 Vecto Opeations: δ ij and ɛ These ae some notes on the use of the antisymmetic symbol ɛ fo expessing coss poducts This is an extemely poweful tool fo manipulating

More information

On the Quasi-inverse of a Non-square Matrix: An Infinite Solution

On the Quasi-inverse of a Non-square Matrix: An Infinite Solution Applied Mathematical Sciences, Vol 11, 2017, no 27, 1337-1351 HIKARI Ltd, wwwm-hikaicom https://doiog/1012988/ams20177273 On the Quasi-invese of a Non-squae Matix: An Infinite Solution Ruben D Codeo J

More information

STUDY OF SOLUTIONS OF LOGARITHMIC ORDER TO HIGHER ORDER LINEAR DIFFERENTIAL-DIFFERENCE EQUATIONS WITH COEFFICIENTS HAVING THE SAME LOGARITHMIC ORDER

STUDY OF SOLUTIONS OF LOGARITHMIC ORDER TO HIGHER ORDER LINEAR DIFFERENTIAL-DIFFERENCE EQUATIONS WITH COEFFICIENTS HAVING THE SAME LOGARITHMIC ORDER UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA doi: 104467/20843828AM170027078 542017, 15 32 STUDY OF SOLUTIONS OF LOGARITHMIC ORDER TO HIGHER ORDER LINEAR DIFFERENTIAL-DIFFERENCE EQUATIONS WITH COEFFICIENTS

More information

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50 woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,

More information

Scattering in Three Dimensions

Scattering in Three Dimensions Scatteing in Thee Dimensions Scatteing expeiments ae an impotant souce of infomation about quantum systems, anging in enegy fom vey low enegy chemical eactions to the highest possible enegies at the LHC.

More information

MATH 415, WEEK 3: Parameter-Dependence and Bifurcations

MATH 415, WEEK 3: Parameter-Dependence and Bifurcations MATH 415, WEEK 3: Paamete-Dependence and Bifucations 1 A Note on Paamete Dependence We should pause to make a bief note about the ole played in the study of dynamical systems by the system s paametes.

More information

New problems in universal algebraic geometry illustrated by boolean equations

New problems in universal algebraic geometry illustrated by boolean equations New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic

More information

GROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS

GROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS Annales Academiæ Scientiaum Fennicæ Mathematica Volumen 32, 2007, 595 599 GROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS Teo Kilpeläinen, Henik Shahgholian and Xiao Zhong

More information

The Substring Search Problem

The Substring Search Problem The Substing Seach Poblem One algoithm which is used in a vaiety of applications is the family of substing seach algoithms. These algoithms allow a use to detemine if, given two chaacte stings, one is

More information

Quasi-Randomness and the Distribution of Copies of a Fixed Graph

Quasi-Randomness and the Distribution of Copies of a Fixed Graph Quasi-Randomness and the Distibution of Copies of a Fixed Gaph Asaf Shapia Abstact We show that if a gaph G has the popety that all subsets of vetices of size n/4 contain the coect numbe of tiangles one

More information

Boundary Layers and Singular Perturbation Lectures 16 and 17 Boundary Layers and Singular Perturbation. x% 0 Ω0æ + Kx% 1 Ω0æ + ` : 0. (9.

Boundary Layers and Singular Perturbation Lectures 16 and 17 Boundary Layers and Singular Perturbation. x% 0 Ω0æ + Kx% 1 Ω0æ + ` : 0. (9. Lectues 16 and 17 Bounday Layes and Singula Petubation A Regula Petubation In some physical poblems, the solution is dependent on a paamete K. When the paamete K is vey small, it is natual to expect that

More information

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Physics 2A Chapter 10 - Moment of Inertia Fall 2018 Physics Chapte 0 - oment of netia Fall 08 The moment of inetia of a otating object is a measue of its otational inetia in the same way that the mass of an object is a measue of its inetia fo linea motion.

More information

Section 8.2 Polar Coordinates

Section 8.2 Polar Coordinates Section 8. Pola Coodinates 467 Section 8. Pola Coodinates The coodinate system we ae most familia with is called the Catesian coodinate system, a ectangula plane divided into fou quadants by the hoizontal

More information

Fractional Zero Forcing via Three-color Forcing Games

Fractional Zero Forcing via Three-color Forcing Games Factional Zeo Focing via Thee-colo Focing Games Leslie Hogben Kevin F. Palmowski David E. Robeson Michael Young May 13, 2015 Abstact An -fold analogue of the positive semidefinite zeo focing pocess that

More information

A proof of the binomial theorem

A proof of the binomial theorem A poof of the binomial theoem If n is a natual numbe, let n! denote the poduct of the numbes,2,3,,n. So! =, 2! = 2 = 2, 3! = 2 3 = 6, 4! = 2 3 4 = 24 and so on. We also let 0! =. If n is a non-negative

More information

Numerical Integration

Numerical Integration MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.

More information

arxiv: v1 [math.nt] 12 May 2017

arxiv: v1 [math.nt] 12 May 2017 SEQUENCES OF CONSECUTIVE HAPPY NUMBERS IN NEGATIVE BASES HELEN G. GRUNDMAN AND PAMELA E. HARRIS axiv:1705.04648v1 [math.nt] 12 May 2017 ABSTRACT. Fo b 2 and e 2, let S e,b : Z Z 0 be the function taking

More information

B. Spherical Wave Propagation

B. Spherical Wave Propagation 11/8/007 Spheical Wave Popagation notes 1/1 B. Spheical Wave Popagation Evey antenna launches a spheical wave, thus its powe density educes as a function of 1, whee is the distance fom the antenna. We

More information

Chapter 3 Optical Systems with Annular Pupils

Chapter 3 Optical Systems with Annular Pupils Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The

More information

arxiv: v1 [math.ca] 12 Mar 2015

arxiv: v1 [math.ca] 12 Mar 2015 axiv:503.0356v [math.ca] 2 Ma 205 AN APPLICATION OF FOURIER ANALYSIS TO RIEMANN SUMS TRISTRAM DE PIRO Abstact. We develop a method fo calculating Riemann sums using Fouie analysis.. Poisson Summation Fomula

More information

AQI: Advanced Quantum Information Lecture 2 (Module 4): Order finding and factoring algorithms February 20, 2013

AQI: Advanced Quantum Information Lecture 2 (Module 4): Order finding and factoring algorithms February 20, 2013 AQI: Advanced Quantum Infomation Lectue 2 (Module 4): Ode finding and factoing algoithms Febuay 20, 203 Lectue: D. Mak Tame (email: m.tame@impeial.ac.uk) Intoduction In the last lectue we looked at the

More information

Vanishing lines in generalized Adams spectral sequences are generic

Vanishing lines in generalized Adams spectral sequences are generic ISSN 364-0380 (on line) 465-3060 (pinted) 55 Geomety & Topology Volume 3 (999) 55 65 Published: 2 July 999 G G G G T T T G T T T G T G T GG TT G G G G GG T T T TT Vanishing lines in genealized Adams spectal

More information

4/18/2005. Statistical Learning Theory

4/18/2005. Statistical Learning Theory Statistical Leaning Theoy Statistical Leaning Theoy A model of supevised leaning consists of: a Envionment - Supplying a vecto x with a fixed but unknown pdf F x (x b Teache. It povides a desied esponse

More information

Semicanonical basis generators of the cluster algebra of type A (1)

Semicanonical basis generators of the cluster algebra of type A (1) Semicanonical basis geneatos of the cluste algeba of type A (1 1 Andei Zelevinsky Depatment of Mathematics Notheasten Univesity, Boston, USA andei@neu.edu Submitted: Jul 7, 006; Accepted: Dec 3, 006; Published:

More information

On the ratio of maximum and minimum degree in maximal intersecting families

On the ratio of maximum and minimum degree in maximal intersecting families On the atio of maximum and minimum degee in maximal intesecting families Zoltán Lóánt Nagy Lale Özkahya Balázs Patkós Máté Vize Mach 6, 013 Abstact To study how balanced o unbalanced a maximal intesecting

More information

Analysis of simple branching trees with TI-92

Analysis of simple branching trees with TI-92 Analysis of simple banching tees with TI-9 Dušan Pagon, Univesity of Maibo, Slovenia Abstact. In the complex plane we stat at the cente of the coodinate system with a vetical segment of the length one

More information

Asymptotically Lacunary Statistical Equivalent Sequence Spaces Defined by Ideal Convergence and an Orlicz Function

Asymptotically Lacunary Statistical Equivalent Sequence Spaces Defined by Ideal Convergence and an Orlicz Function "Science Stays Tue Hee" Jounal of Mathematics and Statistical Science, 335-35 Science Signpost Publishing Asymptotically Lacunay Statistical Equivalent Sequence Spaces Defined by Ideal Convegence and an

More information

Physics 505 Homework No. 9 Solutions S9-1

Physics 505 Homework No. 9 Solutions S9-1 Physics 505 Homewok No 9 s S9-1 1 As pomised, hee is the tick fo summing the matix elements fo the Stak effect fo the gound state of the hydogen atom Recall, we need to calculate the coection to the gound

More information

15 Solving the Laplace equation by Fourier method

15 Solving the Laplace equation by Fourier method 5 Solving the Laplace equation by Fouie method I aleady intoduced two o thee dimensional heat equation, when I deived it, ecall that it taes the fom u t = α 2 u + F, (5.) whee u: [0, ) D R, D R is the

More information

221B Lecture Notes Scattering Theory I

221B Lecture Notes Scattering Theory I Why Scatteing? B Lectue Notes Scatteing Theoy I Scatteing of paticles off taget has been one of the most impotant applications of quantum mechanics. It is pobably the most effective way to study the stuctue

More information

Chem 453/544 Fall /08/03. Exam #1 Solutions

Chem 453/544 Fall /08/03. Exam #1 Solutions Chem 453/544 Fall 3 /8/3 Exam # Solutions. ( points) Use the genealized compessibility diagam povided on the last page to estimate ove what ange of pessues A at oom tempeatue confoms to the ideal gas law

More information

On a generalization of Eulerian numbers

On a generalization of Eulerian numbers Notes on Numbe Theoy and Discete Mathematics Pint ISSN 1310 513, Online ISSN 367 875 Vol, 018, No 1, 16 DOI: 10756/nntdm018116- On a genealization of Euleian numbes Claudio Pita-Ruiz Facultad de Ingenieía,

More information

ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS

ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS L. MICU Hoia Hulubei National Institute fo Physics and Nuclea Engineeing, P.O. Box MG-6, RO-0775 Buchaest-Maguele, Romania, E-mail: lmicu@theoy.nipne.o (Received

More information

Physics 221 Lecture 41 Nonlinear Absorption and Refraction

Physics 221 Lecture 41 Nonlinear Absorption and Refraction Physics 221 Lectue 41 Nonlinea Absoption and Refaction Refeences Meye-Aendt, pp. 97-98. Boyd, Nonlinea Optics, 1.4 Yaiv, Optical Waves in Cystals, p. 22 (Table of cystal symmeties) 1. Intoductoy Remaks.

More information

Functions Defined on Fuzzy Real Numbers According to Zadeh s Extension

Functions Defined on Fuzzy Real Numbers According to Zadeh s Extension Intenational Mathematical Foum, 3, 2008, no. 16, 763-776 Functions Defined on Fuzzy Real Numbes Accoding to Zadeh s Extension Oma A. AbuAaqob, Nabil T. Shawagfeh and Oma A. AbuGhneim 1 Mathematics Depatment,

More information

On the Poisson Approximation to the Negative Hypergeometric Distribution

On the Poisson Approximation to the Negative Hypergeometric Distribution BULLETIN of the Malaysian Mathematical Sciences Society http://mathusmmy/bulletin Bull Malays Math Sci Soc (2) 34(2) (2011), 331 336 On the Poisson Appoximation to the Negative Hypegeometic Distibution

More information

Relating Branching Program Size and. Formula Size over the Full Binary Basis. FB Informatik, LS II, Univ. Dortmund, Dortmund, Germany

Relating Branching Program Size and. Formula Size over the Full Binary Basis. FB Informatik, LS II, Univ. Dortmund, Dortmund, Germany Relating Banching Pogam Size and omula Size ove the ull Binay Basis Matin Saueho y Ingo Wegene y Ralph Wechne z y B Infomatik, LS II, Univ. Dotmund, 44 Dotmund, Gemany z ankfut, Gemany sauehof/wegene@ls.cs.uni-dotmund.de

More information

Let us next consider how one can calculate one- electron elements of the long- range nuclear- electron potential:

Let us next consider how one can calculate one- electron elements of the long- range nuclear- electron potential: Hatee- Fock Long Range In the following assume the patition of 1/ into shot ange and long ance components 1 efc( / a) ef ( / a) + Hee a is a paamete detemining the patitioning between shot ange and long

More information

H.W.GOULD West Virginia University, Morgan town, West Virginia 26506

H.W.GOULD West Virginia University, Morgan town, West Virginia 26506 A F I B O N A C C I F O R M U L A OF LUCAS A N D ITS SUBSEQUENT M A N I F E S T A T I O N S A N D R E D I S C O V E R I E S H.W.GOULD West Viginia Univesity, Mogan town, West Viginia 26506 Almost eveyone

More information

The geometric construction of Ewald sphere and Bragg condition:

The geometric construction of Ewald sphere and Bragg condition: The geometic constuction of Ewald sphee and Bagg condition: The constuction of Ewald sphee must be done such that the Bagg condition is satisfied. This can be done as follows: i) Daw a wave vecto k in

More information