Generalization of the persistent random walk to dimensions greater than 1

Size: px
Start display at page:

Download "Generalization of the persistent random walk to dimensions greater than 1"

Transcription

1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental, Universitat e Barcelona, Diagonal, 647, Barcelona 08028, Spain Receive 16 July 1998 We propose a generalization of the persistent ranom walk for imensions greater than 1. Base on a cubic lattice, the moel is suitable for an arbitrary imension. We stuy the continuum limit an obtain the equation satisfie by the probability ensity function for the position of the ranom walker. An exact solution is obtaine for the projecte motion along an axis. This solution, which is written in terms of the free-space solution of the one-imensional telegrapher s equation, may open a new way to aress the problem of light propagation through thin slabs. S X PACS number s : j, w, r I. INTRODUCTION The persistent ranom walk PRW, first introuce by Fürth 1 an shortly after by Taylor 2, is probably the earliest an simplest generalization of the orinary ranom walk that incorporates some form of momentum in aition to ranom motion. The persistent ranom walk iffers from the orinary ranom walk in that the probabilistic quantity use at each step is the probability of continually moving in a given irection rather than the probability of moving in a given irection regarless of the irection of the preceing step. In this way the PRW introuces some form of momentum, i.e., persistence, into the purely ranom motion. This remarkable feature of the moel is one of the reasons why the PRW has been recently applie to escribe scattering an iffusion in isorere meia 3. For one-imensional lattices an in the continuum or iffusive limit the probability ensity function for the isplacement at time t, p(x,t), satisfies the telegrapher s equation TE 4 : 2 p p 2 2 p v2 2 x 2. As is well known, Eq. 1 has solutions with a finite velocity of propagation given by v 6. This fact has justifie the extensive use of TE as a generalization of the mesoscopic iffusion equations in fiels such as heat propagation 5 an light ispersion in turbi meia 6,7. However, none of the generalizations explore in two an three imensions of persistent ranom walks obeys the TE in the continuum limit On the other han, an besies the recent work of Gooy et al. 10 on two-imensional walks, there have been, to our knowlege, very few attempts to generalize the PRW to imensions higher than 1 in spite of its potential for moeling transport in isorere meia. One of the reasons for the lack of such a generalization is the nonexistence of a unique generalization of the PRW to imensions greater than 1 since several kins of lattices cubic, hexagonal, etc. can be use for the extension. Our main goal in this paper is to propose a generalization of the PRW to higher imensions assuming a 1 cubic lattice, an to obtain the governing equations for the probability ensity function of the process in the continuum limit. Unfortunately, in the continuum limit an for imensions greater than 1, the probability ensity function of the process oes not obey a higher-imensional telegrapher s equation. Nevertheless, in the context of transport in isorere meia, the partial ifferential equation escribing particle concentration i.e., the probability ensity function can suggest new approximations for the transport equation of more realistic moels. This is the case, for example, of light propagation in turbi meia where such an approach becomes extremely useful, especially when photons propagate in constraine geometries such as thin slabs where the Gaussian approximation becomes quite imprecise 12. In aition, the telegrapher s equation has been shown not to furnish better results than the iffusion approximation in two an three imensions 13. As we have mentione, in the one-imensional case the equation satisfie by p(x,t) is the telegrapher s equation. In two imensions, the moel consiere herein was partially analyze by Gooy et al. 10 but the equation for the probability ensity function, p(r,t), of the process was not obtaine. Another goal of this paper is to stuy the projecte motion, along a given irection, of the higher-imensional PRW. This projecte motion is relevant in the stuy of light propagation in turbi meia. Inee, when persistent ranom walks are use as moels for light propagation through slabs, the basic information is containe in the motion projecte along the coorinate orthogonal to the faces of the slab 14. We thus obtain the equation that governs the evolution of the projecte motion an write its solution in terms of the freespace solution of the one-imensional telegrapher s equation. The paper is organize as follows. In Sec. II we set the general analysis for the PRW in a cubic lattice of imension. In Sec. III we obtain the continuum limit an fin the governing equation for the probability ensity function of the resulting process. Section IV is evote to the projecte motion on a given coorinate an we obtain an exact expression of the ensity. Conclusions are rawn in Sec. V. II. GENERAL ANALYSIS We consier a ranom walk in a cubic lattice of arbitrary imension. The istance between nearest lattice points is l X/98/58 6 / /$15.00 PRE The American Physical Society

