arxiv:gr-qc/ v1 26 May 2004

Size: px
Start display at page:

Download "arxiv:gr-qc/ v1 26 May 2004"

Transcription

1 Dirc-like pproch for consistent discretiztions of clssicl constrined theories Cyetno Di Brtolo 1, Rodolfo Gmbini 2, Rfel A. Porto 3 nd Jorge Pullin 4 1. Deprtmento de Físic, Universidd Simón Bolívr, Aptdo , Crcs 1080-A, Venezuel. 2. Instituto de Físic, Fcultd de Ciencis, Iguá 4225, esq. Mtojo, Montevideo, Uruguy. 3. Deprtment of Physics, Crnegie Mellon University, Pittsburgh, PA Deprtment of Physics nd Astronomy, Louisin Stte University, Bton Rouge, LA rxiv:gr-qc/ v1 26 My 2004 We nlyze the cnonicl tretment of clssicl constrined mechnicl systems formulted with discrete time. We prove tht under very generl conditions, it is possible to introduce nonsingulr cnonicl trnsformtions tht preserve the constrint surfce nd the Poisson or Dirc brcket structure. The conditions for the preservtion of the constrints re more stringent thn in the continuous cse nd s consequence some of the continuum constrints become second clss upon discretiztion nd need to be solved by fixing their ssocited Lgrnge multipliers. The guge invrince of the discrete theory is encoded in set of rbitrry functions tht pper in the generting function of the evolution equtions. The resulting scheme is generl enough to ccommodte the tretment of field theories on the lttice. This pper ttempts to clrify nd put on sounder footing discretiztion technique tht hs lredy been used to tret vriety of systems, including Yng Mills theories, BF-theory nd generl reltivity on the lttice. I. INTRODUCTION We hve recently introduced [1, 2] technique for treting the theories tht rise when one discretizes (spce)- time in constrined mechnicl system or continuum field theory. We hve shown tht this technique works for Yng Mills nd BF theories nd implemented it for the grvittionl cse. Previous ttempts to studying systems with discrete time hd concentrted mostly on systems without constrints or with holonomic constrints only (for review with comprehensive reference list see [3].) The ide consists on strting from discretized ction, constructing discrete Lgrnge equtions nd introducing symplectic structure in the discrete spce. The evolution is implemented vi cnonicl trnsformtions nd the consistency of the discrete theory determines in prt the Lgrnge multipliers. In some totlly constrined systems, like generl reltivity, the resulting discrete theories re constrint-free since the constrints re solved for the Lgrnge multipliers. This mkes the quntiztion of the discrete theories considerbly simpler thn the continuum cses. This ws exploited to mke progress in solving the problem of time in quntum grvity [4, 5] nd to implement the Pge Wootters reltionl time [5, 6] nd show tht fundmentl decoherence rises in quntum mechnics from quntum grvity. In this pper, we wnt to ddress in more systemtic wy the issue of the cnonicl formultion of discrete constrined systems. Up to now, most of the nlysis hs been mde on specific exmples, nd cnonicl nlysis, l Dirc[7, 8], is still lcking. In prticulr, the technique relied hevily on defining cnonicl trnsformtion tht ws initilly singulr (nd therefore not true cnonicl trnsformtion) nd showing tht one could eliminte vribles nd end up with true cnonicl trnsformtion. Up to now this ws shown in cse by cse bsis. This pper determines the generl conditions needed for the construction of proper cnonicl trnsformtion by following close nlogue of the Dirc procedure, dpted to the discrete cse. In prticulr, we note tht there re severl wys to proceed tht yield equivlent results, but tht my offer different dvntges for prticulr systems. In section II we will ly out the frmework of how to del with mechnicl systems where time is discrete, including singulr nd non-singulr systems. In section III we will develop clssifiction of constrints into first nd second clss suitble for the discrete context. In section IV we work out specific exmple tht exhibits the detils of the formlism. In section V we discuss n lterntive formultion of the formlism nd we end with conclusions nd discussion. II. MECHANICS WITH DISCRETE TIME We strt by considering continuum theory representing mechnicl system. Its Lgrngin will be denoted by ˆL(q, q ), = 1...M. This setting is generl enough to ccommodte, for instnce, totlly constrined systems. In

2 such cse q will be the derivtive of the cnonicl vribles with respect to the evolution prmeter. It is lso generl enough to include the systems tht result from formulting on discrete spce-time lttice continuum field theory. We discretize the evolution prmeter in intervls (possibly vrying upon evolution) t n+1 t n = ǫ n nd we lbel the generlized coordintes evluted t t n s q n. We define the discretized Lgrngin s 2 where The ction cn then be written s L(n,n+1) L(q n,q n+1 ) ǫ nˆl(q, q ) (1) q = q n nd q q n+1 q n ǫ n (2) S = N L(n,n+1) (3) n=1 If the continuum theory is invrint under reprmeteriztions of the evolution prmeter, one cn show tht the informtion bout the intervls ǫ n my be bsorbed in the Lgrnge multipliers. In the cse of stndrd mechnicl systems it is simpler to use n invrint intervl ǫ n = ǫ. The Lgrnge equtions of motion re obtined by requiring the ction to be sttionry under vritions of the configurtion vribles q fixed t the endpoints of the evolution intervl n = 0,n = N +1, L(n,n+1) q n + L(n 1,n) q n = 0. (4) These equtions define unique evolution if the determinnt 2 L(n,n+1) +1 b q 0. (5) n We will refer to this cse s the nonsingulr cse. When the determinnt vnishes, one will hve to nlyze the sitution differently. Let us strt with the nonsingulr cse. A. Nonsingulr cse In this cse one cn solve the Lgrnge equtions explicitly nd the q n+1 re uniquely given s function of q n nd q n 1. This is the equivlent of the Hessin condition for the non-singulr Lgrngin theories in the continuum. The resulting equtions re second order in the sense tht the q n s re determined provided one knows two previous time levels. One cn introduce first order formultion by introducing cnoniclly conjugte vribles s is usully done when introducing Hmiltonin formultion in the continuum theories. We introduce the following definition of vribles tht we will lter show end up being cnoniclly conjugte moment of the configurtion vribles, p n+1 L(n,n+1) +1 p n L(n 1,n) = L(n,n+1) q n (6) (7) Where we hve used Eq. (4). The equtions (6) nd (7) define cnonicl trnsformtion for the vribles q n,p n to q n+1,p n+1 with the type 1 generting function F 1 = L(qn,q n +1) provided tht condition (5) is fulfilled. Notice tht the evolution scheme is implicit, one cn use the bottom eqution (since we re in the non-singulr cse) to give n expression for q n+1 in terms of q n,p n, which in turn cn be substituted in the top eqution to get n eqution for p n+1 purely in terms of q n,p n. It should be noted tht there re severl other possible choices, when going from the set of equtions (6,7) to n explicit evolution scheme. For exmple, one cn choose to do things in wy tht yields closer nlogy with the

