Multidimensional Stationary Probability Distribution for Interacting Active Particles, Supplemental Material
|
|
- Dwayne Francis
- 5 years ago
- Views:
Transcription
1 Multidimensional Stationary Probability istribution for Interacting Active Particles, Supplemental Material Claudio Maggi 1, a Umberto Marini Bettolo Marconi 2, Nicoletta Gnan 3, and Roberto i Leonardo 1,4 1 ipartimento di Fisica, Università di Roma Sapienza, I-00185, Roma, Italy 2 Scuola di Scienze e Tecnologie, Università di Camerino, Via Madonna delle Carceri, 62032, Camerino, INFN Perugia, Italy 3 CNR-ISC, UOS Sapienza, P.le A. Moro2, I-00185, Roma, Italy and 4 CNR-IPCF, UOS Roma, ipartimento di Fisica Università Sapienza, I Roma, Italy I. PROOF OF THE MAIN RESULT Let us consider the set of stochastic differential equations: ẋ i = xi Φ + η i (1) where i = 1,.., N, xi = / x i and the η i -s are a set of independent Gaussian stochastic processes with zero mean η j = 0 and exponential time-correlation: η i (t)η j (s) = δ ij e t s / The Gaussian colored-noise component η i is necessarily produced by the Ornstein-Uhlenbec process (oob s theorem) 15]: η i = η i + 1/2 Γ i (2) where the Γ i -s are a set of uncorrelated white-noise processes: Γ i (s)γ j (t) = 2δ ij δ(t s) with zero average Γ i = 0. By combining Eq.s (1) and (2) we find 15, 17]: ẍ i + x i Φ + ( δi ) + x ix Φ ẋ = 1/2 Γ i (3) The Unified Colored Noise Approximation (UCNA) is done by neglecting the term ẍ i in Eq. (4). With this simplification we can invert Eq. (4) to obtain ẋ i = M 1 ] i x Φ + 1/2 M 1 ] i Γ (4) where we have defined the matrix M ij = δ ij / + xix j Φ (which is symmetric by construction being xix Φ the Hessian matrix) and its inverse M 1 ] ij. The Foer-Planc equation associated with the (Stratonovich) Langevin equation (4) is given by 18]: a Electronic address: claudio.maggi@roma1.infn.it
2 2 Ω = l,i { } 1 xl M 1 ] li ( xi Φ) Ω + { 2 xl M 1 ( ] li xj M 1 ] ij Ω )} (5) l,i,j where Ω = Ω(x 1,..., x N ) is the probability density function. The stationary solution with vanishing current is found from Eq. (5) by setting: M 1 ] li ( xi Φ)Ω + i M 1 { ] li xj M 1 ] ji Ω } = 0 (6) i,j By multiplying Eq. (6) by M l and summing over l we get: ( x Φ)Ω Ω j xj M 1 ] j = j M 1 ] j xj Ω (7) where we have used the identity l M lm 1 ] li = δ i. We reiterate this step by multiplying by M n and summing over both sides of Eq. (7): Ω M n ( x Φ) + Ω M 1 ] j ( xj M 1 ] n ) = xn Ω (8) j Note that the second term on the l.h.s of Eq. (8) changes sign because xj ( M 1 ] j M n ) = xj (δ jn ) = 0 that gives M 1 ] j ( xj M n ) = ( xj M 1 ] j )M n The second term on l.h.s of Eq. (8) can be further transformed if we use the identity xj M n = xn M j, which stems from the fact that M ij contains the 2 nd -order derivatives of the potential,: M 1 ] j ( xj M 1 ] n ) = j j M 1 ] j ( xn M 1 ] j ) = 1 M x n M where in the last equality M is the determinant of the matrix M and we have used the Jacobi s formula (1/ M ) x M = Tr ( M 1 x M ) 19]. Returning to Eq. (8) we thus have Ω ( ) M n ( x Φ) + 1 M x n M = xn Ω (9) which can be integrated and gives, up to a constant multiplicative (normalization) factor : Ω = exp Φ 2 ] ( xi Φ) 2 M (10) where M is the absolute value of the determinant of M. This last equation which coincides with our main result. i
3 3 II. HAR REPULSIVE POTENTIAL IN 1d: THE GCN CASE For d = 1 the formula (10) reduces to the nown form 15] Ω(x) = exp Φ(x) ] 2 Φ (x) Φ (x) (11) To derive the form taen by Eq. (11) in presence of a steep repulsive potential we consider the fact that, in presence of such a potential, the Ω(x) strongly peas around x = x where Φ (x ) /. This is intuitive since, at the point x, the GCN root mean-squared velocity η 2 = / is balanced by the external (deterministic) velocity component and the particle stops. We further assume that the potential is well approximated by an exponential in the vicinity of x = x Φ λ exp(x x )/λ] (12) which consistently gives Φ (x ) = /. If Ω(x) is peaed around x this must be a maximum which gives the condition 15] Φ (x ) = Φ (x ) (Φ (x ) + 1) 2 If we use Eq.s (12) and (13) into Eq. (11) and we expand it to second order around x we get (13) Ω(x) expa + B(x x ) 2 ] (14) where A = λ + log λ B = 7λ2 + 8λ λ3 2λ 3 + 2λ 2 2 whose integral is Ω(x)dx = 2π exp λ 1 ] ( ) 3/2 + λ 2 7λ 2 (15) + 8λ + 2() 3/2 + 2λ 3 In the limit of λ 0 the potential (12) approaches the potential generated by a hard wall located at x : Φ(x) = if x > x and Φ(x) = 0 if x < x. By taing this limit the integral (15) becomes (note that π/e 1.07): Ω(x)dx = π e (16) Let us now consider one single hard wall of thicness 2x centered at x = 0. We also consider one GCN-driven particle that is free to move around the wall in a 1d space with periodic boundaries located at ±L/2. Using the result given above in this case Eq. (11) becomes Ω = (x x ) + (x + x ) + 1 where (x x ) = 2 δ(x x ), since, in the λ 0 limit, Φ 0 and Ω 1 whenever x > x and Ω 0 if x < x. By integrating Ω we get the normalization constant L/2 L/2 dxω(x) = 2 + L and the normalized probability
4 4 P (x) = (x x ) + (x + x ) + 1 L + 2 (17) where L = L 2x. Eq. (17) can be used to compute the quantity Φ if we recall that Φ (±x ) = ± / and Φ (x) = 0 at any point x ±x. By integrating Eq. (17) we get: Φ = L/2 L/2 dxp (x) Φ 2 (x) = L + 2 (18) III. HAR REPULSIVE POTENTIAL IN 1d: THE RNT CASE To show the similarity between the GCN model and run and tumble (RnT) model in presece of a hard potential we consider one single RnT particle moving in 1d. In the RnT dynamics the noise consists in a two-state Poissonian process happening at rate λ 15]. The RnT noise (also nown as telegraphic noise) is exponentially correlated in time as the GNC, but bimodally distributed switching between two velocities only: η = ±v. The speed v and the characteristic time = (2λ) 1 define the diffusivity = v 2 of the RnT particle in absence of any external potential 6, 11]. We consider the continuity equations that can be written for the RnT model 5]: L(x) = x {L(x) v Φ (x)]} λ 2 L(x) + λ R(x) (19) 2 Ṙ(x) = x {R(x)v Φ (x)]} λ 2 R(x) + λ L(x) (20) 2 where L(x) and R(x) is the fraction of particles with a random velocity component pointing to the left (η = v) and to the right (η = v) respectively. Let us now consider a hard wall as above extending between x and x with periodic boundaries at ±L/2. We consider a stationary solution of Eq.s (19) and (20) of the form L(x) = 2Aδ(x x ) + B/2 (21) R(x) = 2Aδ(x + x ) + B/2 (22) where the constant A represents the fraction of particles going to the left accumulating on the wall at x that is equal to the fraction of particles going to the right accumulating on the wall at x. The constant B in Eq.s (21),(22) represents the fraction of particles that are not in contact with the walls independently on the sign of their velocity. We now focus only on those particle going to the left in the neighbourhood of x (the same conclusions are obtained if we consider those going to the right and x ). Considering the stationary case L = 0 we integrate Eq. (19) it between x and x + ɛ: x +ɛ x +ɛ {L(x + ɛ) v + Φ (x + ɛ)] L(x ) v + Φ (x )]} λ dxl(x) + λ 2 x 2 If we use Eq.s (21),(22) in Eq. (23) we get x dxr(x) = 0 (23) Aλ Φ (0) Φ (ɛ)]b/2 = Aλ + vb/2 = 0 (24) where we have used the fact that Φ (x ) = v. Eq. (24) must be supplemented with the normalization condition L/2 L/2 dx L(x) + R(x)] = 2A + BL = 1 (25) where P (x) = L(x) + R(x) is the overall probability density and L = L x. Using Eq.s (24) and (25) we find B = λ/(v + λl ) and A = v/2(v + λl )] that when expressed as a function of and becomes
5 5 A = B = L L + 2 (26) (27) Using Eq.s (26) and (27) we have P (x) = L(x) + R(x) = (x x ) + (x + x ) + 1 L + 2 (28) (with (x) = 2 δ(x)), which is identical to Eq. (17). IV. RAIALLY SYMMETRIC HAR POTENTIAL Let us consider the case of a spherically symmetric potential for which Eq. (10) becomes Ω(r) = Θ exp Φ(r) ] 2 Φ (r) Φ (r) r + Φ (r) d 1 (29) where Θ is the solid angle in d-dimensions. Let us consider the case d = 2 ( r = x 2 + y 2, Θ = 2π) and a hard potential that is Φ(r) = 0 for r < R and Φ(r) = 0 for r > R. Let us consider one GCN-driven particle that is confined by this circular infinite-well where we expect that Ω(r) peas strongly in correspondence of r = R. We note that Eq. (29) is identical to Eq. (11) multiplied by the term 2π r + Φ (r) when d = 2. We further assume that multiplying Eq. (11) by 2π r + Φ (r) does not change too much the shape of the pea given by the terms where (r) = 2 δ(r). In this way we have that exp Φ ] 2 Φ Φ (r R ) + 1 Ω(r) 2π(r R ) + 1] r + Φ (r) (30) where Φ (R ) = / and Φ (r) = 0 for r < R. The integral of Ω(r) (normalization constant) is computed by noting that, for r < R, Ω(r) = 2πr so that drω(r) = πr 2, while in the neighbourhood of R we get r<r drω(r) = 2π (R + ). Using these results we have from Eq. (30) r R P (r) = 2π(r R ) + 1] r + Φ (r) πr 2 + 2π (R + ) (31) From Eq. (31) we can compute: R Φ = drp (r) Φ 2(R + ) (r) = R R 0 (32) V. COMPARISON OF THE GCN MOEL WITH THE RNT MOEL We briefly illustrate here how the run and tumble (RnT) model compares with the GCN model. The numerical results of the RnT model are compared with the theoretical results obtained by the MUCNASP showing quite unexpected analogies.
6 6 GCN model RnT model (a) (b) (a) (b) (c) (d) (c) (d) FIG. 1: (LEFT) GCN model (a) Probability density function of the position of one GCN-driven particle in presence of the external potential Φ = x 12 (shaded area). Full lines: simulations results, dashed lines: theory, dashed-dotted line: Boltzmann distribution. The more peaed curve corresponds to = 100 µm 2 /s, = 1 s and the less peaed one to = 0.4 µm 2 /s, = 0.1 s (zoomed in the inset) (b) Average value of Φ as a function of for three different values of = 0.1, 0.325, 1 s from top to bottom respectively. Points: simulations results, full lines: theory, dashed lines: theory in the limit of an hard potential, dashed dotted line: equilibrium case with a hard potential. (c) Probability density function of the distance between two GCN-driven particles interacting via the potential Φ = x 12, same legend as Fig. (a). (d) Average value of Φ as a function of for three different values of = 0.1, 0.325, 1 s from top to bottom respectively, same legend as as Fig. (b). (RIGHT) (a) RnT model Probability density function of the position of one RnT particle in presence of the external potential Φ = x 12. Full lines: simulations results, dashed lines: exact theory for the RnT (Eq. (33)). The more peaed curve corresponds to = 100 µm 2 /s, = 1 s and the less peaed one to = 0.4 µm 2 /s, = 0.1 s (b) Average value of Φ as a function of for three different values of = 0.1, 0.325, 1 s from top to bottom respectively. Points: simulations results, full lines: theory (MUCNA), dashed lines: theory in the limit of an hard potential. (c) Probability density function of the distance between two RnT particles interacting via the potential Φ = x 12, same legend as Fig. (a). (d) Average value of Φ as a function of for three different values of = 0.1, 0.325, 1 s from top to bottom respectively, same legend as as Fig. (b). A. One degree of freedom Let us start by the case in which one single particle moves in a 1d space of extension L (periodic boundaries are located at ±L/2) in presence of the external potential Φ = x 12. In this case the probability distribution of the GCN driven particle peas in correspondence of x 1 as shown in Fig. 1(a)(left). This is true also for the RnT model as see Fig. 1(a)(right) where a log-scale has been used for clarity. Note that in this simple case, with one single degree of freedom, the exact P (x) is now for the RnT model 6, 11, 15]: P (x) = A x Φ (x) 2 exp Φ ] (y) dy Φ (y) 2 (33) where A is a normalization constant that absorbs the lower integration limit. Eq. (33) is plotted in Fig. 1(a)(right) as dashed line and it overlaps perfectly the numerical P (x). The average value of the absolute external velocity component Φ can be computed for the GCN case by using the UCNA (Eq. (11)) and is shown as full lines in Fig. 1(b)(left). However we find that this theoretical prediction for the GCN approximates well also the results of the RnT model as shown in Fig. 1(b)(right). This is somehow expected since the expression of the P (x) is almost identical in the two models in the limit of a infinitely steep potential (see Eq.s (17) and (28)) represented by the dashed lines in Fig. 1(b)(left and right).
7 7 GCN model RnT model (a) (b) (a) (b) FIG. 2: (LEFT) GCN model (a) Probability density function of the radial coordinate of one GCN-driven particle in presence of the spherically symmetric external potential Φ = (r R) 12 in 2d (shaded area). Full lines: simulations results, dashed lines: theory, dashed-dotted line: Boltzmann distribution. The more peaed curve corresponds to = 100 µm 2 /s, = 1 s and the less peaed one to = 0.4 µm 2 /s, = 0.1 s (zoomed in the inset) (b) Average value of Φ as a function of for three different values of = 0.1, 0.325, 1 s from top to bottom respectively. Points: simulations results, full lines: theory, dashed lines: theory in the limit of an hard potential, dashed dotted line: equilibrium case with a hard potential. (RIGHT) RnT model (a) Probability density function of the radial coordinate of one RnT particle in presence of the spherically symmetric external potential Φ = (r R) 12 in 2d. Full lines: simulations results, dashed lines: theory. The more peaed curve corresponds to = 100 µm 2 /s, = 1 s and the less peaed one to = 0.4 µm 2 /s, = 0.1 s (b) Average value of Φ as a function of for three different values of = 0.1, 0.325, 1 s from top to bottom respectively. Points: simulations results, full lines: theory (MUCNASP), dashed lines: theory in the limit of an hard potential. B. Two particles in 1d We now consider two particles with coordinates x 1 and x 2 moving in 1d interacting via the pair potential Φ( x) = (x 1 x 2 ) 12 with periodic boundaries at ±L/2. In this case the GCN-driven particles are frequently found in contact leading to the a very peaed P ( x) which is well described by the MUCNASP: Ω( x) = exp Φ( x) 2 ] 2 Φ ( x) Φ ( x) (34) as shown in Fig. 1(c)(left). Eq. (34) allows also to compute the average Φ that is compared with simulation results in Fig. 1(d)(left). ata compare well also with the case of a hard potential obtained by setting 2 in Eq. (18) (dashed lines in Fig. 1(d)(left)) Two RnT particles interacting via the same potential show a more complicated P ( x) that can be obtained numerically but not analytically. This is reported in Fig. 1(c)(right) where P ( x) show a clear secondary pea. This is qualitatively understood as follows: initially the two particle collide with opposite velocities giving rise to the pea at contact, subsequently one of the particles tumble and accelerate for a short distance (being pushed away by the potential Φ = x 12 ), when this becomes negligible the particle has moved at the distance where the 2nd pea is found: in this configuration two particles proceed together in the same direction. We note however that, despite the complicate form of the P ( x), the value of Φ obtained numerically for the RnT model is well approximated by the MUCNASP as shown in Fig. 1(d)(right). This is surprising since the P ( x) of the two models have a quite different form. C. Radially symmetric potential We consider one particle confined in a 2d circular potential well of the form Φ(r) = (r R) 12 where r = x 2 + y 2 and R = 5 µm. The probability density of the GCN-driven particle peas strongly at R R 1 and this is well captured by the MUCNASP (Eq. (29)) as shown in Fig. 2(a)(left). Also the value of the radial external velocity component Φ obtained numerically is well described by the theory (full lines in Fig. 2(b)(left)) and it is close to the hard potential case given by Eq. (32) (dashed lines in Fig. 2(b)(left)). Generalizing the RnT model in 2d we assume that the particle moves at constant speed v and the tumble randomizes completely the swimming direction (see Ref. 27]). The tumble is a Poisson process happening at rate 1 and the axial diffusion coefficient is = v 2 /2 in absence of an external potential. The P (r) of RnT particles interacting with
8 the potential Φ(r) = (r R) 12 is shown in Fig. 2(a)(right) where we can see how the particles accumulate onto the repulsive portion of the potential. This P (r) cannot be predicted theoretically and it is quite different from the P (r) of the GCN-driven particles (see Fig. 2(a)(left)). From the numerical simulations we obtain the quantity Φ that is reported in Fig. 2(b)(right). Again, quite unexpectedly, the values of Φ are close to those found theoretically by the MUCNASP (full lines in Fig. 2(b)(right)) and by the same theory in the hard potential limit (dashed lines in Fig. 2(b)(right)). 8
9 1] R. P. Feynman Statistical Mechanics: A Set Of Lectures, Westview Press, 2nd edition (1998) 2] K. Huang, Statistical Mechanics (John Wiley, New Yor, 1987) 3] M. E. Cates, Rep. Prog. Phys (2012) 4] H. C. Berg, E. coli in Motion (Springer-Verlag, New Yor,2004) 5] M. J. Schnitzer Physical Review E (1993). 6] J. Tailleur and M. E. Cates, EPL (2009) 7] J. Tailleur and M. E. Cates, Phys. Rev. Lett (2008) 8] X. Zheng, B. ten Hagen, A. Kaiser et al. Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys (2013) 9] J. Palacci, S. Sacanna, A. P. Steinberg et al. Science , (2013) 10] J. R. Howse, R. A. Jones Physical review letters 99, (2007) 11] N. Koumais, C. Maggi, R. i Leonardo, Soft Matter, 10, , (2014) 12] C. Maggi, M. Paoluzzi, N. Pellicciotta, A. Lepore, L. Angelani, R. i Leonardo, Phys. Rev. Lett., 113, , (2014) 13] F. O oherty, and J. P. Gleeson. Circuits and Systems II: Express Briefs 54 : (2007) 14] P. Jung and P. Hänggi JOSA B (1988) 15] Hänggi and P. Jung Advances in Chemical Physics 89, (1995) 16] P. Jung and Hänggi Phys. Rev. A 35, 4464 (1987) 17] L. Cao,. Wu and X. Luo Phys. Rev. A 47, 57 (1993) 18] H. Risen, The Foer-Planc Equation: Methods of Solution and Applications (Springer, Berlin, 1984) 19] Bellman, Richard, et al. Introduction to matrix analysis. Vol New Yor: McGraw-Hill, ] Supplemental Material 21] L. Angelani and R. i Leonardo Computer Physics Communications (2011) 22] M. Januszewsi and M. Kostur Computer Physics Communications (2010) 23] V. A. Martinez, R. Besseling, O. A. Croze et al. Biophysical journal (2012) 24] C. Maggi, A. Lepore, J. Solari, A. Rizzo and R. i Leonardo Soft Matter (2013) 25] C. Cammarota and G. Biroli. PNAS (2012) S. Karmaar and G. Parisi PNAS (2013). W. Kob, S. Roldán-Vargas and L. Berthier Nature Physics 8, (2011) 26] L. Angelani, C. Maggi, M. L. Bernardini, A. Rizzo, R. i Leonardo, Phys. Rev. Lett (2011) 27] K. Martens, L. Angelani, R. Leonardo, L. Bocquet, The European Physical Journal E, 35, 1-6, (2012) 28] I.. Vladescu E. J. Marsden, J. Schwarz-Line et al. Phys. Rev. Lett (2014) 9
Velocity distribution in active particles systems, Supplemental Material
Velocity distribution in active particles systems, Supplemental Material Umberto Marini Bettolo Marconi, 1 Claudio Maggi, 2 Nicoletta Gnan, 2 Matteo Paoluzzi, 3 and Roberto Di Leonardo 2 1 Scuola di Scienze
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 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 informationarxiv:chao-dyn/ v1 20 May 1994
TNT 94-4 Stochastic Resonance in Deterministic Chaotic Systems A. Crisanti, M. Falcioni arxiv:chao-dyn/9405012v1 20 May 1994 Dipartimento di Fisica, Università La Sapienza, I-00185 Roma, Italy G. Paladin
More informationPHYSICAL REVIEW LETTERS
PHYSICAL REVIEW LETTERS VOLUME 80 1 JUNE 1998 NUMBER 22 Field-Induced Stabilization of Activation Processes N. G. Stocks* and R. Mannella Dipartimento di Fisica, Università di Pisa, and Istituto Nazionale
More informationNOTE. Application of Contour Dynamics to Systems with Cylindrical Boundaries
JOURNAL OF COMPUTATIONAL PHYSICS 145, 462 468 (1998) ARTICLE NO. CP986024 NOTE Application of Contour Dynamics to Systems with Cylindrical Boundaries 1. INTRODUCTION Contour dynamics (CD) is a widely used
More informationSoft Matter Accepted Manuscript
Accepted Manuscript This is an Accepted Manuscript, which has been through the Royal Society of Chemistry peer review process and has been accepted for publication. Accepted Manuscripts are published online
More informationEffective Electrodiffusion equation for non uniform nanochannels. Abstract
Effective Electrodiffusion equation for non uniform nanochannels Umberto Marini Bettolo Marconi, 1 Simone Melchionna, 2 and Ignacio Pagonabarraga 3 1 Scuola di Scienze e Tecnologie, Università di Camerino,
More informationA Brownian ratchet driven by Coulomb friction
XCIX Congresso Nazionale SIF Trieste, 23 September 2013 A Brownian ratchet driven by Coulomb friction Alberto Petri CNR - ISC Roma Tor Vergata and Roma Sapienza Andrea Gnoli Fergal Dalton Giorgio Pontuale
More informationFrom time series to superstatistics
From time series to superstatistics Christian Beck School of Mathematical Sciences, Queen Mary, University of London, Mile End Road, London E 4NS, United Kingdom Ezechiel G. D. Cohen The Rockefeller University,
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationarxiv: v1 [cond-mat.stat-mech] 26 Jun 2015
EPJ manuscript No. (will be inserted by the editor First-passage time of run-and-tumble particles. Angelani 1a, R. Di eonardo 1 and M. Paoluzzi 1 CNR-IPCF, UOS Roma c/o Dip. di Fisica Università Sapienza,
More informationA Feynman-Kac Path-Integral Implementation for Poisson s Equation
for Poisson s Equation Chi-Ok Hwang and Michael Mascagni Department of Computer Science, Florida State University, 203 Love Building Tallahassee, FL 32306-4530 Abstract. This study presents a Feynman-Kac
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationPHYSICS 219 Homework 2 Due in class, Wednesday May 3. Makeup lectures on Friday May 12 and 19, usual time. Location will be ISB 231 or 235.
PHYSICS 219 Homework 2 Due in class, Wednesday May 3 Note: Makeup lectures on Friday May 12 and 19, usual time. Location will be ISB 231 or 235. No lecture: May 8 (I m away at a meeting) and May 29 (holiday).
More informationarxiv: v2 [cond-mat.stat-mech] 30 Jul 2013
Effective run-and-tumble dynamics of bacteria baths arxiv:135.6475v2 [cond-mat.stat-mech] 3 Jul 213 M. Paoluzzi, R. Di Leonardo, L. Angelani CNR-IPCF, UOS Roma, Dip. Fisica, Università Sapienza, P. le
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Oct 2005
Moment Equations for a Spatially Extended System of Two Competing Species arxiv:cond-mat/0510104v1 [cond-mat.stat-mech] 5 Oct 2005 D. Valenti a, L. Schimansky-Geier b, X. Sailer b, and B. Spagnolo a a
More informationColored extrinsic fluctuations and stochastic gene expression: Supplementary information
Colored extrinsic fluctuations and stochastic gene expression: Supplementary information Vahid Shahrezaei, Julien F. Ollivier, and Peter S. Swain Stochastic simulation algorithm with continuous and discontinuous
More informationPhase Transitions in Nonequilibrium Steady States and Power Laws and Scaling Functions of Surface Growth Processes
Phase Transitions in Nonequilibrium Steady States and Power Laws and Scaling Functions of Surface Growth Processes Term-paper for PHY563 Xianfeng Rui, UIUC Physics Abstract: Three models of surface growth
More informationarxiv: v1 [cond-mat.stat-mech] 15 Sep 2007
Current in a three-dimensional periodic tube with unbiased forces Bao-quan Ai a and Liang-gang Liu b a School of Physics and Telecommunication Engineering, South China Normal University, 56 GuangZhou,
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 informationBoundary Dissipation in a Driven Hard Disk System
Boundary Dissipation in a Driven Hard Disk System P.L. Garrido () and G. Gallavotti (2) () Institute Carlos I for Theoretical and Computational Physics, and Departamento de lectromagnetismo y Física de
More informationQuantifying Intermittent Transport in Cell Cytoplasm
Quantifying Intermittent Transport in Cell Cytoplasm Ecole Normale Supérieure, Mathematics and Biology Department. Paris, France. May 19 th 2009 Cellular Transport Introduction Cellular Transport Intermittent
More informationEffects of colored noise on stochastic resonance in a bistable system subject to multiplicative and additive noise
PHYSICAL REVIEW E, VOLUME 63, 031107 Effects of colored noise on stochastic resonance in a bistable system subject to multiplicative and additive noise Ya Jia, 1, * Xiao-ping Zheng, 1 Xiang-ming Hu, 1
More informationRandom Averaging. Eli Ben-Naim Los Alamos National Laboratory. Paul Krapivsky (Boston University) John Machta (University of Massachusetts)
Random Averaging Eli Ben-Naim Los Alamos National Laboratory Paul Krapivsky (Boston University) John Machta (University of Massachusetts) Talk, papers available from: http://cnls.lanl.gov/~ebn Plan I.
More informationSynchronization of Limit Cycle Oscillators by Telegraph Noise. arxiv: v1 [cond-mat.stat-mech] 5 Aug 2014
Synchronization of Limit Cycle Oscillators by Telegraph Noise Denis S. Goldobin arxiv:148.135v1 [cond-mat.stat-mech] 5 Aug 214 Department of Physics, University of Potsdam, Postfach 61553, D-14415 Potsdam,
More information16. 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 informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 6 Jun 2001
arxiv:cond-mat/0106113v1 [cond-mat.mtrl-sci] 6 Jun 2001 On the power spectrum of magnetization noise G. Durin a,, S. Zapperi b a Istituto Elettrotecnico Nazionale Galileo Ferraris and INFM, Corso M. d
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationTransitions in a Logistic Growth Model Induced by Noise Coupling and Noise Color
Commun. Theor. Phys. (Beijing, China) 46 (2006) pp. 175 182 c International Academic Publishers Vol. 46, No. 1, July 15, 2006 Transitions in a Logistic Growth Model Induced by Noise Coupling and Noise
More informationA priori tests of a stochastic mode reduction strategy
Physica D 170 (2002) 206 252 A priori tests of a stochastic mode reduction strategy A. Majda a,b, I. Timofeyev a,b, E. Vanden-Eijnden a, a Courant Institute of Mathematical Sciences, New Yor University,
More informationarxiv:physics/ v1 [physics.plasm-ph] 5 Nov 2004
Ion Resonance Instability in the ELTRAP electron plasma G. Bettega, 1 F. Cavaliere, 2 M. Cavenago, 3 A. Illiberi, 1 R. Pozzoli, 1 and M. Romé 1 1 INFM Milano Università, INFN Sezione di Milano, Dipartimento
More informationPERSISTENCE AND SURVIVAL IN EQUILIBRIUM STEP FLUCTUATIONS. Chandan Dasgupta
PERSISTENCE AND SURVIVAL IN EQUILIBRIUM STEP FLUCTUATIONS Chandan Dasgupta Centre for Condensed Matter Theory Department of Physics Indian Institute of Science, Bangalore http://www.physics.iisc.ernet.in/~cdgupta
More informationDepartment of Physics PRELIMINARY EXAMINATION 2014 Part II. Long Questions
Department of Physics PRELIMINARY EXAMINATION 2014 Part II. Long Questions Friday May 16th, 2014, 14-17h Examiners: Prof. A. Clerk, Prof. M. Dobbs, Prof. G. Gervais (Chair), Prof. T. Webb, Prof. P. Wiseman
More informationPotential/density pairs and Gauss s law
Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical
More informationRAD229: Midterm Exam 2015/2016 October 19, Minutes. Please do not proceed to the next page until the exam begins.
RAD229: Midterm Exam 2015/2016 October 19, 2015 ---- 75 Minutes Name: Student ID: General Instructions: 1. Write your name legibly on this page. 2. You may use notes including lectures, homework, solutions
More informationActive Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: Hydrodynamics of SP Hard Rods
Active Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: of SP Hard Rods M. Cristina Marchetti Syracuse University Baskaran & MCM, PRE 77 (2008);
More informationSupplementary Material
1 2 3 Topological defects in confined populations of spindle-shaped cells by G. Duclos et al. Supplementary Material 4 5 6 7 8 9 10 11 12 13 Supplementary Note 1: Characteristic time associated with the
More informationPhase separation in systems with absorbing states
EUROPHYSICS LETTERS 1 September 1998 Europhys. Lett., 43 (5), pp. 552-557 (1998) Phase separation in systems with absorbing states M. A. Muñoz 1, U. Marini Bettolo Marconi 1,2 and R. Cafiero 1,3 1 INFM
More informationProbability distributions for the run-and-tumble bacterial dynamics: An analogy to the Lorentz model
Eur. Phys. J. E 35: 84 DOI.4/epje/i-84-y Regular Article THE EUROPEAN PHYSICAL JOURNAL E Probability distributions for the run-and-tumble bacterial dynamics: An analogy to the Lorentz model K. Martens,
More informationAnalytical Dynamics: Lagrange s Equation and its Application A Brief Introduction
Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction D. S. Stutts, Ph.D. Associate Professor of Mechanical Engineering Missouri University of Science and Technology Rolla,
More informationLow-Resistant Band-Passing Noise and Its Dynamical Effects
Commun. Theor. Phys. (Beijing, China) 48 (7) pp. 11 116 c International Academic Publishers Vol. 48, No. 1, July 15, 7 Low-Resistant Band-Passing Noise and Its Dynamical Effects BAI Zhan-Wu Department
More informationGENERATION OF COLORED NOISE
International Journal of Modern Physics C, Vol. 12, No. 6 (2001) 851 855 c World Scientific Publishing Company GENERATION OF COLORED NOISE LORENZ BARTOSCH Institut für Theoretische Physik, Johann Wolfgang
More informationLecture 6: Ideal gas ensembles
Introduction Lecture 6: Ideal gas ensembles A simple, instructive and practical application of the equilibrium ensemble formalisms of the previous lecture concerns an ideal gas. Such a physical system
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 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 informationWeak Localization for the alloy-type 3D Anderson Model
Outline Statement of the result Proof of the Theorem Weak Localization in the alloy-type 3D Anderson Model Zhenwei Cao Math Department Virginia Tech Arizona School of Analysis and Mathematical Physics
More informationNon equilibrium thermodynamic transformations. Giovanni Jona-Lasinio
Non equilibrium thermodynamic transformations Giovanni Jona-Lasinio Kyoto, July 29, 2013 1. PRELIMINARIES 2. RARE FLUCTUATIONS 3. THERMODYNAMIC TRANSFORMATIONS 1. PRELIMINARIES Over the last ten years,
More informationPhysics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017
Physics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
More informationMonopole polarization of C 60 fullerene shell
Monopole polarization of C 6 fullerene shell M. Ya. Amusia 1, and A. S. Baltenkov 3 1 Racah Institute of Physics, the Hebrew University, Jerusalem, 9194 Israel Ioffe Physical-Technical Institute, St. Petersburg,
More information[7] M. Falcioni, E. Marinari, M.L. Paciello, G. Parisi and B. Taglienti, Phys. Lett. B 108
[5] G. Parisi, Statistical Field Theory, Addisson Wesley 1992. [6] U. Wol, Phys. Lett. 228B 3(1989) [7] M. Falcioni, E. Marinari, M.L. Paciello, G. Parisi and B. Taglienti, Phys. Lett. B 108 (1982) 331.
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 informationFree-particle and harmonic-oscillator propagators in two and three dimensions
Free-particle harmonic-oscillator propagators in two three dimensions Ferno Chaos U. Departamento de Física, Escuela Superior de Física y Matemáticas Instituto Politécnico Nacional, México, D. F., MEXICO
More informationLogarithmically modified scaling of temperature structure functions in thermal convection
EUROPHYSICS LETTERS 1 January 2005 Europhys. Lett., 69 (1), pp. 75 80 (2005) DOI: 10.1209/epl/i2004-10373-4 Logarithmically modified scaling of temperature structure functions in thermal convection A.
More informationSTABILITY ANALYSIS OF DYNAMIC SYSTEMS
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY STABILITY ANALYSIS OF DYNAMIC SYSTEMS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationLegendre s Equation. PHYS Southern Illinois University. October 13, 2016
PHYS 500 - Southern Illinois University October 13, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 13, 2016 1 / 10 The Laplacian in Spherical Coordinates The Laplacian is given
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationA Study of the Thermal Properties of a One. Dimensional Lennard-Jones System
A Study of the Thermal Properties of a One Dimensional Lennard-Jones System Abstract In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research,
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationCSCI-6971 Lecture Notes: Monte Carlo integration
CSCI-6971 Lecture otes: Monte Carlo integration Kristopher R. Beevers Department of Computer Science Rensselaer Polytechnic Institute beevek@cs.rpi.edu February 21, 2006 1 Overview Consider the following
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 informationCitation PHYSICAL REVIEW E (2002), 65(6) RightCopyright 2002 American Physical So
TitleStatistical mechanics of two hard d Author(s) Munakata, T; Hu, G Citation PHYSICAL REVIEW E (2002), 65(6) Issue Date 2002-06 URL http://hdl.handle.net/2433/50310 RightCopyright 2002 American Physical
More informationPolymerization and force generation
Polymerization and force generation by Eric Cytrynbaum April 8, 2008 Equilibrium polymer in a box An equilibrium polymer is a polymer has no source of extraneous energy available to it. This does not mean
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationTRANSIENT PHENOMENA IN QUANTUM MECHANICS: DIFFRACTION IN TIME
Paths of Discovery Pontifical Academy of Sciences, Acta 18, Vatican City 2006 www.pas.va/content/dam/accademia/pdf/acta18/acta18-moshinsky.pdf TRANSIENT PHENOMENA IN QUANTUM MECHANICS: DIFFRACTION IN TIME
More informationGeometric projection of stochastic differential equations
Geometric projection of stochastic differential equations John Armstrong (King s College London) Damiano Brigo (Imperial) August 9, 2018 Idea: Projection Idea: Projection Projection gives a method of systematically
More informationLinear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters
Linear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 204 Emo Todorov (UW) AMATH/CSE 579, Winter
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 informationLecture 1. Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev
Lecture Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev Introduction In many cases in nature, like in the Earth s atmosphere, in the interior of stars and planets, one sees the appearance
More informationA scaling limit from Euler to Navier-Stokes equations with random perturbation
A scaling limit from Euler to Navier-Stokes equations with random perturbation Franco Flandoli, Scuola Normale Superiore of Pisa Newton Institute, October 208 Newton Institute, October 208 / Subject of
More informationOn the Asymptotic Convergence. of the Transient and Steady State Fluctuation Theorems. Gary Ayton and Denis J. Evans. Research School Of Chemistry
1 On the Asymptotic Convergence of the Transient and Steady State Fluctuation Theorems. Gary Ayton and Denis J. Evans Research School Of Chemistry Australian National University Canberra, ACT 0200 Australia
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 informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/319/5869/1543/dc1 Supporting Online Material for Synaptic Theory of Working Memory Gianluigi Mongillo, Omri Barak, Misha Tsodyks* *To whom correspondence should be addressed.
More informationDynamics of two coupled particles: comparison of Lennard-Jones and spring forces
Dynamics of two coupled particles: comparison of Lennard-Jones and spring forces Renata Retkute a and James P. Gleeson a a Department of Applied Mathematics, UCC, Cork, Ireland; ABSTRACT The motion of
More informationarxiv: v1 [quant-ph] 15 Dec 2011
Sharp and Infinite Boundaries in the Path Integral Formalism Phillip Dluhy and Asim Gangopadhyaya Loyola University Chicago, Department of Physics, Chicago, IL 666 Abstract arxiv:.3674v [quant-ph 5 Dec
More informationON THE FLUCTUATION-RESPONSE RELATION IN GEOPHYSICAL SYSTEMS
International Journal of Modern Physics B Vol. 23, Nos. 28 & 29 (2009) 5515 5529 c World Scientific Publishing Company ON THE FLUCTUATION-RESPONSE RELATION IN GEOPHYSICAL SYSTEMS GUGLIELMO LACORATA Institute
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationTHREE-BODY INTERACTIONS DRIVE THE TRANSITION TO POLAR ORDER IN A SIMPLE FLOCKING MODEL
THREE-BODY INTERACTIONS DRIVE THE TRANSITION TO POLAR ORDER IN A SIMPLE FLOCKING MODEL Purba Chatterjee and Nigel Goldenfeld Department of Physics University of Illinois at Urbana-Champaign Flocking in
More informationarxiv: v1 [cond-mat.stat-mech] 12 Dec 2018
Active chiral particles under confinement: surface currents and bulk accumulation phenomena arxiv:1812.05134v1 [cond-mat.stat-mech] 12 Dec 2018 1 Introduction Lorenzo Caprini a and Umberto Marini Bettolo
More informationENTIRE FUNCTIONS AND COMPLETENESS PROBLEMS. Lecture 3
ENTIRE FUNCTIONS AND COMPLETENESS PROBLEMS A. POLTORATSKI Lecture 3 A version of the Heisenberg Uncertainty Principle formulated in terms of Harmonic Analysis claims that a non-zero measure (distribution)
More informationarxiv: v1 [cond-mat.stat-mech] 6 Mar 2008
CD2dBS-v2 Convergence dynamics of 2-dimensional isotropic and anisotropic Bak-Sneppen models Burhan Bakar and Ugur Tirnakli Department of Physics, Faculty of Science, Ege University, 35100 Izmir, Turkey
More informationThe Equivalence of Ergodicity and Weak Mixing for Infinitely Divisible Processes1
Journal of Theoretical Probability. Vol. 10, No. 1, 1997 The Equivalence of Ergodicity and Weak Mixing for Infinitely Divisible Processes1 Jan Rosinski2 and Tomasz Zak Received June 20, 1995: revised September
More informationDensity and temperature of fermions and bosons from quantum fluctuations
Density and temperature of fermions and bosons from quantum fluctuations Hua Zheng and Aldo Bonasera 1 1 Laboratori Nazionali del Sud, INFN, via Santa Sofia, 6, 951 Catania, Italy In recent years, the
More informationTable of contents. d 2 y dx 2, As the equation is linear, these quantities can only be involved in the following manner:
M ath 0 1 E S 1 W inter 0 1 0 Last Updated: January, 01 0 Solving Second Order Linear ODEs Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections 4. 4. 7 and
More informationarxiv:nucl-th/ v2 26 Feb 2002
NEGAIVE SPECIFIC HEA IN OU-OF-EQUILIBRIUM NONEXENSIVE SYSEMS A. Rapisarda and V. Latora Dipartimento di Fisica e Astronomia, Università di Catania and INFN Sezione di Catania, Corso Italia 57, I-95129
More informationMOMENTS AND CUMULANTS OF THE SUBORDINATED LEVY PROCESSES
MOMENTS AND CUMULANTS OF THE SUBORDINATED LEVY PROCESSES PENKA MAYSTER ISET de Rades UNIVERSITY OF TUNIS 1 The formula of Faa di Bruno represents the n-th derivative of the composition of two functions
More informationTopic 6: Projected Dynamical Systems
Topic 6: Projected Dynamical Systems John F. Smith Memorial Professor and Director Virtual Center for Supernetworks Isenberg School of Management University of Massachusetts Amherst, Massachusetts 01003
More informationLeast action principle and stochastic motion : a generic derivation of path probability
Journal of Physics: Conference Series PAPER OPEN ACCESS Least action principle and stochastic motion : a generic derivation of path probability To cite this article: Aziz El Kaabouchi and Qiuping A Wang
More informationArticle (peer-reviewed)
Title Author(s) Influence of noise intensity on the spectrum of an oscillator Swain, Rabi Sankar; Gleeson, James P.; Kennedy, Michael Peter Publication date 2005-11 Original citation Type of publication
More informationEXISTENCE OF POSITIVE PERIODIC SOLUTIONS OF DISCRETE MODEL FOR THE INTERACTION OF DEMAND AND SUPPLY. S. H. Saker
Nonlinear Funct. Anal. & Appl. Vol. 10 No. 005 pp. 311 34 EXISTENCE OF POSITIVE PERIODIC SOLUTIONS OF DISCRETE MODEL FOR THE INTERACTION OF DEMAND AND SUPPLY S. H. Saker Abstract. In this paper we derive
More informationScalar electromagnetic integral equations
Scalar electromagnetic integral equations Uday K Khankhoje Abstract This brief note derives the two dimensional scalar electromagnetic integral equation starting from Maxwell s equations, and shows how
More informationProgress in the Search for Continuous Gravitational Waves from Binary Pulsars and a Proposed Approach to Machine-Learning Detection
Progress in the Search for Continuous Gravitational Waves from Binary Pulsars and a Proposed Approach to Machine-Learning Detection Shalma Wegsman Gueron The University of Chicago, Department of Physics
More information8: Source-Sink Problems in 1 Energy Group
8: Source-Sink Problems in 1 Energy Group B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 015 Sept.-Dec. 015 September 1 Contents Solving the 1-group diffusion
More informationAPPLICATION OF A STOCHASTIC REPRESENTATION IN NUMERICAL STUDIES OF THE RELAXATION FROM A METASTABLE STATE
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 7 (1), 83-90 (2001) APPLICATION OF A STOCHASTIC REPRESENTATION IN NUMERICAL STUDIES OF THE RELAXATION FROM A METASTABLE STATE FERDINANDO DE PASQUALE 1, ANTONIO
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationEffects of ion temperature on electrostatic solitary structures in nonthermal plasmas
PHYSICAL REVIEW E VOLUME 55, NUMBER FEBRUARY 1997 Effects of ion temperature on electrostatic solitary structures in nonthermal plasmas A. A. Mamun Department of Physics, Jahangirnagar University, Savar,
More informationSynchronized Queues with Deterministic Arrivals
Synchronized Queues with Deterministic Arrivals Dimitra Pinotsi and Michael A. Zazanis Department of Statistics Athens University of Economics and Business 76 Patission str., Athens 14 34, Greece Abstract
More informationTHE problem of phase noise and its influence on oscillators
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: EXPRESS BRIEFS, VOL. 54, NO. 5, MAY 2007 435 Phase Diffusion Coefficient for Oscillators Perturbed by Colored Noise Fergal O Doherty and James P. Gleeson Abstract
More informationThe correlation between stochastic resonance and the average phase-synchronization time of a bistable system driven by colour-correlated noises
Chin. Phys. B Vol. 19, No. 1 (010) 01050 The correlation between stochastic resonance and the average phase-synchronization time of a bistable system driven by colour-correlated noises Dong Xiao-Juan(
More information