2 PRE 58 GENERALIZATION OF THE PERSISTENT RANDOM Jumps to another lattice point occur after a time interval. At each point, the ranom walker can take 2 ifferent irections. Among all possible events, we are only intereste in three of them: a The forwar scattering, where the ranom walker moves in the same irection as the previous jump, b the backwar scattering, where the walker reverses its previous irection, an c the scattering in the 2( 1) remaining irections. Let us enote by,, an the probabilities of events a, b, an c, respectively. As absorption will not be consiere here, the scattering probabilities satisfy the normalization conition This kin of ranom walk is terme persistent because it introuces a persistent probability an generalizes the oneimensional PRW to arbitrary imensions 3. Let us now set the general equations of this moel. We efine a set of auxiliary probabilities P n ( k) (i 1,...,i ), k 1,...,, where P n ( k) (i 1,...,i ) is the probability that the walker reaches the lattice point (i 1,...,i ) at step n moving along irection k. A similar efinition applies to P n ( k) (i 1,...,i ). Following an analogous reasoning to that of the one-imensional PRW 3, we can see that P n ( k) (i 1,...,i ) obeys the following set of recursive equations: P k n 1 i 1,...,i P k n i 1,...,i k 1,...,i P k n i 1,...,i k 1,...,i P j n i 1,...,i k 1,...,i P n j i 1,...,i k 1,...,i, 3 P k n 1 i 1,...,i P k n i 1,...,i k 1,...,i P k n i 1,...,i k 1,...,i P j n i 1,...,i k 1,...,i P n j i 1,...,i k 1,...,i. 4 Let us briefly explain how these equations can be obtaine. Consier an the origin locate at (0,...,0). Suppose the ranom walker has reache the origin at step n 1 moving along irection 1. Then, it necessarily was at the lattice point ( 1,0,...,0) at step n. The probability that the ranom walker jumps from this point to the origin epens on the arrival irection to the point ( 1,0,...,0). Each irection an there are 2 irections contributes to the total probability P ( 1) n 1 (0,...,0) with a ifferent weight. Inee, the term P ( 1) n ( 1,0,...,0) gives the probability that the ranom walker arrives at ( 1,0,...,0) along irection 1 an keeps going on the same irection. The probability that the ranom walker reverses its arrival irection after reaching ( 1,0,...,0) is P ( 1) n ( 1,0,...,0). The rest of the 2( 1) irections contributes with the terms P ( k) n ( 1,0,...,0), with. Equations 3 an 4 are easily obtaine after generalizing this reasoning to arbitrary points (i 1,...,i ) an irections k. These equations completely characterize our extension of the PRW an they are a convenient starting point for numerical analysis when no further analytical treatment can be mae. III. THE CONTINUUM LIMIT We now procee to the continuum limit. In this situation the length of each step, l, an the time interval between jumps,, both go to zero in such a way that the ranom walker moves at finite velocity v, l v lim. l, 0 5 As in the case of the one-imensional PRW 3,4, we also have to scale the scattering probabilities in the form 1, c, a /2, 6 where is the interval between jumps an is a parameter whose units are T 1. Let us now see what the physical meaning is of the parameters appearing in Eq. 6. We first observe that, as a irect consequence of the scaling 6, the occurrence of collision events is governe by the Poisson law. Inee, in the iscrete case the probability that the ranom walker jumps k times in the same irection is k. Therefore, in the continuum limit the probability that the particle keeps moving in the same irection for a time t k is lim k lim 1 t/ e. 0 0 We thus see the physical meaning of the parameter efine in Eq. 6, since it represents the mean frequency at which the ranom walker changes its irection, that is, is the mean number of scattering events per unit time, or equivalently 1 is the mean time between collision events. As a consequence, since is the probability of reversing irection, we see from Eq. 6 that c/ is the conitional probability of reversing the irection of motion. Analogously a/(2 ) is the probability that there is a turn to an orthogonal irection. Finally, the normalization conition requires that c 1 a. We also note that the angle between irections of motion before an after a collision is necessarily a multiple of /2 7

