More on Cotton Flow

Cotton flow tends to evolve a given initial metric on a three manifold to a conformally flat one. Here we expound upon the earlier work on Cotton flow and study the linearized version of it around a generic initial metric by employing a modified form of the DeTurck trick. We show that the flow around the flat space, as a critical point, reduces to an anisotropic generalization of linearized KdV equation with complex dispersion relations one of which is an unstable mode, rendering the flat space unstable under small perturbations. We also show that Einstein spaces and some conformally flat non-Einstein spaces are linearly unstable. We refine the gradient flow formalism and compute the second variation of the entropy and show that generic critical points are extended Cotton solitons. We study some properties of these solutions and find a Topologically Massive soliton that is built from Cotton and Ricci solitons. In the Lorentzian signature, we also show that the pp-wave metrics are both Cotton and Ricci solitons.


I. INTRODUCTION
Gravity has a smoothing effect on a mass or matter distribution, as is evident from all the sufficiently massive spherical objects in the universe, including our planet which has deviations from surface smoothness with its tallest mountain and deepest point in the ocean amounting to less than 10 −3 of its radius. For more massive compact objects, the deviations from surface smoothness would be even less: For example on a neutron star with 1.4 times the mass of the Sun, the height of a mountain could be around 1 cm (not much to hike!). Black holes are even smoother: They have "no hair". Rotations about an axis can introduce larger deviations from perfect spherical shape, but that is a totally different story. For the earth, this deviation is about 4 × 10 −3 , that is the polar radius is around 25 km smaller than the equatorial one. While gravity's smoothing effect on objects that we can see is clear, what we also know is that space-time, even if devoid of any matter , as a manifold itself, also gravitates due to the non-linearity of the field equations. This raises many questions that fall in the intersection of physics and mathematics. For example, one might wonder if one can study the topological classification of manifolds with the help of gravitational equations of some sort. The answer is affirmative but not obvious and actual execution of such a classification is certainly not straightforward. First, one needs to make sure that there is a connection between the Riemannian structure (the metric, the curvature) and the topology. After all the metric and the curvature are local quantities but global topology a priori should not be determined by the local structure. This cursory look is actually not correct and there is a connection between the topology of a manifold and its Riemannian structure as beautifully exemplified by the formula relating the Euler number of a closed two surface in terms of its total Gaussian curvature: 2π , χ(Σ) =´Σ K √ g d 2 x. This and analogous results in higher dimensions encourage one to study the global topology of manifolds with the help of curvature and the metric. Considering the above mentioned smoothing nature of gravity, and the fact that topological classification boils down to identifying manifolds that can be smoothly deformed to each other, one might perhaps try to let the metric on a manifold evolve in time with Einstein's equation or some modified form of it. While this seems like a feasible idea, it is also well-known that singularities, such as black hole singularity, form in a sufficiently general non-linear theory of gravity. Hence one must somehow deal with singularities or better yet turn them into advantage as is done in Ricci flow with surgery.
The choice of gravity-like equations in order to study the topology of manifolds is an important issue. For example, if one naively takes the Einstein's equations (say in vacuum), then depending on the number of dimensions, one has a very different nature of solutions. Considering three manifolds per se, one can split the four dimensional Einstein's equations into 3 + 1 as space and time hence study the evolution of 3 spaces in time. But it turns out that there are all sorts of solutions, including propagating local modes (gravitons), gravitational waves etc. due to the hyperbolic nature of the equations. Sticking to three dimensional manifolds, the goal is to have a finite dimensional solution space (up to perhaps scaling of the metric and the diffeomorphisms). In the choice of the proper equations, uniformization theorem of two dimensional closed, orientable surfaces, which reduces to finding conformally equivalent metrics of constant Gaussian curvature (usually no rmalized as 1, 0, −1 for the sphere, torus and higher genus tori respectively), gives us an important hint: That is, one must find "constant curvature" metrics on a given manifold. Here by "constant-curvature", we mean either constant scalar curvature or constant Ricci curvature or in the extreme case, constant Riemann curvature. Of course the last one is too restrictive and its solutions are already known; they are the maximally-symmetric spaces. On the other hand, finding constant Ricci curvature metrics boils down to solving Einstein's equations in vacuum (R ij = R n g ij ) for which a general strategy is not known even though there are books compiling plethora of solutions. The question whether a manifold of dimension larger than two admits a constant scalar curvature is called the Yamabe problem [1], which in some sense generalizes the two dimensional uniformization problem to higher dimensions.
All the above discussion suggests that given a differentiable manifold M and an initial metric g 0 on it, we should try to get a "nicer", more symmetric metric by smoothly deforming the initial metric. The way to proceed is clear: As in the case of the usual diffusive scalar heat equation, one must have an equation of the form where the deformation parameter t is "time" (or a parameter related to the energy scale) and it is an external parameter not a dimension on the manifold, while x denotes the local coordinates on the manifold which we shall take to have a positive signature metric. Depending on the property of the manifold one is interested in, a specific E ij tensor can be chosen. For example, to define the Ricci flow, which was introduced by Hamilton [2], one takes the following "volume-normalized " equation where r is the average Ricci scalar Fixed points of the flow are Einstein metrics, therefore this is in some sense a parabolic extension of Einstein's equation, albeit its parabolic nature is not apparent at this non-linear level. Namely with this flow, the metric is deformed in such a diffusive way that the manifold eventually admits a constant Ricci curvature while preserving volume. Hamilton [3] used this equation to give another proof of Poincare's uniformization theorem for two dimensional closed surfaces mentioned above. Over a decade ago, Perelman [4,5], using a modified form of the Ricci flow equation, proved the three dimensional uniformization theorem, then called the Thurston's conjecture [6] with a heroic effort bringing Hamilton's expectation and program to a successful end. Even though the uniformization theorem itself and the equations employed to prove it are so closely related to physics, not much work has been done on Ricci flow in physics (save a couple which we mention below). Let us expound upon how the uniformization theorem and Ricci flow type equations are related to physics. First of all, just about the same time (and actually, a little before) Hamilton introduced the Ricci flow, the equation was found by Friedan [7] as the renormalization group (RG) equations of the couplings in a non-linear two dimensional sigma model. In the non-linear sigma model, the metric is just like any other coupling and runs with energy, hence there is a β function for it and the β function is given as the Ricci tensor plus infinitely many two tensors built from the curvature tensor and its derivatives. Two-dimensional sigma model fields are coordinates on the target space which is a Riemannian manifold. RG flow in the sigma model corresponds to a modified Ricci flow in the target space.
One can also see the relevance of the three dimensional uniformization theorem in general relativity (GR): In the 3 + 1 dimensional formulation of GR, three manifolds evolve in time. Since closed 3 manifolds are classified by the uniformization theorem, in essence, the theorem is at the heart of the question of what the shape of the apparent 3 dimensional spatial universe is. Let us give another example of how these flows can be relevant to physics. Consider a gravity theory (a Euclidean one) defined with certain boundary conditions. There could be many solutions satisfying the boundary conditions. In that case, the question is which solution is the global minimum of the Euclidean action. One can use a Ricci flow type equation to study the stability of the solutions and determine possible transitions between solutions. In fact, this was done in [8] for general relativity in a box with a boundary of S 2 × S 1 with 3 critical points in four dimensions. Their exercise is a demonstration of how Euclidean version of gravity (which could possibly correspond to quantum processes such as black hole decay in the Lorentzian setting) can be studied by geometric flows. Namely, the discrete solution space is extended by the Ricci flow hence an off-shell theory can be constructed. See [9] for a similar analysis using the Ricci-Cotton flow that studies the transition between the two vacua Anti de Sitter (AdS) and warped AdS. See [10], for some other possible applications of Ricci flow and [11] for a modified Ricci flow with additional fields and [12] for some higher derivative flows.
Inspired by the work on Ricci flow, in [13], Cotton tensor (C ij ) was used to define a geometric flow exclusive to three dimensions as where fixed points are locally conformally flat metrics. Note that since the Weyl tensor vanishes identically in three dimensions, a manifold is conformally flat if and only if C ij = 0. Therefore, Cotton flow equation could be used to answer the question: whether a given three manifold M admits a conformally flat metric or not. In this sense pure Cotton flow is orthogonal to the Yamabe flow since Yamabe flow keeps the conformal class intact, while the Cotton flow changes the conformal class. In [13], equation (4) was derived as a gradient flow of an entropy functional and the equation was numerically and analytically solved for homogeneous Thurston's geometries. In particular, it was shown that a deformed initial metric on S 3 flows into the conformally flat round metric on S 3 . In the table below, we compile these homogeneous geometries and compute their Cotton tensors just for the sake of completeness. As it is clear from the  [14]. Analytic solutions show that homogeneous metrics in the N il class tend to a pancake degeneracy and the Solv metrics develop cigar degeneracies.

