arxiv: v2 [nlin.cg] 9 May 2016
|
|
- Violet Francis
- 5 years ago
- Views:
Transcription
1 Phase coexistence phenomena in an extreme case of the misanthrope process with open boundaries Chikashi Arita, Chihiro Matsui Theoretische Physik, Universität des Saarlandes, 664 Saarbrücken, Germany Mathematical Informatics, The University of Tokyo, Tokyo, Japan Abstract arxiv:65.97v [nlin.cg] 9 May 6 The misanthrope process is a class of stochastic interacting particle systems, generalizing the simple exclusion process. It allows each site of the lattice to accommodate more than one particle. We consider a special case of the one dimensional misanthrope process whose probability distribution is completely equivalent to the ordinary simple exclusion process under the periodic boundary condition. By imposing open boundaries, high- and low-density domains can coexist in the system, which we investigate by Monte Carlo simulations. We examine finite-size corrections of density profiles and correlation functions, when the jump rule for particles is symmetric. Moreover, we study properties of delocalized and localized shocks in the case of the totally asymmetric jump rule. The localized shock slowly moves to its stable position in the bulk. Introduction: Exclusion process The exclusion process is one of paradigmatic models of non-equilibrium statistical physics, which is a Markov process of interacting particles on lattices []. et us consider a chain consisting of sites. Each site i of the chain is either empty (τ i = ) or accommodates a particle (τ i = ). A particle stochastically jumps to one of its nearest neighbor sites with the leftward and rightward rates, say, p and qp, if a target site is empty. The stationary distribution of the exclusion process is very simple in the periodic boundary condition where the th site is regarded as i = : all the possible configurations of particles are equally probable. Therefore, the stationary current is given by ( q)pρ( ρ) in the limit. Here ρ denotes the global density, i.e. the number of particles divided by. Open boundary conditions also provide us important perceptions. At the leftmost site (i = ), particles are injected into and extracted from the system with rates α and γ, respectively. Meanwhile we set the injection and extraction rates at the rightmost site (i = ) to be δ and β, respectively. et us consider that the left and right ends are connected to reservoirs of densities ρ and ρ, respectively, by relations α = pρ, β = p( ρ ), γ = pq( ρ ) and δ = pqρ. In the symmetric case q =, the stationary current and the density profile are given by [] J(ρ, ρ ) = p(ρ ρ )/, ρ i = ρ + (ρ ρ )i/. () The diffusion coefficient is simply given by the jump rate p, which is compatible with (). The exact stationary state is given in a matrix product form [], from which the correlation C = τ i τ i+ τ i τ i+ is calculated as C = i( i )(ρ ρ ) /[ ( )] () x( x)(ρ ρ ) / (x = i/, ). (3) Exact solutions and phenomenological arguments have been proposed also in the asymmetric case q []. Misanthrope process with the periodic boundary condition et us consider generalized exclusion processes, where each site can accommodate more than one particle [3, 4, 5]. The jump rates of particles depend on the occupation numbers of both departure and target sites: m n m n+ (rate p mn ), (4) m n m+ n (rate qp nm ). (5)
2 Figure : Misanthrope process with maximal occupancy, connected to density reservoirs at both ends. The arrows represent possible jumps of particles. In the extreme case (p = ), the dashed lines correspond to prohibited jumps. The misanthrope process [6] is a class of generalized exclusion processes, which is one of relevant models for traffic flow[7, 8]. When the maximal occupation number per site is, the misanthrope process is characterized by p = p + p. (6) See Fig. for the illustration of the model. (The open boundary condition will be explained in the next section.) We define ratios a and b of the parameters as p = a p, p = b p. (7) Under the periodic boundary condition, the stationary probability of finding a configuration τ τ is given by the product of single-site weights X m s as [6] P (τ τ ) = X τ X τ, (8) X = Z, X = λ Z, X = aλ Z, Z = + λ + aλ (9) with the fugacity parameter λ. In the limit with global density ρ fixed, λ is specified by the relation ln Z = ρ, which can be easily solved as λ λ λ = ρ/( ρ + R), R = ( 4a)( ρ)ρ. () We shall occasionally use notations X m (ρ) in order to emphasize that they are functions of the density via (). The rightward current of particles is calculated as [7] J (ρ) = m=, n=, X m (ρ)x n (ρ)p mn () = p ρ( ρ) [ ( ρ) + b(ρ + R )/ ] /( + R). () The leftward current and the current are given as J (ρ) = qj (ρ) and J(ρ) = ( q)j (ρ), respectively. The shape of the fundamental diagram changes, depending on the values of the parameters a and b, see Fig.. In this work we are particularly interested in the case where particles are extremely misanthrope, i.e. a =. (This restriction is similar to the repulsion process [9], a special case of the Katz-ebowitz-Spohn model [].) Under the periodic boundary condition, the probabilities of finding configurations are governed by the same master equations as the usual simple exclusion process: when the global density ρ is less than, all the sites are either occupied by one particle or empty. If the target site is occupied, any particle cannot jump due to p =. For ρ >, no empty site appears, and we regard the configurations τ τ as τ τ with τ i = τ i. Then we recover the same master equations as the simple exclusion process. For ρ = only the configuration is found, and J =. In summary, the single-site weights become { ( ρ, ρ, ) (ρ ), (X, X, X ) = (3) (, ρ, ρ ) (ρ > ), and the fundamental diagram of this extreme case consists of the two parabolas J = ( q)p ρ( ρ) ( < ρ ) and J = ( q)p (ρ )( ρ) ( < ρ < ), see Fig.. For ρ >, the dynamics (e.g. the mean-squared displacement) of an individual particle is different from the simple exclusion process.
3 J.5.5 a=.5, b=.5 a=.65, b= a=., b=. a=.5, b= J.5.5 Figure : Fundamental diagrams of the misanthrope process for various values of the parameters a and b with ( q)p =, and the extreme case a = with ( q)p = and ( q)p =. 3 Open boundary condition We consider the misanthrope process on the chain with open ends (Fig. ). At the left boundary, particles are injected and extracted with rates α τ and γ τ, respectively, depending on the occupation number τ of site. In the same way, at the right boundary, we impose injection and extraction rates, δ τ and β τ, respectively, depending on the occupation number τ. To realize the reservoir densities ρ and ρ at virtual sites i =,, we set α τ = p mτ X m (ρ ), β τ = p τm X m (ρ ), (4) m=, γ τ = q m=, m=, p τm X m (ρ ), δ τ = q m=, For the special case ρ = ρ, the stationary state is written in the product form (8). 4 Symmetric case For the symmetric case q =, the current between sites i and i + is given by p mτ X m (ρ ). (5) J i = G(τ i, τ i+ ) ( < i < ), (6) G(τ, σ) = ) p mn (g m (τ)g n (σ) g n (τ)g m (σ). (7) m=, n=, We use the following polynomial representation of Kronecker s delta g m (τ) = (m = τ), (m τ) for m, τ {,, }: g (τ) = g ( τ) = ( τ)( τ), g (τ) = τ( τ). For the boundaries, we have J = G(ρ, τ ), J = G(ρ, τ ), (8) G(ρ, τ) = ) p mn (X m (ρ)g n (τ) X n (ρ)g m (τ). (9) m=, n=, One finds that the current satisfies the gradient condition J i = K i K i+ [3], where K i s are given by { p g (τ i ) + p g (τ i ) ( < i < ), K i = h(ρ i ) (i =, ), () h(ρ) :=p X (ρ) + p X (ρ) () =p ρ [ ( ρ) + b(ρ + R )/ ]/ ( + R). () In the stationary state we have J = J = = J (=: J). Then we achieve Fick s law J(ρ, ρ ) = [ h(ρ ) h(ρ ) ] /, (3) K i =h(ρ ) + [ h(ρ ) h(ρ ) ] i/, (4) 3
4 i p =3, =3 p =9, =4.5 i/ ( x *, ).. p =6 p =3 p = p =9 d ` ( i, ) d /c.. (c) =78 =5 =36 =56 ( i/ x ) * d ( i, ) d /c.. (d) p =9, =7 p =6, =78 p =3, =78 p =, =7 5 5 * ( i/ x ) d /u Figure 3: Density profiles of the symmetric misanthrope process with p =. The values of p and are given in the panel. The discrepancies at site i = x for various values of p. The simulation results and numerical fittings are indicated by circles and solid lines, respectively. (c) Rescaled discrepancies vs. rescaled site for various values of with (p, p ) = (3, ). (d) Rescaled discrepancies vs. rescaled site for various values of p. For all the four sub-figures we have set the boundary densities as (ρ, ρ ) = (, ). For each set of the parameters (p, ) we performed two independent simulation runs. We also averaged the simulation data over 6 t 9 for p =, 3 or 6 t 8 for p = 6, 9. which indicates that the diffusion constant D is given by D(ρ) = d dρ h(ρ) = p This formula can be also expressed in the form ( ρ) + b(ρ + R ). (5) R D(ρ) = χ J (ρ), χ = Rρ( ρ)/( + R), (6) where the factor χ corresponds to the fluctuation of the density χ = λ dρ dλ = m=,, X m(ρ) ( m ρ ). In the special case p = p (b = ), we have h(ρ) = p ρ, D(ρ) = p and the linear profile because of K i = ρ i : ρ i = ρ + (ρ ρ )i/. (7) We expect that the hydrodynamic description is valid in a large system with x = i/ and t = t/ : t ρ(x, t ) = [ D(ρ(x, t )) ] x x ρ(x, t ). (8) In the stationary state, we have leading to an implicit form of a density profile ρ(x): D(ρ(x)) d dx ρ(x) = J(ρ, ρ ), (9) h(ρ(x)) h(ρ ) = xj(ρ, ρ ). (3) Though the formula for the current (3) is correct for any finite, the prediction (3) does not always provide the true analytic formula. However, we expect that the density profile τ x converges to (3) in the limit. When p =, the diffusion coefficient reduces to D(ρ) = p (ρ < ), p (ρ > ). (3) The density profiles in the cases ρ < ρ < and ρ > ρ > are equivalent to the simple exclusion process. An interesting situation is the case ρ > ρ <. From (3), the hydrodynamic prediction (3) of the density profiles becomes piecewise linear: { ρ + x ρ(x) = x ( ρ ) ( < x < x ), ρ + x x ( ρ ) (x (3) < x < ), x =(ρ )p / [ ( ρ )p + (ρ )p ]. (33) One may derive the formula (6) by assuming the factorization (8) with the weights (9) in the bulk sites []. In our case, we did not use any non-rigorous assumption to achieve (5). 4
5 C i x C 5 x=/8 x=/4 x=/.5 i/( ) C i x C 3.5 i/( ) x=/8 x=/4 x=/ C i (c) x=/8 x=/4 x=/ x C i/( ) C i.3 (d) x C 5.5 i/( ) x=/8 x=/4 x=/ Figure 4: Correlation profiles (37) with (m, n) = (,), (,), (,) (c) and (,) (d). The parameters were chosen as (p, p, ρ, ρ ) = (3,,, ) (x = /4). For each, we performed two independent simulation runs. We also averaged the simulation data over 6 t 9. The position x separates the space into high and low density domains, i.e. ρ(x) > for x < x and ρ(x) < for x > x. Under the further restriction p = p, it is completely linear, ρ(x) = ρ + x(ρ ρ ), which is compatible with the true density profile (7). On the other hand, for p p, we observe that the simulation results are deviated from the prediction, in particular, near x = x, see Fig. 3. et us denote the difference between the true density ρ i measured by simulations and the prediction by ρ(i, ) = ρ i ρ(i/). Figure 3 shows the discrepancy at the site i = x vs.. It seems that it exhibits power-law decay ρ(x, ) c d, and we expect that the exponent d is independent of b. The inset shows our numerical estimation of d, depending on the interval of the system size l/5 < < l that we use, e.g. d 4. for l 7. (The error bars are due to changing the value of p.) For the fitting lines in and rescaling of the x- and y-axes of (c) and (d), we use this value. 3 We observe overlap of markers for different values of in Fig. 3 (c), d ρ(x, )/c vs. (x x ) d. Furthermore we can find rescalers u + > and u <, such that the discrepancies of different values of b also overlap, see Fig. 3 (d). Technically we determined u ± via x ± as ρ(x ±, ) d (c) =.5 (x < x < x + ), (34) u ± = (x ± x ) d. (35) This result indicates the existence of a scaling function f ρ(x, ) { c d f ( (u + ) (x x ) d) (x x ), c d f ( (u ) (x x ) d) (x < x ). (36) Now we turn to nearest-neighbor correlation functions C mn i = g m (τ i ) g n (τ i+ ) g m (τ i )g n (τ i+ ) (37) for m, n {,, }, which satisfy m=,, Cmn i = n=,, Cmn i =. et us check whether the asymptotic behaviors are the same as the SSEP case (3). Fig. 4 shows the correlation profiles, multiplied by the system size. It seems that, far from i/( ) x = x, Ci mn s decay being proportional to or faster than the power law, since the lines of different values of are overlapping or almost. On the other hand, in the vicinity of x = x, they decay more slowly than. The insets also support these observations. (Another scenario could be that the finite-size effects are too strong to observe O( )-decay). 5 Totally asymmetric case We consider the totally asymmetric case q =. Since we know the fundamental diagram (Fig. ), the phase diagram (Fig. 5 ) is predictable by the extremal current principal, see the original works [, 3] for details e.g. the phase boundaries. In the extreme case p =, the regions ρ < ρ < and ρ > ρ > are equivalent to the usual simple exclusion process. In the minimal current phase the stationary state is an absorbing state, where each site is occupied by only one particle, and thus the particle current is. Therefore 3 We think that the power-law decay with the exponent slightly bigger than.4 is the most reasonable guess. However the following possibility has not been excluded: the slopes in the logarithmic frame continue to slowly decrease and diverge in the limit. 5
6 bulk = bulk= min bulk = bulk= bulk = bulk= maxi bulk= maxii Figure 5: Phase diagram of the totally asymmetric misanthrope process in the case where J(ρ) has two local maxima J(ρ maxi ) > J(ρ maxii ) (ρ maxi < ρ maxii ). The stationary current is give by J(ρ bulk ), in particular, J(ρ min ) is the local minimum. Phase diagram for p = p = p = (a =, b = ). we focus on the case ρ > > ρ. We also restrict our consideration to the case p = p =, where the two maximal currents have the same value. The phase diagram of the current is simplified as Fig. 5. We are particularly interested in properties of an anti-shock [4]. For the simple exclusion process, the second-class particle [5] is often used in order to microscopically probe the shock position. In this work, we define the shock position S as follows. et us denote by S the rightmost site occupied by two particles, and by S the leftmost empty site. Any configuration is written as τ τ S τ S+ τ (38) with τ i {, } for i < S and τ i {, } for i > S. In particular, S = when there is no site occupied by two particles, and S = when there is no empty site. Then we simply set S = S+S {,, 3,, }. (See [6, 7] for similar microscopic definitions of shocks.) Note that we always have S < S in the case ρ > > ρ. The tag S increases by one, according to jump of a particle at site S under the condition that τ S+ =. On the other hand, the tag S decreases by one, according to jump of a particle at site S under the condition that τ S =. These events shift the shock position leftward or rightward by /. When S S =, a particle at site S jumps to site S, which causes a non-local shift; in this situation, the tags are renewed, i.e. the second rightmost doubly occupied site becomes the new S and the second leftmost empty site becomes the new S. For example, for a configuration shown in the Fig. ; if a particle on the fourth site jumps to the fifth, the shock position changes as (S, S, S ) = (4, 9/, 5) (, 4, 6). The shock is delocalized on the phase transition line < ρ = ρ <.5, as shown in Fig. 6. The densities of the domains i < S and i > S are given by the reservoir densities ρ and ρ, respectively [3]. Figure 6 shows typical density profiles of spatial average r i := l l+ k= l τ k+i with l =, by using snapshots τ τ without time or ensemble average. Under the assumption that the motion of the shock is governed by a random walk with reflective boundaries, the density profile averaged over long time becomes linear, connecting ρ and ρ, see Fig. 6 (c). As shown in Fig. 6 (d), there is a regime where the mean-squared displacement, MSD = (S(t + t ) S(t )), is proportional to time t, supporting the random-walk description of the shock motion. In the inset, we plot the diffusion coefficient D S, which was estimated from simulation data t (S(t + t ) S(t )) in t [T/, T ]. We chose the value of T so as to avoid the saturation of the linearity caused by the reflective boundaries. The line is a guessed form D S = ( ρ)(3 ρ) ρ, which approaches as ρ.5. Actually, on the point ρ = ρ =.5, the location of the shock is restricted to the vicinity of S = /, see Fig. 6,, and (c). Now we examine properties of the shock in the maximal current phase. We denote the rescaled shock position by s = S/. In the sub-phase ρ <.5 ρ <.5, the left and right density domains are ρ eft = ρ and ρ Right =.5, from the general framework [, 3]. The shock velocity is given by v S = J(ρ eft) J(ρ Right ) ρ eft ρ Right, which is negative in this case. Finally the shock reaches the left boundary, in other words s =, see Fig. 7. In the sub-phase ρ >.5 ρ <.5 and its boundaries (dashed lines in Fig. 5 ), the densities of the left and right domains are ρ eft =.5 and ρ Right =.5, and therefore v S =. Indeed, from Fig. 7, it is difficult to see a clear tendency of the shock motion. By changing the time scale as Fig. 7, however, we observe that the shock moves to s =.5 in this sub-phase, and s =.5 on the phase boundary ρ =.5 ρ <.5. The density profiles in Fig. 7 (c) also imply that s =.5 and.5 are the stable positions. In the inset, the deviations from these values are observed, which are expected to be finite-size effect. On the other sub-phase boundary ρ <.5 ρ =.5, s =.75 because of a symmetry. We also measured the first passage time (F P T ), i.e. the 6
7 t =. =.5.5 s =.3 ri =.3 =. =.5.5 i/ i (c).5 i/ =. =.3 = (d) =. =.5 =.4 slope t DS..5 Figure 6: Simulation results on the phase transition line < ρ = ρ <.5. Kymograph of shocks in single simulation runs with = 4 and the initial position s =.5. Density profiles r i for = 4, made from snapshots. (c) Density profiles ρ i for = 3, averaged over time interval 6 t 8 and over two independent runs. (d) Mean-squared displacement vs. time and (inset) diffusion coefficient D S, by averaging over 5 simulation runs. first time when the shock hits a stable position, see Fig. 7 (d). The exponents z (F P T = O( z )) for ρ =.5 and.6 are found to be.4 and., respectively, while z = for ρ =.4 due to v S. We define a shock width by the standard deviation of the shock position S = S S, which seems to be proportional to some power of, i.e. S e, see Fig. 7 (e). The exponent e from fitting is e.74. We also performed fitting for other values of ρ with ρ =., and found the exponent between.7 < e <.8 (not shown here). This result is different from / and /3 observed in the exclusion process with a single-site defect [8]. The probability distribution of S (inset) consists of two curves, according to whether S is an integer or a half-integer. Both of them are well fitted by Gaussian distributions. One may also regard S S as another definition of the shock width, which is instantaneously measured in a given configuration. Figure 7 (f) implies that its average converges to some value as. (The solid line is a fitting curve in the form S S = v w.) This behavior is different from S. Furthermore we expect that the probability distribution of S S exhibits asymptotically exponential decay, see the inset of Fig. 7 (f). 6 Conclusions We investigated the extreme case (p = ) of the misanthrope process by imposing high (ρ > ) and low (ρ < ) densities at the reservoirs. For the symmetric case (q = ), the formula for the current (3) is correct for any system size. On the other hand, the true density profile is, in general, deviated from the hydrodynamic prediction, where the point x separates between high and low density domains. We expect that this discrepancy is a finite-size effect and exhibits a power-law decay at x = x. We also conjecture that there exists a scaling function in the vicinity of x = x. The nearest-neighbor correlations decay more slowly than O( ) in the vicinity of x = x as far as we observed simulations for <. They are proportional to or faster than power-law decays for x x or x x. For the totally asymmetric case (q = ), an anti-shock can exist, and we introduced a microscopic definition of its position. On the transition line ρ + ρ = ( < ρ <.5), where the motion of the shock position is diffusive, we probed the mean-squared displacement and measured the diffusion constant characterizing the motion of the shock. In the sub-phase ρ >.5 ρ <.5 of the maximal current phase, the formula of the shock velocity becomes. However we found that the shock spontaneously reaches a stable position s =.5. On the sub-phase boundaries the stable position becomes s =.5 or.75. The exponent of the shock width was found between.7 and.8. Acknowledgements CM is supported by JSPS Grant-in-Aid No. 5K939. References [] Derrida B., J. Stat. Mech. (7) P73 7
8 5 t 5 7 =.4 =.5 =.6 4 S t 7 4 =.5 =.6 =.4 4 S i (c) =.6 =.5 =.4 h s i i/ (d) Initial Position s =.75 Stable Position.6 s=.5.5 s=.5.4 s= S 7 3 S S.. Prob(S) (e) h S S i (f) Prob( S S ) S S Figure 7: Simulation results in the maximal current phase. The density of the right reservoir is chosen as ρ =.. (a,b) Kymographs of the shock in single simulation runs with = 4 and the initial position S =.75. The dashed line is the theoretical line S = tv S (c) Density profiles ρ i for = 3, and (inset) mean shock positions near the sub-phase boundary with guide for eyes by lines. (d) First passage time. The dashed line is F P T =.75/v S, and thin lines are fitting. (e) Shock width vs. system size for (ρ, ρ ) = (.7,.), and (inset) probability distribution of the shock position for =. (f) Mean instantaneous shock width vs. system size for (ρ, ρ ) = (.7,.), and (inset) probability distribution. Technical details: We performed two independent simulation runs for (c), 5 / runs for (d) or ten runs for (e) and (f). We also averaged over 6 t 8 for (c), (e) and (f). [] Blythe R. A. and Evans M. R., J. Phys. A: Math. Gen. 4, R333 (7) [3] Kipnis C., andim C., and Olla S., Comm. Pure Appl. Math. 47, 475 (994) [4] Arita C., Krapivsky P.., and Mallick K., Phys. Rev. E 9, 58 (4) [5] Matsui C., J. Stat. Phys. 58, 9 (5) [6] Cocozza-Thivent C., Z. Wahrscheinlichkeitstheorie verw. Gebiete 7, 59 (985) [7] Kanai M., Phys. Rev. E 8, 667 () [8] Eymard R., Roussignol M., and Tordeux A., Stoch. Process. Their Appl., 3648 () [9] Krapivsky P.., J. Stat. Mech. (3) P6. [] Katz S., ebowitz J.., and Spohn H., J. Stat. Phys. 34, 497 (984) [] Becker T., Nelissen K., Cleuren B., Partoens B., and Van den Broeck C., Phys. Rev. E 9, 539 (4) [] Popkov V. and Schütz G. M., Europhys. ett. 48, 57 (999) [3] Hager J. S., Krug J., Popkov V., and Schütz G. M., Phys. Rev. E 63, 56 () [4] Belitsky V. and Schütz G. M., J. Phys. A: Math. Theor. 46, 954 (3) [5] Boldrighini C., Cosimi C., Frigio A., and Grasso-Nuñes G.: J. Stat. Phys. 55, 6 (989) [6] Cividini J., Hilhorst H. J., and Appert-Rolland C., J. Phys. A: Math. Theor. 47, (4) [7] de Gier J. and Finn C., J. Stat. Mech. (4) P74. [8] Janowsky S. A. and ebowitz J.., Phys. Rev. A 45, 68 (5) 8
arxiv: v2 [cond-mat.stat-mech] 24 Oct 2007
arxiv:0707434v2 [cond-matstat-mech] 24 Oct 2007 Matrix Product Steady States as Superposition of Product Shock Measures in D Driven Systems F H Jafarpour and S R Masharian 2 Bu-Ali Sina University, Physics
More informationOpen boundary conditions in stochastic transport processes with pair-factorized steady states
Open boundary conditions in stochastic transport processes with pair-factorized steady states Hannes Nagel a, Darka Labavić b, Hildegard Meyer-Ortmanns b, Wolfhard Janke a a Institut für Theoretische Physik,
More informationAvailable online at ScienceDirect. Physics Procedia 57 (2014 ) 77 81
Available online at www.sciencedirect.com ScienceDirect Physics Procedia 57 (204 ) 77 8 27th Annual CSP Workshops on Recent Developments in Computer Simulation Studies in Condensed Matter Physics, CSP
More informationLimiting shapes of Ising droplets, fingers, and corners
Limiting shapes of Ising droplets, fingers, and corners Pavel Krapivsky Boston University Plan and Motivation Evolving limiting shapes in the context of Ising model endowed with T=0 spin-flip dynamics.
More informationTwo-channel totally asymmetric simple exclusion processes
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 37 (24) 997 998 PII: S35-447(4)83426-7 Two-channel totally asymmetric simple exclusion processes Ekaterina
More informationSecond order additive invariants in elementary cellular automata
November 21, 2018 arxiv:nlin/0502037v1 [nlin.cg] 17 Feb 2005 Second order additive invariants in elementary cellular automata Henryk Fukś Department of Mathematics Brock University St. Catharines, Ontario
More informationNon-equilibrium phase transitions
Non-equilibrium phase transitions An Introduction Lecture III Haye Hinrichsen University of Würzburg, Germany March 2006 Third Lecture: Outline 1 Directed Percolation Scaling Theory Langevin Equation 2
More informationarxiv: v1 [cond-mat.stat-mech] 30 Aug 2015
Explosive condensation in symmetric mass transport models arxiv:1508.07516v1 [cond-mat.stat-mech] 30 Aug 2015 Yu-Xi Chau 1, Colm Connaughton 1,2,3, Stefan Grosskinsky 1,2 1 Centre for Complexity Science,
More informationThe Totally Asymmetric Simple Exclusion Process (TASEP) and related models: theory and simulation results
The Totally Asymmetric Simple Exclusion Process (TASEP) and related models: theory and simulation results Chris Daub May 9, 2003 1 Synopsis Simple exclusion processes The Bethe ansatz TASEP Definition
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 informationLarge deviations of the current in a two-dimensional diffusive system
Large deviations of the current in a two-dimensional diffusive system C. Pérez-Espigares, J.J. del Pozo, P.L. Garrido and P.I. Hurtado Departamento de Electromagnetismo y Física de la Materia, and Instituto
More informationEvaporation/Condensation of Ising Droplets
, Elmar Bittner and Wolfhard Janke Institut für Theoretische Physik, Universität Leipzig, Augustusplatz 10/11, D-04109 Leipzig, Germany E-mail: andreas.nussbaumer@itp.uni-leipzig.de Recently Biskup et
More informationD riven diffusive system is a rewarding research topic in recent decades, which exhibits non-vanishing
OPEN SUBJECT AREAS: STATISTICAL PHYSICS NONLINEAR PHENOMENA PHASE TRANSITIONS AND CRITICAL PHENOMENA Bulk induced phase transition in driven diffusive systems Yu-Qing Wang 1, Rui Jiang 1, Anatoly B. Kolomeisky
More informationLocal inhomogeneity in asymmetric simple exclusion processes with extended objects
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 37 (24) 25 23 PII: S35-447(4)69849- Local inhomogeneity in asymmetric simple exclusion processes with
More informationarxiv: v1 [cond-mat.stat-mech] 12 Jun 2007
TOPICAL REVIEW arxiv:0706.678v [cond-mat.stat-mech] 2 Jun 2007 Nonequilibrium Steady States of Matrix Product Form: A Solver s Guide Contents R. A. Blythe and M. R. Evans SUPA, School of Physics, University
More informationPhysics Letters A 375 (2011) Contents lists available at ScienceDirect. Physics Letters A.
Physics Letters A 375 (2011) 318 323 Contents lists available at ScienceDirect Physics Letters A www.elsevier.com/locate/pla Spontaneous symmetry breaking on a mutiple-channel hollow cylinder Ruili Wang
More informationQuasi-Stationary Simulation: the Subcritical Contact Process
Brazilian Journal of Physics, vol. 36, no. 3A, September, 6 685 Quasi-Stationary Simulation: the Subcritical Contact Process Marcelo Martins de Oliveira and Ronald Dickman Departamento de Física, ICEx,
More informationSTATIONARY NON EQULIBRIUM STATES F STATES FROM A MICROSCOPIC AND A MACROSCOPIC POINT OF VIEW
STATIONARY NON EQULIBRIUM STATES FROM A MICROSCOPIC AND A MACROSCOPIC POINT OF VIEW University of L Aquila 1 July 2014 GGI Firenze References L. Bertini; A. De Sole; D. G. ; G. Jona-Lasinio; C. Landim
More informationPhysics of Cellular materials: Filaments
Physics of Cellular materials: Filaments Tom Chou Dept. of Biomathematics, UCLA, Los Angeles, CA 995-766 (Dated: December 6, ) The basic filamentary structures in a cell are reviewed. Their basic structures
More informationCoarsening process in the 2d voter model
Alessandro Tartaglia (LPTHE) Coarsening in the 2d voter model May 8, 2015 1 / 34 Coarsening process in the 2d voter model Alessandro Tartaglia LPTHE, Université Pierre et Marie Curie alessandro.tartaglia91@gmail.com
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 informationarxiv: v2 [cond-mat.stat-mech] 1 May 2018
arxiv:1712.678v2 [cond-mat.stat-mech] 1 May 218 GRADIENT STRUCTURE AND TRANSPORT COEFFICIENTS FOR STRONG PARTICLES DAVIDE GABRIELLI AND P. L. KRAPIVSKY Abstract. We introduce and study a simple and natural
More informationEquilibrium correlations and heat conduction in the Fermi-Pasta-Ulam chain.
Equilibrium correlations and heat conduction in the Fermi-Pasta-Ulam chain. Abhishek Dhar International centre for theoretical sciences TIFR, Bangalore www.icts.res.in Suman G. Das (Raman Research Institute,
More informationMATH 564/STAT 555 Applied Stochastic Processes Homework 2, September 18, 2015 Due September 30, 2015
ID NAME SCORE MATH 56/STAT 555 Applied Stochastic Processes Homework 2, September 8, 205 Due September 30, 205 The generating function of a sequence a n n 0 is defined as As : a ns n for all s 0 for which
More informationSymmetries in large deviations of additive observables: what one gains from forgetting probabilities and turning to the quantum world
Symmetries in large deviations of additive observables: what one gains from forgetting probabilities and turning to the quantum world Marc Cheneau 1, Juan P. Garrahan 2, Frédéric van Wijland 3 Cécile Appert-Rolland
More informationA scalar conservation law with discontinuous flux for supply chains with finite buffers.
A scalar conservation law with discontinuous flux for supply chains with finite buffers. Dieter Armbruster School of Mathematical and Statistical Sciences, Arizona State University & Department of Mechanical
More informationMicroscopic Deterministic Dynamics and Persistence Exponent arxiv:cond-mat/ v1 [cond-mat.stat-mech] 22 Sep 1999
Microscopic Deterministic Dynamics and Persistence Exponent arxiv:cond-mat/9909323v1 [cond-mat.stat-mech] 22 Sep 1999 B. Zheng FB Physik, Universität Halle, 06099 Halle, Germany Abstract Numerically we
More informationMicroscopic structure of travelling wave solutions in aclass of stochastic interacting particle systems
Microscopic structure of travelling wave solutions in aclass of stochastic interacting particle systems KKrebs 1,FHJafarpour 2,3 and G M Schütz 4 1 Physikalisches Institut, Universität Bonn, 53111 Bonn,
More informationL n = l n (π n ) = length of a longest increasing subsequence of π n.
Longest increasing subsequences π n : permutation of 1,2,...,n. L n = l n (π n ) = length of a longest increasing subsequence of π n. Example: π n = (π n (1),..., π n (n)) = (7, 2, 8, 1, 3, 4, 10, 6, 9,
More informationMicrocanonical scaling in small systems arxiv:cond-mat/ v1 [cond-mat.stat-mech] 3 Jun 2004
Microcanonical scaling in small systems arxiv:cond-mat/0406080v1 [cond-mat.stat-mech] 3 Jun 2004 M. Pleimling a,1, H. Behringer a and A. Hüller a a Institut für Theoretische Physik 1, Universität Erlangen-Nürnberg,
More informationNumerical Analysis of 2-D Ising Model. Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011
Numerical Analysis of 2-D Ising Model By Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011 Contents Abstract Acknowledgment Introduction Computational techniques Numerical Analysis
More informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 22 Oct 2004
The Totally Asymmetric Simple Exclusion Process with Langmuir Kinetics A. Parmeggiani,, T. Franosch,2, and E. Frey,2 arxiv:cond-mat/4834v2 [cond-mat.stat-mech] 22 Oct 24 Hahn-Meitner Institut, Abteilung
More informationOn thermodynamics of stationary states of diffusive systems. Giovanni Jona-Lasinio Coauthors: L. Bertini, A. De Sole, D. Gabrielli, C.
On thermodynamics of stationary states of diffusive systems Giovanni Jona-Lasinio Coauthors: L. Bertini, A. De Sole, D. Gabrielli, C. Landim Newton Institute, October 29, 2013 0. PRELIMINARIES 1. RARE
More informationMATH 1040 Objectives List
MATH 1040 Objectives List Textbook: Calculus, Early Transcendentals, 7th edition, James Stewart Students should expect test questions that require synthesis of these objectives. Unit 1 WebAssign problems
More informationarxiv: v2 [cond-mat.stat-mech] 24 Aug 2014
Hyperuniformity of critical absorbing states Daniel Hexner and Dov Levine, Department of Physics, Technion, Haifa, Israel Initiative for the Theoretical Sciences - CUNY Graduate Center 65 Fifth Avenue,
More informationarxiv:cond-mat/ v1 17 Aug 1994
Universality in the One-Dimensional Self-Organized Critical Forest-Fire Model Barbara Drossel, Siegfried Clar, and Franz Schwabl Institut für Theoretische Physik, arxiv:cond-mat/9408046v1 17 Aug 1994 Physik-Department
More informationBLOW-UP OF COMPLEX SOLUTIONS OF THE 3-d NAVIER-STOKES EQUATIONS AND BEHAVIOR OF RELATED REAL SOLUTIONS
BLOW-UP OF COMPLEX SOLUTIONS OF THE 3-d NAVIER-STOKES EQUATIONS AND BEHAVIOR OF RELATED REAL SOLUTIONS BLOW-UP OF COMPLEX SOLUTIONS OF THE 3-d NAVIER-STOKES EQUATIONS AND BEHAVIOR OF RELATED REAL SOLUTIONS
More informationJ. Stat. Mech. (2016)
Journal of Statistical Mechanics: Theory and Experiment PAPER: Classical statistical mechanics, equilibrium and non-equilibrium Emergence of dynamic phases in the presence of different kinds of open boundaries
More informationGeneralized Manna Sandpile Model with Height Restrictions
75 Brazilian Journal of Physics, vol. 36, no. 3A, September, 26 Generalized Manna Sandpile Model with Height Restrictions Wellington Gomes Dantas and Jürgen F. Stilck Instituto de Física, Universidade
More informationShift in the velocity of a front due to a cutoff
PHYSICAL REVIEW E VOLUME 56, NUMBER 3 SEPTEMBER 1997 Shift in the velocity of a front due to a cutoff Eric Brunet* and Bernard Derrida Laboratoire de Physique Statistique, ENS, 24 rue Lhomond, 75005 Paris,
More informationDynamics of Instantaneous Condensation in the ZRP Conditioned on an Atypical Current
Entropy 2013, 15, 5065-5083; doi:10.3390/e15115065 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article Dynamics of Instantaneous Condensation in the ZRP Conditioned on an Atypical Current
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 informationUniversality class of triad dynamics on a triangular lattice
Universality class of triad dynamics on a triangular lattice Filippo Radicchi,* Daniele Vilone, and Hildegard Meyer-Ortmanns School of Engineering and Science, International University Bremen, P. O. Box
More informationarxiv: v1 [nlin.cg] 27 Mar 2019
arxiv:1903.11319v1 [nlin.cg] 27 Mar 2019 Velocity control for improving flow through a bottleneck 1. Introduction Hiroki Yamamoto 1, Daichi Yanagisawa 1 and Katsuhiro Nishinari 2 1 Department of Aeronautics
More informationUniversal scaling behavior of directed percolation and the pair contact process in an external field
Journal of Physics A 35, 005, (00) Universal scaling behavior of directed percolation and the pair contact process in an external field S. Lübeck, and R. D. Willmann,3 Weizmann Institute, Department of
More informationAPERITIFS. Chapter Diffusion
Chapter 1 APERITIFS Broadly speaking, non-equilibrium statistical physics describes the time-dependent evolution of many-particle systems. The individual particles are elemental interacting entities which,
More informationarxiv: v2 [cond-mat.dis-nn] 24 Oct 2013
arxiv:137.2446v2 [cond-mat.dis-nn] 24 Oct 213 Extinction transition in stochastic population dynamics in a random, convective environment Róbert Juhász Institute for Solid State Physics and Optics, Wigner
More informationarxiv: v1 [nlin.ps] 9 May 2015
Scaling properties of generalized two-dimensional Kuramoto-Sivashinsky equations V. Juknevičius Institute of Theoretical Physics and Astronomy, Vilnius University, A. Goštauto 2, LT-008 Vilnius, Lithuania
More informationarxiv: v2 [physics.soc-ph] 29 Sep 2014
Universal flow-density relation of single-file bicycle, pedestrian and car motion J. Zhang, W. Mehner, S. Holl, and M. Boltes Jülich Supercomputing Centre, Forschungszentrum Jülich GmbH, 52425 Jülich,
More informationAgeing properties of three-dimensional pure and site-diluted Ising ferromagnets
Journal of Physics: Conference Series OPEN ACCESS Ageing properties of three-dimensional pure and site-diluted Ising ferromagnets To cite this article: Vladimir V Prudnikov et al 2014 J. Phys.: Conf. Ser.
More informationBethe Ansatz Solutions and Excitation Gap of the Attractive Bose-Hubbard Model
Journal of the Korean Physical Society, Vol. 39, No., August 001, pp. 03 08 Bethe Ansatz Solutions and Excitation Gap of the Attractive Bose-Hubbard Model Deok-Sun Lee and Doochul Kim School of Physics,
More informationLatent voter model on random regular graphs
Latent voter model on random regular graphs Shirshendu Chatterjee Cornell University (visiting Duke U.) Work in progress with Rick Durrett April 25, 2011 Outline Definition of voter model and duality with
More informationTraffic by multiple species of molecular motors
PHYSICAL REVIEW E 8, 41928 29 Traffic by multiple species of molecular motors Yan Chai, Stefan Klumpp, Melanie J. I. Müller, and Reinhard Lipowsky Max Planck Institute of Colloids and Interfaces, Science
More informationPhase Diagram, Density, And Current Density Profiles Of Totally Asymmetric Simple Exclusion Process For A Junction With Two Entrances and Two Exits
Jurnal Materi dan Pembelajaran Fisika (JMPF) 38 Phase Diagram, Density, And Current Density Profiles Of Totally Asymmetric Simple Exclusion Process For A Junction With Two Entrances and Two Exits M. Za
More informationarxiv: v1 [cond-mat.stat-mech] 20 Sep 2010
JYFL Finite-size effects in dynamics of zero-range processes Janne Juntunen 1, Otto Pulkkinen 2, and Juha Merikoski 1 1 Deparment of Physics, University of Jyväskylä P.O. Box 35, FI-40014 Jyväskylä, Finland
More informationarxiv: v1 [physics.comp-ph] 14 Nov 2014
Variation of the critical percolation threshold in the Achlioptas processes Paraskevas Giazitzidis, 1 Isak Avramov, 2 and Panos Argyrakis 1 1 Department of Physics, University of Thessaloniki, 54124 Thessaloniki,
More informationLecture 6 non-equilibrium statistical mechanics (part 2) 1 Derivations of hydrodynamic behaviour
Lecture 6 non-equilibrium statistical mechanics (part 2) I will outline two areas in which there has been a lot of work in the past 2-3 decades 1 Derivations of hydrodynamic behaviour 1.1 Independent random
More informationarxiv: v1 [cond-mat.stat-mech] 23 Feb 2017
First passage properties in a crowded environment Vincent Tejedor 1 1 Cour des comptes, 13, rue Cambon, 75 001 Paris, France February 27, 2017 arxiv:1702.07394v1 [cond-mat.stat-mech] 23 Feb 2017 Introduction
More informationv n,t n
THE DYNAMICAL STRUCTURE FACTOR AND CRITICAL BEHAVIOR OF A TRAFFIC FLOW MODEL 61 L. ROTERS, S. L UBECK, and K. D. USADEL Theoretische Physik, Gerhard-Mercator-Universitat, 4748 Duisburg, Deutschland, E-mail:
More informationMetropolis Monte Carlo simulation of the Ising Model
Metropolis Monte Carlo simulation of the Ising Model Krishna Shrinivas (CH10B026) Swaroop Ramaswamy (CH10B068) May 10, 2013 Modelling and Simulation of Particulate Processes (CH5012) Introduction The Ising
More informationInconsistencies in Steady State. Thermodynamics
Inconsistencies in Steady State Thermodynamics Ronald Dickman and Ricardo Motai Universidade Federal de Minas Gerais OUTLINE Nonequilibrium steady states (NESS) Steady state thermodynamics (Oono-Paniconi,
More informationDiffusion constant for the repton model of gel electrophoresis
PHYSICAL REVIEW E VOLUME 56, NUMBER 3 SEPTEMBER 1997 Diffusion constant for the repton model of gel electrophoresis M. E. J. Newman Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, New Mexico 87501 G.
More informationTASEP on a ring in sub-relaxation time scale
TASEP on a ring in sub-relaxation time scale Jinho Baik and Zhipeng Liu October 27, 2016 Abstract Interacting particle systems in the KPZ universality class on a ring of size L with OL number of particles
More informationSegregation in a noninteracting binary mixture
Segregation in a noninteracting binary mixture Filip Krzyżewski and Magdalena A. Załuska-Kotur,2 Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland 2 Faculty of
More informationEffect of Diffusing Disorder on an. Absorbing-State Phase Transition
Effect of Diffusing Disorder on an Absorbing-State Phase Transition Ronald Dickman Universidade Federal de Minas Gerais, Brazil Support: CNPq & Fapemig, Brazil OUTLINE Introduction: absorbing-state phase
More informationExact results for deterministic cellular automata traffic models
June 28, 2017 arxiv:comp-gas/9902001v2 2 Nov 1999 Exact results for deterministic cellular automata traffic models Henryk Fukś The Fields Institute for Research in Mathematical Sciences Toronto, Ontario
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 informationDuality to determinants for ASEP
IPS talk Page 1 Duality to determinants for ASEP Ivan Corwin (Clay Math Institute, MIT and MSR) IPS talk Page 2 Mystery: Almost all exactly solvable models in KPZ universality class fit into Macdonald
More informationarxiv: v2 [cond-mat.stat-mech] 11 Aug 2015
arxiv:4.83v2 [cond-mat.stat-mech] Aug 25 Effect of diffusion in simple discontinuous absorbing transition models Salete Pianegonda Instituto de Física, Universidade Federal do Rio Grande do Sul, Caixa
More informationTitle of communication, titles not fitting in one line will break automatically
Title of communication titles not fitting in one line will break automatically First Author Second Author 2 Department University City Country 2 Other Institute City Country Abstract If you want to add
More informationCover times of random searches
Cover times of random searches by Chupeau et al. NATURE PHYSICS www.nature.com/naturephysics 1 I. DEFINITIONS AND DERIVATION OF THE COVER TIME DISTRIBUTION We consider a random walker moving on a network
More informationUniversal examples. Chapter The Bernoulli process
Chapter 1 Universal examples 1.1 The Bernoulli process First description: Bernoulli random variables Y i for i = 1, 2, 3,... independent with P [Y i = 1] = p and P [Y i = ] = 1 p. Second description: Binomial
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Jul 2005
Driven lattice gas with nearest-neighbor exclusion: shear-like drive arxiv:cond-mat/57114v1 [cond-mat.stat-mech] 5 Jul 25 Fabricio Q. Potiguar and Ronald Dickman Departamento de Física, ICEx, Universidade
More informationarxiv:cond-mat/ v2 7 Dec 1999
Molecular dynamics study of a classical two-dimensional electron system: Positional and orientational orders arxiv:cond-mat/9906213v2 7 Dec 1999 Satoru Muto 1, Hideo Aoki Department of Physics, University
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Jun 1997
SISSA Ref. 39/97/EP arxiv:cond-mat/9706114v1 [cond-mat.stat-mech] 11 Jun 1997 First-order phase transitions in one-dimensional steady states Peter F. Arndt, Thomas Heinzel and Vladimir Rittenberg SISSA,
More informationCompetition between multiple totally asymmetric simple exclusion processes for a finite pool of resources
Iowa State University From the SelectedWorks of Beate Schmittmann September 28, 29 Competition between multiple ally asymmetric simple exclusion processes for a finite pool of resources L. Jonathan Cook,
More informationPhase Transition & Approximate Partition Function In Ising Model and Percolation In Two Dimension: Specifically For Square Lattices
IOSR Journal of Applied Physics (IOSR-JAP) ISS: 2278-4861. Volume 2, Issue 3 (ov. - Dec. 2012), PP 31-37 Phase Transition & Approximate Partition Function In Ising Model and Percolation In Two Dimension:
More informationCHINO VALLEY UNIFIED SCHOOL DISTRICT INSTRUCTIONAL GUIDE ALGEBRA II
CHINO VALLEY UNIFIED SCHOOL DISTRICT INSTRUCTIONAL GUIDE ALGEBRA II Course Number 5116 Department Mathematics Qualification Guidelines Successful completion of both semesters of Algebra 1 or Algebra 1
More informationA STUDY ON CONDENSATION IN ZERO RANGE PROCESSES
J. KSIAM Vol.22, No.3, 37 54, 208 http://dx.doi.org/0.294/jksiam.208.22.37 A STUDY ON CONDENSATION IN ZERO RANGE PROCESSES CHEOL-UNG PARK AND INTAE JEON 2 XINAPSE, SOUTH KOREA E-mail address: cheol.ung.park0@gmail.com
More informationAnalytic solution of the Domain Wall initial state. Jacopo Viti. ECT & IIP (UFRN), Natal, Brazil
Analytic solution of the Domain Wall initial state Jacopo Viti ECT & IIP (UFRN), Natal, Brazil Joint work with M. Collura (Oxford Un.) and A. De Luca (Oxford Un.) based on 1707.06218 Background: Non-equilibrium
More informationarxiv: v1 [cond-mat.stat-mech] 1 Feb 2010
Asymmetric Exclusion Process in a System of Interacting Brownian Particles José Eduardo de Oliveira Rodrigues and Ronald Dickman Departamento de Física, ICEx, Universidade Federal de Minas Gerais, Caixa
More informationSupporting Information
Supporting Information A: Calculation of radial distribution functions To get an effective propagator in one dimension, we first transform 1) into spherical coordinates: x a = ρ sin θ cos φ, y = ρ sin
More informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 22 Jan 1998
1 arxiv:cond-mat/9811v [cond-mat.stat-mech] Jan 1998 Investigation of the dynamical structure factor of the Nagel-Schreckenberg traffic flow model S. Lübeck, L. Roters, and K. D. Usadel Theoretische Physik,
More informationScaling exponents for certain 1+1 dimensional directed polymers
Scaling exponents for certain 1+1 dimensional directed polymers Timo Seppäläinen Department of Mathematics University of Wisconsin-Madison 2010 Scaling for a polymer 1/29 1 Introduction 2 KPZ 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 informationSteady-state properties of a totally asymmetric exclusion process with periodic structure
Steady-state properties of a totally asymmetric exclusion process with periodic structure Greg Lakatos, 1, * Tom Chou, 2, and Anatoly Kolomeisky 3, 1 Department of Physics, University of British Columbia,
More informationMATH-1420 Review Concepts (Haugen)
MATH-40 Review Concepts (Haugen) Unit : Equations, Inequalities, Functions, and Graphs Rational Expressions Determine the domain of a rational expression Simplify rational expressions -factor and then
More informationKinetic Monte Carlo. Heiko Rieger. Theoretical Physics Saarland University Saarbrücken, Germany
Kinetic Monte Carlo Heiko Rieger Theoretical Physics Saarland University Saarbrücken, Germany DPG school on Efficient Algorithms in Computational Physics, 10.-14.9.2012, Bad Honnef Intro Kinetic Monte
More informationarxiv: v1 [cond-mat.stat-mech] 1 Jan 2019
Occupation time statistics of a gas of interacting diffusing particles Tal Agranov, 1, P. L. Krapivsky,, and Baruch Meerson 1, 1 Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 91904,
More informationUMUC MATH-107 Final Exam Information
UMUC MATH-07 Final Exam Information What should you know for the final exam? Here are some highlights of textbook material you should study in preparation for the final exam. Review this material from
More informationOPTIMIZATION BY SIMULATED ANNEALING: A NECESSARY AND SUFFICIENT CONDITION FOR CONVERGENCE. Bruce Hajek* University of Illinois at Champaign-Urbana
OPTIMIZATION BY SIMULATED ANNEALING: A NECESSARY AND SUFFICIENT CONDITION FOR CONVERGENCE Bruce Hajek* University of Illinois at Champaign-Urbana A Monte Carlo optimization technique called "simulated
More informationarxiv: v1 [cond-mat.str-el] 22 Sep 2016
arxiv:1609.07191v1 [cond-mat.str-el] 22 Sep 2016 Simulated floating zone method Ryo Ozawa, Yasuyuki Kato, and Yukitoshi Motome Department of Applied Physics, University of Tokyo, Tokyo 113-8656, Japan
More informationSPA07 University of Illinois at Urbana-Champaign August 6 10, 2007
Integral Formulas for the Asymmetric Simple Exclusion Process a Craig A. Tracy & Harold Widom SPA07 University of Illinois at Urbana-Champaign August 6 10, 2007 a arxiv:0704.2633, supported in part by
More informationMarkov chains of nonlinear Markov processes and an application to a winner-takes-all model for social conformity
LETTER TO THE EDITOR Markov chains of nonlinear Markov processes and an application to a winner-takes-all model for social conformity T.D. Frank Center for the Ecological Study of Perception and Action,
More informationPolymer dynamics in repton model at large fields
JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 16 22 APRIL 2004 Polymer dynamics in repton model at large fields Anatoly B. Kolomeisky Department of Chemistry, Rice University, Houston, Texas 77005-1892
More informationarxiv: v1 [nlin.cd] 22 Feb 2011
Generalising the logistic map through the q-product arxiv:1102.4609v1 [nlin.cd] 22 Feb 2011 R W S Pessoa, E P Borges Escola Politécnica, Universidade Federal da Bahia, Rua Aristides Novis 2, Salvador,
More informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 15 Jul 2004
Avalanche Behavior in an Absorbing State Oslo Model Kim Christensen, Nicholas R. Moloney, and Ole Peters Physics of Geological Processes, University of Oslo, PO Box 148, Blindern, N-316 Oslo, Norway Permanent:
More informationFluctuations of TASEP on a ring in relaxation time scale
Fluctuations of TASEP on a ring in relaxation time scale Jinho Baik and Zhipeng Liu March 1, 2017 Abstract We consider the totally asymmetric simple exclusion process on a ring with flat and step initial
More informationModule 16. Diffusion in solids II. Lecture 16. Diffusion in solids II
Module 16 Diffusion in solids II Lecture 16 Diffusion in solids II 1 NPTEL Phase II : IIT Kharagpur : Prof. R. N. Ghosh, Dept of Metallurgical and Materials Engineering Keywords: Micro mechanisms of diffusion,
More informationPhysics 212: Statistical mechanics II Lecture XI
Physics 212: Statistical mechanics II Lecture XI The main result of the last lecture was a calculation of the averaged magnetization in mean-field theory in Fourier space when the spin at the origin is
More information