3 6994 BOGUÑÁ, PORRÀ, AND MASOLIVER PRE 58 p r,t p k r,t p k r,t. 10 For 1, we alreay explaine that p(x,t) evolves accoring to a TE, Eq. 1. We now present the partial ifferential equation escribing the evolution of this ensity when 2. We show in Appenix A that this equation reas t t b t a 1 t b 1 t a 2 t a v 2 2 p r,t 0, 11 FIG. 1. A realization of the continuum moel in three imensions escribe in the paper. Note that the ranom walker moves along irecctions that are parallel to the axes. recall that the iscrete moel has been built on a cubic lattice. Therefore, the only possible velocities of the ranom walker are the 2 values vî k (,...,), where î k is the unit vector in the x k coorinate irection. Figure 1 shows a realization of the moel in three imensions. In orer to get the iffusive limit of Eqs. 3 an 4, we efine t n an x k i k l,,...,, an let l 0 an 0 with the conition that v efine by Eq. 5 is constant. In this limit, probabilities P ( ) n (i 1,...,i ) become P k n 1 i 1,...,i j,...,i p k r,t p k r,t O 2, P n k i 1,...,i j 1,...,i p k r,t l p k r,t x j O l 2, 8 where p ( k) (r,t)r is the joint probability ensity function for the position r (x 1,...,x ) an the velocity of the ranom walker at time t. All the information about the continuous moel is containe in the following set of equations: p k r,t p k r,t v p k r,t p k r,t cp k r,t 1 2 a p j r,t p j r,t, v p k r,t p k r,t cp k r,t 1 2 a p j r,t p j r,t, which are the result of applying the continuum limit 8 to the general recursive equations 3 an 4. In many applications the most interesting quantity is the probability ensity p(r,t) for the position inepenent of the velocity, that is, 9 where 2 is the Laplacian of the spatial coorinates, b c, an is an operator incluing all spatial partial erivatives of fourth orer or greater. For 2 an 3, the operator has the following expressions: 3 v 6 6 x 2 y 2 z 2 2 v 4 4 x 2 y 2, 12 v 4 t b t 2a x 2 y 2 x 2 z 2 y 2 z We have thus obtaine the 2 orer partial ifferential equation that satisfies the probability ensity function of the persistent ranom walk in higher imensions. This equation oes not show spherical symmetry because of the velocities allowe by the moel. Nevertheless, using the so-calle ominant balance technique 18 one can easily see that the behavior of the probability ensity function at long times, an for positions r sufficiently far away from the moving bounary, is given by the lowest-orer partial erivatives. This central limit approximation transforms Eq. 11 into the iffusion equation t p r,t D 2 p r,t, where the iffusion constant D is D v2 b Note that the iffusion constant is in agreement with the result that follows from transport theory D vl t /, where l t v/b is the transport mean free path 15. IV. THE PROJECTED MOTION Let us now stuy the projecte motion of the PRW on a given axis. The probability ensity function of this motion is given by p x,t p r,t x2 x. This marginal ensity obeys a simpler equation than Eq. 11 mainly because the integration of p(r,t) over the coorinates x 2 x is zero. We show in Appenix B that the equation for p(x,t) reas t t b t a p x,t v 2 t a 2 x 2p x,t. 16