Thurston Geometries Geometry Metric
Curvature Tensors Cotton Tensors In this current work, we shall study several pertinent points in Cotton flow, the most important being the linear (in)stability of flat space. We also show the instability of Einstein spaces and some conformally flat spaces. We shall also linearize the Cotton flow equations about a generic background and apply DeTurck's technique [15] to get a slightly better equation, yet we shall not be able to prove short-time existence for an arbitrary initial metric. Nevertheless, ensuing discussions along this line will need our linearized equations. We also refine the earlier gradient flow formulation of the Cotton flow and define new objects which we shall call extended Cotton solitons that are critical points of the entropy functional. We study two examples of Cotton and Ricci solitons and define Topologically Massive solitons. There has been several other interesting works in the literature regarding the Cotton flow (or its more complicated cousin Ricci-Cotton flow). For example, see [16][17][18] for Cotton-like flow in the Horava-Lifshitz gravity context. See also [19][20][21] where Cotton solitons were found. The lay-out of the paper is as follows: In Sec.II, we linearize the Cotton flow equation about a generic background and identify the symmetries of the equation. In Sec.III, we fix the diffeomorphisms and scalings of the metric via the DeTurck trick to simplify the equation and compute the principal symbol of the relevant operator. In Sec.IV, we study the linearized stability of the flat space which is a critical point of the flow and find an unstable mode. In Sec.V, we refine the gradient flow formulation of the Cotton flow, and expand the entropy functional for a generic fixed point of the flow up to second order in perturbation theory, and show that Einstein spaces as critical points are unstable. In Sec VI, we give details of the properties of extended Cotton solitons and show that there are no non-trivial compact solutions. We give examples of Cotton and Ricci solitons as well as Topologically Massive solitons. In the Appendices A,B,C, we expound upon some of our calculations and in Appendix D, we show how the Cotton tensor flow under arbitarary flows and specifically the Ricci flow.