3 stndrd Hmiltonin description in the continuum by introducing type two cnonicl trnsformtions. To do this, we choose to invert eqution (6) for qn+1 q n+1 (qb n,pb n+1 ), which is possible only if 2 L(n,n+1) +1 qb 0. (8) n+1 We cn now introduce Legendre trnsform nd define 3 F 2 (q n,p n+1 ) q n+1 p n+1 L(q n,q n+1 ). (9) From here it is immedite to obtin q b n+1 = F 2 p b n+1 F 2 q b n = L(n,n+1) q b n (10) = p b n, (11) where we hve used the Lgrnge eqution in order to obtin the lst equlity. We esily recognize from here tht F 2 behves s type two generting function of the cnonicl trnsformtion connecting level n with level n+1. We my define now sort of type 2 Hmiltonin (in the sense tht it depends on p n+1 nd q n, H 2 (p n+1,q n ) given by, H 2 (q n,p n+1 ) p n+1(q n+1 q n) L(q n,q n+1 ) = F 2 (q n,p n+1 ) p n+1q n. (12) which leds to the discrete Hmilton equtions, q b n+1 = qb n + H 2 p b n+1 (13) p b n = pb n+1 + H 2 b. (14) It should be noted tht lthough this formultion hs degree of nlogy with the trditionl Hmiltonin formultion, there re significnt differences due to the fct tht the conjugte vribles live t different time slices. It would not be possible therefore to use this formultion to ttempt to construct Schrödinger eqution strting from the bove Hmiltonin. Provided tht the cnonicl mp defined by F 2 is invertible we end up with discrete evolution implemented by cnonicl trnsformtion. It cn be esily seen by using the Legendre trnsform tht F 2 will be invertible if nd only if 2 F 2 (n,n+1) q np b n+1 = 2 L(n,n+1) q n+1 b 2 L(n,n+1) +1 qb n (15) Thus, in order to hve cnonicl trnsformtions generted by type two functions, the Hessin condition in the continuum time mechnics leds to two independent conditions in the discrete theory given by Eqs. (5) nd (8). Notice however tht (8) is not necessry for introducing symplectic structure. It is cler tht when one builds cnonicl discrete theory there re four possibilities depending on which pir of vribles one chooses to construct the generting functionl of the cnonicl trnsformtion, either q n,q n+1, q n,p n+1, p n,q n+1, p n,p n+1. In this subsection we considered only two, but the others cn be esily generlized from the discussion here. B. The singulr cse Let us consider s before the generic discrete Lgrngin L(n,n+1) = L(qn,q n+1 ) with = 1...M. It leds to the equtions we lredy discussed, p n+1 = L(n,n+1) +1 p n = L(n,n+1) (16) (17)

4 If 2 L(n,n+1) +1 qb vnishes nd its rnk is K, the system is singulr nd hs M K constrints of the form n 4 Φ A (q n,p n ) = 0 (18) tht result from Eq(17), nd M K constrints Ψ A (q n+1,p n+1 ) = 0 (19) resulting from Eq(16). The evolution of the configurtion vribles from level n to n+1 is given by solving for q n+1 the Eqs. (17). As the system is singulr the evolution depends of (M K) rbitrry functions V A q n+1 = f (q b n,p b n,v A ). (20) We shll follow closely the stndrd Dirc cnonicl procedure of continuum mechnics. In this cse, the nlysis of constrint system goes trough two steps. The first step consists in the definition of set of evolution equtions tht wekly preserve the constrints nd the Poisson symplectic structure. To do tht one defines the totl Hmiltonin H T = H 0 +V α φ α where φ α re the primry constrints, nd the V α re prtilly determined in order to preserve ll the constrints of the system. Even though some of the V s my be rbitrry functions, once they re specified, the evolution generted by H T preserves the Poisson brckets nd the dynmicl evolution is consistent t the clssicl level with the constrint structure. The second step is only required to quntize the system nd consists in the identifiction of the first nd second clss constrints nd the introduction of the Dirc brckets tht enforce strongly the second clss constrints. As we shll see, the sme procedure my be followed in the discrete cse. The min difference is the implementtion of the cnonicl trnsformtion tht is not generted by Hmiltonin but by cnonicl trnsformtion of type 2, 3, or 4. Let us strt by completing the evolution equtions. We need to dd to Eq(20) n eqution for p n+1 p n+1 = L(q n,q n+1 ) q n+1 We now impose the preservtion of the constrints Φ A (q n+1,p n+1 ) = Φ A q n+1 =f (q b n,pb n,v A ). (21) ( f n, L(q n,f ) ) +1 = 0. (22) Furthermore, we need to impose the n + 1 level constrints t level n, Ψ A (qn,p n) = 0 nd impose the consistency conditions ( Ψ A f n, L(q n,f ) ) +1 = 0. (23) Three different cses my occur: ) Eqs.(22,23) vnish utomticlly, nd therefore we re not led to new conditions. b) They led to inconsistencies, nd the dynmicl system is inconsistent. c) New secondry constrints C(q n,p n) pper or/nd some of the rbitrry functions V A re determined, tht is V A = V A (q n,p n,vα ) with α = 1...R (N K), nd v α rbitrry functions. The process is repeted until consistency is chieved. Tht is, until the consistency conditions re utomticlly stisfied without further constrints nd conditions for V. SubstitutingV A in(20)nd(21)wegettheevolutionequtionsthtpreservelltheconstrints: primry,secondry, tertiry, nd so on, q n+1 = f (q b n,p b n,v A (q,p,v)) = f (q b n,p b n,v α ) (24) p n+1 = g (q b n,pb n,vα ) (25)