4 PRE 58 GENERALIZATION OF THE PERSISTENT RANDOM This equation is a thir-orer hyperbolic partial ifferential equation with finite velocity of propagation v. The reuction in the orer from 2 for Eq. 11 to 3 for Eq. 16 is ue to the fact that the projecte motion is couple to the motion along any other irection by a first-orer equation see Appenix B. Equation 16 may provie a better approximation to transport problems than the iffusion equation, especially in thin slabs where the effect of taking into account the spee of propagation is crucial 16. Note that Eq. 16 possesses a richer structure than the telegrapher s equation. Moreover, we have been able to erive Eq. 16 from a microscopic moel we recall that any attempt to erive TE from a microscopic moel of transport fails in higher imensions. Nevertheless, Eq. 16 shares with TE at least two important features: i a finite velocity of propagation, an ii a similar asymptotic behavior of the moments see below. An equation similar to Eq. 16 was use several years ago in the context of heat propagation in rigi solis 17. It is possible to exactly solve Eq. 16 in the Fourier- Laplace space for the isotropic initial conitions p k r,t p k r,t 1 2 r, 17 which provie the three initial conitions require for solving Eq. 16, p x,t 0 x, 2 p x,t 2 p x,t 0, 0 0 v2 x. 18 In Appenix C we show how to erive these conitions from Eq. 17. It is possible to solve Eq. 16 in the Fourier- Laplace space. To this en we efine the joint transform pˆ,s e st e i x p x,t x t. 0 Then the transformation of Eq. 16 an the use of Eq. 18 lea to the solution pˆ,s s a s b v / s s a s b v 2 2 s a. 19 We observe that this case is equivalent to a three-state continuous-time ranom walk where the particle is moving to the right or left with velocity v ( v) or is at rest, which correspons to the motion in an orthogonal irection. The transition matrix of this three-state walk is 3 c/ 1 a/ p i j 0 c/ 0 1 a/, 20 a/ 2 a/ 2 1 a/ where i 1 correspons to the state with velocity v, i 2 with velocity v, an i 3 to the state without isplacement along the x irection. Let us now stuy the moments m n (t), n 1,2,3,..., of the istribution p(x,t). In terms of the characteristic function, p (,t), the moments are given by m n (t) i n pˆ (,t) 0. Due to the isotropic initial conitions, Eq. 18, all o moments vanish, m 2n 1 (t) 0. Then it follows from Eq. 19 that the Laplace transform of m 2n (t) is given by mˆ 2n s 2n! s a s s a The leaing behavior for small s is v n 1 s s b 2 n. 21 mˆ 2n s v2n 2n! b n n s 1 n s 0 an the leaing behavior for large s is mˆ 2n s v2n 2n! s 1 2n s. Then, by the Tauberian theorems 3,18, we get the asymptotic behaviors an m 2n 2n! n! v2 n b t 22 m 2n v2n t2n Note that the moments as t are ientical to those of orinary iffusion with D v 2 /b in agreement with the result of the preceing section, which means that p(x,t) behaves in this limit as p x,t 1 4 Dt exp x2 /4Dt. 24 On the other han, when t 0, the motion behaves as if it were eterministic. In fact, taking into account the isotropic initial conitions, we can easily show that the ensity p(x,t) has the following expansion at short times: p x,t 1 2 x vt x vt x 0, 25 which immeiately leas to the result 23 for the moments. Therefore, the behavior of moments is similar to that of the solution of the telegrapher s equation. Inee, we see from Eqs. 22 an 23 that if t, we have orinary iffusion where, for instance, the secon moment goes as t), while if t 0 the behavior is eterministic where the secon moment goes as t 2 ). We have not been able to invert Eq. 19 in general. However, the solution in real space can be written for the two-

5 6996 BOGUÑÁ, PORRÀ, AND MASOLIVER PRE 58 imensional case 2 an when the backscattering probability vanishes, c 0. In effect, now a b, an the expression for pˆ (,s) reas pˆ,s 1 2 s 1 2 pˆ te,s 1 2 s where s 2 pˆ te,s s s 2 v 2 2 pˆ te,s, is the Fourier-Laplace transform of the free-space solution of the one-imensional telegrapher s equation for isotropic initial conitions. Inee, the Fourier-Laplace transform of Eq. 1 along with the isotropic initial conitions p te (x,0) (x) an t p te (x,0) 0 leas to Eq. 27 6,19. Therefore, we see from Eq. 26 that the sought-after expression reas p x,t 1 2 e x 1 2 p te x,t 2 e 0 t e p te x,, where p te (x,t) is 4,19 p te x,t e 2v e x /v 2v x /v I 0 2 x 2 /v 2 t 2 x 2 /v I 1 2 x 2 /v 2, where I n (z) are moifie Bessel functions. Surprisingly, p te (x, ) is also the solution of a three-state ranom walk with transition matrix 20 when the initial conitions have an equal probability of moving to the left or to the right, an zero probability to be at rest. This fact clarifies the origin of solution 28. Inee, the initial conitions for that solution were 1/2 probability to be at rest an 1/4 to move in either x irection. The term in p(x,t) accounts for the probability that the ranom walker keeps moving along the y axis for all perio t. The secon term is the contribution of ranom walkers that are not at rest initially, accoring to the interpretation for p te (x, ) given above. Finally, particles at rest at t 0 that begin to move at some time between (0,t) can be consiere as a source of particles that evolve as p te (x, ). The thir term is therefore the convolution of the source function ( e ) with p te (x, ). In Fig. 2, we plot the solution for ifferent times without terms. When t is less than 1, the -function terms account for most of p(x,t). As time grows, the contribution of these terms ecays exponentially an the solution converges to the Gaussian istribution, Eq. 24. V. CONCLUSIONS FIG. 2. Probability ensity function p(x,t) for the motion projecte along the x irection at ifferent times as a function of the position x. Soli line: 2xt 1; ashe line: 2 5; ot-ashe line: The generalization of the persistent ranom walk to imensions higher than 1 oes not evolve accoring to a telegrapher s equation as it oes in one imension. We have explicitly obtaine the 2 orer partial ifferential equation governing the evolution of the probability ensity function. Therefore, higher-imensional persistent ranom walks in cubic lattices cannot be consiere as microscopic moels to erive TE in imensions greater than 1. Moreover, even the motion projecte along an axis, besies being one-imensional, oes not evolve accoring to a TE. We have foun that the probability ensity of the projection obeys a thir-orer partial ifferential equation. The orer of the equation is completely consistent with the fact that the projecte motion is equivalent to a three-state ranom walk. The evolution equation for the projecte motion may be very useful when consiere as an approximation to the full transport equation for moels of light propagation because Eq. 16 may overcome the problems of the iffusive approximation when light propagates through thin slabs. In this sense we have shown that the evolution equation for the projecte motion correctly links the expecte short-time behavior eterministic with the long-time behavior iffusive of these moels. ACKNOWLEDGMENTS This work has been supporte in part by Dirección General e Investigación CientíficayTécnica uner Contract No. PB , an by Societat Catalana e Física Institut Estuis Catalans. APPENDIX A: DERIVATION OF EQ. 11 Let us efine the following functions: U k r,t p k r,t p k r,t, V k r,t p k r,t p k r,t, A1 for,...,. The probability ensity of the position becomes