II. LINEARIZATION OF THE COTTON FLOW ABOUT AN ARBITRARY CURVED BACKGROUND
Cotton flow on a 3 manifold is defined as 1 where we have scaled t to set κ = 1 and from now on we suppress the arguments of the tensors. The Cotton tensor is given as with ǫ ikl being the anti-symmetric tensor density defined as ǫ 123 = +1. [In Appendix A, we compute how some other tensors and scalars built out of tensors flow under Cotton flow.] In what follows we will use the anti-symmetric tensor C ij is symmetric and covariantly conserved ∇ i C ij = 0, and traceless g ij C ij = 0. These constraints tell us that out of 6 components of the (symmetric) metric g ij , only 2 can be determined by the Cotton flow equation (5), hence a gauge fixing must be done. To understand the local nature of the Cotton flow about a flat background, equation (5) was linearized for flat metricsḡ ij = δ ij in [13]. In this section, we will generalize this result for a generic backgroundḡ ij to understand the linearized version of the flow. The results are relevant to the existence and uniqueness of the flow as well as the linear stability of the critical points, namely whether the critical points are saddle points, minima or maxima of the action from which Cotton flow is derived. For the existence and uniqueness of the flow, unfortunately we have nothing to say since the equation is a third order PDE and necessary mathematical technology does not seem to exist yet, but stability issue will be studied in the next section.
To linearize the Cotton tensor about an arbitrary backgroundḡ ij , it is useful to recast it in an explicitly symmetric form as [22] where G j l = R j l − 1 2 δ j l R is the Einstein's tensor. Assuming g ij =ḡ ij + h ij , and h ij is small compared to the background metricḡ ij which at this stage arbitrary (namely not a critical point of the flow), then one obtains the linearized Cotton tensor as whereS ij =R ij − 1 3ḡ ijR is the traceless Ricci tensor and¯ =ḡ ij∇ i∇j and h ≡ḡ ij h ij . Note that "L" refers to the linearization and all the barred quantities are taken with respect to the background metric, which also raises and lowers the indices of linearized quantities. One must be careful with the linearization of the up and down indices, for example Derivation of (9) is somewhat tedious, hence we relegate the details into Appendix B. This is the main equation that should be used in the study of the stability of various critical points of the flow. Let us note several properties of (9). First of all one obtains the contractionḡ ij (C ij ) L = −h ijC ij , as expected from the linearized tracelessness condition (g ij C ij ) L = 0. Secondly, the linearized version of the divergence-free condition of the Cotton tensor reads at this order as which is satisfied by (9) which can be shown after a rather long computation. Let us also understand the symmetries of the linearized Cotton tensor (9). Under infinitesimal diffeomorphisms and scalings of the metric as with λ(x) a scalar function and ζ i a vector field, (9) transforms as where the first three terms come from the Lie derivative L ζ (C ij ) L as expected. The important point is that (13) vanishes at the fixed points (C ij = 0) of the flow, hence diffeomorphism and scalings apparently make the critical point a saddle point. But this is actually a red-herring coming from the symmetry of the theory and so must be moded out or gauged out to study genuine flows. This property of the linearized flow equation can also be reached in the non-linear flow equation.
Diffeomorphism are easy to understand (they are similar to the Ricci flow case), therefore let us concentrate on the scalings of the metric. Consider changing the t-parameter to s → s(t), and define a new metric on the manifold as such that g ij (t) satisfies the Cotton flow equation Then, lets check the flow satisfied byg ij (s): Since the Cotton tensor transforms under conformal scalings asC ij (g) = ϕ − 1 2 C ij (g), (namely √ g C i j is conformally invariant) one has the following flow equation Therefore by parameterizing the time variable one can produce the scaling term. Next having observed the symmetry of the flow, we can introduce the DeTurck trick [15] which was used to show that the Ricci flow is parabolic and hence diffusive. Without this trick, parabolic nature of the Ricci flow for an arbitrary initial metric is a highly cumbersome task to show [2].
Since the Cotton tensor is of third order, one should not expect to apply the theory of elliptic differential operators. But, as we shall see, a modified version of DeTurck trick removes the zero modes of the relevant operator at the critical point and hence modes out the symmetries of the flow, yielding a somewhat simpler equation.
Before closing this section let us briefly note how the Cotton flow equation should be linearized in the first order formalism where instead of the metric one uses the dreibein and the spin connection. For this purpose, let e a and w a b denote the dreibein and spin-connection 1-form respectively. Then the curvature 2-form reads from which the Ricci 1-form and the scalar curvature can be computed via the inner product.
Then the Cotton 2-form reads with the Schouten 1-form given as Y a = (Ric) a − 1 4 Re a . With these definitions, the first order form of the Cotton flow (which was employed in [13]) to discuss the flow of homogeneous quantities) reads where * denotes the Hodge dual. Assumingē a satisfies following the notation of [23], let us expand around the critical point of the flow as e a ≡ē a + ϕ a bē b .
Hereῑ bē a = δ a b . The linearized flow equation reads To compute the linearized Cotton 2-form, let us first find the linearized spin connection whereD a ≡ῑ aD and the background spin connection satisfyingDē a = dē a +w a b ∧ē b = 0. One obtains the curvature 2-form as and the Ricci 1-form as whereD 2 =D aD a , from which the linearized Schouten and Cotton forms can be computed.

