arxiv:cond-mat/ v1 [cond-mat.stat-mech] 19 Jan 2006
|
|
- Thomasina Logan
- 6 years ago
- Views:
Transcription
1 Statistical Properties of Functionals of the Paths of a Particle Diffusing in a One-Dimensional Random Potential arxiv:cond-mat/6455v cond-mat.stat-mech] 9 Jan 6 Sanjib Sabhapandit,, Satya N. Majumdar, and Alain Comtet, Laboratoire de Physique Théorique et Modèles Statistiques, Université Paris-Sud, Bâtiment, 945 Orsay Cedex, France Université Pierre et Marie Curie - Paris 6, Institut Henri Poincaré, rue Pierre et Marie Curie, Paris, F-755, France Dated: June 6, 7 We present a formalism for obtaining the statistical properties of functionals and inverse functionals of the paths of a particle diffusing in a one-dimensional quenched random potential. We demonstrate the implementation of the formalism in two specific examples: where the functional corresponds to the local time spent by the particle around the origin and where the functional corresponds to the occupation time spent by the particle on the positive side of the origin, within an observation time window of size t. We compute the disorder average distributions of the local time, the inverse local time, the occupation time and the inverse occupation time, and show that in many cases disorder modifies the behavior drastically. PACS numbers: 5.4.-a,.5.-r, g I. INTRODUCTION The statistical properties of functionals of a one dimensional Brownian motion have been extensively studied and have found numerous applications in diverse fields ranging from probability theory,, 3], finance 4, 5, 6], mesocopic physics 7], computer science 8], and in understanding weather records 9]. The position xτ of a one dimensional Brownian motion evolves with time τ via the Langevin equation dx dτ = ητ, starting from x = x, where ητ is a thermal Gaussian white noise with mean ητ = and a correlator ητητ = δτ τ. A functional T is simply the integral up to time t T = t V xτ dτ, where V x is a prescribed non-negative function whose choice depends on the specific observable of interest. For a fixed initial position x of the Brownian motion and a fixed observation time t, the value of T varies from one history or realization of the Brownian path {xτ} to another see Fig. and a natural question is: what is the probability density function pdf PT t, x? Following the path integral methods devised by Feynman ], Kac showed, ] that the calculation of the pdf PT t, x can essentially be reduced to a quantum mechanics problem, namely solving a single particle Shrödinger equation in an external potential V x. This formalism is known in the literature as the celebrated Feynman-Kac formula. Subsequently, this method has been widely used to calculate the pdf of T with different choices of V x as demanded by specific applications. This has been reviewed recently in Ref. 8]. In particular, two most popular applications correspond respectively to the choices, V x = δx a and V x = θx, where δx is the Dirac s delta function and θx is the Heaviside step function. In the former T t t t 3 t 4 R R R 3 R 4 FIG. : Schematic plots of T defined by Eq. as a function of t, corresponding to four different realizations of the paths {xτ}, for τ t] denoted by R, R, R 3 and R 4 respectively. For fixed t t or t or t 3 or t 4, shown by vertical dashed lines, T takes different value for different realizations. On the other hand for a fixed T horizontal dashed line corresponding t is different for different realizations: t for R, t for R, t 3 for R 3 and t 4 for R 4. case, the corresponding functional Ta = t δxτ adτ has the following physical meaning: Tada is just the time spent by the particle in the vicinity of the point a in space, i.e., in the region a, a + da], out of the total observation time t. Note that, by definition, Tada = t. The functional Ta is known as the local time density in the literature. In the second case V x = θx, the functional T = t θxτdτ measures the time spent by the particle on the positive side of the origin out of the total time t and is known as the occupation time. The probability distribution of the occupation time was originally computed by Lévy ], t
2 T T/t PT t, dt = π arcsin, and is known as the arcsine law of Lévy. Since then, the local and the occupation times for pure diffusion, have been studied extensively by mathematicians, 3, 4, 5, 6, 7]. Recently, the study of the occupation time has seen a revival in physics literature and has been used in understanding the dynamics out of equilibrium in coarsening systems 8, 9], ergodicity properties in anomalously diffusive processes, ], in renewal processes ], in models related to spin glasses 3], in simple models of blinking quantum dots 4], and also in the context of persistence 5, 6]. Local and Occupation times have been also studied in the context of stochastic ergodicity breaking 7], first passage time 8], diffusion controlled reactions activated by catalytic sites 9] and diffusion on graphs 3, 3]. In polymer science, a long flexible polymer of length t is often modeled by a Brownian path up to time t. In this context, the local time at a position r is proportional to the concentration of monomers at r in a polymer of length t. A natural and important question is how to generalize the Feynman-Kac formalism to calculate the statistical properties of the functionals of the type in Eq. when xτ is not just a pure diffusion process, but it represents the position of a particle in an external random medium. While various properties of diffusion in random media have been widely studied in the past 3, 33, 34, 35], the study of statistical properties of functionals in random media is yet to receive its much deserved attention. In this paper we undertake this task. More precisely, we are interested in calculating the pdf PT t, x of a functional T as in Eq. where xτ now evolves via the Langevin equation dx dτ = F xτ + ητ 3 where ητ represents the thermal noise as in Eq. and Fx = du/dx represents the external force, the derivative of the potential Ux, felt by the particle. Most generally, the external potential consists of two parts, Ux = U d x + U r x, a deterministic part U d x and a random part U r x. The random part of the potential U r x is quenched in the sense that it does not change during the time evolution of the particle, but fluctuates from one sample to another according to some prescribed probability distribution. Consequently, the pdf PT t, x will also fluctuate from one sample of the random potential to another and the goal is to compute the disorder averaged pdf PT t, x where the... denote the average over the distribution of the random potential. A popular model for the random potential is the celebrated Sinai model 36], where various disorder averaged physical quantities can be computed analytically 3, 33, 37, 38, 39, 4, 4, 4], and yet the results exhibit rich and nontrivial behaviors and also capture many of the qualitative behaviors of more complex realistic disordered systems. The Sinai model assumes that U r x = Bx where Bx represents a Brownian motion in space, i.e., db dx = ξx 4 where ξx is a Gaussian white noise with mean ξx = and a correlator ξxξx = δx x. The constant represents the strength of the random potential. In this paper, we first present a generalization of the Feynman-Kac formalism to calculate the pdf PT t, x in presence of an arbitrary external potential Ux. To obtain explicit results using this formalism, we next assume that the random part of the potential is as in the Sinai model, i.e., our external potential is of the form Ux = U d x + Bx, where Bx is a Brownian motion in space and U d x is the non-random deterministic part of the potential. It turns out that the asymptotic behavior of the disorder averaged pdf PT t, x, quite generically, has three different qualitative behavior depending on the curvature of the deterministic potential U d x, i.e., whether U d x has a convex concave upward shape representing a stable potential i.e., attractive force towards the origin, a concave concave downward shape representing unstable potential a repulsive force away from the origin or just flat indicating the absence of any external potential. To facilitate explicit calculation, we model the deterministic potential simply by, U d x = µ x, so that µ < represents a stable potential, µ > represents an unstable potential and µ = represents a flat potential. This specific choice facilitates explicit calculation, but the results are qualitatively similar if one chooses another form of this potential. Thus, in our model, we will consider the external potential as Ux = µ x + Bx 5 where Bx = x ξx dx is the trajectory of a Brownian motion in space see Fig.. Note that the case µ = corresponds to the pure Sinai model. Figure shows typical realization of potentials for µ =, µ > and µ <. The corresponding force in Eq. 3 is simply given by Fx = µ signx + ξx. 6 We will demonstrate how to calculate explicitly, using our generalized Feynman-Kac formalism, the disorder averaged pdf PT t, x when the external potential is of the form given by Eq. 5. Despite the simplicity of the choice of the external potential, a variety of rich and interesting behavior can be obtained by tuning the parameter µ/, as shown in this paper. We will present detailed results for the two functionals, namely for the local time and the occupation time corresponding to the choices V x = δx and V x = θx respectively in Eq.. Also, to keep the discussion simple, we will present our final results for x = corresponding to the particle starting at the origin. However, our method is not limited only to this specific choice. Some of these results were briefly announced in a previous letter 43]. In addition, in this paper we also introduce the notion of inverse functional, which is defined as follows. If V x in Eq. is non-negative, then for each path {xτ}, T is a non-decreasing function of t, which we formally denote by T = g t {xτ}, x. Therefore for a given realization of path {xτ} and given T there is a unique value of t see
3 3 Ux µ = Ux µ > x Ux µ < x x FIG. : A classical particle represented by diffusing in a typical realization of the potential Ux = µ x + Bx, where Bx represents the trajectory of a Brownian motion in space with B =. The three figures are for µ =, µ > and µ < respectively. The dash lines show the potential for =. Fig., which we formally write as the inverse of the functional g 49] t = g T {xτ}, x. 7 This inverse time t physically means the observation time that is required for any given path {xτ} in order to produce a prescribed value of T. Of course, for the same value T, for a different path {xτ}, the value of t will be different. Thus, t is a random variable for a fixed T, which takes different values for different realizations of paths and we would like to compute its pdf, which we denote by It T, x and by definition It T, dt =. Clearly, this pdf will also differ from sample to sample of the external potential in Eq. 5 and our goal is to obtain the disorder averaged distribution It T, x. In this paper, we present detailed results for It T, again for the two choices of V x = δx and V x = θx corresponding to the local time and the occupation time respectively. The inverse local and occupation times are useful for experimentalists as they provide an estimate on the required measurement time. For example, in the context of polymers, the inverse local time is the typical length of a polymer required to obtain a certain monomer concentration. The rest of the paper is organized as follows. In Sec. II, we present our general approach for computing the pdf PT t, x of the functional T defined by Eq. for a given t, and the pdf It T, x of the inverse functional defined by Eq. 7 for a given T, for a given sample of the random potential, for arbitrary starting position of the particle x = x and for arbitrary but non-negative V x. After this section we consider the specific examples of local time and occupation time by setting V x = δx and V x = θx respectively. We will use different notations for the pdfs in the two examples to avoid any misunderstanding. In the first example, where T is the local time, we denote the pdf of the local time PT t, for a given t by P loc T t, and the pdf of the inverse local time It T, for a given T by I loc t T. In the second example, where T is the occupation time we denote the pdf of the occupation time PT t, for a given t by P occ T t and the pdf of the inverse occupation time It T, for a given T by I occ t T. While our final goal is to obtain the disorder averaged distributions P loc T t, I loc t T, P occ T t and I occ t T, it is however, instructive to study the pure case first, before tackling the problem with disorder which is obviously harder. In the same spirit, we have presented the detailed discussions on the local time, inverse local time, occupation time and inverse occupation time for the pure case = in Secs. III, V, VII and IX respectively, before computing their disorder average in Secs. IV, VI, VIII and X respectively. Sec. XI contains some concluding remarks. Some of the details are relegated to the appendixes. The results are summarized in Tables I, II, and III. II. GENERAL APPROACH In this section we will show how to compute the pdfs PT t, x and It T, x for arbitrary non-negative V x and arbitrary starting position x = x, for each realization of random force Fx, by using a backward Fokker-Planck equation approach. In the following discussion we will denote the functional defined in Eq. by g t {xτ}, x, and use T as the value of the functional for a given path {xτ}, for τ t]. Since V x is considered to be non-negative, T defined by Eq. has only positive support. Therefore, a natural step is to introduce the Laplace transform of the pdf PT t, x with respect to T Q p x, t = PT t, xe pt dt = e pg t {xτ},x = { exp p t x=x V xt ] dt } x=x, 8 where x=x denotes the average over all paths that start at the position x = x and propagate up to time t. Our aim is to derive a backward Fokker-Planck equation for Q p x, t with respect to the initial position x = x. We consider a particle starting from the initial position x, evolves via Eq. 3 up to time t + t. Now we split the time interval, t + t] into two parts: an infinitesimal interval, t], over which the particle experiences an infinitesimal displacement x from its initial position x and the remaining interval t, t + t] in which the particle evolves
4 4 from a starting position x + x. Since x = x, one gets t V xt ] dt = V x t + O t ]. Therefore using Eq. 8, and splitting the integral over t into the above two time intervals we obtain { } t+ t Q p x, t + t = exp p V xt ] dt x=x = e Q pv x t p x + x, t, 9 x where x denotes the average with respect to all possible displacements x. Now, in the limit t, integrating Eq. 3 one gets, x = Fx t + t ητdτ + O t ]. Hence, using the zero mean and the uncorrelated properties of the noise we get x x lim = Fx and lim =, t t t t Therefore, by Taylor expanding Eq. 9 for small x, and taking the limit t, one arrives at the backward Fokker- Planck equation, Q p t = Q p x + Fx Q p x pv xq p, with the initial condition Q p x, =, which is easily checked by Eq. 8. The advantage of the above equation over the usual Feynman-Kac formalism, ] is that, in the later case one has a forward Fokker-Planck equation spatial derivative with respect to the final position where after obtaining the solution of the differential equation, one has to again perform an additional step of integration over the final position. In contrast, Eq. involves the spatial derivatives with respect to the initial position of the particle, and hence no additional step of integration over the final position is required. The standard practice of attacking the partial differential equations of above type is by using the method of Laplace transform. We define the Laplace transform of Q p x, t with respect to t ux = = Q p x, te αt dt dt e αt dt e pt PT t, x, 3 where for notational convenience, we have suppressed the α and p dependence of ux. Now by taking Laplace transform of Eq. with respect to t we obtain the ordinary differential equation u x + Fxu x α + pv x]ux =, 4 where u x = du/dx. The appropriate boundary conditions ux ± are to be derived from the observation that if the particle starts at x ± it will never cross the origin in finite time. Note that Eq. 4 is valid for each sample of the quenched random force Fx. Thus in principle, from the solution ux one obtains PT t, x by inverting the double Laplace transform in Eq. 3 for each sample of quenched random potential, and then takes the average over the disorder. Our next goal in this section, is to show how to compute the pdf It T, x for a given sample of the quenched random force Fx. It turns out that It T, x is related to the pdf PT t, x in their Laplace space as shown below. By definition we have, It T, x = δ t g T {xτ}, x x. 5 However, it is elementary that for each realization of path {xτ} δ t g T {xτ}, x = δ T g t {xτ}, x dt dt, 6 where dt/dt is the usual Jacobian of the transformation, which is simply dt/dt as both T and t have only positive support. It immediately follows from the above two equations that It T, x = δ T g dt t {xτ}, x. 7 dt x Therefore, Laplace transform of It T, x with respect to T reads dt e pt It T, x = e t {xτ},x dg pg dt x = p t Q px, t, 8 where Q p x, t is given by Eq. 8. Now taking a further Laplace transform in Eq. 8 with respect to t, it is straightforward to obtain dt e αt dt e pt It T, x = αux. 9 p Thus, we have established via Eqs. 8 and 9, the relationships between the Laplace transforms of the pdf of the functional T defined by Eq. and the pdf of the inverse functional defined by Eq. 7. Hence, again in principle, from the solution ux of the ordinary differential equation 4, one obtains It T, x by inverting the double Laplace transform in Eq. 9 for each sample of quenched random potential, and then takes the average over the disorder. Note that putting α = in Eq. 9 and inverting the Laplace transform with respect to p immediately verifies the normalization condition It T, xdt =. In the rest of the paper, we will demonstrate how to implement this formalism for the particular examples of the local time corresponding to the choice V x = δx and the occupation time corresponding to the choice V x = θx. Since in these examples we consider the starting position of the particle to be the origin, we need to only find the solution u of
5 5 the differential equation 4. In each example, we will consider the pure cases = first, which help us anticipate the general features of the results in the disordered case > studied later. III. LOCAL TIME WITHOUT DISORDER = In this case V x = δx, corresponds to T in Eq. being the local time in the vicinity of the origin and PT t, in Eq. 3 being P loc T t the pdf of the local time T for a given observation time window of size t and the starting position of the particle x =. For our purpose we only need the solution u of the differential equation 4, which corresponds to the starting position of the particle being the origin. However, to obtain u we have to solve Eq. 4 in the entire region of x with the boundary conditions ux ± which are derived from the following observation. If the initial position x ±, the particle can not reach the origin in finite time, which means that the local time T =. Therefore, by substituting PT t, x ± δt in Eq. 3 one obtains the boundary conditions ux ± = α. We have to obtain the solutions ux = u + x for x > and ux = u x for x < by solving Eq. 4 separately in the respective two regions, u ± x + Fxu ± x αu ±x =, with the boundary conditions u + x = /α and u x = /α, and then matching the two solutions u + x and u x at x =. The matching conditions are u + = u = u, u + u = pu. The first condition follows from the fact that the solution must be continuous at x = and the second one is derived by integrating Eq. 4 across x =. By making a constant shift u ± x = /α+a ± y ± x, from Eq. one finds that y ± x satisfy the p-independent homogeneous equation y ±x + Fxy ±x αy ± x = 3 with the boundary conditions y + x and y x. The constants A ± are determined by the matching conditions given in Eq., which can be rewritten as A + y + = A y = u α, A + y + A y = pu. 4b 4a Eliminating the constants A ± from Eq. 4, we obtain the Laplace transform u, defined by Eq. 3 with PT t, P loc T t, as u = = where λα is simply given by λα = z z + dt e αt dt e pt P loc T t λα αp + λα]. 5 with z ± x = y ±x y ± x. 6 Note that putting p = in Eq. 5 and then inverting the Laplace transform with respect to α readily verifies the normalization condition P loc T tdt =. 7 Since λα is independent of p, inverting the Laplace transform in Eq. 5 with respect to p yields Gα = = λα α dt e αt P loc T t exp λαt], 8 which is valid for any arbitrary force Fx. In the following subsections we will consider qualitatively different types of deterministic potentials to derive more explicit results. A. Flat potential We first consider the simple Brownian motion without any external force, Fx =. In this case the solutions of Eq. 3 are obtained as y ± x = y ± exp x ] α. 9 Using the solutions in Eq. 6 one gets λα = α and hence the Laplace transform Gα in Eq. 8 becomes Gα = dt e αt P loc T t = e αt. 3 α Now inverting the Laplace transform with respect to α, one finds that the distribution of the local time is Gaussian for all T and t, P loc T t = πt exp T B. Unstable potential t ]. 3 Now we consider the case of a Brownian particle moving in an unstable potential Ux such that Ux ±. The corresponding repulsive force Fx drives the particle
6 6 eventually either to + or to. The pdf of the local time P loc T t in the case of an unstable potential tends to a steady distribution P loc T as t, which can be computed explicitly. To see this consider the function Gα in Eq. 8. By making a change of variable τ = αt, it follows from Eq. 8, Gα = α dτ P loc T τ. 3 α Assuming P loc T t = P loc T, we find form the above equation that Gα P loc T/α as α. Comparing this behavior with Eq. 8 gives P loc T = λexp λt], 33 provided λ is a finite positive number. Thus generically for all repulsive force Fx, the local time distribution has a universal Poisson distribution in the limit t. The only dependence on the precise form of the force Fx is through the rate constant λ. The rate constant λ can be expressed in terms of the force Fx in a more explicit manner. Putting α = in Eq. 3 and solving the resulting equation with the boundary conditions y + x and y x we get x y + x = y + ψ ydy, x >, 34 ψ ydy x y x = y ψ ydy, x <, 35 ψ ydy where ψy = exp y Fxdx]. Substituting these results in Eq. 6 gives the rate constant as ] λ = ψ ydy ψ ydy Let us now consider a simple example where the potential Ux = µ x with µ >, corresponding to the repulsive force Fx = µ signx from the origin. In this case ψy = exp µ y ] and hence from Eq. 36 we get λ = µ. C. Stable potential We now turn our attention to the complementary situation when the potential Ux is stable, i.e., Ux ±. In this case the force Fx is attractive towards the origin so that the system eventually reaches a well defined stationary state. The stationary probability distribution px for the position of the particle is given by the Gibbs measure px = e Ux, 37 Z where Ux = x Fx dx and Z is the partition function, Z = e Ux dx. 38 In this case the Laplace transform Gα of the pdf of the local time P loc T t is still given by Eq. 8. However, unlike the unstable potential in the previous section, the distribution P loc T t does not approach a steady state as t. Instead it has a rather different asymptotic behavior. To deduce this asymptotic behavior, let us first consider the average local time T = t δxt ] dt. For large t, the average δxt ] approaches its stationary value δxt ] p, where p = /Z from Eq. 37. Hence as t the ratio T/t approaches the limit T t Z, 39 where Z is given by Eq. 38. Thus for large t, the average local time scales linearly with time t, which indicates that the natural scaling limit in this case is when t, T but keeping the ratio r = T/t fixed. We will see that in this scaling limit the local time distribution P loc T t tends to the following asymptotic form ] T P loc T t exp tφ, 4 t where Φr is a large deviation function. To compute the large deviation function we first substitute this presumed asymptotic form of P loc T t given by Eq. 4 in the Laplace transform Gα = e αt P loc T tdt and then make a change of variable r = T/t in the integration. The resulting integral can be evaluated in the large T limit by the method of steepest descent, which gives Gα exp TWα] where Wα = min r {α + Φr}/r]. Comparing this result with Eq. 8 gives min r α + Φr r ] = λα, 4 where λα is given by Eq. 6. Thus λα is just the Legendre transform of Φr. Inversion of this transform gives the exact large deviation function Φr = max α α + rλα], 4 with λα given by Eq. 6. This is a general result valid for any confining potential Ux. We will now explicitly compute the large deviation Φr for the particular potential given by Eq. 5 with µ < and =. Substituting the corresponding force Fx = µ signx in Eq. 3 and solving the resulting differential equations with the boundary conditions y + x and y x we get y ± x = y ± exp µ + ] µ + α x. 43 Substituting these results in Eq. 6 we get λα = µ + µ + α. From Eq. 4 one then gets the large deviation function Φr = r µ. 44
7 7 It turns out that for this particular form of the force Fx = µ signx, the Laplace transform in Eq. 8 can be inverted to get the local time distribution P loc T t exactly for all T and t. The calculations are presented in appendix A. We find that in the large t limit, the distribution reduces to the asymptotic form P loc T t exp t πt ] T t µ, 45 near the mean T = µ t, which verifies the result obtained above by the large deviation function calculation. In fact, the limiting Gaussian form of the distribution of the local time near its mean value is quite generic for any stable potentials where the system eventually becomes ergodic and is just the manifestation of the central limit theorem. From the definition, T T = t {δxt ] δxt ] } dt, it follows that when T T, the random variables δxt ] δxt ] at different times t become only weakly correlated. Then in the limit when t is much larger than the correlation time between these variables, one expects the central limit theorem to hold which predicts a Gaussian form for T near its mean value T. IV. LOCAL TIME WITH DISORDER > So far we have considered the case where the random part of the potential was not present. In this section we will study the effect of the randomness by adding a random part to the potential. In particular, we will consider the diffusive motion of the particle when the force Fx is given by Eq. 6 with >. Equation 8 still remains valid for each realization of the force Fx, i.e., for each realization of {ξx}. Our aim is to compute the average of the pdf of the local time P loc T t over the noise history {ξx}. From Eq. 8, one needs to know the distribution of λα = z z + ]/, which is now a random variable since Fx is random. The variables z + and z are independent of each other and therefore their joint probability distribution factorizes to the individual distributions. The calculations of these distributions are presented in appendix B. Using the distributions of z ± from Eqs. B9 and B respectively with a ± = α, one gets with qt = exp λαt] = ] qt, 46 q w µ/ exp { w + T + α }] dw w 47 α + T = α µ/ + T µ/ K µ/, 48 where K ν x is the modified Bessel function of order ν 44] and K ν x = K ν x. Averaging Eq. 8 over disorder we finally get the exact formula dt e αt P loc T t = d q αq T ]. 49 dt However, it is not an easy task to invert the Laplace transform to get the exact distribution P loc T t for all T and t. In the following subsections we will extract the asymptotic behaviors of P loc T t, for the three cases, when the deterministic part of the potential is: i flat corresponding to µ =, ii unstable corresponding to µ > and iii stable corresponding to µ <. A. Flat potential µ = Sinai model We first consider a particle diffusing in the continuous Sinai potential, i.e., µ = in Eq. 5. Our aim is to find out how this random potential modifies the behavior of the local time. In this case substituting qt and q from Eq. 48 with µ = in Eq. 49 we get the Laplace transform of the disorder averaged local time distribution as dt e αt P loc T t = αk α/ T K α + T. 5 We will now consider the interesting limit where both t and T are large, but the ratio y = T/t is kept fixed. This corresponds to taking the limit α with αt = s keeping fixed. In this limit α K log α 5 and α + T s K K. 5 Therefore substituting t = s/αy and T = s/α in Eq. 5, in the limit α we get dy e s/y s αy P loc s α ] s = αy 4 log α s s K. 53 The above equation suggests that, in the limit t and T, while their ratio T/t is kept fixed, P loc T t should have the scaling form P loc T t = t log t f T/t. 54
8 ½ ݵ 8 Now substituting the above form in Eq. 53 and making the change of variable ỹ = /y, we obtain after straightforward simplification ] dỹ e sỹ f /ỹ s ỹ = 4K. 55 Note that the right hand side of the above equation is simply the Laplace transform of the function f /ỹ/ỹ. Therefore by using the identity e a /4ω ω e sω dω = K a s 56 and the convolution property of Laplace transform, we can invert the Laplace transform in Eq. 55 with respect to s. Inverting the Laplace transform and after simplification we finally get / dx f /ỹ = ỹ x x exp ỹx x Therefore the scaling function f y is simply given by f y = y / dx x x exp By substituting x x = /z gives f y = y 4 y x x ]. 57 ]. 58 dz exp y zz 4 z, 59 where the integral can be evaluated exactly 44], which finally gives the scaling function in Eq. 54 as f y = y e y/ K y/. 6 However the scaling given by Eq. 54 breaks down for very small y very small T when y. The scaling function is displayed in Fig. 3. In the large y limit, using the asymptotic the behavior K ν x π/xe x from Eq. 6 we find that f y π y 3/ e y/. B. Unstable potential µ > In this case the behavior in the limit t can be obtained by either setting α = in the integral in Eq. 47 or taking the α limit in K ν. in Eq. 48, which gives qt Γµ/ µ/ + T µ/, 6 where Γx is the Gamma function 45]. Substituting qt and q in Eq. 49 and inverting the Laplace transform with respect to α gives P loc T t = µ + T µ/, Ý FIG. 3: The scaling function f y in Eq. 54. The solid line is plotted by using Eq. 6, and the dash line is plotted by using the limiting form f y π y 3/ e y/ as y. i.e., in the limit t, the distribution P loc T t tends to a steady state distribution P loc T for all T. The disorder averaged local time distribution has a broad power law distribution even though for each sample the local time has a narrow exponential distribution see Eq.33 in Sec. III B. This indicates wide sample to sample fluctuations and lack of self-averaging. C. Stable potential µ < In this case substituting qt and q from Eq. 48 in Eq. 49 and denoting ν = µ / we get dt e αt P loc T t = αk ν α/ ] α + T + T ν/ K ν. 63 T We consider the scaling limit where both t and T are large, but their ratio y = T/t is kept fixed. This corresponds to taking the limit α with keeping αt = s fixed, which gives the following limiting forms + T s α, 64 α K ν Γν ν, 65 α K ν α + T K ν s. 66
9 9 Substituting the above limits on the right hand side of Eq. 63 and making change of variables t = s/αy and T = s/α on the left hand side, it is straightforward to get dy e s/y s s αy P ] s loc α = αy in the limit α. P loc T t 4 ν Γ ν s s ν/ K ν s ], 67 This suggests the limiting form for P loc T t t f T/t, 68 fy ւ ν < ւ ν > in the scaling limit t and T with a fixed ratio y = T/t. To compute the scaling function we substitute the above scaling form in Eq. 67, and make the change of variable ỹ = /y. Then Eq. 67 simplifies to the Laplace transform dỹ e sỹ f /ỹ ỹ ] = 4 ν Γ ν s ν/ K ν s ], 69 which can be inverted with respect to s, by using the identity a ν e a /4ω ω ν+ e sω dω = s ν/ K ν a s 7 and the convolution property of the Laplace transform. After simplification, the inverse Laplace transform gives / f y = yν dx ν Γ ν x ν+ x ν+ ] y exp. 7 x x By making a change of variable x x = /z in the above integral, it can be presented in the form f y = yν ν Γ ν 4 z ν / z 4 / exp y z dz, 7 which now can be expressed in more elegant forms 45] as ] π y ν f y = e y/ Γ U/, + ν, y/, ν 73 where Ua, b, x is the Confluent Hypergeometric Function of the Second Kind also known as Kummer s function of the second kind 45], which has the following limiting behaviors: U/, + ν, x Γν π x ν for small x, 74a U/, + ν, x x for large x. 74b y/ FIG. 4: The scaling function f y in Eq. 68, plotted by using Eq. 73. ν = µ /. The scaling function f y is displayed in Fig. 4. Using the limiting behaviors from Eq. 74, one finds that the scaling function decays as f y y 4ν 3/ e y/ for large y. For small y, the scaling function behaves as f y y ν, which increases with y for ν >, however, diverges when y for ν <, a behavior qualitatively similar to the Sinai case see Fig. 3. For ν <, the disorder wins over the strength of the stable potential. In that situation when the particle gets trapped in the wells of the random potential, the weak external deterministic force often can not lift it out of the well and send towards the origin. Therefore, the scaling function f y carries very large weight near y = which corresponds to very small local time T for a given observation time t. Note that, for the particular value ν = /, the scaling function has a simple form f y = /πy exp y/. V. INVERSE LOCAL TIME WITHOUT DISORDER = The inverse local time means how long one has to observe the particle until the total time spent in the infinitesimal neighborhood of the origin is T. The double Laplace transform of the pdf of the inverse local time is obtained by simply putting x = in Eq. 9. Corresponding u in Eq. 9, which is nothing but the double Laplace transform of the pdf of local time, has already been evaluated in Sec. III and is given by Eq. 5. Substituting u and replacing It T, with the pdf of the inverse local time I loc t T, after straightforward simplification, for x = Eq. 9 reads dt e αt dt e pt I loc t T = p + λα, 75 where λα is given by Eq. 6, which depends on the force Fx through Eq. 3. Inverting the Laplace transform with
10 T I loct T STABLE POTENTIAL µ < FLAT POTENTIAL µ = it decays exponentially I loc t T exp µ t/. On the other hand when the potential is unstable, µ >, as we see from Eq. 79, the distribution I loc t T is not normalized to unity. In this case the particle escapes to ± with probability e µt and Eq. 78 gives the distribution only for those events where the particle does not escape to ±. Therefore for µ >, it is appropriate to represent the full normalized distribution as I loc t T = T ] exp T + µt πt 3 t + e µt δt. 8 UNSTABLE POTENTIAL µ > ւ Note that the second term does not show up in the Laplace transform of I loc t T with respect to t t/t VI. INVERSE LOCAL TIME WITH DISORDER > FIG. 5: The pdfs of the inverse local time for stable µ = /, flat µ = and unstable µ = / potentials, plotted using Eq. 78 and T =. respect to p gives the general formula dt e αt I loc t T = exp λαt], 76 valid for arbitrary force Fx, a result known in the mathematics literature 7, 46]. We first consider the pure case where the force given by Eq. 6 with =. Substituting solutions of Eq. 3 for Fx = µ signx in Eq. 6 we obtain λα = µ + µ + α. Now using this λα in Eq. 76 and making a change of the parameter α = β µ / we get ] dt e βt e µt/ I loc t T = e µt e βt, 77 where the right hand side is simply the Laplace transform of e µt/ I loc t T with respect to t. The Laplace transform can be inverted to obtain the exact pdf of the inverse local time I loc t T = T exp πt 3 T + µt t ]. 78 with the normalization condition { for µ, I loc t Tdt = e µ+ µ T = e µt for µ >, 79 which is simply obtained by putting α = in Eq. 76. As we infer from Eq. 78, although in the limit t the inverse local time distribution I loc t T exp T /t is independent of µ, for large t it depends on the nature of the potential, as shown in Fig 5. While in the absence of any force, i.e., µ = the inverse local time distribution has a power-law tail I loc t T t 3/, for the stable potential, i.e., µ <, In this section, we switch on the disorder by considering > in the force given by Eq. 6. In the presence of disorder, taking the disorder average of Eq. 76 gives dt e αt I loc t T = exp λαt], 8 with λα = z z + ]/, where z + and z are independent random variables, whose distributions are given by Eqs. B9 and B respectively with a ± = α. The object exp λαt] on the right hand side of Eq. 8, has already been evaluated in Sec. IV, which is given by Eq. 46. In the following subsections we will determine the behavior of I loc t T in the scaling limit t, T, while keeping their ratio x = t/t fixed, for the three qualitatively different cases: i µ =, ii µ > and iii µ <. A. Flat potential µ = Sinai model Following the analysis of Sec. IV A, in the limit α with keeping αt = s fixed, exp λαt] 4 s log α K. 8 Therefore, substituting t = xs/α and T = s/α, in the limit α, Eq. 8 reads s sx dxe sx α I loc s ] = 4 s α α log α K. 83 This suggest the scaling form I loc t T = T log T g t/t, 84 in the limit t, T but keeping x = t/t fixed. Substituting this scaling form in Eq. 83, after straightforward
11 ½ ܵ FIG. 6: The scaling function g y in Eq. 84. The solid line is plotted by using Eq. 86, and the dash line is plotted by using the limiting forms: g x πx / exp /x for small x and g x logx/x for large x. simplification one obtains Ü dxe sx g x = 4K s. 85 Now direct comparison of the above equation with Eq. 55 gives g x = f /x/x, where f x is given by Eq. 6. Substituting f /x one obtains the scaling function g x as g x = x e /x K /x, 86 which is displayed in Fig. 6. The scaling function increases as g x πx / exp /x for small x and decays as g x logx/x at large x. B. Unstable potential µ > In this case the right hand side of Eq. 8 is given by α + T/ exp λαt] = +T ν K ν K ν α/, 87 with ν = µ/. Putting α = in the above equation gives the normalization condition, I loc t Tdt = + T ν, which implies that for the unstable potential, where the force is repulsive from the origin, the particle escapes to ± with probability + T ν, and the disorder averaged pdf I loc t T obtained by inverting the Laplace transform in Eq. 8 represents only those events where the particle does not escape to ±. Now in the limit of α with αt = s keeping fixed, one gets α K ν Γν α + T K ν Therefore Eq. 87 becomes exp λαt] T ν/ s ν, 88 s K ν. 89 ] ν 4T s 3 ν Γ s ν/ K ν. 9 ν In the corresponding limit T, t, but keeping their ratio x = t/t fixed, using the scaling form I loc t T = T ν+ g t/t 9 in Eq. 8 one finally arrives at the Laplace transform ] e sx 4 s g xdx = 3 ν Γ s ν/ K ν. ν 9 The Laplace transform can be inverted with respect to s to obtain the scaling function g x and in fact the inversion has already been done in Sec. IV C. Comparing the above equation with Eq. 69 readily gives g x = ν f /x/x where f x is given by Eq. 73. Substituting f /x gives g x = π ] e /x ν Γ ν U/, + ν, /x, 93 x ν+ where Ua, b, x is the Confluent Hypergeometric Function of the Second Kind, whose its small and large x behaviors are given in Eq. 74. The scaling function g x is displayed in Fig. 7. The scaling function increases as g x exp /x for small x and eventually decreases for large x as g x /x ν. In particular, for ν = / it has a very simple form g x = /π 3 x 3/ exp /x. C. Stable potential µ < Following the analysis of Sec. IV C, in the limit α, keeping αt = s fixed one gets exp λαt] = ] 4 s ν Γ s ν/ K ν, 94 ν with ν = µ /. On the other hand, in the corresponding limit T, t, but keeping t/t = x fixed, using the scaling form I loc t T = T g 3t/T, 95
12 ν gx ν = 3 4 x ւ ν = FIG. 7: The scaling functions g x in Eq. 9 plotted by using Eq. 93. ν = µ /. one gets dt e αt I loc t T = dxe sx g 3 x, 96 with s = αt. Therefore, in this scaling limit Eq. 8 becomes ] dxe sx 4 s g 3 x = ν Γ s ν/ K ν. ν 97 Now comparing the above equation with Eq. 9 one gets g 3 x = ν g x, 98 where the scaling function g x is given by Eq. 93 and displayed in Fig. 7. While I loc t T has the same scaling function up to a multiplicative factor of ν for both stable and unstable potential, the physical behaviors, however, are quite different. For the stable potential, I loc t T is normalized to unity. Note that the scaling function g 3 x becomes narrower as one increases ν, as expected since the particle becomes more localized near the origin. For the unstable potential, on the other hand, the weight of I loc t T decreases as T ν, as one increases ν, as expected since when the repulsive force from the origin becomes stronger, the particle escapes to ± with a higher probability. VII. OCCUPATION TIME WITHOUT DISORDER = In this case V x = θx, corresponds to T in Eq. being the occupation time in the region x > and PT t, in Eq. 3 being P occ T t the pdf of the occupation time for a given observation time window of size t and the initial position of the particle x =. Again as before, we need to solve the differential equation 4 for x > and x < separately and then match the solutions at x =. The matching condition for the slope of the solutions is obtained by integrating Eq. 4 across x =. Thus the matching conditions are u + = u = u, and u + = u, 99 where u ± x satisfy the following differential equations: u +x + Fxu +x α + pu + x =, for x >, and u x + Fxu x αu x =, for x <. The boundary conditions of u ± x when x ± are obtained from the fact that if the starting position goes to ±, the particle will never cross the origin in finite time: PT t, x = δt T and PT t, x = δt, and hence from Eq. 3 u + = α + p, and u = α. Writing u + x = /α + p + B + y + x and u x = /α + B y x, we obtain the homogeneous differential equations for y ± x as y + x + Fxy + x α + py +x =, 3 for x >, and y x + Fxy x αy x =, 4 for x <, with the boundary conditions y + x = and y x =. The constants B ± are determined by the matching conditions given in Eq. 99, which can be rewritten as α + p + B +y + = α + B y = u, B + y + = B y. 5a 5b Eliminating the constants from Eq. 5, we obtain the double Laplace transform of the pdf of the occupation time where u = = l α, p α t dt e αt dt e pt P occ T t + l α, p α + p, 6 ] z l α, p =, 7 z z + ] z + l α, p =, 8 z z +
13 3 and z ± x = y ± x/y ±x. Note that l α, p + l α, p =. 9 Putting p =, in Eq. 6 gives u = /α, and hence inverting the Laplace transform with respect to α readily verifies the normalization t P occ T tdt =. For any symmetric deterministic potential the distribution of the occupation time is symmetric about its mean T = t/, i.e., P occ T t = P occ t T t. Then, it follows from this symmetry that l α + p, p = l α, p. In other words, the double-integral in Eq. 6 remains invariant under the following simultaneous replacements: α + p α and α α + p. Thus under these replacements one must have z + z and vice versa, which also implies that z + = z for p =. Equivalently, l α, = l α, = /, which also directly follows from Eqs. 9 and. Therefore if one splits the distribution function into two parts: P occ T t = R L T t + R R T t such that dt e αt t dt e αt t dt e pt R L T t = l α, p, α dt e pt R R T t = l α, p α + p, 3 then it follows from the above discussion that R L t T t = R R T t. This symmetry of the distribution will come handy later. Moreover, putting p = and inverting the Laplace transforms with respect to α gives the normalization for each part separately t R L T tdt = t R R T tdt =. 4 As an example, we first consider the pure case: = in the force given by Eq. 6. For Fx = µ signx, the solutions of Eqs. 3 and 4 are obtained as y + x = y + exp for x >, and, y x = y exp µ + µ + α + p ] x, 5 µ + ] µ + α x, 6 for x <. These give the expressions for z ± = y ±/y ± as z + = µ + ] µ + α + p, z = µ + ] 7 µ + α. In the following subsections we will consider the three different cases: i µ =, ii µ > and iii µ <. A. Flat potential µ = For µ =, using z + = α + p and z = α from Eq. 6 we get t dt e αt dt e pt P occ T t = αα + p. 8 Inverting the double Laplace transform with respect to p and then with respect to α finally reproduces the well known Lévy s arcsine law ] for the pdf of the occupation time of an ordinary Brownian motion, P occ T t = π, < T < t. 9 Tt T The distribution P occ T t diverges on both ends T = and T = t, which indicates that the Brownian particle tends to stay on one side of the origin. B. Unstable potential µ > Since for µ >, the force is repulsive from the origin x =, one would expect the occupation time distribution to be convex concave upward, with minimum at T = t/. Now in the limit of large window size t, the part of the distribution P occ T t to the left of the midpoint T = t/ approaches R L T t, as the midpoint itself goes to. By making a change of variable z = αt, it follows from Eq. z/α dz e z dt e pt R L T z/α = l α, p. Now the large t limit of R L T t can be obtained by taking α in the above equation, where one realizes that R L T t approaches a steady t independent distribution, R L T t R L T, whose Laplace transform is given by dt e pt R L T = l, p, where l, p is obtained from Eq. 7, by using z ± from Eq. 7, which gives l, p = µ 3µ + µ + p. The above Laplace transform can be inverted with respect to p, which gives R L T = µ e µ T/ 3µ exp πt 9µ T ] 3µ erfc T. 3 with the normalization R L TdT = l, = /.
14 4 The limiting behavior of the distribution is given by R L T µ πt, 4 for small T and decays exponentially for large T, R L T 9µ e µ T/. 5 π T 3/ C. Stable potential µ < As we discussed earlier in Sec. III C in the context of the local time, for generic stable potential Ux the system eventually becomes ergodic at large t and hence the average θxt] approaches its stationary value θxt] Z + /Z, where Z = e Ux dx is the equilibrium partition function and Z + = e Ux dx is the restricted partition function. Therefore, for large t the average occupation time T = t θxt ] dt scales linearly with t Z+ T t. 6 Z Note that when the potential Ux is symmetric about zero, the average occupation time T = t/ for all t. From the definition, T T = t {θxt ] θxt ] } dt, it follows that when T T, the random variables θxt ] θxt ] at different times t become only weekly correlated. Then in the limit when t is much larger than the correlation time between these variables, one expects the central limit theorem to hold, which predicts a Gaussian form for the distribution of the occupation time T near the mean value T, ] T T P occ T t exp, 7 where the variance = T T can be obtained from the Laplace transform of the moments T n e αt dt = n n u p n, 8 p= with u given by Eq. 6. For the particular attractive force Fx = µ signx, using z ± from Eq. 7 in Eq. 6 and taking derivatives with respect to p in we get u p u p = p= α, 9 = p= α 3 + 4µ α + O. 3 α Therefore inverting the Laplace transform in Eq. 8 with respect to α immediately gives T = t/ for all t, and T = t /4 + t/4µ for large t which gives = t/4µ. VIII. OCCUPATION TIME WITH DISORDER > Now we consider the occupation time when the disorder is switched on: > in Eq. 6. Our aim is to calculate the disorder averaged P occ T t. As one realizes from Eqs. 6, 7 and 8, to calculate P occ T t one needs the distribution of z + and z, which are given by Eqs. B9 and B with a + = α + p and a = α respectively. In the following subsections, we will consider the three cases: i µ =, ii µ > and iii µ <. A. Flat potential µ = Sinai model We first consider the diffusive motion of a particle in a continuous Sinai potential, where the potential itself is a Brownian motion in space. In the limit of large window size t the left half of the disorder averaged pdf of the occupation time R L T t for T t/ is obtained by taking the disorder average in Eq.. The right half of the distribution for t/ T t is just the symmetric reflection of the left part. The detailed calculations for R L T t are presented in appendix C. We find that R L T t has a large t behavior R L T t RT, 3 log t where the function RT is independent of t. The limiting behaviors of RT are given by RT, 3 πt as T and for large T. RT T, 33 B. Unstable potential µ > For µ >, we find that disorder does not change the asymptotic behavior of the distribution for the pure case qualitatively. The calculations are presented in appendix D. We find that in the limit t the left half of the disorder averaged occupation time distribution tends to a t independent form R L T t = R L T. 34 In fact the small T limit of R L T remains same as in the pure case R L T µ πt. 35 For large T, the distribution R L T still decays exponentially R L T e bt, 36
15 Ó Ýµ 5 where the decay coefficient b is, however, different from the pure case see Eq. D7. C. Stable potential µ <.5 ½ This particular situation, where one finds the interplay between two competing processes, is a very interesting one. On one hand, as we discussed in Sec. VII C, the stable potential in the absence of the disordered potential makes the system ergodic in the large t limit, and as a result the pdf of the occupation time is peaked at = t/ and decays fast away from it. On the other hand, as we discussed in Sec. VIII A, without any underlying deterministic potential the disorder averaged pdf of the occupation time is convex concave upward with a minimum at T = t/ and diverges at the both ends T and T t. Therefore, if both the stable potential and disordered potential are included, as their relative strength ν = µ / is varied, one expects a phase transition at some critical value ν c where the system looses ergodicity. In the scaling limit where both t and T, but their ratio y = T/t is kept fixed, we find that the disorder averaged pdf of the occupation time has a scaling form P occ T t = t f ot/t. 37 The calculation of the scaling function f o y is presented in appendix E, where we find the Beta law f o y = Bν, ν y y]ν, y, 38 where ν = µ / and Bν, ν is the Beta function 44]. Now if one tunes the parameter ν by varying either µ or the disorder strength, the distribution P occ T t exhibits a phase transition in the ergodicity of the particle position at ν c = Fig. 8. For ν < ν c the distribution f o y in Eq. 37 is convex with a minimum at y = / and diverges at the two ends y =,. This means that particle tends to stay on one side of the origin such that T is close to either or t. In other words the paths with small number of zero crossings carry more weight than the ones that cross many times. For ν > ν c the scenario is exactly opposite, where f o y is maximum at the mean value y = / indicating that particle tends to spend equal times on both sides of the origin x =, such that paths with large number of zero crossings, for which T is closer to t/ carry larger weight. Similar phase transition in the ergodicity properties of a stochastic process as one changes a parameter, was first noted in the context of diffusion equation 9], and later found for a class of Gaussian Markov processes ] and in simple models of coarsening 47, 48]. A very interesting observation about Eq. 38 is that for ν = /, the result is same as Lévy s result for the onedimensional Brownian motion given by Eq. 9. It seems as if the attractive force cancels the effect of disorder exactly at ν = /. However, this is no more true in the context of the local time..5 ½.5 FIG. 8: The scaling functions f ox in Eq. 37 plotted by using Eq. 38. IX. INVERSE OCCUPATION TIME WITHOUT DISORDER = In this case It T, in Eq. 9 is replaced with I occ t T, which is the distribution of the time t needed to observe the particle with a starting position x =, until the total amount of time spent on the positive side x > is T. Corresponding u in Eq. 9 with x =, which is the double Laplace transform of the pdf of the occupation time, has already been evaluated in Sec. VII, which is given by Eq. 6. Substituting u in Eq. 9 gives dt e pt T Ý dt e αt I occ t T = l α, p α + p, 39 where l α, p is given by Eq. 8. Comparing the above equation with Eq. 3, one can infer that I occ t T and R R T t have the same functional form, i.e., I occ t T = R R T t and especially for symmetric deterministic potential I occ t T = R R T t = R L t T t. It is useful to present the above equation in the following form, dz e z dτ e ατ I occ τ + z z β = l β, α β, β 4 which has been obtained by substituting p = β α in Eq. 39 and subsequently making the change of variables βt = z and τ = t T. On the right hand side, we have substituted l α, β α = l β, α β, using Eq. and l α, p is given by Eq. 7. Now by taking the limit β in Eq. 4, one obtains the large T behavior of I occ t T. For the pure case, = in Eq. 6, we have already obtained z ± in Sec. VII, which are given by Eq. 7 and hence we can evaluate l α, p and l α, p by using
16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationWeak Ergodicity Breaking WCHAOS 2011
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University Bel, Burov, Korabel, Margolin, Rebenshtok WCHAOS 211 Outline Single molecule experiments exhibit weak ergodicity breaking. Blinking quantum dots,
More informationWeak Ergodicity Breaking
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University 215 Ergodicity Ergodicity: time averages = ensemble averages. x = lim t t x(τ)dτ t. x = xp eq (x)dx. Ergodicity out of equilibrium δ 2 (, t) = t
More informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationMaximal distance travelled by N vicious walkers till their survival arxiv: v1 [cond-mat.stat-mech] 17 Feb 2014
Noname manuscript No. will be inserted by the editor Anupam Kundu Satya N. Majumdar Grégory Schehr Maximal distance travelled by N vicious walkers till their survival arxiv:4.3964v [cond-mat.stat-mech]
More informationAnomalous Lévy diffusion: From the flight of an albatross to optical lattices. Eric Lutz Abteilung für Quantenphysik, Universität Ulm
Anomalous Lévy diffusion: From the flight of an albatross to optical lattices Eric Lutz Abteilung für Quantenphysik, Universität Ulm Outline 1 Lévy distributions Broad distributions Central limit theorem
More informationKolmogorov Equations and Markov Processes
Kolmogorov Equations and Markov Processes May 3, 013 1 Transition measures and functions Consider a stochastic process {X(t)} t 0 whose state space is a product of intervals contained in R n. We define
More informationCHAPTER V. Brownian motion. V.1 Langevin dynamics
CHAPTER V Brownian motion In this chapter, we study the very general paradigm provided by Brownian motion. Originally, this motion is that a heavy particle, called Brownian particle, immersed in a fluid
More informationTheory of fractional Lévy diffusion of cold atoms in optical lattices
Theory of fractional Lévy diffusion of cold atoms in optical lattices, Erez Aghion, David Kessler Bar-Ilan Univ. PRL, 108 230602 (2012) PRX, 4 011022 (2014) Fractional Calculus, Leibniz (1695) L Hospital:
More informationUncertainty quantification and systemic risk
Uncertainty quantification and systemic risk Josselin Garnier (Université Paris Diderot) with George Papanicolaou and Tzu-Wei Yang (Stanford University) February 3, 2016 Modeling systemic risk We consider
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 23 Feb 1998
Multiplicative processes and power laws arxiv:cond-mat/978231v2 [cond-mat.stat-mech] 23 Feb 1998 Didier Sornette 1,2 1 Laboratoire de Physique de la Matière Condensée, CNRS UMR6632 Université des Sciences,
More informationStationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos
N Cufaro Petroni: CYCLOTRONS 2007 Giardini Naxos, 1 5 October, 2007 1 Stationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos CYCLOTRONS 2007 Giardini Naxos, 1 5 October Nicola
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationSmoluchowski Diffusion Equation
Chapter 4 Smoluchowski Diffusion Equation Contents 4. Derivation of the Smoluchoswki Diffusion Equation for Potential Fields 64 4.2 One-DimensionalDiffusoninaLinearPotential... 67 4.2. Diffusion in an
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 1 Aug 2001
Financial Market Dynamics arxiv:cond-mat/0108017v1 [cond-mat.stat-mech] 1 Aug 2001 Fredrick Michael and M.D. Johnson Department of Physics, University of Central Florida, Orlando, FL 32816-2385 (May 23,
More informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
More informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 2 Dec 2004
Airy Distribution Function: From the Area Under a Brownian Excursion to the Maximal Height of Fluctuating Interfaces Satya N. Majumdar 1, and Alain Comtet,3 1 aboratoire de Physique Théorique (UMR C515
More informationODE Homework Solutions of Linear Homogeneous Equations; the Wronskian
ODE Homework 3 3.. Solutions of Linear Homogeneous Equations; the Wronskian 1. Verify that the functions y 1 (t = e t and y (t = te t are solutions of the differential equation y y + y = 0 Do they constitute
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationPerturbed Gaussian Copula
Perturbed Gaussian Copula Jean-Pierre Fouque Xianwen Zhou August 8, 6 Abstract Gaussian copula is by far the most popular copula used in the financial industry in default dependency modeling. However,
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationWeak Ergodicity Breaking. Manchester 2016
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University Burov, Froemberg, Garini, Metzler PCCP 16 (44), 24128 (2014) Akimoto, Saito Manchester 2016 Outline Experiments: anomalous diffusion of single molecules
More informationMath 308 Exam I Practice Problems
Math 308 Exam I Practice Problems This review should not be used as your sole source of preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More informationNew Physical Principle for Monte-Carlo simulations
EJTP 6, No. 21 (2009) 9 20 Electronic Journal of Theoretical Physics New Physical Principle for Monte-Carlo simulations Michail Zak Jet Propulsion Laboratory California Institute of Technology, Advance
More informationFokker-Planck Equation with Detailed Balance
Appendix E Fokker-Planck Equation with Detailed Balance A stochastic process is simply a function of two variables, one is the time, the other is a stochastic variable X, defined by specifying: a: the
More informationThe Kramers problem and first passage times.
Chapter 8 The Kramers problem and first passage times. The Kramers problem is to find the rate at which a Brownian particle escapes from a potential well over a potential barrier. One method of attack
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationBrownian motion and the Central Limit Theorem
Brownian motion and the Central Limit Theorem Amir Bar January 4, 3 Based on Shang-Keng Ma, Statistical Mechanics, sections.,.7 and the course s notes section 6. Introduction In this tutorial we shall
More information1 Elementary probability
1 Elementary probability Problem 1.1 (*) A coin is thrown several times. Find the probability, that at the n-th experiment: (a) Head appears for the first time (b) Head and Tail have appeared equal number
More informationVIII.B Equilibrium Dynamics of a Field
VIII.B Equilibrium Dynamics of a Field The next step is to generalize the Langevin formalism to a collection of degrees of freedom, most conveniently described by a continuous field. Let us consider the
More informationF n = F n 1 + F n 2. F(z) = n= z z 2. (A.3)
Appendix A MATTERS OF TECHNIQUE A.1 Transform Methods Laplace transforms for continuum systems Generating function for discrete systems This method is demonstrated for the Fibonacci sequence F n = F n
More informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationAn efficient approach to stochastic optimal control. Bert Kappen SNN Radboud University Nijmegen the Netherlands
An efficient approach to stochastic optimal control Bert Kappen SNN Radboud University Nijmegen the Netherlands Bert Kappen Examples of control tasks Motor control Bert Kappen Pascal workshop, 27-29 May
More information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationPhysics 116C The Distribution of the Sum of Random Variables
Physics 116C The Distribution of the Sum of Random Variables Peter Young (Dated: December 2, 2013) Consider a random variable with a probability distribution P(x). The distribution is normalized, i.e.
More informationStatistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
More informationA new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009
A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationDiffusion in a Logarithmic Potential Anomalies, Aging & Ergodicity Breaking
Diffusion in a Logarithmic Potential Anomalies, Aging & Ergodicity Breaking David Kessler Bar-Ilan Univ. E. Barkai (Bar-Ilan) A. Dechant (Augsburg) E. Lutz (Augsburg) PRL, 105 120602 (2010) arxiv:1105.5496
More information9 More on the 1D Heat Equation
9 More on the D Heat Equation 9. Heat equation on the line with sources: Duhamel s principle Theorem: Consider the Cauchy problem = D 2 u + F (x, t), on x t x 2 u(x, ) = f(x) for x < () where f
More informationI forgot to mention last time: in the Ito formula for two standard processes, putting
I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy
More informationSemiclassical analysis and a new result for Poisson - Lévy excursion measures
E l e c t r o n i c J o u r n a l o f P r o b a b i l i t y Vol 3 8, Paper no 45, pages 83 36 Journal URL http://wwwmathwashingtonedu/~ejpecp/ Semiclassical analysis and a new result for Poisson - Lévy
More informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
More informationarxiv: v1 [cond-mat.stat-mech] 31 Mar 2008
Real Roots of Random Polynomials and Zero Crossing Properties of Diffusion Equation. arxiv:83.4396v [cond-mat.stat-mech] 3 Mar 28 Grégory Schehr and Satya N. Majumdar 2 Laboratoire de Physique Théorique
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More informationMATH3203 Lecture 1 Mathematical Modelling and ODEs
MATH3203 Lecture 1 Mathematical Modelling and ODEs Dion Weatherley Earth Systems Science Computational Centre, University of Queensland February 27, 2006 Abstract Contents 1 Mathematical Modelling 2 1.1
More informationENGI 9420 Lecture Notes 1 - ODEs Page 1.01
ENGI 940 Lecture Notes - ODEs Page.0. Ordinary Differential Equations An equation involving a function of one independent variable and the derivative(s) of that function is an ordinary differential equation
More informationUNDERSTANDING BOLTZMANN S ANALYSIS VIA. Contents SOLVABLE MODELS
UNDERSTANDING BOLTZMANN S ANALYSIS VIA Contents SOLVABLE MODELS 1 Kac ring model 2 1.1 Microstates............................ 3 1.2 Macrostates............................ 6 1.3 Boltzmann s entropy.......................
More informationBasic Concepts and Tools in Statistical Physics
Chapter 1 Basic Concepts and Tools in Statistical Physics 1.1 Introduction Statistical mechanics provides general methods to study properties of systems composed of a large number of particles. It establishes
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationLANGEVIN THEORY OF BROWNIAN MOTION. Contents. 1 Langevin theory. 1 Langevin theory 1. 2 The Ornstein-Uhlenbeck process 8
Contents LANGEVIN THEORY OF BROWNIAN MOTION 1 Langevin theory 1 2 The Ornstein-Uhlenbeck process 8 1 Langevin theory Einstein (as well as Smoluchowski) was well aware that the theory of Brownian motion
More informationHW2 Solutions. MATH 20D Fall 2013 Prof: Sun Hui TA: Zezhou Zhang (David) October 14, Checklist: Section 2.6: 1, 3, 6, 8, 10, 15, [20, 22]
HW2 Solutions MATH 20D Fall 2013 Prof: Sun Hui TA: Zezhou Zhang (David) October 14, 2013 Checklist: Section 2.6: 1, 3, 6, 8, 10, 15, [20, 22] Section 3.1: 1, 2, 3, 9, 16, 18, 20, 23 Section 3.2: 1, 2,
More informationPath integrals for classical Markov processes
Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field
More informationLANGEVIN EQUATION AND THERMODYNAMICS
LANGEVIN EQUATION AND THERMODYNAMICS RELATING STOCHASTIC DYNAMICS WITH THERMODYNAMIC LAWS November 10, 2017 1 / 20 MOTIVATION There are at least three levels of description of classical dynamics: thermodynamic,
More informationReal roots of random polynomials and zero crossing properties of diffusion equation
Real roots of random polynomials and zero crossing properties of diffusion equation Grégory Schehr Laboratoire de Physique Théorique et Modèles Statistiques Orsay, Université Paris XI G. S., S. N. Majumdar,
More informationDynamics of a tagged monomer: Effects of elastic pinning and harmonic absorption. Shamik Gupta
: Effects of elastic pinning and harmonic absorption Laboratoire de Physique Théorique et Modèles Statistiques, Université Paris-Sud, France Joint work with Alberto Rosso Christophe Texier Ref.: Phys.
More informationLet's transfer our results for conditional probability for events into conditional probabilities for random variables.
Kolmogorov/Smoluchowski equation approach to Brownian motion Tuesday, February 12, 2013 1:53 PM Readings: Gardiner, Secs. 1.2, 3.8.1, 3.8.2 Einstein Homework 1 due February 22. Conditional probability
More informatione (x y)2 /4kt φ(y) dy, for t > 0. (4)
Math 24A October 26, 2 Viktor Grigoryan Heat equation: interpretation of the solution Last time we considered the IVP for the heat equation on the whole line { ut ku xx = ( < x
More information1. Stochastic Process
HETERGENEITY IN QUANTITATIVE MACROECONOMICS @ TSE OCTOBER 17, 216 STOCHASTIC CALCULUS BASICS SANG YOON (TIM) LEE Very simple notes (need to add references). It is NOT meant to be a substitute for a real
More information1 2 predators competing for 1 prey
1 2 predators competing for 1 prey I consider here the equations for two predator species competing for 1 prey species The equations of the system are H (t) = rh(1 H K ) a1hp1 1+a a 2HP 2 1T 1H 1 + a 2
More information5 Applying the Fokker-Planck equation
5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p
More informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More informationHandbook of Stochastic Methods
C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical
More informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
More informationThe Laplace Transform
C H A P T E R 6 The Laplace Transform Many practical engineering problems involve mechanical or electrical systems acted on by discontinuous or impulsive forcing terms. For such problems the methods described
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationMath Applied Differential Equations
Math 256 - Applied Differential Equations Notes Existence and Uniqueness The following theorem gives sufficient conditions for the existence and uniqueness of a solution to the IVP for first order nonlinear
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More information1 Review of di erential calculus
Review of di erential calculus This chapter presents the main elements of di erential calculus needed in probability theory. Often, students taking a course on probability theory have problems with concepts
More informationNonconstant Coefficients
Chapter 7 Nonconstant Coefficients We return to second-order linear ODEs, but with nonconstant coefficients. That is, we consider (7.1) y + p(t)y + q(t)y = 0, with not both p(t) and q(t) constant. The
More informationMath 308 Exam I Practice Problems
Math 308 Exam I Practice Problems This review should not be used as your sole source for preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationDynamical Systems. August 13, 2013
Dynamical Systems Joshua Wilde, revised by Isabel Tecu, Takeshi Suzuki and María José Boccardi August 13, 2013 Dynamical Systems are systems, described by one or more equations, that evolve over time.
More information6. Brownian Motion. Q(A) = P [ ω : x(, ω) A )
6. Brownian Motion. stochastic process can be thought of in one of many equivalent ways. We can begin with an underlying probability space (Ω, Σ, P) and a real valued stochastic process can be defined
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More informationlaplace s method for ordinary differential equations
Physics 24 Spring 217 laplace s method for ordinary differential equations lecture notes, spring semester 217 http://www.phys.uconn.edu/ rozman/ourses/p24_17s/ Last modified: May 19, 217 It is possible
More informationIntroduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!
Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationBMIR Lecture Series on Probability and Statistics Fall 2015 Discrete RVs
Lecture #7 BMIR Lecture Series on Probability and Statistics Fall 2015 Department of Biomedical Engineering and Environmental Sciences National Tsing Hua University 7.1 Function of Single Variable Theorem
More informationLarge deviations of the top eigenvalue of random matrices and applications in statistical physics
Large deviations of the top eigenvalue of random matrices and applications in statistical physics Grégory Schehr LPTMS, CNRS-Université Paris-Sud XI Journées de Physique Statistique Paris, January 29-30,
More informationErrata for Stochastic Calculus for Finance II Continuous-Time Models September 2006
1 Errata for Stochastic Calculus for Finance II Continuous-Time Models September 6 Page 6, lines 1, 3 and 7 from bottom. eplace A n,m by S n,m. Page 1, line 1. After Borel measurable. insert the sentence
More informationMath 215/255 Final Exam (Dec 2005)
Exam (Dec 2005) Last Student #: First name: Signature: Circle your section #: Burggraf=0, Peterson=02, Khadra=03, Burghelea=04, Li=05 I have read and understood the instructions below: Please sign: Instructions:.
More informationStochastic continuity equation and related processes
Stochastic continuity equation and related processes Gabriele Bassi c Armando Bazzani a Helmut Mais b Giorgio Turchetti a a Dept. of Physics Univ. of Bologna, INFN Sezione di Bologna, ITALY b DESY, Hamburg,
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationarxiv: v1 [math.pr] 18 Nov 2018
Integral Equation Approach to Stationary Stochastic Counting Process with Independent Increments arxiv:1811.7262v1 [math.pr] 18 Nov 218 Enzhi Li Suning R & D Center, Palo Alto, USA November 2, 218 Abstract:
More informationStochastic (intermittent) Spikes and Strong Noise Limit of SDEs.
Stochastic (intermittent) Spikes and Strong Noise Limit of SDEs. D. Bernard in collaboration with M. Bauer and (in part) A. Tilloy. IPAM-UCLA, Jan 2017. 1 Strong Noise Limit of (some) SDEs Stochastic differential
More informationMath 216 Final Exam 24 April, 2017
Math 216 Final Exam 24 April, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationExample 4.1 Let X be a random variable and f(t) a given function of time. Then. Y (t) = f(t)x. Y (t) = X sin(ωt + δ)
Chapter 4 Stochastic Processes 4. Definition In the previous chapter we studied random variables as functions on a sample space X(ω), ω Ω, without regard to how these might depend on parameters. We now
More informationHJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011
Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance
More informationLecture 4: Introduction to stochastic processes and stochastic calculus
Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London
More information