6 PRE 58 GENERALIZATION OF THE PERSISTENT RANDOM an the set of Eqs. 9 reas U k r,t p r,t U k r,t, v V k r,t au r,t ap r,t, k A2 where x 1 was enote by x an U k (x,t) an V k (x,t) correspon to the functions U k (r,t) an V k (r,t) integrate over the coorinates x 2 x. The first equation of system B1 implies that the quantity q x,t U k x,t k 2 V k r,t v U k r,t bv r,t, A3 k where b c. If we efine the operator D t, where is a constant, then the equation for U k (r,t) becomes D b ap r,t D a U k r,t v 2 2 U k r,t,,...,. 2 A4 For 1, this equation coincies with the telegrapher s equation because, in this case, a 0 an U 1 (x,t) p(x,t). For 2, the efinition of the operator S k D b D a v 2 2 xk further simplifies Eq. A4 into S k U k r,t ad b p r,t. A5 Finally, applying i 1 S i to both sies of Eq. A2 an using Eq. A5, we obtain the equation for p(r,t), i 1 S i p r,t ad b i k S i p r,t. A6 The expansion of the prouct i 1 S i in powers of D b D a leas to Eq. 11. APPENDIX B: DERIVATION OF EQ. 16 The integration of the system A3 over the coorinates x 2 x yiels satisfies the first-orer equation q x,t aq x,t a 1 U 1 x,t q x,t. B3 This equation along with system B2 leas to a thir-orer partial ifferential equation for p x,t U 1 x,t q x,t. Inee, from Eqs. B1 an B2 we get an t D b p x,t v 2 x 2 U 1 x,t D a U 1 x,t D a p x,t. B4 B5 Finally, the combination of Eqs. B4 an B5 immeiately leas to Eq. 16. APPENDIX C: INITIAL CONDITIONS The integration of the isotropic initial conitions 17 over the coorinates x 2 x gives the initial conitions for the marginal probability p k x,t p k x,t 1 2 x. The first initial conition follows irectly from them, p x,t 0 p k x,t 0 p k x,t 0 x. U k x,t au k x,t ap x,t, B1 Taking into account that V k (x,t 0) 0, we see from Eq. A3 that for k 2,..., an U 1 x,t V k x,t bv k x,t, v V 1 x,t au x 1 x,t ap x,t, V 1 x,t v U 1 x,t bv x 1 x,t, B2 U k x,t a 1 0 x a x 0, an therefore the secon initial conition reas p x,t 0. 0 C1 C2 After taking the erivative with respect to time in the equa-