III. DETURCK TRICK IN THE COTTON FLOW
By moving the¯ operator to the left in the second term and playing with the relevant indices, we have deliberately arranged (9) in such a way that we can gauge away some terms by extending the Cotton flow as As in the Ricci flow case we can take the symmetry breaking vector as but to simplify our equation further by removing the third term and its partner in (9) we will take it as whereΓ α lλ is a fixed connection on the manifold which we shall take it to be coming from the initial metricḡ ij . There is a subtle issue here which can be seen for the flat (R ij ) critical points: One might mistakingly think that the first piece in (31) would be sufficient to remove all the zero modes but it turns out that it only removes the diffeomorphisms, but not the scalings δ λ h ij = λ(x)δ ij . Linearization of (31) yields Since one has , the third term in (9) and its symmetric partner are canceled with this gauge-fixing, DeTurck-Cotton flow becomes where the reaction terms do not contribute to the principal symbol of the operator on the left-hand side. Of course, even after this simplification one cannot say much about the diffusive nature of the equation, since the principal symbol of the relevant operator vanishes as we show here. Let us calculate the principal symbol of the Cotton tensor as an operator on the metric using the above linearized form. We have with ϕ 2 = g ij ϕ i ϕ j and ϕ i is a non-vanishing vector field. Then the inner product vanishes This class of the operators are known as sub-elliptic operators when the operator is of second order. The classification for third order operators has not been done to the best of our knowledge, so we have not much to conclude from (35). As we noted, the critical points of the flow satisfyC ij = 0, that is, they are the locally conformally flat metrics. Among these critical points, one clearly has Einstein spacesS ij = 0 as a subclass, including the flat space. In 3 dimensions, since the Riemann tensor can be written in terms of the Ricci tensor and scalar curvature, Ricci flat metrics are Riemann flat. In the next section we study the stability of flat metrics as a subclass of critical points of the flow.