5 5 Initil vlues need to be restricted by the n = 0 level constrints, Φ A (q0,p 0 ) = Ψ A(q0,p 0 ) = 0. (26) In order to hve complete nlogy with the continuum cse, we still need to nlyze under wht conditions this evolution lso preserves the Poisson brcket structure. As in the continuum cse we ssume tht the rbitrry functions v hve been fixed. Three different cses my rise depending on if one chooses to implement things in terms of cnonicl trnsformtion of type 2, 3 or 4: Cse I) The eqution(24) is invertible for qn: tht is, q n+1 b 0 nd therefore one cn write qn = h (qn+1 b,pb n ). (27) Notice tht, under these hypotheses, there re no pseudo-constrints of the form G(q n+1,p n) = 0. We cll these pseudo-constrints becuse they involve vribles t different instnts of time. We my define type 3 generting functions of cnonicl trnsformtions, Then we hve F 3 (qn+1 b,pb n ) = [pb n qb n +L(q n,q n+1 )] q b n =hb. (28) F 3 (q n+1,p n) p b n = q b n = h b (q n+1,p n), (29) nd F 3 (q n+1,p n ) q b n+1 = L(q n,q n+1 ) q b n+1 q n=h = p b n+1. (30) Notice tht in the lst equlity there re lso contributions coming from the dependence on q n+1 of the level n vriblesq n, but these contributionscncelbecuse ofthe definition ofthe cnoniclmoment. The informtionbout the moment is completely encoded in the evolution equtions qn+1 = h nd the constrints. As the first eqution is equivlent to (24) one ends up recovering the fundmentl evolution equtions s cnonicl trnsformtion generted by F 3. Furthermore, F 3 (qn+1,p n ) p b n q = b n (31) due to the fct tht we hve ssumed tht q b n+1 0, nd consequently F 3 is non singulr generting function nd therefore the resulting cnonicl trnsformtion preserves the Poisson brcket structure. Cse II) The eqution(25) is invertible for p n : tht is p n+1 p b 0 nd therefore one cn write, n p n = g (q b n,p b n+1). (32) Notice tht, under these hypotheses, there re no pseudo-constrints of the form G(qn,p n+1 ) = 0. By substituting (32) in (24), one gets which llows to introduce type 2 generting function qn+1 = k (qn b,pb n+1 ), (33) F 2 (q b n,p b n+1) = [p b n+1q b n+1 L(q n,q n+1 )] q b n+1 =k b (34) One cn now esily check tht this generting function reproduces the evolution equtions (24) nd (25) nd defines non singulr cnonicl trnsformtion tht preserves the Poisson Brckets in the evolution. Cse III) Even when the system hs pseudo-constrints of the form G(qn+1,p n ) = 0 nd F(q n,p n+1 ) = 0 one my be ble to find cnonicl trnsformtion provided tht the system does not hve pseudo-constrints of the form F(p n,p n+1 ) = 0

6 6 In fct, by using Eq.(25) one cn invert for qn = l (p b n,pb n+1 ), (35) nd substituting in (24) get qn+1 = m (p b n,pb n+1 ). (36) A generting function of type 4 tht does the sme job tht the two previous ones my now be introduced, q F 4 (p b n,pb n+1 ) = [ pb n+1 qb n+1 +pb n qb n +L(q b n+1 n,q n+1 )] =mb (37) qn b=lb All of the discrete systems tht hve been treted up to now in the literture my be nlyzed by following this cnonicl procedure, llowing to preserve the constrints nd the Poisson brcket structure. Lter on we will show n exmple of system of this type in order to nlyze how this procedure works in concrete cse. It should be noted tht there my exist mechnicl systems tht do not fll into ny of the bove clssifictions. For instnce, system could hve pseudo-constrints of ll the types listed bove. In such cses one will need to develop further techniques to tret them. For instnce one could introduce cnonicl trnsformtions of given type for some of the vribles nd of different type for other vribles. This would require further study nd it does not pper necessry for the systems tht hve been nlyzed up to present. III. CLASSIFICATION OF THE CONSTRAINTS At this point we hve set of constrints primry, secondry, tertiry, etc, of the form χ Z (q,p ) = 0 with Z = 1...A with A the totl number of independent constrints, tht re preserved under the evolution given by the equtions(24,25) provided prt or ll the rbitrry functions V re conveniently fixed. As in the continuum cse it is convenient to introduce the notion of first nd second clss constrints, in order to quntize the theory. A constrint is of first clss if it commutes with ll the constrints, if tht is not the cse it is of second clss. As in the continuum cse one cn define first clss functions of the cnonicl vribles f(q,p) tht re not necessrily constrints. Such function will be first clss if it commutes with ll the constrints. Second clss constrints my be imposed strongly by introducing Dirc brckets. As the evolution equtions preserve the Poisson structure, they will preserve de Dirc structure becuse Dirc brckets re defined in terms of Poisson brckets. One ends up with theory with set of evolution equtions tht preserve the symplectic structure of the system, nd therefore my be quntized by describing the evolution in terms of unitry opertors. In the discrete cse there is not strightforwrd reltion between the number of first nd second clss constrints nd the number of phse spce degrees of freedom. This is due to the fct tht now the evolution of the constrints is not directly relted with their Poisson brckets with totl Hmiltonin. Thus, the fct tht constrint does not commute with others is not esily relted with the determintion of n rbitrry function. It is still very esy to determine the number of phse spce degrees of freedom. In fct this number is given by two times the number of configurtion vribles minus the totl number of constrints minus the number of rbitrry functions v. IV. AN EXAMPLE OF CONSTRAINED SYSTEM WITH SECOND CLASS CONSTRAINTS To illustrte the techniques elborted bove, we would like to discuss model tht is simple, yet ddresses in non-trivil wy the min points we discussed. This exmple hd been treted using d-hoc techniques in [4]. The model consists of prmeterized free prticle in two dimensionl spce-time under the influence of liner potentil. The discrete Lgrngin is given by, L n L(q n,π n,n n,q n+1,π n+1,n n+1 ) = π n(q n+1 q n) N n [π 0 n (π1 n) 2 +αq 1 n]. (38) We hve chosen first order formultion for the prticle. However, this Lgrnginis of the type we considered in this pper, one simply needs to consider ll vribles, q,π,n s configurtion vribles. The system is clerly singulr since the π s nd N only pper t level n (or in the continuum Lgrngin, their time derivtives re bsent). When

7 7 considered s Type I generting function, the bove Lgrngin leds to the equtions nd p π,n = L n π n p q,n = L n q n p π,n+1 = L n π n+1 = 0, (39) p q,n+1 = L n +1 = πn, (40) p N,n+1 = L n N n+1 = 0, (41) = (q n+1 q n)+π 1 nn n δ 1 +N n δ 0, (42) = π n +δ 1 αn n, (43) p N,n = L n N n = π 0 n (π1 n) 2 +αq 1 n. (44) One cn esily recognize tht the system hs six constrints: three t the n+1 level, nd three t the n level. They re: Therefore the evolution depends on three rbitrry functions V N,n,V π,n, ψ1 p π,n+1 = 0, (45) ψ 2 p N,n+1 = 0, (46) Φ 1 p q,n (πn +δ1αn n ), (47) Φ 2 p N,n [π 0 n (π1 n) 2 +αq 1 n]. (48) qn+1 = qn +πnn 1 n δ1 +N n δ0 p π,n (49) πn+1 = π n +V π,n (50) N n+1 = N n +V N,n. (51) The preservtion of the ψ constrints from level n to level n+1 is utomticlly ensured from (23). Now we impose the preservtion of the Φ constrints upon evolution. Let us begin with Φ 0 1 : which tking into ccount the constrint Φ 0 1n implies V 0 π,n = 0. For the Φ 2 one gets the eqution tht tking into ccount the constrint Φ 2n implies tht where ǫ = ±1. Finlly we hve tht fter imposing the constrint t level n leds to Φ 0 1n+1 p0 q,n+1 π0 n+1 = p0 q,n π0 n V0 π,n = 0, (52) Φ 2n+1 = πn 0 +α(qn 1 +πnn 1 n )+(πn 1 +Vn π1 ) 2 /2 = 0, (53) V 1 π,n = π 1 n +ǫ (π 1 n) 2 2π 1 nn n α (54) Φ 1 1n+1 = p1 q,n N nα π 1 n V1 π,n α(n n +V N,n ) = 0, (55) V N,n = 1 α V 1 π,n N n. (56)