7 6998 BOGUÑÁ, PORRÀ, AND MASOLIVER PRE 58 tion for U k (r,t) of system A3 an using the equation for V k (r,t) in the same system, we get 2 U k r,t 2 v 2 2 U k r,t 2 ap r,t. bv V k r,t au r,t k The sum over all k an the integration over coorinates x 2 x of this equation gives the thir initial conition 2 p x,t 2 0 v2 x. 1 R. Furth, Ann. Phys. Leipzig 53, G. I. Taylor, Proceeing of the Lonon Mathematical Society 20, G. H. Weiss, Aspects an Applications of the Ranom Walk North-Hollan, Amsteram, S. Golstein, Q. J. Mech. Appl. Math. 4, D. D. Joseph an L. Preziosi, Rev. Mo. Phys. 61, ; 62, J. Masoliver an G. H. Weiss, Eur. J. Phys. 17, A. Ishimaru, Wave Propagation an Scattering in Ranom Meia Acaemic, New York, 1978, Vol. I.; Appl. Opt. 28, J. Masoliver, J. M. Porrà, an G. H. Weiss, Physica A 182, J. Masoliver, J. M. Porrà, an G. H. Weiss, Physica A 193, S. Gooy an L. S. García-Colín, Phys. Rev. E 55, J. M. Porrà, J. Masoliver, an G. H. Weiss, Physica A 218, B. B. Das, Feng Liu, an R. R. Alfano, Rep. Prog. Phys. 60, J. M. Porrà, J. Masoliver, an G. H. Weiss, Phys. Rev. E 55, D. J. Durian an J. Runick, J. Opt. Soc. Am. A 14, J. J. Duerstant an W. R. Martin, Transport Theory J. Wiley, New York, K. M. Yoo, Feng Liu, an R. R. Alfano, Phys. Rev. Lett. 64, M. E. Gurtin an A. C. Pipkin, Arch. Ration. Mech. Anal. 31, J. Masoliver, Phys. Rev. A 45, J. Masoliver, J. M. Porrà, an G. H. Weiss, Phys. Rev. A 45,

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The effect of nonvertical shear on turbulence in a stably stratified medium The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Average value of position for the anharmonic oscillator: Classical versus quantum results

Average value of position for the anharmonic oscillator: Classical versus quantum results verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive

More information

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Acute sets in Euclidean spaces

Acute sets in Euclidean spaces Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of

More information

Applications of the Wronskian to ordinary linear differential equations

Applications of the Wronskian to ordinary linear differential equations Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.

More information

arxiv:nlin/ v1 [nlin.cd] 21 Mar 2002

arxiv:nlin/ v1 [nlin.cd] 21 Mar 2002 Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Topic 2.3: The Geometry of Derivatives of Vector Functions

Topic 2.3: The Geometry of Derivatives of Vector Functions BSU Math 275 Notes Topic 2.3: The Geometry of Derivatives of Vector Functions Textbook Sections: 13.2 From the Toolbox (what you nee from previous classes): Be able to compute erivatives scalar-value functions

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

P. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom

P. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom Time-harmonic torsional waves in a composite cyliner with an imperfect interface J. R. Berger a) Division of Engineering, Colorao School of Mines, Golen, Colorao 80401 P. A. Martin b) Department of Mathematics,

More information

model considered before, but the prey obey logistic growth in the absence of predators. In

model considered before, but the prey obey logistic growth in the absence of predators. In 5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an

More information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the

More information

The continuity equation

The continuity equation Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Conservation laws a simple application to the telegraph equation

Conservation laws a simple application to the telegraph equation J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

Extinction, σ/area. Energy (ev) D = 20 nm. t = 1.5 t 0. t = t 0

Extinction, σ/area. Energy (ev) D = 20 nm. t = 1.5 t 0. t = t 0 Extinction, σ/area 1.5 1.0 t = t 0 t = 0.7 t 0 t = t 0 t = 1.3 t 0 t = 1.5 t 0 0.7 0.9 1.1 Energy (ev) = 20 nm t 1.3 Supplementary Figure 1: Plasmon epenence on isk thickness. We show classical calculations

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Lower Bounds for the Smoothed Number of Pareto optimal Solutions Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.

More information

WUCHEN LI AND STANLEY OSHER

WUCHEN LI AND STANLEY OSHER CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Calculus in the AP Physics C Course The Derivative

Calculus in the AP Physics C Course The Derivative Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Introduction to Markov Processes

Introduction to Markov Processes Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav

More information

4. Important theorems in quantum mechanics