IV. STABILITY ANALYSIS OF THE FLAT SPACE
Let us now study the stability of the flat space in Cotton flow. This problem is interesting on its own since flat space is a fixed point of the flow, but it is also relevant to understand the stability of any manifold under small wave-length (or high momentum) perturbations. Here what we mean by small wave-length perturbation is that λ |R ij | ≪ 1, with |R ij | a typical curvature of the manifold in some coordinates. Namely for these sufficiently high energy modes, any space is locally flat. One expects that these modes die out, or diffuse as time passes if the critical point is a stable point under linear perturbations. But it will turn out that there is a potentially unstable mode in flat space at the linearized level.
For the flat background, the linear Cotton flow equation with the DeTurck terms becomes 2 where ∂ 2 is the usual Laplacian in flat space. In a broad sense, this equation is a tensor version of the linearized Korteweg-de Vries (KdV) equation albeit an anisotropic one, in the sense that the flow of h ij with given i and j is determined by h mn with m and n not necessarily equal to i and j. This anisotropy is extremely crucial for the diffusive or non-diffusive nature of the equation since the sign of h ij is not necessarily positive. Therefore, in the stability analysis, the sign of the amplitude of the perturbation is also relevant. So to see the nature of this equation, we must first diagonalize it. Given the initial perturbation at t = 0 as h ij (0, x) with the asymptotic fall off as lim |x|→∞ h ij (0, x) → 0, we can construct "plane-wave" type solutions or Fourier modes from which we can construct the general solution that satisfies the asymptotic boundary condition.
Let us insert the plane-wave ansatz into (36) where ξ ij (p) are (complex) polarizations which only depend on the mode (or momentum p) and w = w(p) is the dispersion relation to be determined below. This insertion leads to a consistency condition on the polarizations where p = |p| and the trace of the polarization tensor is ξ = δ ij ξ ij . We will now play with this equation to find w = w(p) and the constraints on the polarization tensor. [Note that these constraints come from the linearized Bianchi identity and the tracelessness of the Cotton tensor.] Taking the trace of (38), one obtains and contracting (38) with p i p j yields To get the possible dispersion relations, let us assume that ξ = 0 and w = ip 2 then from (39) p i p j ξ ij = 0 is obtained. But (40) gives ξ = 0 so this is a contradiction which can be resolved by assuming w = ip 2 and ξ = 0 and p i p j ξ ij = 0. This necessarily says that ξ ij = ξ T T ij , that is ξ ij is a transverse traceless tensor. Going back to (38), one has This equation can be diagonalized once the left-hand side is plugged to the ξ T T ml in the right-hand side and contraction of η's are carried out, to get 3 Since we cannot take ξ T T ij to be zero (otherwise h ij ≡ 0), we have which contradicts with our earlier assumption w = ip 2 . Hence we conclude that the assumed dispersion relation is not correct: w = ip 2 . Clearly the same reasoning works for the other possible dispersion relations w = −2ip 2 coming from (40). Namely, two apparently possible dispersion relation in (39) and (40) are do not yield possible solutions of (38). Hence we must assume w = ip 2 and w = −2ip 2 . Then, assuming ξ = 0 and p i p j ξ ij = 0 equations (39) and (40) yield which leads to another possible dispersion relation and the identity ξ = p i p j ξ ij . Now, let us evaluate whether this mode is viable or not. For this purpose decompose ξ ij as follows where A(p) is arbitrary at this stage. Using ξ = p i p j ξ ij in (47) gives ξ = A(p)p 2 , and so one arrives at Substituting (48) into (38) yields 3 Alternatively, for T T -modes, by taking one more derivative of (36), one obtains whose Fourier transform yields (44).
Inserting w = −ip 2 , leads to the conclusion that for these modes, the T T polarizations vanish: ξ T T ij = 0. Therefore this is a valid dispersion relation for the following modes For these modes h ij (t, x) = ξ ij e ip.x−p 2 t , hence as time passes, these modes decay namely they are diffusive modes.
Let us now consider the final possibility, that is ξ T T ij = 0, then ξ = 0, p i p j ξ ij = 0, hence we have the complex cubic dispersion relation With this dispersion relation, we still need to make sure that there is a non-trivial solution of equation (41). We can recast that equation for both signs of the dispersion relation as a matrix equation with ξ = (ξ T T ij ) as a 3× 3 symmetric matrix and B as a 3× 3 anti-symmetric matrix whose elements are given as Equation (52)