8 8 Thus, the evolution equtions for the configurtion vribles re q n+1 = q n +π1 n N nδ 1 +N nδ 0 (57) π 0 n+1 = π 0 n, (58) π 1 n+1 = ǫ (π 1 n )2 2π 1 n N nα, (59) N n+1 = 1 α [π1 n ǫ (π 1 n) 2 2π 1 nn n α]. (60) We re now redy to define n invertible cnonicl trnsformtion with the help of type 3 generting function. Notice tht these evolution equtions re invertible for q n, nd therefore we re in the cse I. The inverse is given by set of equtions of the form q n (q n+1 ), explicitly given by, q 1 n = q1 n+1 π1 n+1 N n+1 α 2 (N n+1) 2, (61) π qn 0 = qn+1 0 n+1 1 N n+1 + α 2 (N n+1) 2 πn+1 1 +αn, (62) n+1 π 0 n = π 0 n+1, (63) π 1 n = π1 n+1 +αn n+1, (64) N n = Reclling tht F 3 is given by eqution (28) we obtin, π 1 n+1n n+1 + α 2 (N n+1) 2 π 1 n+1 +αn n+1. (65) F 3 = p 0 π,n π0 n+1 +p1 π,n (αn n+1 +πn+1 1 )+p0 q,n q0 n N [ n+1 α 2 Nn αN n+1πn (π1 n+1 )2] N n+1(αn n+1 +2π 1 n+1 ) 2(αN n+1 +π 1 n+1 ) [ p 0 q,n p N,n (π1 n+1 )2 +αq 1 n+1 ] ( +p 1 q,n qn+1 1 α ) 2 N2 n+1 N n+1πn+1 1. (66) One cn check tht Eqs. (61) to (65) re esily recovered by tking the prtil derivtive with respect to p n. By differentiting with respect to qn+1 one gets p 0 q,n+1 = p 0 q,n, p 1 q,n+1 = p 1 q,n α 2 A nn n+1 (αn n+1 +2π 1 n+1), (68) (67) p 0 π,n+1 = p0 π,n (69) p 1 π,n+1 = p1 π,n + N n+1 2 [ 3αNn+1 2p 1 q,n +4π1 n+1 +2B na n A n (π 1 n+1 +B na n )(αn n+1 +2π 1 n+1 )], (70) p N,n+1 = 3 2 α2 Nn+1 2 +αp1 π,n p1 q,n π1 n+1 +(π1 n+1 )2 +B n +αn n+1 [ p 1 q,n +3πn ] 2 B na 2 n(αn n+1 +2πn+1) 1, (71) where we hve introduced A n (αn n+1 +π 1 n+1 ) 1, (72) B n p N,n p 0 q,n 1 2 (π1 n+1 )2 αq 1 n+1. (73) By substituting in Eqs.(67) to (71) the vribles q n+1 nd using the level n constrints one gets Eqs. (39) to (41), which is cnonicl trnsformtion tht reproduces, on the constrint surfce, the evolution equtions of the discrete prticle. Wht remins to be done is to identify the second clss constrints nd impose them strongly. The complete set of six constrints of this model ψ nd Φ re second clss nd llow to solve for π nd N nd eliminte completely these vribles nd their complex conjugtes P π nd P N.

9 One cn proceed in two different wys. The first lterntive is to strt by observing tht P π nd P N vnish strongly, nd then to solve for π nd N in terms of the n+1 level vribles, leding to π 0 n = p 0 q,n+1, (74) π 1 n = p 1 q,n+1, (75) N n = C n+1 αp 1, (76) q,n+1 where C n+1 = p 0 q,n+1 +(p1 q,n+1 )2 /2+αqn+1 1. The relevnt evolution equtions re obtined from (57,67,68),nd re given by qn+1 0 = qn 0 +N n = qn 0 + C n+1 αp 1, (77) q,n+1 qn+1 1 = q1 n +N nπn 1 = q1 n + C n+1 α, (78) p 0 q,n = p0 q,n+1, (79) p 1 q,n = p1 q,n+1 +αn n = p 1 q,n+1 + C n+1 p 1, (80) q,n+1 nd we recover the evolution equtions obtined in [4]. The second lterntive consists in solving for π nd N in terms of the n level vribles, leding to π 0 n = p0 q,n, (81) πn 2(p 1 = ǫ 0 q,n +αq1 n ), (82) N n = 1 ( ) p 1 q,n ǫ 2(p α 0 q,n +αqn) 1, (83) nd from here, computing the evolution equtions by using (57,67,68). The two methods yield evolution schemes of different functionl form since one propgtes forwrd in time nd the other bckwrd. The inequivlence in the functionl form stems from the fct tht the discretiztion of the time derivtives chosen in the Lgrngin is not centered. It should be emphsized tht if one strts from given initil dt nd propgtes forwrds with the first system of equtions nd then bckwrds using the second, one will return to the sme initil dt. Notice tht we hve 6 second clss constrints, nd the initil number of phse spce degrees of freedom ws 10. By noticing tht there re no rbitrry functions left, one is left with 4 degrees of freedom on the constrint surfce. The continuum model hd two degrees of freedom. The procedure we hve followed here is completely generl nd my be simplified when one is treting specific cses. For instnce, s it hppens in the continuum theory [8] it is sometimes possible to implement the cnonicl nlysis by first solving the constrints for the unphysicl degrees of freedom N,π,p N,p π nd then introducing generting functionl on the physicl degrees of freedom by following the procedure of the previous sections. In this prticulr cse it is esy to show tht, for p q1 0, q n+1 b 0, where q (q 1,q 0 ), nd thus it is possible to construct n F 3 generting functionl. In ll the models treted up to now in the literture the unphysicl degrees of freedom were eliminted before obtining the cnonicl trnsformtion for the evolution of the physicl degrees of freedom. In order to keep the nlysis generl in simple model, here we hve kept ll the vribles involved in this pproch. 9 V. TREATMENT IN TERMS OF TYPE II GENERATING FUNCTIONS Up to now, we hve tken s the strting point singulr Type I generting function given by F 1 = L(q n,q n+1). It is interesting to nlyze how singulr systems my be described in terms of other types of generting functions. Let us ssume tht neither condition (8) nor condition (5) re fulfilled by the discrete system nd, therefore, we re in the singulr cse. Then, it my be immeditely seen, by using (16),nd tking into ccount tht the determinnt tht ppers in (8) now vnishes, tht singulr systems with first order Lgrngins in the continuum led to the presence of pseudo-constrints φ α (q n,p n+1 ) = 0, in the theory. Recll tht these pseudo-constrints re similr to the constrints tht rise in the continuum cse, but mixing configurtion vribles t one level with moment in the

