{'fraction_non_alphanumeric': 0.05812599832146628, 'fraction_numerical': 0.01851801716436094, 'mean_word_length': 3.7320010248526776, 'pattern_counts': {'":': 0, '<': 20, '<?xml version=': 0, '>': 21, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 15, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Stability of gradient Kaehler-Ricci solitons

We study stability of non-compact gradient Kaehler-Ricci flow solitons with positive holomorphic bisectional curvature. Our main result is that any compactly supported perturbation and appropriately decaying perturbations of the Kaehler potential of the soliton will converge to the original soliton under Kaehler-Ricci flow as time tends to infinity. To obtain this result, we construct appropriate barriers and introduce an L^p-norm that decays for these barriers with non-negative Ricci curvature.


Introduction
In [16], Hamilton introduced the Ricci flow to find Einstein metrics on a compact Riemannian manifold by evolving an existing metric by its negative Ricci curvature. By solving a complex Monge-Ampère equation Yau [22] found Einstein metrics on compact Kähler manifolds. The results of Yau were later re-established by Cao [6], using the Kähler-Ricci flow d dt g i = −R i (1.1) to find Einstein metrics on a compact Kähler manifold. On C n this equation corresponds to the evolution equation d dt U = log det(U i ) (1.2) for the Kähler potential U of the Kähler metric g i . (For the notations and conventions used in this paper, we refer to Section 2.) In general Ricci and Kähler-Ricci flow may develop singularities well before converging to an Einstein metric, and blow-up analysis gives rise to complete non-compact solitons for these flow equations [4,15]. Thus an understanding of the flow on non-compact manifolds is essential. The general theory for the Ricci flow on non-compact manifolds was established in a series of papers by Shi [18] where in particular, parabolic maximum principles on non-compact manifolds are established. For a survey article concerning the Ricci flow and singularity analysis in particular see [7], a list of solitons for the Kähler-Ricci flow can be found in [11]. Stability questions for the Ricci flow have been  [13,14,21]. Questions of stability and uniqueness in a limit as time tends to infinity of (Kähler-)Ricci flow on compact manifolds have also been studied in many papers, e. g. [16]. All these references represent only a small selection of many articles concerning this subject. Another area of research is to rule out special types of solitons after a blow-up by using geometric restrictions and to prove uniqueness of those solitons.
In this paper, we focus on the non-compact complete gradient Kähler-Ricci solitons found in [5]. They are rotationally symmetric with positive holomorphic bisectional curvature and their existence has been proved by solving an ordinary differential equation. It turns out that these solutions are unique (up to scaling and dilatations) in the class of rotationally symmetric gradient solitons with positive holomorphic bisectional curvature. To learn more about these solitons, it is desirable to know, whether they are stable under appropriate perturbations. In this paper we answer this question in the affirmative. We will show that the gradient solitons in [5] on C n are stable under appropriately decaying perturbations of the Kähler potential.
A further question in this direction is, whether there exist other solitons without rotational symmetry. We wish to remark that both questions also seem to be unsolved for the corresponding problems concerning strictly convex non-compact hypersurfaces that arise as translating solitons for the mean curvature flow.
We wish to give a heuristic argument why we impose the condition that the perturbation should decay at infinity. In our situation, we don't expect to get stronger results than for the standard heat equation on R n . In this case, however, we can take a bounded perturbation of the stationary solution u = 0 that satisfies for t = 0 u((x 1 , x 2 , . . . , x n ), t) → 1 as x 1 → ∞, −1 as x 1 → −∞, (1.3) uniformly in (x 2 , . . . , x n ). It follows directly from the heat kernel representation of a solution that (1.3) remains true during the evolution, i. e. for t > 0. Similarly, we expect that for general bounded perturbations of the potential in our equation, the oscillation of a perturbation will not tend to zero. Of course in such special cases as the one above one can show that the solution flattens out on compact sets. This explains why it is natural to have decay assumptions in our stability theorem. To simplify its formulation, we give the following Definition 1.1. A function u 0 : C n → R is called a C-potential, if it is rotationally symmetric, (u 0 ) i is a Kähler metric with positive holomorphic bisectional curvature, and gives rise to a gradient Kähler-Ricci soliton.
We refer to Theorem 1 in [5], where C-potentials are shown to exist and to be unique up to scaling and holomorphic transformations. In particular, using Notation 2.3, it is shown that C-potentials are characterized by the Equations (5.1) and (5.2).
Theorem 1.2. Let u 0 be a C-potential in complex dimension n ≥ 2 andũ a smooth perturbation such that (1)ũ i defines a complete Kähler metric on C n equivalent to (u 0 ) i with bounded curvature.
Then with u as initial condition, (2.2) has a long time smooth solution converging to 0 as time tends to infinity.
In the special case of a compactly supported perturbation, we get Corollary 1.3. Let u 0 be a C-potential in complex dimension n ≥ 2 andũ a smooth perturbation such thatũ i defines a complete Kähler metric on C n equivalent to (u 0 ) i and u =ũ − u 0 is compactly supported. Then with u as initial condition, (2.2) has a long time smooth solution converging to 0 as time tends to infinity.
A geometric interpretation of the decay condition (1.4) is a follows. If the eigenvalues of the perturbed metric minus the soliton metricũ i −(u 0 ) i with respect to the soliton metric (u 0 ) i decay like (2 log |x|) −2−α and u tends to 0 at infinity, then (1.4) is fulfilled for an appropriate value of K. This is obtained by integrating u radially. Note that the barriers introduced in (5.3) have the same decay in terms of the metric.
To give an overview over the method used here to prove Theorem 1.2, we describe our proof in words. It is convenient to transform our flow equation such that the gradient solitons, we are interested in, become stationary solutions of an equation, namely (2.1). In Section 4, we obtain smooth longtime existence and get uniform estimates for the perturbation of a C-potential. Then we want to apply the maximum principle to deduce that the oscillation of the perturbation is strictly decreasing in time (or zero). Due to the noncompactness of C n , however, we have to make sure that the supremum is attained somewhere. We don't know how to prove this directly. Instead, we enclose our perturbation from above and from below by radially symmetric barriers that decay at infinity (as a function of r = |x|) and correspond to Kähler metrics with positive holomorphic bisectional curvature. During the evolution, the upper barrier stays positive, monotone in r, and rotationally symmetric. This ensures that it attains its maximum at the origin and we can apply the strong maximum principle to deduce that the oscillation, as a function of time t, is strictly decreasing. However, this is not enough to show that the barrier converges to zero. It might happen (and seems to be an interesting question, for which equations it actually happens), that the fact that the perturbation tends to zero at infinity is destroyed during the evolution as t → ∞. This would imply that the perturbation would converge to a positive constant as t → ∞. For the standard heat equation, however, it is quite easy to exclude this phenomenon as the L 2 -norm of a solution is nonincreasing, so it remains finite during the evolution, at least for H 1,2 initial data. The argument extends to any smooth solution. In our situation, we can find a quantity (Lemma 6.3) that is equivalent to the intrinsic L p -norm of the perturbation and is also decaying. This property relies heavily on the construction of special barriers with positive holomorphic bisectional curvature (this is preserved during the evolution [18]). As the total (intrinsic) volume of our soliton is infinite, this excludes the possibility that the perturbation tends to a positive constant. Similar considerations apply to the lower barrier and thus the original perturbation, enclosed in between these two barriers during the evolution, converges to zero. Once C 0 -convergence is established, smooth convergence follows immediately from our a priori estimates.
In Section 2, we introduce the notations that we will use throughout the paper and transform the evolution equation (1.2) to other coordinate systems. We explain in Section 3, why we have shorttime existence of solutions, and prove uniform a priori estimates that guarantee longtime existence in Section 4. We sketch how to construct barriers in Section 5 and refer to the Appendices B and C for details. In the proof of Theorem 1.2, we use the evident Lemma 6.2 and give it's proof in Appendix A. Finally, Section 6 contains the proof, that our barrier converges to zero, the crux of the proof of Theorem 1.2.
The authors wish to thank Jürgen Jost, Shing-Tung Yau, the Alexander von Humboldt foundation (Feodor Lynen Research Fellowship), Harvard University, Cambridge, MA, U.S.A., and the Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany, for their support. We thank Gerhard Huisken for his interest in the paper and the advice to reformulate the decay condition geometrically.

Preliminaries.
Notation 2.1. We use indices to denote partial derivatives, If for a function u : C n → R, the matrix (u i ) is positive definite, we call u a Kähler potential. Then (u i ) is a Kähler metric and we denote its inverse by u i . Lower case Latin indices range from 1 to n. We use the Einstein summation convention with a special convention for Latin capitals, e. g. zu.
We will use z i and (z) to denote standard flat coordinates on C n . Sometimes, it will be appropriate to use standard Euclidean coordinates x i . The Laplace operator with respect to the metric (u i ) is defined by In the estimates that follow, we will use c to denote a fixed positive constant that does not depend on time, but may change its value from line to line. Indices preceded by a comma, e. g. u ,ik , indicate covariant differentiation with respect to the background metric (u 0 ) i introduced in (2.2). As usually, we use R i to denote the Ricci tensor, Rm for the Riemannian curvature tensor, · to denote a (pointwise) norm with respect to the induced metric, and ∇ g to indicate covariant differentiation with respect to the metric g i . We do not use different notations for an the initial value u : C n → R and for the corresponding solution to Kähler-Ricci flow u : Notation 2.2. For a Kähler metric (g i ), we obtain Christoffel symbols as follows Γ i jk = g il ∂g jl ∂z k , Γīk = Γ i jk , and other components are identically zero. Covariant differentiation is defined by Xc ,ab =Xc ,ba , X c,ab =X c,ba − R abcd gd e X e , Xc ,ab =Xc ,ba + R abdc g dē Xē.
In these formulae, the (holomorphic) Riemannian curvature tensor appears, which is defined by Contracting with respect to the metric yields the (holomorphic) Ricci tensor Finally, a Kähler manifold has positive biholomorphic sectional curvature, if R iīj > 0.
Note, that we do not sum here.
Notation 2.3. It will be convenient for the proofs to switch between different evolution equations. If we assume in the following that a function U fulfills (1.2), we will also assume that the function u is obtained from U by applying the transformations leading to (2.2), so it solves (2.2). Similarly, we assume thatũ fulfills (2.1). When we consider a rotationally symmetric solution of (2.1) depending on the variable s = log |z| 2 , we denote this byû. Analogous notations are used for other functions solving (1.2). Given a C-potential u 0 , we will denote by U 0 the corresponding solution to (1.2).

2.2.
Transformations. We will now fix a C-potential u 0 . Here and in the rest of the paper, all evolution equations are defined on C n .
For further considerations, it will be convenient, to change coordinates such that the evolution of a C-potential in time is as simple as possible. This can be obtained by introducing where U is a solution to (1.2). We immediately get d dtũ = 1 2 U I z I e 1 2 t + d dt U +nt and det(ũ i ) = e nt det(U i ). Using (1.2), we see that the evolution equation forũ is given by where we write again z i and z for e 1 2 t z i and e 1 2 t z, respectively, i. e. we evaluateũ at z i , (z) , t . The initial value is clearly unchanged,ũ(x, 0) = U (x, 0). We now show that the C-potentials introduced in [5] are in fact stationary for (2.1). These potentials are characterized in [5] by radially symmetric functionsû(s) in the variable s = log |z| 2 for which the following conditions are satisfied. Forû ′ (s) ≡ ϕ(s) it is required that ϕ(s) → 0 for s → −∞ and ϕ fulfills (when normalized appropriately) the ordinary differential equation By differentiatingû log |z| 2 =ũ as it is seen that a C-potential is a stationary solution to (2.1). To show stability of (2.1) at u 0 , it will be convenient to writeũ = u 0 + u. Asũ and u 0 solve (2.1) we get directly from the definition of u its evolution equation The advantage of this evolution equation is that u 0 is time-independent and it will turn out that it allows to consider functions u(x, t) that are uniformly bounded in time.
2.3. Hölder spaces. We now define the parabolic and elliptic Hölder spaces of a non-compact Kähler manifold (M, g i ). We will use these spaces to apply Schauder estimates in proving a priori estimates for (2.2). These are parabolic versions of the elliptic spaces defined in [8,19,20].
where k 1 > 0 and k 2 < ∞ are independent of i. Also, for all a, b, we have g * ab i,k+α ≤ C where the constant C is independent of a, b, and i and · i,k+α denotes the standard Hölder norm on V i . A manifold is said to have bounded geometry of order infinity, if it has bounded geometry of order k + α for every k ∈ N and 0 < α < 1.
Consider now the following norm defined for any smooth function u on is the standard parabolic Hölder norm on V i × [0, T ) in the standard coordinates on V i . We can now state the following definition. When to time interval [0, T ) is understood from the context, we will simply use C k+α, k ). Also, one can define the elliptic Hölder spaces C k+α in an obvious way.

Short Time Existence
We now establish the following general short time existence result. Lemma 3.1. Let (M, g i ) be a complete non-compact Kähler manifold such that Rm ≤ c 0 and f : M → R is a smooth potential of the Ricci tensor, i. e. R i = −f i . Then for some T > 0 depending only on c 0 , the following initial value problem has a smooth solution u(x, t) for t ∈ (0, T ]. Moreover, for any t ∈ (0, T ], the Kähler metric g i (x)+ u i (x, t) is equivalent to g i and has bounded geometry of order ∞ and f (x) This is the Kähler potential version of the following theorem of Shi [18].
for t ∈ (0, T ]. Moreover, for all t ∈ (0, T ],g i (x, t) is a complete Kähler metric on M equivalent to g i and we have the following estimates for the covariant derivatives of the curvature tensor ofg i (x, t).
Proof of Lemma 3.1. Under the hypothesis of the lemma, Theorem 3.2 guarantees a short time solutiong i to the Kähler-Ricci flow (3.2). Using this solution, we solve the following ordinary differential equation for a smooth function u(x, t). It is then straight forward to verify that we must haveg i (x, t) = g i (x) + u i and thus u(x, t) is a smooth solution to (3.1). The details of this verification can be found in [9]. To complete the proof of Lemma 3.1, we need to show that for any t ∈ (0, T ] the Kähler metricg i (x, t) = g i (x) + u i (x, t) has bounded geometry of order ∞. In [20] the authors prove that on a non-compact Kähler manifold, one has bounded geometry of order 2+α provided one has bounded curvature, scalar curvature, and gradient of scalar curvature. Their proof can in fact be extended to show that one has bounded geometry of infinite order provided one has all covariant derivatives of curvature bounded. Thus sinceg i (x, t) has all covariant derivatives of its curvature bounded, Theorem 3.2, we see thatg i (x, t) in fact has bounded geometry of order infinity. This completes the proof of Lemma 3.1.

A Priori Estimates And Longtime Existence
In this section we prove a priori estimates for solutions of (2.2). We follow the approach first used by Cao [6] which is to adapt the elliptic estimates proved by Yau [22] and Aubin [3] for the elliptic complex Monge-Ampère equation to the parabolic case. As we may transform (2.2) to an equation of the form (3.1), we get short time existence. We may assume that we have Choosing T smaller if necessary, we may also assume that v(·, t) gives rise to a complete Kähler metric uniformly equivalent to (u 0 ) i .
Proof. This follows directly from the maximum principle in [10].
Proof. We estimate Then interpolation gives the claimed inequality. We differentiate (1.2) and obtain d 2 where R V and R U 0 are the scalar curvatures of the metrics V i and (U 0 ) i , respectively. As V and U 0 give rise to solutions of Kähler-Ricci flow, the corresponding scalar curvatures are uniformly bounded [18].
We are not able to prove gradient estimates directly. Instead, we have Proof. We apply L p -estimates [12,Thm. 9.11] to 0 < n + ∆ u 0 v and obtain spatial H 2,p -bounds for v. Then the Sobolev imbedding theorem implies the result. Note that we only used 0 < n + ∆ u 0 v and the C 0 -bound.

Second Order Estimates. Consider the quantity
where ∆v denotes the Laplacian of v with respect to (u 0 ) i (x) and the constant k ≫ 1 is to be chosen later. Clearly a bound on |A| implies a bound on |∆v|. We will bound A from above using the maximum principle. The bound from below will follow directly from some simple inequalities.
Proof. We compute and Note that z i , z, and the Riemannian curvature tensor induced by u 0 are bounded with respect to the metric ((u 0 ) i ), see Remark C.3. This allows to estimate the terms obtained by interchanging the order of covariant differentiation.
We combine (4.2) and (4.3), and get in Ω Here we have used that at a fixed point, we can always choose holomorphic coordinates such that (u 0 ) lk = δ lk and v ik = 0 for i = k, so we get and deduce the second inequality in (4.4). In such coordinates, the terms 1 + v iī are positive for each i and are simply the eigenvalues of the tensor (u 0 ) i + v i with respect to metric (u 0 ) i = δ i . Finally, in passing to the last line of (4.4), we have used Lemma 4.3.
We are now in a position to prove the following Proof. For the proof of the upper bound for ∆v, we will only consider those (x, t) ∈ C n × [0, T ] such that sup with inverseṽ i . We use the geometric-arithmetic means inequality and which is proved easily as we may assume that 1 = λ 1 ≤ λ 2 ≤ . . . ≤ λ n , to obtain We apply Lemma 4.4 and estimate dA dt Fixing k ≫ 1 so large that k −c is bounded below by some positive constant, the maximum principle can be applied to the evolution equation (1 + v iī ) is bounded below by a positive constant. The function v gives rise to a Kähler metric. So all the factors are positive. As we have seen that v iī is uniformly bounded above for each i, the lower bound follows.
Corollary 4.6. There is a constant K 1 > 0 depending only on K 0 and K 2 such that v(·, t) 1 < K 1 for all t ∈ [0, T ].
Moreover, the equivalence factor depends only on K 0 , K 1 and K 2 .
Proof. This follows from the proof of Lemma 4.5.

Higher Order Estimates and Long Time Existence. Consider the quantities
where the covariant differentiation is with respect to (u 0 ) i andṽ i represents the inverse of the time dependent metricṽ i (x, t) = (u 0 ) i (x) + v i (x, t). By the previous section, this norm is equivalent to that using (u 0 ) i . Lemma 4.8. There are constants K 3 , K 4 > 0 depending only on K 0 , K 1 , K 2 such that |Q 3 | C 0 < K 3 and |Q 4 | C 0 < K 4 for all t ∈ [0, T ].
Proof. The above estimates are known in the special case that v is a solution to (3.1) and have appeared in several places in various equivalent forms. We describe some of these briefly. Calabi first estimated |Q 3 | C 0 for the elliptic Monge-Ampère equation on a compact manifold. This estimate was later used by Aubin [3] and by Yau [22] in proving the Calabi conjecture. Calabi's estimate was applied directly by Cao [6] to (3.1) and later by Shi [18] to (1.1). In [18] Shi goes further to estimate an appropriate second derivative of the solution to (1.1) and observes that this is equivalent to estimating the curvature tensor of the evolving metric. An equivalent estimate can be found in [9] where |Q 4 | C 0 is estimated for (3.1).
In our case that v is a solution to (2.2), we point out that it is straight forward to adapt the arguments of the authors cited above to our case.
Notice that while the estimates in Lemma 4.8 follow rather painlessly from the corresponding estimates for (3.1), such is not the case for our laplacian estimate in Lemma 4.5. The difference is that in Lemma 4.8 we have already estimated all derivatives of lower order and second derivatives of the form v i , while in the case of Lemma 4.5 do not have gradient estimates a priori. This does not cause a problem in (3.1) while in our case it does.
Let v be a solution to (2.2) and let C k+α, k 2 + α 2 be the Hölder spaces on C n relative to the metricṽ i . Then for every k, v is bounded in Proof. We prove the respective result for Hölder spaces with respect to the background metric. The corresponding results in these Hölder spaces imply the claimed estimates. Consider an arbitrary coordinate neighborhood V β with coordinates z i as in Definition 2.4. Differentiating (2.2) with respect to z i in this coordinates and rearranging terms gives In what follows, all bounds stated will be independent of β and t. It is readily seen that our estimates from above provide us with a C α, α 2 bound for the coefficients and terms of (4.5). We may then apply standard parabolic Schauder estimates to obtain a C 2+α,1+ α 2 bound for v i (x, t) in an interior domain of V β . A standard bootstrapping argument [6] combined with the fact that the metricṽ i has bounded geometry of order ∞ then allows us to obtain a C k+α, k 2 + α 2 bound on v(x, t) in V β for all k. The lemma now follows readily from Definition 2.4.
uniformly over all t and the curvature ofṽ i (x, t) remains bounded on C n independent of t.
Proof. By our a priori estimates, it is straight forward to see that the curvature of the metricṽ i (x, t) stays uniformly bounded on [0, T ]. Corollary 4.7 implies that the metrics stay uniformly equivalent. Thus to prove the corollary it suffices to prove the assertion of long time existence. Moreover, long time existence for v follows from long time existence for (1.2) with initial conditionṽ(x, 0) = u 0 (x) + v(x, 0). Begin by assuming that T is the maximal time up to which we have a smooth solution. Choosing a time T ′ < T arbitrarily close to T and applying Lemma 3.1 to the metricṽ i (x, T ′ ), we may extendṽ past T as a solution to (1.2) thus arriving at a contradiction and thus proving the corollary.

Barrier Construction
Before we can construct a barrier, we have to determine the precise asymptotic behavior of our soliton. According to [5], we may assume that the function ϕ(s) =û ′ (s) fulfilling where s = log |z| 2 , gives rise to our soliton. The second condition is required to obtain a smooth solution at the origin. We derive in Appendix B the following expansions for ϕ and its derivatives at infinity In the following, we construct barriers in the case n ≥ 2. Now, we assume that our perturbation u(x, 0) of the initial value is such that |u(x, 0)| ≤ K · min 1, s −α , where s = 2 log |x|, 0 < α < 1.
For our barrier we make the ansatz with ϕ as above, that corresponds to the barrier Here ψ is a smooth monotone function such that Assume from now on that R ≥ 1 2 . It is straight-forward to check thatb(s) lies above/below our perturbed initial value.
To prove that for R ≫ 1 fixed sufficiently large is again a technical calculation, we refer to Appendix C.
Note that it is essentially the integrability condition for ϕ b and not (5.4) that determines the possible exponents in the decay condition.
For n = 1, our method does not seem to work. In this case, ϕ(s) is even explicitly known to be log (1 + e s ), but seems to exclude such a barrier construction. For results concerning longtime behavior of solutions to Ricci flow in the corresponding real dimension 2, we refer to [14,21].

Convergence to Zero
In this section we prove Theorem 1.2. We use the radially symmetric decaying barriers constructed in Section 5 to enclose our initial perturbation from above and from below. By smoothly evolving our barriers and perturbed initial value using (2.2) for all time, the maximum principle of [10] implies that our perturbation will converge to zero provided such is true of our barriers. In particular, the perturbed soliton converges back to the original soliton as t → ∞. We will only show that the upper barrier converges back to the original soliton. Studying the behavior during Kähler-Ricci flow is simpler for the barriers as they are rotationally symmetric and decaying in |z|.
Lemma 6.1. Let b be the upper (lower) barrier constructed in Section 5. Then (2.2) with initial condition b has a long time smooth solution, which we also denote by b, which converges to zero as t → ∞ in the C 0 norm.
The proof is divided into several steps. We sketch the proof for the case of the upper barrier and note that the case of the lower barrier is similar. Part of the argument is a modification of the convergence proof in [17]. We first show that the condition that b initially decays monotonely in |z| is preserved for all time, so we get especially b(0, t) ≥ b(x, t) for all (x, t) ∈ C n × [0, ∞). We do this in Lemma 6.2. The strong maximum principle then guarantees that sup C n b is strictly decreasing in t. In fact, we claim that b must converge to a constant. This can be seen as follows. In view of our a priori estimates, we can find for every sequence t n → ∞ a subsequence, again denoted by t n , such that the maps converge locally uniformly in any C k -norm to a smooth function b ∞ (x, t) satisfying the evolution equation (2.2) everywhere in C n × R. Moreover, since the oscillation of b decreases strictly in time by the strong maximum principle [1], it must converge to some nonnegative constant. In other words, the limit solution b ∞ (x, t) has nonnegative oscillation which is constant in time. But it is easy to see that the rotational symmetry and decay condition on b(x, t) also holds for b ∞ (x, t) and thus by the strong maximum principle, the oscillation of b ∞ (x, t) cannot be a positive constant. Thus b ∞ (x, t) is constant in space. The monotonicity of b(0, t) shows that this constant is independent of the chosen subsequence and hence b actually converges to a constant. In Corollary 6.4 we show that during the evolution the L p -norm, for some p ≥ 2, of b is dominated by its value at t = 0. We compute the L p -norm with respect to an evolving volume form which stays uniformly equivalent to the volume form for the initial soliton metric, thus the integral of any positive constant over C n with respect to this volume form for fixed t is infinite. By Remark C.4, the L p -norm is finite for t = 0 provided p > n+1 α with α as in Section 5. So b has to converge uniformly to zero on compact subsets of C n as t → ∞. Note that the monotonicity in |z| is preserved during the evolution and when we extract subsequences. Moreover, b(0, t) is decreasing in t. So b(0, t) has to converge to zero and it follows that the perturbed soliton converges back to the original soliton in C n . It is quite evident that this lemma is true. Thus, we defer it's proof to Appendix A. Lemma 6.3. Let u 0 be a C-potential and b a barrier as constructed in Section 5. Then there exists a metric a i , uniformly equivalent to (u 0 ) i and (u 0 ) i + b i , such that for b evolving according to (2.2), we have C n |b(t)| p det(a i (t)) ≡ I p (t) ≤ I p (0) for p ≥ 2 and t ≥ 0.
Proof. Interpolating between the two determinants in (2.2) and using upper indices to denote inverses, we get Now we define (a i ) to be the inverse of a i . By definition, (a i ) is uniformly equivalent to (u 0 ) i and (u 0 +b) i as these two metrics stay uniformly equivalent during the evolution. For showing the definiteness of terms like d dt (u 0 ) i and d dt ((u 0 ) i +b i ), it will be convenient to substitute so that we almost come back to the original evolution equation (1.2) (1.2). As the metric B i has positive holomorphic bisectional curvature for t = 0 (Appendix C), this is preserved during the evolution [18], so the Ricci curvature also stays positive definite. From (1.1), we obtain that in the sense of matrices. The second inequality follows by noting that (u 0 ) i = (U 0 ) i and thus also corresponds to a solution to the (3.2) with positive holomorphic bisectional curvature. The chain rule and the transformation formula for integrals imply that as det(a i ) = e −nt det(a i ) and the volume elements differ by a factor e nt .
Here (a i ) is the inverse of Note that and we get from (6.1) d dt As (a i ) is obtained by taking the inverse of this matrix, integrating, and taking the inverse once more, the definiteness for the time derivative is inverted twice, so d dt a i ≤ 0. Finally, (a i ) is positive definite, so it follows that It is not obvious, whether I p (t) is differentiable with respect to t or not. Therefore, we define for radii R > 0 In order to compute d dt I p,R (t), we have to compute the evolution equation Using (6.2) and (6.3) To estimate further, we denote by g the real metric corresponding to (a i ), see e. g. [18], and obtain in real coordinates We apply the divergence theorem and use ν to denote the exterior unit normal to B R with respect to the metric g which coincides with x |x| up to a positive factor Here we used suggestive invariant notation. We apply Lemma 6.2 to see that the boundary integral is non-positive and get I p,R (t 1 ) ≥ I p,R (t 2 ) for 0 ≤ t 1 ≤ t 2 . Finally, we let R → ∞ and obtain the claimed inequality.
Corollary 6.4. Let u 0 be a C-soliton and b : C n → R the barrier constructed in Section 5. Assume that p ≥ 2 is chosen such that the L p -norm is finite for t = 0. Then the L p -norm of b stays uniformly bounded when b evolves by Kähler-Ricci flow (2.2) where the constant depends only on the uniform equivalence of the metrics (u 0 ) i and (u 0 ) i + b i that is guaranteed during the evolution.
Proof of Lemma 6.1 and Theorem 1.2. Lemmata 6.2 and 6.4 together with the arguments at the beginning of the section complete the proof of Lemma 6.1 and thus of Theorem 1.2.

Appendix A. Preserving Monotonicity
Proof of Lemma 6.2. It is clear that the rotational symmetry is preserved during the evolution.
If b is not a monotone decaying function of |z| for all t > 0, we choose 0 ≤ T < ∞ maximal such that b is monotone decaying in |z| for t ∈ [0, T ]. Note that b is clearly monotone on a relatively closed subset in time. Our lemma follows if we can show that b stays monotone for a while after T .
To simplify notation, we note, that applying (the independently proven) First, we consider b on C n \ B R (0) for R ≫ 1. The radius R depends only on the fact, that certain coefficients in the ordinary differential equation are not too far from the corresponding values in the asymptotic expansion for the soliton. So the value of R depends only on b(·, 0) and our initial soliton as the initial soliton and the perturbed soliton stay uniformly equivalent during the evolution. We have b(R, T ) > 0 as otherwise the strong maximum principle would imply b(·, T ) ≡ 0, so b(·, t) ≡ 0 for t > T , contradicting the maximality of T . Due to the uniformly bounded geometry during the evolution, there exists T * > T such that b(R, t) Note that both u 0 andb = u 0 + b solve (2.1). We consider u 0 and b as functions of s = log |z| 2 and t and use u 0 and b to indicate that. Equation (2.1) implies that Considering the difference of these two evolution equations gives As u ′ 0 = ϕ in the notation of Section B, we see that b fulfills a parabolic equation of the form where α, β, α −1 ∈ L ∞ ((log R 2 , ∞)) for R ≫ 1 fixed appropriately. As b(R, t) − 1 2 b(R, T ) ≥ 1 c > 0 for T ≤ t ≤ T * , we can extend α, β, and b from [log R 2 , ∞] × [T, T * ] to R × [T, T * ] as in the case with boundary in [2] and apply the result of this paper to see that for h ∈ 0, 1 2 b(R, T ) , #{r ≥ R : b(r, t) = h} = 1 for fixed t ∈ (T, T * ]. This implies monotonicity for r ≥ R. It remains to prove that monotonicity is preserved for b > 1 2 b(R, t). Similarly as above, we can fix a radius R * > R and T * > T such that b(R * , t) < 1 2 b(R, T ) for T ≤ t ≤ T * . Fix ε > 0 and assume that for t 0 ∈ [T, T * ], there exist 0 ≤ r 1 < r 2 < R * such that b(r 2 , t 0 ) ≥ b(r 1 , t 0 ) + ε and t 0 is chosen minimal with this property. b(r, t 0 ) tends to zero as r → ∞. Choose r 3 > r 2 minimal such that b(r 3 , t 0 ) = 1 2 (b(r 1 , t 0 ) + b(r 2 , t 0 )) and r 0 < r 1 maximal such that b(r 0 , t 0 ) = 1 2 (b(r 1 , t 0 ) + b(r 2 , t 0 )) (if such an r 0 exists). Set Ω := B r 3 \ B r 0 if r 0 with this property exists, otherwise Ω := B r 3 . From our assumptions, we get that is strictly smaller than ε for T ≤ t < t 0 and equals ε for t = t 0 . Note that for t close to t 0 , b is close to 1 2 (b(r 1 , t 0 ) + b(r 2 , t 0 )) on ∂Ω. So b| ∂Ω does not "contribute" to the oscillation for t close to t 0 and we get a contradiction to the strong maximum principle as a positive oscillation has to be strictly decreasing in time (Huisken, see e. g. [1]).
As ε was arbitrary, we see that monotonicity is preserved in B R * (0) for T ≤ t ≤ T * , so monotonicity is preserved everywhere for T ≤ t ≤ min{T * , T * } and our lemma follows. and ϕ → 0 for s → −∞ has the following asymptotic behavior at infinity ϕ =ns − log n n s n−1 + (n − 1) log n n s n−1 ns + (n − 1) 1 ns + 1 2 (n − 1) log 2 n n s n−1 n 2 s 2 − (n − 1)(n − 2) log n n s n−1 n 2 s 2 − 1 2 (n − 1)(3n − 5) 1 n 2 s 2 + 1 3 (n − 1) log 3 n n s n−1 (ns) 3 − 1 2 (n − 1)(3n − 5) log 2 n n s n−1 (ns) 3 + (n − 1) n 2 − 6n + 7 log n n s n−1 (ns) 3 + 1 6 (n − 1) 11n 2 − 46n + 47 1 (ns) 3 + o We wish to emphasize that for the application we have in mind, we don't need the high precision of (B.2) explicitly. But as we are not only aiming for the asymptotic expansion for ϕ, but also for ϕ ′ , ϕ ′′ and ϕ ′′′ , we have to compute the expansion for ϕ with high precision, as we have to use (B.1) and derivatives of this equation to determine derivatives of ϕ iteratively. Obviously, derivatives of the expansion of a function do not necessarily have to coincide with expansions of the derivatives. In our situation, however, these two operations commute. This is essentially due to the fact that ϕ satisfies (B.1).
To get (C.1), we study s −1−α R α ψ in detail. When we choose R sufficiently large, |s −1−α R α | becomes arbitrarily small for s ≥ R. For s ≤ R, however, ψ s R vanishes. Thus (C.1) follows for s ≥ R when R is sufficiently large and is true for s < R by the calculations in [5].
Equations (C.2), (C.3), and (C.4) are proved similarly. Note, however, that the term s −1−α R α−2 ψ ′′ is estimated by choosing R large, as s −1−α decays slower as the "leading" term n−1 s 2 as a function of s. This works as ψ ′′ is zero outside R ≤ s ≤ 2R. The same arguments can also be applied to ψ ′′′ . Thus for ϕ b , ϕ ′ b , ϕ ′′ b , and ϕ ′′′ b , the additional terms with a factor K can all be absorbed in the original error terms for R ≫ 1 fixed sufficiently large. We wish to stress, that the sign of ϕ ′′′ b , as s → ∞, is important to get (C.5). For this reason, we had to do all the approximations in Section B up to such a high precision.
Remark C.3. The expression for the Riemannian curvature tensor for a radially symmetric Kähler potential in [5] and the expansions of ϕ and ϕ b at infinity imply Rm ≤ c for the Kähler metrics corresponding to ϕ and ϕ b , respectively. Moreover, the vector fields z i and (z) have finite length with respect to these metrics.
Remark C.4. It follows from Lemma B.1 that the L p -norm and the uniformly equivalent quantity considered in Corollary 6.4 and Lemma 6.3, respectively, are finite for t = 0, if p ≥ 2 is chosen so large that ∞ e (log r) n−1−αp 1 r dr < ∞. (C.6) Choose p such that αp > n+1, with 1 > α > 0 as in (5.3). Introducing a new variable for log r, we see, that the integral in (C.6) is finite as ∞ 1 ρ −2 dρ < ∞.