Equation (52) yields
Normalizing ξ 11 = 1, ξ 12 = ∓i depending on the sign of the dispersion w = ± i 2 p 2 |p|. Therefore we have non-trivial solutions for both signs of the dispersion relation.
Since the equation (36) is linear, then the most general solution reads as We can recast this in terms of the initial perturbation as where This integral can be evaluated for the minus sign, in terms of hyper-geometric and Bessel-kind functions. On the other hand, the positive sign mode is problematic. It is clear that the negative sign modes decay in time and hence these perturbations are diffused. As for the positive sign modes, there is a linear instability. But bear in mind that this analysis has been at the perturbative level and the perturbation theory breaks down for these growing modes hence a strict conclusion cannot be made with regard to the non-linear instability of flat space in Cotton flow.
In fact, even though we have not been able to prove stability at a non-linear level, we conjecture that these growing modes could be tamed at the non-linear level. If this turns out not to be the case, then, since the above analysis is also valid for any manifold for high momentum fluctuations that see the manifold locally flat, any critical point of the flow would be unstable, leading to the conclusion that the theory has no minima but just saddle points. This would be highly unexpected.

V. COTTON FLOW AS A GRADIENT FLOW AND THE ENTROPY FUNCTIONAL
In [13] a gradient flow formulation of Cotton flow was given which we refine it here to better understand the role of symmetries and the Cotton solitons and also study the issue of stability using the entropy functional. The Chern-Simons action [25] on a manifold without a boundary has the first variation under an arbitrary change of the metric as Observe that for δg ij = ∇ i V j + ∇ j V i + λ(x)g ij , the first variation vanishes, δF = 0, owing to the fact that Cotton tensor is traceless and covariantly conserved. Hence such a variation gives us the symmetries of the theory: which are diffeomorphisms and arbitrary scalings of the metric. To get a steepest descend, we can choose the following flow which leads to To get a gradient flow we must require dF dt ≥ 0 where the bound is saturated for the critical points that satisfy the following equation which we shall call it as the extended Cotton solitons. [In [13], the case of constant λ metrics were called Cotton solitons]. It is clear that (61) vanishes at the critical points, since by definition δg ij = 0 there. Requiring the volume density to be a conserved quantity of the flow ∂ t √ g = 0, yields, from (60), the constraint for all metrics g ij that are solutions of the flow.
To understand the stability of the fixed points, let us expand the Entropy functional of the Cotton flow up to the second order in δg ij for a fixed time t. Then from (59), one obtains whereḡ ij refers to the critical points that satisfy (62). Making use of (9) one arrives at where we have dropped some boundary terms and made use of the Lichnerowicz operator action on a symmetric two tensor as which yields in three dimensions To ensure the stability of the critical points one must have For a generic conformally flat fixed point, or for Cotton soliton backgrounds, even though it is possible that δ 2 F is positive, we have not been able to show this. Instability of Einstein spaces S ij = 0, can be seen quite easily if one decomposes the perturbation as which after integration by parts reduces the action (65) to whereR is constant. The first derivative on the perturbation leads to a linear instability. This is consistent with our earlier flat space analysis. Hence Einstein spaces, including the flat space, critical points of the flow are linearly unstable. This is almost evident but let us see it more explicitly for the S 3 metric in the following coordinates Assuming second variation of the entropy about S 3 for this perturbation yields which is not positive in general. Thus, S 3 is a saddle point not a minimum of the flow. In [13], it was shown that as a critical point S 3 is stable under homogeneous deformations.