10 next level. As in the continuum cse, one cn then introduce the pseudo-constrint surfce S φ which we shll ssume hs well defined functions, ( φ α =... φ ) φ...,......, (84) where vector τ in the tngent spce of S φ, T (q,p) S φ, is such tht Let us introduce now s before the type-2 Hmiltonin (12), p n+1 τ φ α = 0 α. (85) 10 H 2 p n+1 (q n+1 q n ) L(q n,q n+1 ),. (86) For which it is esy to see, using the pseudo-constrints, tht it is function of q n,p n+1. Let us consider now n infinitesiml vrition of H 2 long the pseudo-constrint surfce, [ ( (qn+1 qn)dp n+1 + p n+1 L ) ( dqn+1 p n+1 + L ) ] dqn dh 2 = = q n+1 [ ( (qn+1 qn)dp n+1 p n+1 + L ) ] dqn, (87) where in the lst step we used (6). H 2 is well defined in S φ but cn be extended to the whole phse spce s it is done in the cse of constrints in the continuum theory. In order to obtin the cnonicl eqution of motion we strt from the identity, dh 2 (q n,p n+1 ) = ( H2 q dqn + H ) 2 n p dp n+1, (88) n+1 nd evluting it for n infinitesiml displcement (dq,dp) in T (q,p) S φ, we cn use (87) nd the Lgrngin equtions of motion (7) to obtin, [( p n+1 p n + H ( 2 )dq n + qn q n+1 + H ) ] 2 p dp n+1 = 0. (89) n+1 Since (dq,dp) is n rbitrrytngent vector to S φ, following (85) we obtin tht the coefficient must be proportionl to the grdient. Introducing therefore the Lgrnge multipliers λ s the proportionlity fctors, we will end with the set of equtions, q n+1 = q n + H 2 p n+1 p n+1 = p n H 2 q n φ α +λ α n p n+1 λ α n q n, (90) φ α, (91) φ α (q n,p n+1 ) = 0. (92) This system is similr to the continuum set of equtions with the inclusion of pseudo-constrints, which hve the sme functionl form s the continuum constrints, but involve configurtion vribles t level n nd moment t level n + 1, insted of n level vribles. From this initil evolution equtions one my follow procedure similr to the one developed in the previous section in order to study the consistency of the (primry) constrints nd pseudo-constrints. VI. DISCUSSION AND CONCLUSIONS We hve provided cnonicl procedure for the introduction of preserved symplectic structure in discrete constrined systems. The nlogy with Dirc s procedure in the continuum is quite remrkble. It is possible to define notion of discrete evolution tht wekly preserves constrints nd Poisson brckets. The distinction between first

11 nd second clss constrints is still useful nd when second clss constrints re imposed strongly the resulting Dirc brckets re preserved. A feture of the discretized theories is tht they my hve smller number of first clss constrints, nd consequently more degrees of freedom thn the continuum counterprts. The extr degrees of freedom come from the fct tht the discrete theories my not necessrily hve the sme symmetries s the continuum theories. For instnce, in the cse of homogeneous cosmologies studied in [9] the extr pir of phse spce degrees of freedom re ssocited with the fct tht in the discrete theory different choices of refinements in the discretiztion in time correspond to different solutions in the discrete theory tht nevertheless pproximte the sme solution in the continuum theory. This is only first step for complete understnding of the dynmics of discrete guge systems. The reltion between the number of constrints of first nd second clss nd the number of degrees of freedom, nd the connection between the first clss constrints nd the guge invrince of the discrete dynmicl system need to be further studied. Moreover, s discussed in the body of the rticle, if one wishes to consider more pthologicl systems thn the ones considered here, more elborte cnonicl trnsformtions my need to be introduced. The issue of the continuum limit is well understood in the non singulr cse, where there is n externl step prmeter tht controls the pproximtion. However, it needs further study in the cse of singulr systems, prticulrly in the cse of totlly constrined systems where the step of the pproximtion is encoded in the dditionl degrees of freedom of the discrete theory. This issue hs been studied in severl models[4, 9] but complete chrcteriztion of the possible behviors is still lcking. A similr comment pplies to the role of sptil discretiztions when one is considering lttice field theories. 11 VII. ACKNOWLEDGMENTS This pper ws motivted by series of questions by Krel Kuchř bout our previous work. This work ws supported by grnt nsf-phy nd funds from the Horce Herne Jr. Lbortory for Theoreticl Physics nd the Abdus Slm Interntionl Center for Theoreticl Physics. [1] C. Di Brtolo, R. Gmbini, J. Pullin, Clss. Qun. Grv. 19, 5475 (2002). [2] R. Gmbini nd J. Pullin, Phys. Rev. Lett. 90, , (2003). [3] J. Mrsden, M. West, Act Numeric 10, 357 (2001). [4] R. Gmbini, R.A. Porto nd J. Pullin, in Recent developments in grvity, K. Kokkots, N. Stergiouls, World Scientific, Singpore, (2003) [rxiv:gr-qc/ ]. [5] R. Gmbini, R. A. Porto nd J. Pullin, New Journl of Physics 6, 45 (2004). [6] R. Gmbini, R. A. Porto nd J. Pullin, Clss. Qun. Grv. 21 L51 (2004). [7] P. A. M. Dirc, Lectures on Quntum Mechnics, Dover (2001) [8] M. Henneux nd C. Teitelboim, Quntiztion of Guge Systems, Princeton University Press (1994). [9] R. Gmbini nd J. Pullin, Clss. Qun. Grv. 20, 3341 (2003).

Week 10: Line Integrals

Week 10: Line Integrals Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.

More information

221B Lecture Notes WKB Method

221B Lecture Notes WKB Method Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using

More information

EXTENDED BRST SYMMETRIES IN THE GAUGE FIELD THEORY

EXTENDED BRST SYMMETRIES IN THE GAUGE FIELD THEORY Romnin Reports in hysics olume 53 Nos.. 9 4 EXTENDED BRST SYETRIES IN TE GUGE FIELD TEORY UREL BBLEN RDU CONSTNTINESCU CREN IONESCU Deprtment of Theoreticl hysics University of Criov 3.I. Cuz Criov Romni

More information

221A Lecture Notes WKB Method

221A Lecture Notes WKB Method A Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using ψ x, t = e

More information

Symmetries of Dynamically Equivalent Theories

Symmetries of Dynamically Equivalent Theories 13 Brzilin Journl of Physics, vol. 36, no. 1B, Mrch, 006 Symmetries of Dynmiclly Equivlent Theories D. M. Gitmn nd I. V. Tyutin Institute of Physics, University of São Pulo, Brzil nd Lebedev Physics Institute,

More information

MAC-solutions of the nonexistent solutions of mathematical physics

MAC-solutions of the nonexistent solutions of mathematical physics Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE

More information

1.9 C 2 inner variations

1.9 C 2 inner variations 46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for

More information

1 1D heat and wave equations on a finite interval