4. Important theorems in quantum mechanics TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012 arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Convergence of Random Walks

Convergence of Random Walks Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

Some Examples. Uniform motion. Poisson processes on the real line

Some Examples. Uniform motion. Poisson processes on the real line Some Examples Our immeiate goal is to see some examples of Lévy processes, an/or infinitely-ivisible laws on. Uniform motion Choose an fix a nonranom an efine X := for all (1) Then, {X } is a [nonranom]

More information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (II)

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (II) Laplace s Equation in Cylinrical Coorinates an Bessel s Equation (II Qualitative properties of Bessel functions of first an secon kin In the last lecture we foun the expression for the general solution

More information

Applications of First Order Equations

Applications of First Order Equations Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more

More information

Infinite Body Centered Cubic Network of Identical Resistors. Abstract. and any other lattice site ( n

Infinite Body Centered Cubic Network of Identical Resistors. Abstract. and any other lattice site ( n Infinite Boy Centere Cubic Network of Ientical esistors J. H. Asa,*, A. A. Diab,. S. Hijjawi, an J. M. Khalifeh Dep. of Physics- Tabuk University- P. O Box 7- Tabuk 79- Saui Arabia General Stuies Department-

More information

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie

More information

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Bair, faculty.uml.eu/cbair University of Massachusetts Lowell 1. Pre-Einstein Relativity - Einstein i not invent the concept of relativity,

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

Momentum and Energy. Chapter Conservation Principles

Momentum and Energy. Chapter Conservation Principles Chapter 2 Momentum an Energy In this chapter we present some funamental results of continuum mechanics. The formulation is base on the principles of conservation of mass, momentum, angular momentum, an

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

More information

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into

More information

Power law behavior in diffraction from fractal surfaces

Power law behavior in diffraction from fractal surfaces Surface Science 409 (1998 L703 L708 Surface Science Letters Power law behavior in iffraction from fractal surfaces Y.-P. Zhao *, C.-F. Cheng 1, G.-C. Wang, T.-M. Lu Department of Physics, Applie Physics

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010 NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles

More information

Chapter 2 Governing Equations

Chapter 2 Governing Equations Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

arxiv: v1 [math-ph] 2 May 2016

arxiv: v1 [math-ph] 2 May 2016 NONLINEAR HEAT CONDUCTION EQUATIONS WITH MEMORY: PHYSICAL MEANING AND ANALYTICAL RESULTS PIETRO ARTALE HARRIS 1 AND ROBERTO GARRA arxiv:165.576v1 math-ph] May 16 Abstract. We stuy nonlinear heat conuction

More information

We G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies

We G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll

More information

arxiv: v1 [math-ph] 5 May 2014

arxiv: v1 [math-ph] 5 May 2014 DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic

More information

Monotonicity for excited random walk in high dimensions

Monotonicity for excited random walk in high dimensions Monotonicity for excite ranom walk in high imensions Remco van er Hofsta Mark Holmes March, 2009 Abstract We prove that the rift θ, β) for excite ranom walk in imension is monotone in the excitement parameter

More information

The Ehrenfest Theorems

The Ehrenfest Theorems The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent

More information

A. Incorrect! The letter t does not appear in the expression of the given integral

A. Incorrect! The letter t does not appear in the expression of the given integral AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question

More information

3 The variational formulation of elliptic PDEs

3 The variational formulation of elliptic PDEs Chapter 3 The variational formulation of elliptic PDEs We now begin the theoretical stuy of elliptic partial ifferential equations an bounary value problems. We will focus on one approach, which is calle

More information

Tractability results for weighted Banach spaces of smooth functions

Tractability results for weighted Banach spaces of smooth functions Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March

More information

Armenian Transformation Equations For Relativity

Armenian Transformation Equations For Relativity Armenian Transformation Equations For Relativity Robert Nazaryan an Haik Nazaryan 00th Anniversary of the Special Relativity Theory This Research one in Armenia 968-988, Translate from the Armenian Manuscript

More information

On the number of isolated eigenvalues of a pair of particles in a quantum wire

On the number of isolated eigenvalues of a pair of particles in a quantum wire On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

Markov Chains in Continuous Time

Markov Chains in Continuous Time Chapter 23 Markov Chains in Continuous Time Previously we looke at Markov chains, where the transitions betweenstatesoccurreatspecifietime- steps. That it, we mae time (a continuous variable) avance in

More information