VI. EXTENDED COTTON SOLITONS
In the previous section, we have shown that generic fixed points of Cotton flow are extended Cotton solitons that satisfy where λ(x) is not necessarily constant as was assumed in [13]. As was shown in [19], if λ(x) is constant and the manifold is compact without a boundary, then there are no non-trivial Cotton soliton solutions, namely λ = 0, C ij = 0, and V i is a Killing vector For extended Cotton solitons, let us prove a similar theorem.
Theorem: On a compact Riemannian manifold without boundary, all solutions of (73) have L V g = −λg, namely, V is a conformal Killing vector and C ij = 0: There are no non-trivial extended Cotton solitons.
Proof: The square of the Lie derivative of the metric reads Trace and divergence of (73) can be computed as Making use of these in (75), one arrives at Integrating this over the manifold, once the boundary terms are dropped, yieldŝ which can be written asˆd where we made use of (76). Since we are in a Riemannian manifold with a positive metric, from (80) it follows that V i is a conformal Killing vector with the conformal factor λ(x). Then from (73) it follows that C ij = 0. So on a compact Riemannian manifold without a boundary, there are no non-trivial Cotton solitons besides the conformally flat ones. This proves the theorem. Of course for the Lorentzian signature and the non-compact manifolds, the theorem does not follow. Let us study some properties of extended Cotton solitons. Taking the divergence of (77) and making use of (76) one arrives at which can be recast as is a "conserved current". For some manifolds one can define the following conserved total charge. Suppose two space Σ foliate the three manifold M and letn i be the normal to the surface Σ. And let γ ij denotes the induced metric on the Σ. Then is a conserved total charge of the manifold, if the following condition is satisfied. Say we separate one of the coordinates as "r" and Σ as the constant r surface. Then for large r, one has for i = r. This leads to While, at the moment, we do know explicit examples of extended Cotton solitons, (one could perhaps find these using the method of [26] since these are less restrictive solutions compared to the Cotton solitons with constant λ). In fact let us give an example in the Riemannian setting with a constant λ [19]. For constant λ, there are many in the Lorentzian setting [19][20][21]. Consider the following metric with ℓ > 0. This metric has the following properties Also, the components of the Cotton tensor can be computed as This metric is a Cotton soliton for the scale factor λ = − 16 ℓ 3 with the components of the vector field given as where c i 's are constants. For this metric, one can check that which is the λ = constant version of (82).
It is quite interesting that (86) is also a Ricci soliton (R ij + ∇ i X j + ∇ j X i + Λg ij = 0 ) with the scale factor Λ = 6 ℓ 2 and the components of the vector field Using this fact we can define a new object that we shall call Topologically Massive solitons that solve which is satisfied by (86) with the following choices Here µ is the topological mass parameter.