1 1D heat and wave equations on a finite interval 1 1D het nd wve equtions on finite intervl In this section we consider generl method of seprtion of vribles nd its pplictions to solving het eqution nd wve eqution on finite intervl ( 1, 2. Since by trnsltion

More information

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004 Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when

More information

Here we study square linear systems and properties of their coefficient matrices as they relate to the solution set of the linear system.

Here we study square linear systems and properties of their coefficient matrices as they relate to the solution set of the linear system. Section 24 Nonsingulr Liner Systems Here we study squre liner systems nd properties of their coefficient mtrices s they relte to the solution set of the liner system Let A be n n Then we know from previous

More information

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011 Clssicl Mechnics From Moleculr to Con/nuum Physics I WS 11/12 Emilino Ippoli/ October, 2011 Wednesdy, October 12, 2011 Review Mthemtics... Physics Bsic thermodynmics Temperture, idel gs, kinetic gs theory,

More information

PHYS 4390: GENERAL RELATIVITY LECTURE 6: TENSOR CALCULUS

PHYS 4390: GENERAL RELATIVITY LECTURE 6: TENSOR CALCULUS PHYS 4390: GENERAL RELATIVITY LECTURE 6: TENSOR CALCULUS To strt on tensor clculus, we need to define differentition on mnifold.a good question to sk is if the prtil derivtive of tensor tensor on mnifold?

More information

7.2 The Definite Integral

7.2 The Definite Integral 7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where

More information

NOTES ON HILBERT SPACE

NOTES ON HILBERT SPACE NOTES ON HILBERT SPACE 1 DEFINITION: by Prof C-I Tn Deprtment of Physics Brown University A Hilbert spce is n inner product spce which, s metric spce, is complete We will not present n exhustive mthemticl

More information

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties

More information

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz

More information

Continuous Quantum Systems

Continuous Quantum Systems Chpter 8 Continuous Quntum Systems 8.1 The wvefunction So fr, we hve been tlking bout finite dimensionl Hilbert spces: if our system hs k qubits, then our Hilbert spce hs n dimensions, nd is equivlent

More information

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1 Tor Kjellsson Stockholm University Chpter 5 5. Strting with the following informtion: R = m r + m r m + m, r = r r we wnt to derive: µ = m m m + m r = R + µ m r, r = R µ m r 3 = µ m R + r, = µ m R r. 4

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct

More information

The Basic Functional 2 1

The Basic Functional 2 1 2 The Bsic Functionl 2 1 Chpter 2: THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 2.1 Introduction..................... 2 3 2.2 The First Vrition.................. 2 3 2.3 The Euler Eqution..................

More information

Recitation 3: More Applications of the Derivative

Recitation 3: More Applications of the Derivative Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech

More information

Numerical integration

Numerical integration 2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter

More information

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah 1. Born-Oppenheimer pprox.- energy surfces 2. Men-field (Hrtree-Fock) theory- orbitls 3. Pros nd cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usully does HF-how? 6. Bsis sets nd nottions 7. MPn,

More information

1.3 The Lemma of DuBois-Reymond

1.3 The Lemma of DuBois-Reymond 28 CHAPTER 1. INDIRECT METHODS 1.3 The Lemm of DuBois-Reymond We needed extr regulrity to integrte by prts nd obtin the Euler- Lgrnge eqution. The following result shows tht, t lest sometimes, the extr

More information

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique? XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk bout solving systems of liner equtions. These re problems tht give couple of equtions with couple of unknowns, like: 6 2 3 7 4

More information

Variational problems of some second order Lagrangians given by Pfaff forms

Variational problems of some second order Lagrangians given by Pfaff forms Vritionl problems of some second order Lgrngins given by Pfff forms P. Popescu M. Popescu Abstrct. In this pper we study the dynmics of some second order Lgrngins tht come from Pfff forms i.e. differentil

More information

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a). The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples

More information

Bernoulli Numbers Jeff Morton

Bernoulli Numbers Jeff Morton Bernoulli Numbers Jeff Morton. We re interested in the opertor e t k d k t k, which is to sy k tk. Applying this to some function f E to get e t f d k k tk d k f f + d k k tk dk f, we note tht since f

More information

Recitation 3: Applications of the Derivative. 1 Higher-Order Derivatives and their Applications

Recitation 3: Applications of the Derivative. 1 Higher-Order Derivatives and their Applications Mth 1c TA: Pdric Brtlett Recittion 3: Applictions of the Derivtive Week 3 Cltech 013 1 Higher-Order Derivtives nd their Applictions Another thing we could wnt to do with the derivtive, motivted by wht

More information

4. Calculus of Variations

4. Calculus of Variations 4. Clculus of Vritions Introduction - Typicl Problems The clculus of vritions generlises the theory of mxim nd minim. Exmple (): Shortest distnce between two points. On given surfce (e.g. plne), nd the

More information

Stuff You Need to Know From Calculus

Stuff You Need to Know From Calculus Stuff You Need to Know From Clculus For the first time in the semester, the stuff we re doing is finlly going to look like clculus (with vector slnt, of course). This mens tht in order to succeed, you

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS. THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem

More information

The Regulated and Riemann Integrals

The Regulated and Riemann Integrals Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue

More information

Separation of Variables in Linear PDE

Separation of Variables in Linear PDE Seprtion of Vribles in Liner PDE Now we pply the theory of Hilbert spces to liner differentil equtions with prtil derivtives (PDE). We strt with prticulr exmple, the one-dimensionl (1D) wve eqution 2 u

More information

SYDE 112, LECTURES 3 & 4: The Fundamental Theorem of Calculus

SYDE 112, LECTURES 3 & 4: The Fundamental Theorem of Calculus SYDE 112, LECTURES & 4: The Fundmentl Theorem of Clculus So fr we hve introduced two new concepts in this course: ntidifferentition nd Riemnn sums. It turns out tht these quntities re relted, but it is

More information

Name Solutions to Test 3 November 8, 2017

Name Solutions to Test 3 November 8, 2017 Nme Solutions to Test 3 November 8, 07 This test consists of three prts. Plese note tht in prts II nd III, you cn skip one question of those offered. Some possibly useful formuls cn be found below. Brrier

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

Math Lecture 23

Math Lecture 23 Mth 8 - Lecture 3 Dyln Zwick Fll 3 In our lst lecture we delt with solutions to the system: x = Ax where A is n n n mtrix with n distinct eigenvlues. As promised, tody we will del with the question of

More information

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018 Physics 201 Lb 3: Mesurement of Erth s locl grvittionl field I Dt Acquisition nd Preliminry Anlysis Dr. Timothy C. Blck Summer I, 2018 Theoreticl Discussion Grvity is one of the four known fundmentl forces.

More information

Chapter 3 Polynomials

Chapter 3 Polynomials Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling

More information

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008 Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of

More information

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1.2. Linear Variable Coefficient Equations. y + b ! = a y + b  Remark: The case b = 0 and a non-constant can be solved with the same idea as above. 1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt

