arxiv: v1 [math.ap] 28 Aug 2018

Transcription

1 ASYMPTOTICS FOR SCALED KRAMERS-SMOLCHOWSKI EQATIONS IN SEVERAL DIMENSIONS WITH GENERAL POTENTIALS arxiv: v1 [math.ap] 8 Aug 18 INSK SEO AND PEYAM TABRIZIAN Abstract. In this paper, we generalize the results of Evans and Tabrizian [3], by deriving asymptotics for the time-rescaled Kramers-Smoluchowski equations, in the case of a general non-symmetric potential function with multiple wells. The asymptotic limit is described by a system of reaction-diffusion equations whose coefficients are determined by the Kramers constants at the saddle points of the potential function and the Hessians of the potential function at global minima. 1. Introduction In this paper, we consider the following Kramers-Smoluchowski equation { (ρ t a x ρ ) = div [ Dρ + 1 ρ DΦ ] in [, T ], ρ = ρ on {t = }, (1.1) where > is a scaling parameter, ρ = ρ (x, ξ, t) is the chemical density, and Φ = Φ(ξ) is a smooth potential function on with multiple wells. This PDE models a simple chemical reaction on the atomic level. For more information on the chemical background, consult [13, 16] and the references therein. Our primary concern is the limiting behavior of ρ when tends to. In this paper, we show that the asymptotic limit of ρ satisfies a system of reaction-diffusion equations. See Theorem.1 for the rigorous formulation of this result. The one-dimensional case d = 1 has already been investigated in [13, 14, 6, 1]. In those works, Φ is assumed to be an even potential function with two wells, and the limit of ρ is derived using tools such as Γ-convergence [13, 14], a Raleigh-type dissipation functional [6], and a Wasserstein gradient flow [1]. We refer to [3, 16] for more detailed survey of the history of the one-dimensional problem. In [3], Evans and Tabrizian developed a new and direct approach for this problem, based on a clever test function that satisfies an elliptic PDE, as well as using capacity estimates from []. The techniques in [3] are robust enough to be generalized in higher dimensions, where Φ is a double-well potential on. The limitation, however, is that it only works for the case where Φ is symmetric. In this paper, we remove the symmetry-assumption and further allow Φ to have more than two wells. In that case, our analysis becomes more delicate, and requires a generalized version of variational principle in [3], which is Theorem 4. of the current paper. 1 Mathematics Subject Classification. 35K15, 35K57, 35J. Key words and phrases. Kramers-Smoluchowski equation, system of reaction-diffusion equations, Morse potential, scaling limit. 1

2 INSK SEO AND PEYAM TABRIZIAN We note that a similar result to the current paper has recently been derived by [1] using tools from semiclassical analysis. We would like to emphasize that the tools developed in Theorem 5.3 are also useful for analyzing metastable random processes, which are processes with multiple stable equilibria. It has been noted in [7, 8, 11, 15] that, by investigating the inhomogeneous version of our main theorem (Theorem 5.3), one can obtain a complete analysis of the metastability of such processes. In addition, this method turns out to be extremely effective in the investigation of metastable diffusions. Two recent papers [11, 15] obtained scaling limits of metastable diffusions known as small random perturbations of dynamical systems. Although such a scaling limit has already been developed for a wide class of metastable Markov chains, it has not been previously known for metastable diffusions. Our paper is organized as follows: In Section, we introduce the detailed model and our assumptions on Φ, as well as the main result of this paper. In Section 3, we derive some preliminary estimates, in Section 4 we state and prove the generalized variational principle mentioned above, and in Section 5 we construct the auxiliary test function. Finally, Section 6 contains the proof of our main result.. Model and Main Result.1. Potential Φ. Let Φ : R be a smooth potential function with multiple minima. In this section, we state our assumptions on Φ, and introduce some notation about the structure of its valleys. First, we assume that Φ(ξ) grows to + as ξ. Furthermore, suppose Φ has exponentially tight level sets, meaning that for all a there exists a constant C(a) > such that e Φ(ξ)/ dξ C(a)e a/ (.1) {ξ:φ(ξ) a} for all (, 1). Note that (.1) is achieved if Φ grows at least linearly as ξ. Moreover, as observed in [, Assumption H.1], (.1) is also valid if lim inf ξ Φ(ξ) = lim inf [ Φ(ξ) Φ(ξ)] =. ξ Now we introduce the inter-valley structure corresponding to the potential function Φ. We refer to Figure 1 for the illustration of the definitions below. We will assume that Φ has finitely many critical points and achieves minimum at several points. This feature can be characterized more precisely by first defining the valleys of Φ. Fix H R and let S = {σ 1, σ,, σ L } be the set of saddle points of Φ with height H, i.e., Φ(σ) = H. Denote by W 1, W,, W K the connected components/valleys of the set {ξ : Φ(ξ) < H}. Assume that W 1 W W K is connected (here A is the closure of the set A). The minimum of Φ on the valley W i, 1 i K, is achieved at m i W i and we suppose that Φ(m 1 ) = Φ(m ) = = Φ(m K ) = h so that valleys W 1, W,, W K have the same depth H h. Hence, m 1,, m K are minima of Φ. Let S i,j = W i W j S ; 1 i j K (.)

3 ASYMPTOTICS FOR SCALED KS EQATION WITH GENERAL POTENTIAL 3 Figure 1. (Left) An example of the potential function Φ with four valleys, i.e., K = 4. (Right) Visualization of the inter-valley structure corresponding to Φ. be the set of saddle points between valleys W i and W j. We select small enough η (, H h) so that there is no critical point ξ of Φ such that Φ(ξ) (H η, H). Fix such η and define = {ξ W i : Φ(ξ) < H η} ; 1 i K. (.3) Then, the set, 1 i K, is connected. Define ( K ) c =. (.4) Finally, we assume that, for each saddle point σ S, the Hessian (Dξ Φ)(σ) has one negative eigenvalue λ σ and (d 1) positive eigenvalues, and for each minimum m i, 1 i K, the Hessian (Dξ Φ)(m i) is non-degenerate... Kramers-Smoluchowski equation. We now describe the scaled Kramers- Smoluchowski equation. Define = 1 e (H h)/ and σ (ξ) = Z 1 e Φ(ξ)/, (.5) where the normalizing factor Z is defined by Z = e Φ(ξ)/ dξ (.6) so that σ dξ = 1. Note that Z R m < because of (.1). Let be a bounded, smooth domain in R n for some n N and let ρ ν = D xρ ν be the outward normal derivative along the boundary. Let a : R n R be a

4 4 INSK SEO AND PEYAM TABRIZIAN smooth and bounded function such that a( ) a > for some constant a. Fix T > and consider the equation [ (ρ t a x ρ ) = div ξ Dξ ρ + 1 ρ D ξ Φ ] in [, T ], ρ ν = on R d [, T ], (.7) ρ = ρ on {t = }. For 1 i K, we write µ i = 1, µ = µ i, and a i = a(m i ). (.8) det(dξ Φ)(m i) For σ S, denote by λ σ the unique negative eigenvalue of the matrix Dξ Φ(σ), and define the Kramers constant at σ by λ σ κ σ =. π det(dξ Φ)(σ) Recall S i,j from (.) and define κ i,j = σ S i,j κ σ ; 1 i j K. (.9) For convenience we set κ i,i = for all 1 i K. Define the rate constants by r i,j = κ i,j ; 1 i j K. (.1) µ i Now we explain our assumptions on the initial data. Consider the normalized initial data u (x, ξ) = ρ (x, ξ) σ. (ξ) Then, we assume that u is bounded on R, is differentiable with respect to x and ξ, and satisfies ( u + D x u + 1τ D ξ u ) σ dxdξ <. (.11) Finally, assume that, for smooth functions α1,, αk following convergence as tends to : : R, we have the u (x, ξ) µ µ i α i locally uniformly in W i ; 1 i K. nder this set of assumptions, we are now ready to state the main result of our paper: Theorem.1. For all t [, T ], we have, in the sense of Remark., ρ (x, ξ, t) dξ α i (x, t) δ mi as, (.1) where the smooth functions α 1,, α K on [, T ] solve the system of linear reaction-diffusion equations given by t α i a i α i = K j=1 (r j,iα j r i,j α i ) in [, T ] α i ν = on [, T ] (.13) α i = αi on t =

5 ASYMPTOTICS FOR SCALED KS EQATION WITH GENERAL POTENTIAL 5 for all 1 i K. Remark.. The weak convergence (.1) means that for all f = f(x, ξ, t) C( A [, T ]), ρ (x, ξ, t)f(x, ξ, t)dξdxdt = α i (x, t)f(x, m i, t)dxdt. lim [,T ] A [,T ] A.3. Graph structure of valleys and an associated Markov chain. The main result described above is closely related to a Markov chain on a graph whose vertices are the valleys of potential Φ. More precisely, denote by V = {1,,, K} the set of vertices, in such a way that i V corresponds to the valley. Moreover, two vertices i, j V are connected by an edge if and only if W i W j φ, or equivalently κ i,j. Denote by the G the resulting graph. Since we have assumed that the set W 1 W W K is connected, the graph G is a connected graph. Let {X t : t } be a Markov chain on V where the jump rate from i V to j s r i,j (cf. (.1)). Since r i,j = if κ i,j =, X t becomes a Markov chain on G. Define µ i =: µ i µ for 1 i K and µ := ( µ 1,, µ K ). Then, observe that the probability measure µ on s the invariant measure for the Markov chain X t, and furthermore, the Markov chain is reversible with respect to µ in the sense that µ i r i,j = µ j r j,i for all i j. The generator L of this Markov chain can be regarded as a linear operator on R K. More precisely, for b = (b 1,, b K ) R K, the ith component of Lb R K is given by (Lb) i = r i,j (b j b i ). j=1 Remark.3. Assume that a so that α i, 1 i K, is a function of time only. Then, define α i (t) = α i (t)/ µ i, and let α(t) = ( α 1 (t),, α K (t)) R K. Then, we can deduce from (.13) that d α i dt (t) = K j=1 r i,j ( α j (t) α i (t)) = (L α(t)) i. Therefore, α(t) = (α 1 (t),, α K (t)) is the marginal density of the Markov chain X t with respect to the invariant measure µ, whose starting (possibly deterministic) measure is (α 1,, α K ). 3. Preliminary Estimates In this section, we state and prove estimates. Denote by o (1) the term vanishing as. Lemma 3.1. We have that e Φ(ξ)/ dξ = [1 + o (1)] e h/ (π) d/ µ i ; 1 i K, (3.1) V i e Φ(ξ)/ dξ = o (1)e h/ d/, and (3.) Z = [1 + o (1)] e h/ (π) d/ µ. (3.3)

6 6 INSK SEO AND PEYAM TABRIZIAN Proof. The proof of (3.1) is an easy consequence of Laplace s method. The estimate (3.) is a direct consequence of (.1). Finally, (3.3) follows immediately from (3.1) and (3.) because of the definition of Z (cf. (.6)). Lemma 3.. For A, suppose that there exists c > such that Φ(ξ) H + c for all ξ A. Then, σ dξ = o (1). Proof. By (3.3), A σ = [1 + o (1)] Hence, the lemma immediately follows from (.1). (π) d/ µ e(h Φ)/. (3.4) Now we establish several compactness estimates similar to [3, Section 3]. Let u (x, ξ, t) := ρ (x, ξ, t) σ (ξ) Then, by (.7), the u satisfies u t a x u = 1 [ ] σ σ div ξ D ξ u. (3.5) The next lemma is an energy estimate that is similar to that of [3, Lemma 3.1]. However, instead of skipping the proof, we refer the readers to the Appendix, since the notation here is more involved than [3]. Lemma 3.3. For some constant C >, we have the bound and the energy estimate ( sup t T. u C on [, T ] (3.6) u + D x u + τ 1 D ξ u ) σ dxdξ T + u t σ dxdξdt C. (3.7) Define T = (, T ). We next develop some pre-compactness results similar to [3, Lemmas 3. and 3.3]. Again, proofs can be found in the Appendix since they are more involved. Lemma 3.4. There exist a sequence { n } n=1 of positive real numbers converging to and functions α 1, α,, α K H 1 ( T ) that satisfy the following: (1) For all 1 i K, we have that, as n, ρ n (x, ξ, t) dxdξ α i (x, t) weakly in L ( T ) and (3.8) ρ n (x, ξ, t) dξ. (3.9) sup t T

7 ASYMPTOTICS FOR SCALED KS EQATION WITH GENERAL POTENTIAL 7 () For all 1 i K, we have that, as n, t ρ n (x, ξ, t) dξ t α i (x, t) weakly in L ( T ), (3.1) V i D x ρ n (x, ξ, t) dξ D x α i (x, t) weakly in L ( T ), (3.11) V i t ρ n (x, ξ, t) dxdξ strongly in L (, T ), and (3.1) sup t T D x ρ n (x, ξ, t) dxdξ. (3.13) (3) For all t [, T ], for all 1 i K, and almost every x, we have that, as, u (x, ξ, t) α i(x, t) µ i for almost every ξ. (3.14) 4. A variational Problem Throughout the rest of the paper, elements of R K are denoted by bold lowercase letters such as a =(a 1,, a K ), and subsets of R K are denoted by bold capital letters like A and B. Define D : R K R by D(b) = 1 µ κ i,j (b j b i ) ; b R K. (4.1) i,j=1 Note that D(b) = implies b 1 = b = = b K since the graph G is connected. Remark 4.1. The function D is the so-called Dirichlet form associated with the generator L defined in Section.3. More precisely, we can write D(b) = µ i b i ( Lb) i. For b = (b 1,, b K ) R K, define { F b = ψ H 1 ( ) : ψ } Vi b i for all 1 i K. (4.) In the current and the next section, we only consider functions on, that is only depending on ξ and independent of the variable x. Hence, for a function φ : R, the notations Dφ and φ are used to represent D ξ φ and ξ φ, respectively. Then the following result is a generalization of [, Theorem 3.1]. Theorem 4.. For any b R K, we have that Rd σ inf Dψ dξ = [1 + o (1)] D(b). (4.3) ϕ F b Proof. By (3.4) and definition of D( ) we can rewrite the identity (4.3) as Rd inf e Φ/ Dψ H/ (π)d/ dξ = [1 + o (1)] e ϕ F b κ i,j (b j b i ). (4.4) i,j=1

8 8 INSK SEO AND PEYAM TABRIZIAN Denote by ϕ b the minimizer of the left-hand-side. Then, ϕ b solves the following Euler-Lagrange equation: [ ] div e Φ/ Dϕ b = on and ϕ b = b i on for all 1 i K. For 1 i d, write e i = (,,, 1,,, ) the ith standard basis vector of. Then, by linearity and uniqueness of the Euler-Lagrange equation, it follows that Therefore, we can write e Φ/ Dϕ b dξ = ϕ b = d b i ϕ e i. (4.5) b i e Φ/ Dϕ e i dξ + b i b j e Φ/ D(ϕ e i + ϕ e j ) dξ (4.6) 1 i j K In [, Theorem 3.1], it is shown that e Φ/ Dϕ e i dξ = [1 + o (1)] e H/ (π) d/ κ i,l. (4.7) and that, for i j, e Φ/ D(ϕ e i + ϕ e j ) dξ = [1 + o (1)] e H/ (π) d/ 1 l K : l i,j l=1 (κ i,l + κ j,l ). (4.8) By (4.7) and (4.8), we have that e Φ/ Dϕ e i Dϕ e j dξ = [1 + o (1)] e H/ (π) d/ κ i,j. (4.9) We can complete the proof by combining (4.6), (4.7) and (4.9). 5. Construction of the Test Function 5.1. Preliminaries. Let M be the symmetric K K matrix defined by { 1 K µ l=1 M ij = κ i,l if i = j 1 µ κ ; 1 i, j K, i,j if i j so that Define two subsets of R K by D(x) = x T Mx. (5.1) N = {x R K : x 1 = x = = x K }, R = {x R K : x 1 + x + + x K = }. Lemma 5.1. The null-space and range of the matrix M are N and R respectively.

9 ASYMPTOTICS FOR SCALED KS EQATION WITH GENERAL POTENTIAL 9 Proof. Suppose that Mb =. Then, by (5.1) we have D(b) =. Hence, b N as we observed in the line following (4.1). On the other hand, any b N satisfies Mb =. Hence, the null-space of M is N. Since the dimension of the null-space is 1, that of the range of M must be (K 1) dimensional. Since K (Mb) i = for all b R K, the range of M is a subset of R. Since dim(r) = K 1, we can conclude that R is the range of M. For c R, write M 1 c = {b R K : Mb = c} Then, for b M 1 c, we can write M 1 c = N + b. Hence, we can observe that D( ) is a constant function on M 1 c. Now define a function D c : R K R by Then, by (5.1), for b M 1 c, D c (x) = D(x) c x. (5.) D c (b) = b Mb c b = b Mb = D(b), (5.3) and hence the function D c ( ) restricted to M 1 c is a constant function as well. Let us denote that constant by D c (M 1 c), with slight abuse of notation. Lemma 5.. Fix c R. Then, D c (M 1 c) is the minimum of D c ( ). Furthermore, if D c (x) D c (M 1 c) + δ for some δ >, then there exists b M 1 c such that x b C δ for some constant C > not depending on δ. Proof. Let b M 1 c. Then, since M is symmetric, it is easy to observe that for any x, D c (x) = D c (b) + D(t). where t := x b. The first part of the lemma follows since D is a non-negative function. As for the second part, we must have D(t), and hence, by continuity of D and the fact that the nullspace of D is N, we can find t R such that t i t C δ for all 1 i K, where t i is such that t = (t 1,, t K ). Then, b := b + (t, t,, t) b + N = M 1 c fulfills the requirement of the second part of the lemma. 5.. Test function. Denote by χ A ( ) the indicator function of the set A. We emphasize that the following construction of the test function ψ is the main ingredient in the proof of Theorem.1, and contains most of technical difficulties of the problem. Theorem 5.3. Fix a non-zero vector c R and b M 1 c. Then, for each >, there exists a function ψ W,p loc (Rd ) L ( ) for all p [1, ) that satisfies the equation ( ) σ div Dψ c i = χ, (5.4) and the uniform energy estimate and finally, lim sup <<1 sup 1 i K σ Dψ < (5.5) sup ψ (ξ) b i =. (5.6) ξ

10 1 INSK SEO AND PEYAM TABRIZIAN The proof of this theorem is divided into several lemmas. We start by simplifying the problem and by introducing relevant notions before starting these lemmas. By linearity, it suffices to prove the theorem for c = e i e j for some i j. Therefore, without loss of generality, we assume that c = e 1 e = (1, 1,,, ), so that c 1 = 1, c = 1, c i = for i 3. For φ Hloc 1 (Rd ), define a functional I by I[φ] = 1 Rd σ Dφ dξ 1 V 1 V 1 φ dξ + 1 V V φ dξ, (5.7) and let φ be a minimizer of I[φ] on Hloc 1 (Rd ). Then the Euler-Lagrange equation for φ is (5.4) for c = e 1 e R K, and moreover φ W,p loc (Rd ) for all p [1, ). For 1 i K, define λ,i = 1 φ dξ. Since I[φ ] = I[φ + c] for all c R, we can assume without loss of generality that λ,1 = λ, := λ. Note that λ since otherwise we can replace φ with φ. Let µ,i := sup ξ Vi φ (ξ) and define µ := max{µ,1, µ, }. (5.8) Then we can assume that sup φ (ξ) = µ, ξ since otherwise, φ = Λ(φ ) gives a lower value of I, where µ if s (, µ ) Λ(s) = s if s [ µ, µ ] µ if s (µ, ) With the simplification and notations above, we now start the proof of the Theorem 5.3. The first step is the following lemma. Lemma 5.4. We have that σ Dφ dξ = λ. Proof. First observe that, since φ is a minimizer, I[φ ] I[], and so Rd σ Dφ dξ 4 λ. (5.9) For R >, let B R := {ξ : ξ R} and let ζ R : [, 1] be a smooth cutoff function with compact support such that ζ R 1 on B R, and Dζ R 1. Let us select R large enough so that B R for all 1 i K and Φ(ξ) > H + 1 on (B R ) c Then, multiplying (5.4) by ζ R φ and integrating by parts, we obtain Rd σ Dφ ζ ξ = λ Rd σ φ Dφ Dζ ξ. (5.1) Because Dζ R 1, the square of the last term is bounded by (Rd σ ) ( ) (BR)c Dφ σ dξ dξ. (5.11)

11 ASYMPTOTICS FOR SCALED KS EQATION WITH GENERAL POTENTIAL 11 Note that by the assumption Φ(ξ) > H + 1 on (B R ) c and by Lemma 3., the last integral converges to as R Hence, by our priori bound (5.9), the expression in (5.11) vanishes as R. Hence, the proof is completed by letting R in (5.1). Recall the definition of from (.3). Let us take < η < η and let Ṽi, 1 i K, be the connected component of {ξ W i : Φ(ξ) < H η } containing. Then, we can obtain the following L -estimate for φ λ, i on the extended valley Ṽi, for all 1 i K. Lemma 5.5. For all 1 i K, it holds that Proof. Define φ λ,i L (Ṽi) = o (1)λ 1. λ,i = 1 φ dξ ; 1 i K. Ṽi Ṽ i By using Poincaré s inequality, as well as (3.4) and Lemma 5.4, we get that for all 1 i K, Ṽ i φ λ,i dξ C Hence, we can derive Ṽ i Dφ dξ C (d/) 1 e η Ṽi σ λ,i λ 1,i = (φ λ,i ) dξ C φ [Ṽ λ,i dξ i Now, combining (5.1) and (5.13) completes the proof of lemma. Dφ dξ C (d/) 1 e η λ. (5.1) ] 1 = o (1) λ 1. (5.13) The next step is to enhance the previous L -estimate on extended valley Ṽi to the L -estimated on the original valley. The proof is based on the elliptic estimate on φ, and on a bootstrapping argument. Let us fix p (d, ) from now on, and regard p just as a constant. Lemma 5.6. For all 1 i K, it holds that, φ λ,i L () o (1) ( 1 + λ 1 1 p ). Proof. We fix 1 i K. On the set W i, the function φ satisfies ( ) σ div Dφ = c i χ. This equation can be rewritten as (φ λ,i ) = 1 div [(φ λ,i ) DΦ] + 1 (φ λ,i ) Φ + σ c i χ,

12 1 INSK SEO AND PEYAM TABRIZIAN and, therefore, standard regularity estimates for elliptic PDE (cf. 8.17]) imply [5, Theorem φ λ,i L () C φ λ,i L (Ṽi) + C φ λ,i L p (Ṽi) + C σ L () since p (d, ). By Lemma 5.5 and Hölder s inequality, (similar argument to [3, (3.31)]) we obtain [ ] φ λ,i L (V o i) (1) λ 1 + C φ λ,i 1 p λ 1 p + 1. (5.14) L (Ṽi) Recall the definition of µ from (5.8) and write µ = µ,k where k is either 1 or. Then, we obtain φ λ,i L (Ṽi) µ = µ,k λ,k + φ λ,k L (V k ) = λ + φ λ,k L (V k ). (5.15) By inserting this result into (5.14) with i = k, we derive [ φ λ,k L (V k ) o (1) 1 + φ λ,k 1 p L (V k ) λ 1 p Therefore, by Hölder s inequality, we conclude that φ λ,k L (V k ) By inserting (5.16) into (5.15), we obtain, φ λ,i L (Ṽi) o (1) ( 1 + λ 1 1 p λ + o (1) ( 1 + λ 1 1 p + λ 1 1 p + λ 1 ]. Finally, the proof of lemma is completed by inserting this into (5.14). ). (5.16) ). (5.17) In view of the previous lemma, it is important to prove that λ is bounded by a constant for small enough. Indeed, we are able to establish more than this, as in the following lemma. The following lemma is the most renovative part of the current paper. Lemma 5.7. We have that, lim λ = 1 D c(b) = 1 D(b). Proof. Recall b M 1 c from the statement of theorem. Hence, by (5.3), the second identity of the lemma is obvious. Now we focus on the first identity of the lemma. Let ϕ b be the minimizer of the variational problem on the left-hand-side of (4.3). Since I[φ ] = λ by Lemma 5.4 and since φ is the minimizer of I[φ], by Theorem 4., we obtain λ = I[φ ] I[ϕ b] = 1 Rd σ Dϕ τ b dξ c i b i = 1 D c(b) + o (1). (5.18) Therefore, we get lim inf λ 1 D c(b) = 1 D(b) >, (5.19) where the last inequality is strict since c. This proves the half of the first identity.,

13 ASYMPTOTICS FOR SCALED KS EQATION WITH GENERAL POTENTIAL 13 Figure. Illustration of the piecewise linear function Γ(x): In this example, B 1,3 and B 4,5 are good sets, since (β 1, β, β 3) and (β 4, β 5) form cluster in the sense defined in (5.1), and accordingly, we have j 1 = 3 and j = 5. In (5.) we show that the slope of each piece is bounded by 1 + o L(1). We now have to prove the reverse inequality, namely, lim sup λ 1 D c(b). (5.) This is the crux of the proof. Let β 1 β β K be the enumeration of the numbers λ,1 λ λ,,,k λ. Strictly speaking, β i = β i,, but we will ignore the dependency on for the time being. Fix δ > and L >. We now introduce an auxiliary function Γ(x). We refer to Figure for the visualization of the construction below. For 1 i j K, we say that B i,j = {β k : i k j} is a good set if { β k+1 β k L k δ for i k j 1, β k+1 β k > L k (5.1) δ for k = i 1 and k = j, where β := and β K+1 :=. Enumerate all good sets by B j+1,j 1, B j1+1,j,, B jm 1 +1,j M, where j := and j M := K. For 1 k M, define I k = [ β jk 1 +1 δ, β jk + δ, ]. We now define a piecewise linear function Γ = Γ,δ,L : R R by: x + δ if x (, β 1 δ) β jk 1 +1 if x I k, 1 k M Γ(x) = β jk β j k +1 β jk 1 +1 β jk +1 β jk δ (x (β if x (β j k + δ)) jk + δ, β jk +1 δ), 1 k M 1 x (β K + δ) + β jm 1 +1 if x (β K + δ, ).

14 14 INSK SEO AND PEYAM TABRIZIAN By construction, the function Γ is continuous. Now we estimate the slope of Γ on the interval (β jk + δ, β jk +1 δ). By (5.1), β jk +1 β jk 1 +1 β jk +1 β jk δ = β jk 1 j k+1 β jk + i=j k 1 +1 (β i+1 β i ) β jk +1 β jk δ < Lj k δ + j k 1 i=j k 1 +1 Li δ L j k δ δ Hence, since L >, we conclude that < Ljk + KL jk 1 L j k. Γ 1 + K + L = 1 + o L(1), (5.) where o L (1) is a term vanishing when L and independent of and δ. Finally, observe that Γ is constant on intervals of the form [(λ,i /λ ) δ, (λ,i /λ ) + δ], 1 i K, and denote that constant by γ i = γ,δ,l,i. Then, we obtain, γ i λ,i λ max (β j k β jk 1 +1) 1 k M max j k 1 1 k M i=j k 1 +1 L i δ KL K δ. (5.3) Define g := (γ 1, γ,, γ K ). Then, since g c = γ 1 γ and since λ,1 = λ, = λ, we deduce from (5.3) that Now, define g c KL K δ. (5.4) φ = φ,δ,l := ( ) φ Γ. Then, in view of Lemma 5.6, (5.19), and (5.3), we have that φ F g (cf. (4.)) for all sufficiently small. Assume from now on that is small enough so that this condition is valid. (Hence, we should send before taking any limit for δ or L.) Then, by Theorem 4. and (5.4), 1 Rd σ D φ dξ 1 + o (1) D(g) = 1 + o (1) λ D(λ g) = 1 + o [ ] (1) 1 λ D c(λ g) + λ (g c) 1 + o [ ] (1) 1 λ D c(λ g) + λ ( KL K δ). On the other hand, by (5.) and by Lemma 5.4, 1 Rd σ D φ dξ 1 + o L(1) λ Combining those two inequalities, we obtain 1 + o L (1) λ 1 + o (1) λ Rd σ λ Dφ dξ = 1 + o L(1) λ. [ ] 1 D c(λ g) + λ ( KL K δ) We select small enough so that the o (1) term is greater than 1. Then, we can re-organize the previous inequality as [ (1 KL K δ) + o (1) + o L (1) ] λ 1 D c(λ g) 1 D c(b), (5.5).

15 ASYMPTOTICS FOR SCALED KS EQATION WITH GENERAL POTENTIAL 15 since b = M 1 c is the minimizer of D c by Lemma 5.. Now take L large enough so that o L (1) < 1/3 and then take δ small enough so that KL K δ < 1/3. In this way, the quantity in brackets converges to a positive number as. Then, by taking lim sup in the previous inequality, we obtain [ 1 KL K δ o L (1) ] lim sup λ 1 D c(b). Finally, send δ and then L to conclude (5.). This finished the proof. As a direct consequence of the previous lemma, we obtain the following boundedness results. Corollary 5.8. There exist constants C > and > such that, for all (, ), λ C and λ,i C ; 1 i K. Proof. The first inequality is immediate from Lemma 5.7. By this inequality and Lemma 5.6, we can conclude that µ λ = o (1). This implies the second inequality of the corollary since λ,i µ λ + o (1). Now we arrive the last ingredient for the proof of Theorem 5.3 Lemma 5.9. Define l = (λ,1, λ,,, λ,k ). Then, we have that, lim D c(l ) = D c (b). Proof. By (5.3) and the boundedness of λ,i obtained in the previous corollary, we have D c (λ g) D c (l ) C(L K δ + L K δ) (5.6) for some constant C >. Combining this bound and the first inequality of (5.5) yields D c (l ) [ 1 KL K δ + o L (1) + o (1) ] λ + C(L K δ + L K δ). Thus, by Lemma 5.7 lim sup D c (l ) [ 1 KL K δ o L (1) ] D c (b) + C(L K δ + L K δ) By letting δ and then L, we deduce lim sup D c (l ) D c (b). (5.7) On the other hand, we know by Lemma 5. that D c (l ) D c (b) and thus, lim inf D c(l ) D c (b). (5.8) Hence, by (5.7) and (5.8), we can finish the proof of lemma. Now we arrived at the final stage of the proof. Proof of Theorem 5.3. By Lemma 5.9, we can write D c (l ) = D c (b) + o (1). By the second part of Lemma 5., there exists t R such that (λ i, t ) b i = o (1) for all 1 i K. (5.9) We now define ψ := φ t, and claim that ψ satisfies all the requirements of the theorem.

16 16 INSK SEO AND PEYAM TABRIZIAN First, the condition (5.4) holds for φ since this function is chosen as a minimizer if the functional I defined in (5.7). Hence, the condition (5.4) is also valid for ψ as well since t is merely a constant so that Dφ = Dψ. Second, we can deduce (5.5) for φ by combining Lemma 5.4 and Corollary 5.8. By the same reason as above, the condition (5.5) also holds for ψ as well. Finally, the condition (5.6) is immediately follows from (5.9). 6. Proof of Theorem.1 The proof of Theorem.1 is similar to the proof of [3, Theorem 3.7] and we present it below for sake of completeness. Proof of Theorem.1. Define A 1 = {ξ : Φ(ξ) H + 1} and A = {ξ : Φ(ξ) H + }. Let ζ : [, 1] be a smooth cutoff function such that ζ 1 on A 1 and ζ on (A ) c. Fix b R K \ N and let c = Mb( ). Then, denote by ψ the function in Theorem 5.3 with c R and b M 1 c Let f = f(x, t) C ( T ) be a smooth test function. Multiplying (3.5) by ζfψ and integrating by parts, we obtain T A ζfψ ρ t dξdxdt + T = T ψ ζa D x f D x ρ dξdxdt A A σ fd ξ u D ξ (ζψ ) dξdxdt. (6.1) Now we consider three integrals in (6.1) separately. Write A = V 1 V K. Then, the first integral of (6.1) can be split into T T ζfψ ρ t dξdxdt + ζfψ ρ t dξdxdt. (6.) A Since ζ 1 on A, by (3.1) and by (5.6), we get T A ζfψ ρ t dξdxdt = [1 + o (1)] T A \A T fb i ρ t dξdxdt f b i t α i dxdt. (6.3) The second term of (6.) becomes negligible since, by (3.1), T T ζfψ ρ t dξdt C ρ t dξdxdt. (6.4) A \A Now we consider the second integral of (6.1). Similarly, we split it into an integral on A and A \ A respectively. Then, by (3.13), it is easy to verify that the integral on A \ A vanishes as, while by (3.11) the integral on A converges to T as. D x f T a i b i D x α i dxdt = f a i b i x α i dxdt (6.5)

17 ASYMPTOTICS FOR SCALED KS EQATION WITH GENERAL POTENTIAL 17 Finally, integrating by parts again, the last term in (6.1) becomes T A σ f [ ψ D ξ u + u D ξ ψ ] D ξ ζ dξdxdt T [ ] σ + fζu div ξ D ξ ψ dξdxdt. (6.6) A We claim that the first term is negligible. To this end, since D ξ ζ on A 1, by (3.6), by (5.5), the square of this integral is bounded by T σ T σ [ C dξdxdt D ξ u + D ξ ψ ] dξdxdt. A c 1 A c 1 In the last expression, the first integral vanishes as by Lemma 3., and the second integral is bounded because of (3.7) and (5.5). This proves the claim. On the other hand, by Theorem 5.3 and (3.14), the second integral of (6.6) converges to T fc i α i µ i dxdt = T f α i r i,j (b j b i ) dxdt. (6.7) i,j=1 By combining (6.1), (6.3), (6.5) and (6.7), we obtain T f b i ( t α i a i x α i ) dxdt = T f α i r i,j (b j b i ) dxdt. i,j=1 Now we select b = e i for 1 i K. Then, the previous identity implies that This completes the proof. t α i a i x α i = (r j,i α j r i,j α i ). j=1 Appendix: Proof of Lemmas 3.3 and 3.4 Proof of Lemma 3.3. First of all, (3.6) follows from the assumption u C and from the maximum principle applied to (3.5) As for the energy estimate (3.7), multiplying (3.5) by u and integrating over [, t] (where t is fixed), we get t σ u t u dξdxdt t σ a ( x u ) u dξdxdt = 1 t sing σ u tu = 1 σ t u, the first term becomes t σ u t u dξdxdt = 1 σ u (x, ξ, t) dξdx 1 div ξ [σ D ξ u ] u dξdxdt. σ u dξdx.

18 18 INSK SEO AND PEYAM TABRIZIAN Applying the divergence theorem with respect to x (note that there are no boundary terms since u ν = on Rd (, t)), the second term becomes t t σ a ( x u ) u dξdxdt = σ a D x u dξdxdt. Lastly, integrating by parts with respect to ξ (there are again no boundary terms), the term on the right becomes 1 t div ξ [σ D ξ u ] u dξdxdt = 1 t σ D ξ u dξdxdt. Putting everything together, we obtain 1 σ u (x, ξ, t) dξdx + t σ (a D x u + 1 ) D ξ u dξdxdt = 1 σ u dξdx. By our assumption (.11) on the initial data u, the right-hand-side is bounded, and therefore, taking the supremum over t [, T ], we obtain T ( σ u dξdx + σ a D x u + 1 ) D ξ u dξdxdt C. sup t T This gives us one part of our desired estimate; to obtain the other part, multiply (3.5) by u t and integrate to obtain t σ u t dξdxdt t = 1 t a σ ( x u ) u t dξdxdt div ξ [σ D ξ u ] u t dξdxdt. (6.8) The first term stays as it is; as for the second integral, integrating by parts with respect to x and using D x u D x u t = 1 t D x u, we can deduce that it equals to 1 a σ D x u (x, ξ, t) dxdξ 1 a σ D x u dxdξ. Similarly, for the last term of (6.8), integrating by parts with respect to ξ, we can rewrite it as 1 σ D ξ u (x, ξ, t) dxdξ + 1 σ D ξ u τ dxdξ. Putting everything together, we get t σ u t dξdxdt + 1 ( σ a D x u (x, ξ, t) + 1 ) D ξ u (x, ξ, t) dxdξ = 1 ( σ a D x u + 1 ) D ξ u τ dxdξ.

19 ASYMPTOTICS FOR SCALED KS EQATION WITH GENERAL POTENTIAL 19 Now again, by our assumption (.11) on the initial condition, the right-hand side is bounded, and finally taking the supremum over t, we obtain ( sup σ a D x u (x, ξ, t) + 1 ) D ξ u (x, ξ, t) dxdξ t T T + σ u t dξdxdt C, which, combined with the above and the fact that a a >, gives our desired estimate. Proof of Lemma 3.4. Writing ρ = u σ = ((u ) σ ) 1 (σ ) 1, we get ( ( ) ( ) sup ρ dxdξ) sup u σ dxdξ σ dξ. t T t T This follows because the first term on the right-hand-side is bounded by Lemma 3.3, and because the second term on the right-hand-side goes to by (3.) and (3.3). Hence (3.9) follows. Similarly, we deduce sup t T ( D x ρ dxdξ) sup t T from which (3.13) follows, as well as T ( T ( ρ t dxdξ) dt ( from which (3.1) follows. Now define αi = α i (x, t) by αi(x, t) := ρ (x, ξ, t) dξ = D ) ( ) D x u σ dxdξ σ dξ ) ( ) u t σ dxdξ σ dξ dt u (x, ξ, t) σ dξ. In the same way as above, but this time using that σ dξ σ = 1, we get T T T α i dxdt α i,t dxdt T D x α i dxdt T ( T ) ( ) u σ dxdξ σ dξ dt C, ) ( ) u t σ dxdξ σ dξ dt C, and ( ) ( ) D x u σ dxdξ σ dξ dt C. ( Hence for each i, {αi } is bounded in H1 ( [, T ]), a reflexive Banach space, and so by weak compactness, we can extract a subsequence { n } n=1 with n as n, such that, for some limit functions α i = α i (x, t), we have α n i α i weakly in H 1 ( [, T ]) as n. The results (3.8), (3.1), (3.11) then follow by construction.

20 INSK SEO AND PEYAM TABRIZIAN Finally, for (3.14), notice that D ξ u dxdξ ( σ ) 1 ( ) 1 D ξ u dxdξ Vi σ dξ ( C To show that the last integral converges to as, we have σ dξ = 1 H h e [1 + o (1)] e h d (π) µ e Φ dξ 1 e H ) 1 σ dξ [1 + o (1)] (π) d µ e H η = 1 e η [1 + o (1)] (π) d µ Vi, where the first identity follows from the definitions of, σ, and Z, and (3.3). Since η > the last term converges to as. Therefore, it follows that, on, u u i a.e. for some function u i = u i (x, t). But using ρ = σ u, integrating with respect to ξ on and using σ dξ = ˆµ i, we finally obtain u i = αi ˆµ i. Acknowledgments. The authors wish to thank Professor Lawrence C. Evans for his fruitful discussions. The research of I. Seo is supported by the National Research Foundation of Korea NRF grant funded by the Korean government MSIT (Project 18R1C1B66896). References [1] S. Arnrich, A. Mielke, M. A. Peletier, G. Savare, and M. Veneroni: Passage to the limit in a Wasserstein gradient flow: From diffusion to reaction, Calc. Var. Partial Differential Equations, 44 (1), pp [] A. Bovier, M. Eckhoff, V. Gayrard, and M. Klein: Metastability in reversible diusion processes. I. Sharp asymptotics for capacities and exit times, J. Eur. Math. Soc., 6 (4), pp [3] C. Evans and P. Tabrizian: Asymptotics for scaled Kramers-Smoluchoswski equations. Siam J. Math. Anal. 48 (16), pp [4] M. I. Freidlin, A. D. Wentzell: Random Perturbations. In: Random Perturbations of Dynamical Systems. Grundlehren der mathematischen Wissenschaften 6. Springer, New York, NY, 1998 [5] D. Gilbarg and N. Trudinger, Elliptic Partial Differential Equations of Second Order, nd ed. Springer, [6] M. Herrmann and B. Niethammer: Kramers formula for chemical reactions in the context of Wasserstein gradient flows, Commun. Math. Sci., 9 (11), pp [7] C. Landim: Metastable Markov chains. arxiv: (18) [8] C. Landim: Variational Formulae for the capacity induced by second-order elliptic differential operators. Proc. Int. Cong. of Math. 18 Rio de Janeiro, Vol. (18) pp [9] C. Landim, M. Mariani and I. Seo: A Dirichlet and a Thomson principle for non-selfadjoint elliptic operators, Metastability in non-reversible diffusion processes. (17). submitted. [1] C. Landim, R. Misturini, K. Tsunoda: Metastability of reversible random walks in potential field. J. Stat. Phys. 16, (15) [11] C. Landim and I. Seo: Metastability of one-dimensional, non-reversible diffusions with periodic boundary conditions. (17). submitted. [1] L. Michel and M. Zworski: A semiclassical approach to the Kramers-Smoluchowski equation. (17). [13] M. A. Peletier, G. Savare, and M. Veneroni: From diffusion to reaction via Γ convergence, SIAM J. Math. Anal., 4 (1), pp [14] M. A. Peletier, G. Savare, and M. Veneroni: Chemical reactions as Γ limit of diffusion, SIAM Rev., 54 (1), pp

21 ASYMPTOTICS FOR SCALED KS EQATION WITH GENERAL POTENTIAL 1 [15] F. Rezakhanlou and I. Seo: Scaling Limit of Metastable Diffusion Processes. Preprint. (18). [16] P. R. Tabrizian: Asymptotic PDE models for Chemical Reactions and Diffusions. Ph.D. Thesis. niversity of California, Berkeley. (16). Insuk Seo, Seoul National niversity, Department of Mathematical Sciences and Research Institute of Mathematics, GwanAk Ro 1, Seoul 886, Republic of Korea. insuk.seo@snu.ac.kr Peyam Tabrizian, niversity of California, Irvine, Department of Mathematics, 34 Rowland Hall, Irvine, CA 9697, SA. ptabrizi@uci.edu