P P -Wave Solution of Cotton Flow
In the Lorentzian setting, we can show that the pp-wave metric, (see [27] for its properties), where u and v are null coordinates. The components of the metric can be written as g µν = η µν + 2U λ µ λ ν with the λ µ vector satisfying the following properties: With these definitions, the Ricci tensor and the Cotton tensor can be computed as The pp-wave metric (94) satisfies the Cotton soliton equation with where c i 's are arbitrary constants and f i (u)'s are arbitrary functions. Interestingly, this Cotton soliton is also a gradient Cotton soliton [19] namely One can also show that the pp-wave metric is also a Ricci soliton with the following vector field X i and the metric function U given as Finally, one can show that (94) is a gradient Ricci soliton with

VII. CONCLUSIONS
We have studied the linearized Cotton flow equation about a generic background with the help of a modified form of the DeTurck trick that removes the zero modes of the relevant operator in the flow equation. Then we specifically studied the perturbations about the critical points of the flow. We found that flat space and Einstein spaces, as critical points of the flow, have unstable modes at the perturbative level, making these spaces saddle points rather than minima. We have also shown that (in Appendix C) certain conformally flat fixed points that are not Einstein metrics, also have unstable modes with a dispersion relation cubic in the Fourier momentum. We have also supported our arguments by computing the second variation of the entropy functional about the critical points and show that the second variation is negative for Einstein metrics making them saddle points. Finding the true minima (if there is any) of the Cotton flow equations is an open problem.
We have refined the gradient flow formulation of Cotton flow and gave a definition of extended Cotton solitons and worked out some properties of these critical points. We have also given an example of Topologically Massive solitons that is constructed from a solution to both Ricci and Cotton soliton equations. Our work could be relevant to both mathematics, as discussed in the Introduction and physics. For example it would be interesting to see the connection between the Cotton flow and various holographic theories that have the pure Chern-Simons action in the bulk of a three manifold [28].

VIII. APPENDIX A: EVOLUTION OF SOME QUANTITIES TENSOR UNDER THE COTTON FLOW
Under the Cotton flow, let us note how some tensors and scalars evolve, 1. The Cotton tensor itself evolves as 2. Square of the Cotton tensor evolves as 3. The Riemann tensor evolves as 4. Square of the Riemann tensor evolves as 5. Time independence of Bianchi identity ∂ t (∇ i C ij ) = 0 leads to Note that at the critical points, the flow stops for all the above quantities.

IX. APPENDIX B: LINEARIZATION OF THE COTTON TENSOR AROUND A GENERIC BACKGROUND
Let us give some details of the calculations leading to equation (9). By perturbing the metric about an arbitrary background as g ij ≡ḡ ij + h ij , one obtains From now on we will drop the O(h 2 ) terms. The background metric is compatible with the background covariant derivative∇ kḡij = 0, hence the linearized Christoffel connection reads Our task is to rewrite (106) in such a way that one can see the terms that can be related to diffeomorphisms of the background metric. For this purpose, let us first evaluate the A term in (106). Note that the linearization of the Einstein tensor yields Where (G ij ) L = (R ij ) L − 1 2ḡ ij R L − 1 2 h ijR . By using the explicit forms of the Linearized Ricci tensor and Ricci scalar one can write the A part as (110) Let us now calculate each term of (110) separately: the C term can be written as where we have used the three dimensional identitȳ