More information

Time Optimal Control of the Brockett Integrator

Time Optimal Control of the Brockett Integrator Milno (Itly) August 8 - September, 011 Time Optiml Control of the Brockett Integrtor S. Sinh Deprtment of Mthemtics, IIT Bomby, Mumbi, Indi (emil : sunnysphs4891@gmil.com) Abstrct: The Brockett integrtor

More information

arxiv:gr-qc/ v1 14 Mar 2000

arxiv:gr-qc/ v1 14 Mar 2000 The binry blck-hole dynmics t the third post-newtonin order in the orbitl motion Piotr Jrnowski Institute of Theoreticl Physics, University of Bi lystok, Lipow 1, 15-2 Bi lystok, Polnd Gerhrd Schäfer Theoretisch-Physiklisches

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon Phys463.nb 49 7 Energy Bnds Ref: textbook, Chpter 7 Q: Why re there insultors nd conductors? Q: Wht will hppen when n electron moves in crystl? In the previous chpter, we discussed free electron gses,

More information

Calculus of Variations

Calculus of Variations Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function

More information

Lecture XVII. Vector functions, vector and scalar fields Definition 1 A vector-valued function is a map associating vectors to real numbers, that is

Lecture XVII. Vector functions, vector and scalar fields Definition 1 A vector-valued function is a map associating vectors to real numbers, that is Lecture XVII Abstrct We introduce the concepts of vector functions, sclr nd vector fields nd stress their relevnce in pplied sciences. We study curves in three-dimensionl Eucliden spce nd introduce the

More information

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 UNIFORM CONVERGENCE Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 Suppose f n : Ω R or f n : Ω C is sequence of rel or complex functions, nd f n f s n in some sense. Furthermore,

More information

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below. Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we

More information

Math 124A October 04, 2011

Math 124A October 04, 2011 Mth 4A October 04, 0 Viktor Grigoryn 4 Vibrtions nd het flow In this lecture we will derive the wve nd het equtions from physicl principles. These re second order constnt coefficient liner PEs, which model

More information

arxiv:hep-th/ v1 18 Jan 2002

arxiv:hep-th/ v1 18 Jan 2002 Hmiltoniztion of theories with degenerte coordintes. D.M. Gitmn nd I.V. Tyutin Mrch 10, 2018 rxiv:hep-th/0201143v1 18 Jn 2002 Abstrct We consider clss of Lgrngin theories where prt of the coordintes does

More information

Thin-Plate Splines. Contents

Thin-Plate Splines. Contents Thin-Plte Splines Dvid Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Cretive Commons Attribution 4.0 Interntionl License. To view copy of this

More information

Matrices, Moments and Quadrature, cont d

Matrices, Moments and Quadrature, cont d Jim Lmbers MAT 285 Summer Session 2015-16 Lecture 2 Notes Mtrices, Moments nd Qudrture, cont d We hve described how Jcobi mtrices cn be used to compute nodes nd weights for Gussin qudrture rules for generl

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy

More information

Review of Gaussian Quadrature method

Review of Gaussian Quadrature method Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge

More information

Math 270A: Numerical Linear Algebra

Math 270A: Numerical Linear Algebra Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner

More information

The solutions of the single electron Hamiltonian were shown to be Bloch wave of the form: ( ) ( ) ikr

The solutions of the single electron Hamiltonian were shown to be Bloch wave of the form: ( ) ( ) ikr Lecture #1 Progrm 1. Bloch solutions. Reciprocl spce 3. Alternte derivtion of Bloch s theorem 4. Trnsforming the serch for egenfunctions nd eigenvlues from solving PDE to finding the e-vectors nd e-vlues

More information

Chapter 14. Matrix Representations of Linear Transformations

Chapter 14. Matrix Representations of Linear Transformations Chpter 4 Mtrix Representtions of Liner Trnsformtions When considering the Het Stte Evolution, we found tht we could describe this process using multipliction by mtrix. This ws nice becuse computers cn

More information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1 Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution

More information

Frame-like gauge invariant formulation for mixed symmetry fermionic fields

Frame-like gauge invariant formulation for mixed symmetry fermionic fields Frme-like guge invrint formultion for mixed symmetry fermionic fields rxiv:0904.0549v1 [hep-th] 3 Apr 2009 Yu. M. Zinoviev Institute for High Energy Physics Protvino, Moscow Region, 142280, Russi Abstrct

More information

13: Diffusion in 2 Energy Groups

13: Diffusion in 2 Energy Groups 3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups

More information

Conservation Law. Chapter Goal. 5.2 Theory

Conservation Law. Chapter Goal. 5.2 Theory Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very

More information

4 The dynamical FRW universe

4 The dynamical FRW universe 4 The dynmicl FRW universe 4.1 The Einstein equtions Einstein s equtions G µν = T µν (7) relte the expnsion rte (t) to energy distribution in the universe. On the left hnd side is the Einstein tensor which

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

Introduction to the Calculus of Variations

Introduction to the Calculus of Variations Introduction to the Clculus of Vritions Jim Fischer Mrch 20, 1999 Abstrct This is self-contined pper which introduces fundmentl problem in the clculus of vritions, the problem of finding extreme vlues

More information

Green s functions. f(t) =

Green s functions. f(t) = Consider the 2nd order liner inhomogeneous ODE Green s functions d 2 u 2 + k(t)du + p(t)u(t) = f(t). Of course, in prctice we ll only del with the two prticulr types of 2nd order ODEs we discussed lst

More information

Chapter 3. Vector Spaces

Chapter 3. Vector Spaces 3.4 Liner Trnsformtions 1 Chpter 3. Vector Spces 3.4 Liner Trnsformtions Note. We hve lredy studied liner trnsformtions from R n into R m. Now we look t liner trnsformtions from one generl vector spce

More information

Generalizations of the Basic Functional

Generalizations of the Basic Functional 3 Generliztions of the Bsic Functionl 3 1 Chpter 3: GENERALIZATIONS OF THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 3.1 Functionls with Higher Order Derivtives.......... 3 3 3.2 Severl Dependent Vribles...............

More information

Quadratic Forms. Quadratic Forms

Quadratic Forms. Quadratic Forms Qudrtic Forms Recll the Simon & Blume excerpt from n erlier lecture which sid tht the min tsk of clculus is to pproximte nonliner functions with liner functions. It s ctully more ccurte to sy tht we pproximte

More information

Theoretical foundations of Gaussian quadrature

Theoretical foundations of Gaussian quadrature Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of

More information

Quantum Physics II (8.05) Fall 2013 Assignment 2

Quantum Physics II (8.05) Fall 2013 Assignment 2 Quntum Physics II (8.05) Fll 2013 Assignment 2 Msschusetts Institute of Technology Physics Deprtment Due Fridy September 20, 2013 September 13, 2013 3:00 pm Suggested Reding Continued from lst week: 1.

More information

In Section 5.3 we considered initial value problems for the linear second order equation. y.a/ C ˇy 0.a/ D k 1 (13.1.4)

In Section 5.3 we considered initial value problems for the linear second order equation. y.a/ C ˇy 0.a/ D k 1 (13.1.4) 678 Chpter 13 Boundry Vlue Problems for Second Order Ordinry Differentil Equtions 13.1 TWO-POINT BOUNDARY VALUE PROBLEMS In Section 5.3 we considered initil vlue problems for the liner second order eqution

More information

5.5 The Substitution Rule

5.5 The Substitution Rule 5.5 The Substitution Rule Given the usefulness of the Fundmentl Theorem, we wnt some helpful methods for finding ntiderivtives. At the moment, if n nti-derivtive is not esily recognizble, then we re in

More information

N 0 completions on partial matrices

N 0 completions on partial matrices N 0 completions on prtil mtrices C. Jordán C. Mendes Arújo Jun R. Torregros Instituto de Mtemátic Multidisciplinr / Centro de Mtemátic Universidd Politécnic de Vlenci / Universidde do Minho Cmino de Ver

More information

Conducting Ellipsoid and Circular Disk

Conducting Ellipsoid and Circular Disk 1 Problem Conducting Ellipsoid nd Circulr Disk Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 (September 1, 00) Show tht the surfce chrge density σ on conducting ellipsoid,

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

MatFys. Week 2, Nov , 2005, revised Nov. 23

MatFys. Week 2, Nov , 2005, revised Nov. 23 MtFys Week 2, Nov. 21-27, 2005, revised Nov. 23 Lectures This week s lectures will be bsed on Ch.3 of the text book, VIA. Mondy Nov. 21 The fundmentls of the clculus of vritions in Eucliden spce nd its

More information

Math 8 Winter 2015 Applications of Integration

Math 8 Winter 2015 Applications of Integration Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

More information

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1 MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further

More information

p-adic Egyptian Fractions

p-adic Egyptian Fractions p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction

More information

Best Approximation. Chapter The General Case

Best Approximation. Chapter The General Case Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given

More information

SOME PROPERTIES OF CHEBYSHEV SYSTEMS

SOME PROPERTIES OF CHEBYSHEV SYSTEMS SOME PROPERTIES OF CHEBYSHEV SYSTEMS RICHARD A. ZALIK Abstrct. We study Chebyshev systems defined on n intervl, whose constituent functions re either complex or rel vlued, nd focus on problems tht my hve

More information

INTERNATIONAL CENTRE FOR THEORETICAL PHYSICS THE ALGEBRAIC APPROACH TO THE SCATTERING PROBLEM ABSTRACT

INTERNATIONAL CENTRE FOR THEORETICAL PHYSICS THE ALGEBRAIC APPROACH TO THE SCATTERING PROBLEM ABSTRACT IC/69/7 INTERNAL REPORT (Limited distribution) INTERNATIONAL ATOMIC ENERGY AGENCY INTERNATIONAL CENTRE FOR THEORETICAL PHYSICS THE ALGEBRAIC APPROACH TO THE SCATTERING PROBLEM Lot. IXARQ * Institute of

More information

Math 115 ( ) Yum-Tong Siu 1. Lagrange Multipliers and Variational Problems with Constraints. F (x,y,y )dx

Math 115 ( ) Yum-Tong Siu 1. Lagrange Multipliers and Variational Problems with Constraints. F (x,y,y )dx Mth 5 2006-2007) Yum-Tong Siu Lgrnge Multipliers nd Vritionl Problems with Constrints Integrl Constrints. Consider the vritionl problem of finding the extremls for the functionl J[y] = F x,y,y )dx with

More information

Do the one-dimensional kinetic energy and momentum operators commute? If not, what operator does their commutator represent?

Do the one-dimensional kinetic energy and momentum operators commute? If not, what operator does their commutator represent? 1 Problem 1 Do the one-dimensionl kinetic energy nd momentum opertors commute? If not, wht opertor does their commuttor represent? KE ˆ h m d ˆP i h d 1.1 Solution This question requires clculting the

More information

A Convergence Theorem for the Improper Riemann Integral of Banach Space-valued Functions

A Convergence Theorem for the Improper Riemann Integral of Banach Space-valued Functions Interntionl Journl of Mthemticl Anlysis Vol. 8, 2014, no. 50, 2451-2460 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/10.12988/ijm.2014.49294 A Convergence Theorem for the Improper Riemnn Integrl of Bnch

More information

Notes on length and conformal metrics

Notes on length and conformal metrics Notes on length nd conforml metrics We recll how to mesure the Eucliden distnce of n rc in the plne. Let α : [, b] R 2 be smooth (C ) rc. Tht is α(t) (x(t), y(t)) where x(t) nd y(t) re smooth rel vlued

More information

The Active Universe. 1 Active Motion

The Active Universe. 1 Active Motion The Active Universe Alexnder Glück, Helmuth Hüffel, Sš Ilijić, Gerld Kelnhofer Fculty of Physics, University of Vienn helmuth.hueffel@univie.c.t Deprtment of Physics, FER, University of Zgreb ss.ilijic@fer.hr

More information

Review of Riemann Integral

Review of Riemann Integral 1 Review of Riemnn Integrl In this chpter we review the definition of Riemnn integrl of bounded function f : [, b] R, nd point out its limittions so s to be convinced of the necessity of more generl integrl.

More information

1 Online Learning and Regret Minimization

1 Online Learning and Regret Minimization 2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in

More information

Some Methods in the Calculus of Variations

Some Methods in the Calculus of Variations CHAPTER 6 Some Methods in the Clculus of Vritions 6-. If we use the vried function ( α, ) α sin( ) + () Then d α cos ( ) () d Thus, the totl length of the pth is d S + d d α cos ( ) + α cos ( ) d Setting

More information

Lecture 13 - Linking E, ϕ, and ρ

Lecture 13 - Linking E, ϕ, and ρ Lecture 13 - Linking E, ϕ, nd ρ A Puzzle... Inner-Surfce Chrge Density A positive point chrge q is locted off-center inside neutrl conducting sphericl shell. We know from Guss s lw tht the totl chrge on

More information

Solution to HW 4, Ma 1c Prac 2016

Solution to HW 4, Ma 1c Prac 2016 Solution to HW 4 M c Prc 6 Remrk: every function ppering in this homework set is sufficiently nice t lest C following the jrgon from the textbook we cn pply ll kinds of theorems from the textbook without

More information

1 The Lagrange interpolation formula

1 The Lagrange interpolation formula Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x

More information

Chapter 0. What is the Lebesgue integral about?

Chapter 0. What is the Lebesgue integral about? Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous

More information

Sufficient condition on noise correlations for scalable quantum computing

Sufficient condition on noise correlations for scalable quantum computing Sufficient condition on noise correltions for sclble quntum computing John Presill, 2 Februry 202 Is quntum computing sclble? The ccurcy threshold theorem for quntum computtion estblishes tht sclbility

More information