Holomorphic differentials, thermostats and Anosov flows

We introduce a new family of thermostat flows on the unit tangent bundle of an oriented Riemannian $2$-manifold. Suitably reparametrised, these flows include the geodesic flow of metrics of negative Gauss curvature and the geodesic flow induced by the Hilbert metric on the quotient surface of divisible convex sets. We show that the family of flows can be parametrised in terms of certain weighted holomorphic differentials and investigate their properties. In particular, we prove that they admit a dominated splitting and we identify special cases in which the flows are Anosov. In the latter case, we study when they admit an invariant measure in the Lebesgue class and the regularity of the weak foliations.


Introduction
We introduce a new family of flows on the unit tangent bundle SM of a closed oriented Riemannian 2-manifold (M, g) of negative Euler characteristic. The flows are (generalised) thermostat flows and are generated by C ∞ vector fields of the form F := X + (a − V θ)V , where X, V denote the geodesic and vertical vector fields on SM , θ is a 1-form on Mthought of as a real-valued function on SM -and a represents a differential A of degree m 2 on M . The triple (g, A, θ) determining the flow is subject to the equations (1.1) K g = −1 + δ g θ + (m − 1)|A| 2 g and ∂A = where i = √ −1 and where K g denotes the Gauss-curvature, δ g the co-differential and g the Hodge-star with respect to g and the orientation. The case m = 3 of these equations appeared previously in [31 31] (assuming θ is closed), where it is related to certain torsionfree connections on T M which admit an interpretation as Lagrangian minimal surfaces. Here we prove that our flows admit a dominated splitting and moreover, that this family of flows admits a parametrisation in terms of holomorphic data. Indeed, we show that a triple (g, A, θ) satisfying the equations (1.1 1.1) determines a holomorphic line bundle structure on the smooth complex line bundle L m := Λ 2 (T M ) (m−1)/2 ⊗ C, so that the "weighted differential" P = (det g) −(m−1)/4 ⊗ A is a holomorphic section of L m ⊗ K m M and such that a certain negative curvature condition holds. Here K M denotes the canonical bundle of (M, g). Conversely, given a closed hyperbolic Riemann surface (M, [g]), a holomorphic line bundle structure on L m and a holomorphic section P of L m ⊗ K m M satisfying a certain negative curvature condition, we construct a triple (g, A, θ) solving (1.1 1.1) and hence one of our flows, by using the uniformisation theorem and by solving an algebraic equation only.
In [41 41], Wojtkowski introduced W-flows by suitably reparametrising the geodesics of a Weyl connection (or conformal connection). We show that the case where A vanishes identically corresponds to W-flows associated to conformal connections on the tangent bundle of a surface that have negative definite symmetrised Ricci curvature. In particular, we recover [41 41,Theorem 5.2], by showing that the flow associated to a triple (g, 0, θ) solving (1.1 1.1) is Anosov. This is achieved by providing sufficiency conditions for a general thermostat flow to admit a dominated splitting and to have the Anosov property, see Proposition 3.5 3.5 and Theorem 3.7 3.7.
We then turn to the case where θ vanishes identically, so that A is holomorphic, hence we have (1.2) K g = −1 + (m − 1)|A| 2 g and ∂A = 0.
Note that applying standard quasi-linear elliptic PDE techniques we obtain a unique solution g to (1. 2) admit an interpretation as coupled vortex equations, see in particular [ 10 10, §5]. The case m = 2 was considered in [33 33] in the context of Anosov thermostats admitting smooth weak bundles (see Section 6 6 for more details). In the case m = 3, the first equation is known as Wang's equation in the affine sphere literature. In [38 38], Wang related its solution to a complete hyperbolic affine 2-sphere in R 3 , in particular g is known as the Blaschke metric. Moreover, for m = 3, a pair (g, A) on M solving (1.2 1.2) defines a properly convex projective structure on M and hence turns M into a properly convex projective surface, see [25 25] and [29 29]. The universal cover Ω of a properly convex projective surface of negative Euler characteristic is a strictly and properly convex domain in the projective plane RP 2 which admits a cocompact action by a group Γ of projective transformations. Consequently, we obtain a (two-dimensional) divisible convex set.
Since Ω is convex, it is equipped with the Hilbert metric and moreover, the Hilbert metric descends to define a Finsler metric on the quotient surface M Ω/Γ, see in particular [21 21] for a nice survey of these ideas. We observe that the geodesic flow of the Finsler metric is a C 1 reparametrisation of the flow we associate to the pair (g, A). Benoist has shown [3 3] that if (Ω, Γ) is a divisible convex set (not necessarily two-dimensional), then the geodesic flow of the Finsler metric F induced on Ω/Γ -henceforth just called the Hilbert geodesic flow -is Anosov if and only if Ω is strictly convex. Since the Anosov property is invariant under reparametrisation, we may ask if the thermostat flow associated to a pair (g, A) solving (1.2 1.2) is Anosov for all m 2. This is indeed the case, we obtain: Theorem 5.1. Let (g, A) be a pair satisfying the coupled vortex equations∂A = 0 and K g = −1 + (m − 1)|A| 2 g . Then the associated thermostat flow is Anosov.
The hyperbolicity properties of thermostats satisfying (1.1 1.1) are not apparent. To expose them, we first conjugate the derivative cocycle to another one in which we can see the effect of equations (1.1 1.1). This conjugation requires a careful choice of gauge, but once that is established, standard methods using quadratic forms give rise to a dominated splitting. To upgrade this dominated splitting to hyperbolic as in the case of Theorem 5.1 5.1 requires an additional ingredient in the form of Lemma 5.2 5.2 below which asserts that K g < 0; this gives control on the potentially problematic size of A.
In the same way as geodesic flows are paradigms of conservative systems, thermostats may be seen as paradigms of dissipative systems. The special case of Gaussian thermostats (a = 0) has provided interesting models in nonequilibrium statistical mechanics [11 11, 12 12, 35 35]. The next theorem shows that Anosov thermostat flows determined by the coupled vortex equations are indeed dissipative except when A = 0.
Theorem 5.5. Let (g, A) be a pair satisfying the coupled vortex equations∂A = 0 and K g = −1 + (m − 1)|A| 2 g . Then the associated thermostat flow preserves an absolutely continuous measure if and only if A vanishes identically.
We remark that due to a theorem of Ghys [13 13] Anosov thermostat flows are Hölder orbit equivalent to the geodesic flow of (any) negatively curved metric of M and hence transitive (to be precise, [13 13] establishes a topological equivalence and the Hölder orbit equivalence follows from [20 20,Theorem 19.1.5]).
In [3 3], Benoist also observes that the regularity of the weak foliations of the Hilbert geodesic flow coincides with the regularity of the boundary of the divisible convex set (Ω, Γ). By a result of Benzécri [5 5], the boundary has regularity C 2 if and only if Ω is an ellipsoid, in which case the induced Finsler metric is Riemannian and hyperbolic. Hence one might speculate that if a solution to the coupled vortex equations (1.2 1.2) gives rise to an Anosov flow having a weak foliation of regularity C 2 , then A vanishes identically. While we cannot prove this in general, we use Theorem 5.5 5.5 to resolve the odd case: Theorem 7.1. Suppose an Anosov thermostat given by the coupled vortex equations has a weak foliation of class C 2 and m is odd. Then A vanishes identically.
The orbits of our flow -when projected to the surface M -define what is known as a path geometry on M , that is, a prescription of a path on M for every direction in each tangent space. In the case where A vanishes identically the paths are the geodesics of a hyperbolic metric and in the case where m = 3 the paths are the geodesics of a properly convex projective structure. In both cases, the path geometry is flat, by which we mean it is locally equivalent to the path geometry of great circles on the 2-sphere. In the final section of the article we show: Theorem 8.3. Let (g, A) be a pair satisfying the coupled vortex equations∂A = 0 and K g = −1 + (m − 1)|A| 2 g . Then the path geometry defined by the thermostat associated to (g, A) is flat if and only if m = 3 or A vanishes identically.
Holomorphic differentials appear naturally in higher Teichmüller theory and here we briefly provide some context for our results while referring the reader to the recent survey [39 39] [26 26] was recently able to construct a canonical parametrisation of H(M, G) in the case where G is PSL(3, R), PSp(4, R) or the split form G 2,0 of the exceptional group G 2 (see also [25 25] and [29 29] for the case G = PSL(3, R)). More precisely, Labourie obtains a mapping class group equivariant identification of H(M, G) with the fibre bundle over Teichmüller space whose fibre at J is H 0 (M,  [25 25, 26 26] and Loftin [29 29] yields the aforementioned description of properly convex projective structures in terms of pairs ([g], A) with A a holomorphic cubic differential.
Acknowledgements. The authors are grateful to Nigel Hitchin, Rafael Potrie and Andy Sanders for helpful conversations and the anonymous referee for her/his careful reading and many useful suggestions. GPP was partially funded by EPSRC grant EP/M023842/1 and TM was partially funded by the priority programme SPP 2026 "Geometry at Infinity" of DFG.

Preliminaries on general thermostats
Let M be a closed oriented surface equipped with a Riemannian metric g, SM its unit circle bundle and π : SM → M the canonical projection. The latter is in fact a principal SO(2)-bundle and we let V be the infinitesimal generator of the action of SO (2).
Given a unit vector v ∈ T x M , we will denote by Jv the unique unit vector orthogonal to v such that {v, Jv} is an oriented basis of T x M . There are two semibasic 1-forms ω 1 and ω 2 on SM , which are defined by the formulas: The form ω 1 is the canonical contact form of SM whose Reeb vector field is the geodesic vector field X.
A basic theorem in 2-dimensional Riemannian geometry asserts that there exists a unique 1-form ψ on SM -the Levi-Civita connection form of g -such that ψ(V ) = 1 and where K g denotes the Gaussian curvature of g. In fact, the form ψ is given by where Z : (−ε, ε) → SM is any curve with Z(0) = (x, v),Ż(0) = ξ and DZ dt is the covariant derivative of Z along the curve π • Z.
For later use it is convenient to introduce the vector field H uniquely defined by the conditions ω 2 (H) = 1 and ω 1 (H) = ψ(H) = 0. The vector fields X, H, V are dual to ω 1 , ω 2 , ψ and as a consequence of (2. 3) also imply that the vector fields X, H and V preserve the volume form ω 1 ∧ dω 1 and hence the Liouville measure. Note that the flow of H is given by t is the geodesic flow of g.
Let λ be an arbitrary smooth function on SM . For several of the results that we will describe below, we will not need λ to be a special polynomial in the velocities. We consider a (generalised) thermostat flow on (M, g), that is, a flow φ defined by It is easy to check that F := X + λV is the generating vector field of φ. Now let Θ := −ω 1 ∧ dω 1 = ω 1 ∧ ω 2 ∧ ψ. This volume form generates the Liouville measure dµ of SM .
Lemma 2.1. We have: where Ω a is the area form of M , we see that 2.1. Jacobi equations. It is easy to derive the ODEs governing the behaviour of dφ t using the bracket relations above. Given ξ ∈ T (x,v) SM (the initial conditions), if we write 2.2. Quotient cocycle. We consider the rank two quotient vector bundle E = T SM/RF . We use the notation [ξ] with ξ ∈ T SM for the elements of E. Note that dφ t descends to the quotient to define a mapping satisfying ρ t • ρ s = ρ t+s for all t, s ∈ R. The basis of vector fields (F, H, V ) on SM defines a vector bundle isomorphism T SM SM × R 3 and consequently an identification E SM × R 2 . Therefore, for each t ∈ R, we obtain a unique map Ψ t : SM → GL(2, R) defined by the rule for all (x, v) ∈ SM and t, s ∈ R, and hence defines an GL(2, R)-valued cocycle on SM with respect to the R-action defined by φ. Explicitly, Ψ t is the matrix whose action on R 2 is given by Note that for thermostats the 2-plane bundle spanned by H and V is in general not invariant under dφ t .
2.3. Infinitesimal generators and conjugate cocycles. Given a cocycle Ψ t : SM ×R → GL(2, R) we define its infinitesimal generator B : SM → gl(2, R) as The cocycle Ψ t can be recovered from B as the unique solution to For the case of thermostats, it is immediate to check that Given a smooth map P : SM → GL(2, R) (a gauge) we can define a new cocycle by conjugation asΨ . It is easy to check that the infinitesimal generatorB ofΨ t is related to B by (2.12)B = P −1 BP + P −1 F P.

Dominated splitting and hyperbolicity for thermostats
We are interested in the questions: when is this cocycle hyperbolic? When does it have a dominated splitting? We start with some definitions.
, and constants C > 0 and 0 < ζ < 1 < η such that for all t > 0 we have We also say: The cocycle Ψ t is said to have a dominated splitting if there is a continuous ρ-invariant splitting E = E u ⊕ E s , and constants C > 0 and 0 < τ < 1 such that for all Obviously hyperbolicity implies dominated splitting. It also implies that there are no conjugate points [9 9]. Moreover the cocycle Ψ t is hyperbolic if and only if the thermostat flow φ is Anosov (cf. for instance [40 40,Proposition 5.1] where it is proved that the subbundles E s,u of E lift to subbundles of T SM to give the usual definition of Anosov flow). We shall say that φ has a dominated splitting if Ψ t has a dominated splitting (this is the adequate notion of dominated splittings for flows, see e.g. [1 1, Definition 1]). For the case of flows on 3manifolds, as it is our case, the existence of a dominated splitting can produce hyperbolicity if one has additional information on the closed orbits. Indeed [1 1, Theorem B] implies that if all closed orbits of φ are hyperbolic saddles, then SM = Λ ∪ T where Λ is a hyperbolic invariant set and T consists of finitely many normally hyperbolic irrational tori.
A very convenient way to establish the aforementioned properties for cocycles is to use quadratic forms as in [28 28, 41 41, 42 42] In what follows it will be helpful to understand how the spaces E u,s are constructed using Q. This is explained in detail in [42 42,Proposition 4.1], so here we just give a brief summary adapted to our situation. We let L + (x, v) denote the set of all 1-dimensional subspaces W such that Q (x,v) is positive on W . The condition on the quadratic form Q ensures that Ψ t acts as a contraction on L + and hence there is a unique point of intersection All our quadratic forms Q below will have the property that Q(0, b) = 0 (using the identification E SM × R 2 ) and hence we can construct E u (and E s ) simply by applying the procedure (3.1 3.1) to the vertical subspace R(0, 1), that is, Let us put these ideas to use. Define K = κ + F V λ.
Proof. We let (a, b) denote the standard coordinates on R 2 . Using the identification E SM × R 2 we define a quadratic form on E by the rule Then where y is the unique solution ofÿ for (a, b) = 0 and such that Q (x,v) (a, b) = 0. Then Proposition 3.4 3.4 immediately implies that Ψ t has a dominated splitting. We can upgrade that to hyperbolic as follows. If we let z :=ẏ − V (λ)y, then the quadratic form is just zy. By the construction of the subspaces E s,u (cf. (3.1 3.1)) we see that E s,u do not contain neither z = 0, nor y = 0. Hence there exist continuous functions r s,u : This gives exponential growth for z and hence the desired exponential growth for Ψ t on E u . Arguing in a similar way with E s , we deduce that Ψ t is hyperbolic.
In fact we can generalise this further as follows.
Theorem 3.7. Let p : SM → R be a smooth function such that Then φ has a dominated splitting. If in addition κ p + (V λ) 2 4 < 0, then the flow is Anosov.
Remark 3.8. Let us see the main issue with upgrading the last theorem to "hyperbolic" as in the proof of Proposition 3.5 3.5. Certainly we get continuous (Hölder in fact) functions r s,u .
To be definite consider the case of E u and initial conditions (y(0),ẏ(0)) ∈ E u . Thenẏ = r u y and z = (r u −p)y with r u −p > 0 as before.
To get exponential growth we either need: and it is not clear how to get any of these conditions in this generality. In the special cases above p = 0 or p = V λ, we do get one of these conditions. In all these cases the function r = r u,s satisfies the Riccati equation which is easily derived using the invariance of E s,u and the Jacobi equationÿ −V (λ)ẏ +κy = 0. Observe that h := r − p satisfies the Riccati equation Using (3.2 3.2) we can also give a construction of functions r u,s at the level of the Riccati equation as follows. Fix (x, v) and consider for each R > 0, the unique solution u R to the Riccati equation . These limiting solutions exist whenever the cocycle Ψ t has no conjugate points [2 2]. It is easy to check that in all the cases we consider below, the cocycle Ψ t is free of conjugate points.
Remark 3.9. This remark attempts to clarify the role of the function p in terms of conjugate cocycles and infinitesimal generators as in Subsection 2.3 2.3. As we have already pointed out, the infinitesimal generator B for a thermostat is given by Consider a gauge transformation P : SM → GL(2, R) given by P = 1 0 p 1 .
A calculation using (2.12 2.12) shows that the conjugate cocyleΨ t via P has infinitesimal generator given byB The cocycles Ψ t andΨ t share the same dominated splitting/hyperbolicity properties by virtue of being conjugate, but the form ofB exposes clearly the origins of these properties via κ p < 0 (cf. [42 42,Introduction]). The trace of both matrices, which is −V λ (minus divergence of F ), indicates the dissipative nature of thermostats.

Applications
We consider now some special choices of λ. To this end let θ be a 1-form on M which we may equivalently think of as a function θ : SM → R satisfying V V θ = −θ. For later use we record that the co-differential of θ and its Hodge-star satisfy Moreover, let A be a differential of degree m on M with m 2. By this we mean a section of the m-th tensorial power of the canonical bundle K M of (M, g). Likewise, we may equivalently think of a differential A of degree m on M as a real-valued function a : SM → R satisfying V V a = −m 2 a, explicitly, we obtain The thermostat flows we investigate are of the form λ = a − V θ. We will see next that they admit a dominated splitting provided a natural pair of equations is satisfied by the triple (g, A, θ). In order to derive these equations we first need a Lemma.
Remark 4.2. Note that applying V we see that (4.4 4.4) is equivalent to Proof of Lemma 4.1 4.1. We use the complex notationã = V a/m + ia and ω = ω 1 + iω 2 . Since V V a = −m 2 a, we compute that there exist unique complex-valued functionsã andã so that dã =ã ω +ã ω + imãψ. In particular, we have π * (∂A) =ã ω ⊗ ω m . Since We also have π * (θ − i g θ) = (θ + iV θ) ω. Hence (4.   . We could also consider the case λ = a − V θ with a representing a differential of degree m = 1, that is, a (1,0)-form. We exclude this case since it corresponds to the case where A vanishes identically by defining θ = V a and considering λ = −V (θ − θ) = λ. Flows defined by λ = −V θ = were studied previously under the name W -flows as they arise naturally by reparametrising the geodesics of a Weyl connection, see [41 41]. In particular in [41 41,Theorem 5.2] it is proved that W -flows are Anosov provided K g −δ g θ < 0. A simple computation gives that K = K g −δ g θ hence we recover [41 41, Theorem 5.2] by applying Proposition 3.5 3.5. In particular, we see that if A is a holomorphic 1-form and g satisfies K g < 0, then the associated thermostat flow is Anosov.
We now want to apply Theorem 3.7 3.7 to the case λ = a − V θ for some good choice of p.
Lemma 4.5. Suppose λ = a − V θ and take p = V a/m + θ. Then κ p ≡ −1 if and only if the following two equations are identically satisfied Proof. Taking p = V a/m + θ gives where we have used (4.   We also observe: Proposition 4.7. Consider a pair (g, A) with A holomorphic and K g < 0. Then the associated thermostat flow has a dominated splitting. Moreover, for m = 2, the flow is Anosov.
Proof. The fact that there is a dominated splitting follows from κ p < 0. For m = 2 we note that  Before we prove Lemma 4.8 4.8, we first recall some basic facts about conformal connections. Let us fix a Riemannian metric g ∈ [g]. It follows from Koszul's identity that the [g]conformal connections are of the form where θ ∈ Ω 1 (M ), g ∇ denotes the Levi-Civita connection of g and θ the g-dual vector field of θ. Moreover, for u ∈ C ∞ (M ), we have [6 6, Theorem 1.159] exp(2u)g ∇ = g ∇ − g ⊗ g ∇u + du ⊗ Id + Id ⊗ du from which one easily computes (exp(2u)g,θ+du) ∇ = (g,θ) ∇.
Proof of Lemma 4. (i) the triples (g, A, θ) consisting of a Riemannian metric g, a differential A of degree m and a 1-form θ such that (ii) the triples ([g], ∂ Lm , P ) consisting of a conformal structure [g], a holomorphic line bundle structure ∂ Lm on L m and a holomorphic section P of L m ⊗ K m M having the property that the symmetric part of the Ricci curvature of the conformal connection associated to ∂ Lm plus (1 − m)P is negative definite.
Proof. Suppose (g, A, θ) is a triple satisfying We equip L m with the holomorphic line bundle structure induced by the conformal connection (g,θ) ∇. Define P := (det g) −(m−1)/4 ⊗ A, then P is a holomorphic section of L m ⊗ K m M . Indeed, we compute In addition, we observe that the symmetric part of the Ricci curvature of (g,θ) ∇ satisfies Sym Ric (g,θ) ∇ + (1 − m)P = K g − δ g θ + (1 − m)|A| 2 g g = −g which is obviously negative definite. Clearly, the just described map from the first set of triples into the second set of triples is injective.
Conversely, suppose L m is equipped with a holomorphic line bundle structure ∂ Lm and let P be a holomorphic section of L m ⊗ K m M . Assume furthermore that the symmetric part of the Ricci curvature of the conformal connection associated to ∂ Lm plus (1 − m)P is negative definite. We will next use these data to construct a triple (g, A, θ) solving the above equations. Let g 0 ∈ [g] denote the hyperbolic metric in the conformal equivalence class and define Note that (det g 0 ) (m−1)/4 is a non-vanishing section of L −1 m and hence A 0 is a section of K m M . Since P is holomorphic it follows that there exists a unique 1-form θ 0 on M such that Now make the Ansatz g = e 2u g 0 for u ∈ C ∞ (M ) and A = (det g) (m−1)/4 ⊗ P = A 0 e u(m−1) . Then where now all norms and operators are with respect to g 0 . This simplifies to become an algebraic equation for u e 2u − (m − 1)|A 0 | 2 = 1 + δθ 0 .
Clearly, this equation uniquely determines u provided 1 + δθ 0 + (m − 1)|A 0 | 2 is positive. Note that this happens if and only if is negative definite, but (g 0 ,θ 0 ) ∇ is just the conformal connection induced by ∂ Lm . Finally, by construction, the triple associated to (g, A, θ) is ([g], ∂ Lm , P ).

The case of holomorphic differentials
We have seen that a triple (g, A, θ) solving (4.6 4.6) and (4.7 4.7) yields a holomorphic section of L m ⊗ K m M with respect to some appropriate holomorphic line bundle structure on L m . We now restrict to the case where the differential A is already holomorphic so that we obtain the coupled vortex equations K g = −1 + (m − 1)|A| 2 g and ∂A = 0.

Anosov flows.
It is possible to upgrade Corollary 4.6 4.6 in the case where A is holomorphic as follows: Let (g, A) be a pair satisfying the coupled vortex equations∂A = 0 and K g = −1 + (m − 1)|A| 2 g . Then the associated thermostat flow is Anosov. Proof. We already know that there is a dominated splitting, so taking into account Remark 3.8 3.8, the strategy will be to show that r u > 0 and r s < 0. We will do this using the following lemma.
Lemma 5.2. Let (g, A) be a pair satisfying the coupled vortex equations∂A = 0 and Proof. The proof is quite similar to the proof of [4 4, Proposition 3.3], the reader may also compare with [10 10, Theorem 5.1]. The claim is obviously correct if A vanishes identically, hence we assume this not to be the case. We first prove the inequality K g 0. As before let g 0 denote the hyperbolic metric in the conformal equivalence class of g and write g = e 2u g 0 for u ∈ C ∞ (M ). Using where we write α = |A| 2 g 0 . The inequality K g 0 is equivalent to and is clearly satisfied at the points where A vanishes. Therefore, taking the logarithm of (5.3 5.3), we see that K g 0 follows from the non-negativity of the smooth function As M is compact, the Gauss curvature K g attains its maximum at some point x 0 and moreover x 0 ∈ M • . Consequently, the function f attains its infimum at x 0 . A straightforward calculation gives ∆ log α = −2m, where we use that A is holomorphic. At the minimum x 0 of f we thus obtain where we have used (5.4 5.4). It follows that f (x 0 ) 0 and hence f 0 on all of M • . This shows that K g 0. It order to prove K g < 0, we first remark that the function f − 1 + e −f is non-negative on M • . Consequently, (5.5 5.5) gives where ∆ g = e −2u ∆ denotes the Laplacian with respect to g. In particular, it follows that for every point x ∈ M • there exists a constant c > 0, an x-neighbourhood U x and a flat metric g 0 on U x which lies in the conformal equivalence of g, so that on U x . Therefore, by applying the strong maximum principle [16 16,Theorem 3.5] to the operator ∆ g 0 − c, it follows that if f vanishes at some point in U x , then it vanishes on all of U x and consequently on M • . Since A is holomorphic, its zeros are isolated and hence M • is dense in M . Since K g is continuous we conclude that if K g vanishes at some point on M , then it vanishes identically on M , but this possibility is excluded by the Gauss-Bonnet theorem. Since α 0 we have G(x, u) −1 + e 2u and hence G(x, u) < 0 for u < 0. On the other hand, for u > sup x∈M Since ∂G ∂u (x, u) = 2α(x)(m − 1) 2 e −2(m−1)u + 2e 2u > 0 standard quasi-linear elliptic PDE methods (see for instance [37 37, Proposition 1.9]) imply that (5.2 5.2) has a unique smooth solution u for every smooth non-negative function α. Consequently, for every holomorphic differential A on (M, [g]) we obtain a unique solution (g, A) to the coupled vortex equations K g = −1 + (m − 1)|A| 2 g and ∂A = 0. We now show that r u > 0 (the proof that r s < 0 is similar). Set h = r u − V (a)/m. Then h satisfies Given (x, v) ∈ SM , consider for each R > 0, the unique solution h R to the Riccati equation where w := f R + B + 2 and q := − 2 − B + 1. Observe that q 0 by our definitions of c and . We can solve the inhomogeneous linear equation (5.7 5.7) and use that q 0 to derive f R (t) 0 and thus h R (t) . By taking limits, and using (5.6 5.6), we obtain Remark 5.4. As we have seen, Corollary 4.6 4.6 asserts that given a triple (g, A, θ) satisfying (4.6 4.6) and (4.7 4.7), the associated thermostat flow has a dominated splitting. When θ = 0, Theorem 5.1 5.1 tells us that we can do better and in fact the thermostat flow is Anosov. At the "other end", that is, when A = 0, we also know by Proposition 3.5 3.5 that the thermostat flow is also Anosov (in this case K = K g − δ g θ = −1). These two "ends" are Anosov for different reasons, connected with the discussion in Remark 3.8 3.8. In the case θ = 0, as we have just seen, one uses that r u > 0, that is, the first case in (3.3 3.3). In the case A = 0, we use the second case in (3.3 3.3). It is conceivable that the thermostat flow is always Anosov for any triple (g, A, θ) satisfying (4.6 4.6) and (4.7 4.7), but at the time of writing it is not at all clear how to prove this. It should be noted that for the special case of the geodesic flow it is well known that a dominated splitting must be Anosov. We can see this fairly quickly using quadratic forms as follows. Suppose r u,s : SM → R are two continuous functions such that Xr u,s + [r s,u ] 2 + K g = 0 and r u − r s = 0 everywhere. Define Then a calculation showsQ = (ẏ − r u y) 2 + (ẏ − r s y) 2 > 0 unless y =ẏ = 0. Hence by Proposition 3.4 3.4 the geodesic flow is Anosov.

Dissipation and volume.
We will now prove the following result stated in the introduction. We write the volume form as e −u Θ for some real-valued function u on SM . Thus, using (2.6 2.6), we obtain Hence the claim follows by showing that if u solves F u = V a, then a vanishes identically. In order to show this we use the following L 2 identity proved in [22 22,Equation (5)] which is in turn an extension of an identity in [36 36] for geodesic flows. The identity holds for arbitrary thermostats F = X + λV . If we let H c := H + cV where c : SM → R is any smooth function then where u is any smooth function. All norms and inner products are L 2 with respect to the volume form Θ. In our case λ = a and a calculation shows that if we pick c = V (a)/m, then , hence for this choice of c, (5.8 5.8) simplifies to If F u = V a, then V F u = −m 2 a and we compute using that X and H preserve Θ and that XV a − mHa = 0: where the last equation is obtained using that Xu = V a − aV u and c = V (a)/m. Inserting this back into (5.9 5.9), we see that the equality obtained can only hold if V a and hence a vanishes identically.
6. The cases m = 2 and m = 3 In this section we consider the special cases of m = 2, 3 and their peculiarities. These flows have appeared in different contexts and for different reasons and in this section we explain these features.
6.1. The case m = 2. Consider a pair (g, A) where A is a quadratic differential with ∂A = 0 and K g = −1 + |A| 2 g . By Theorem 5.1 5.1, the associated thermostat flow is Anosov. These flows have the distinctive feature that their weak bundles are of class C ∞ . Indeed for this case p = V (a)/2, κ p = −1 and equation (3.4 3.4) From this we clearly see that r u,s = ±1 + V (a)/2 and hence the weak bundles RF ⊕ R(H + r s,u V ) are smooth. This class of thermostats flows was first considered in [33 33], where the coupled vortex equations for m = 2 were derived assuming that the weak foliations were smooth. Theorem 4.6 in [ 15 15] asserts that a smooth Anosov flow on a closed 3-manifold with weak stable and unstable foliations of class C 1,1 , is smoothly orbit equivalent to a suspension or to a quasi-fuchsian flow as described in [14 14,Théorème B]. (In our case, since we are working with circles bundles the latter alternative holds.) A quasi-fuchsian flow ψ depends on a pair of points ([g 1 ], [g 2 ]) in Teichmüller space, has smooth weak stable foliation C ∞ -conjugate to the weak stable foliation of the constant curvature metric g 1 and smooth weak unstable foliation C ∞ -conjugate to the weak unstable foliation of the constant curvature metric g 2 . Moreover, ψ preserves a volume form if and only if [g 1 ] = [g 2 ]. The analogous result on the thermostat side is provided by Theorem 5.5 5.5 which asserts that the thermostat flow preserves a volume form iff A = 0. It is an interesting question (first raised in [33 33]) to decide if the thermostat flows originating from the coupled vortex equations∂A = 0, K g = −1 + |A| 2 g describe all possible quasi-fuchsian flows ψ.
6.2. The case m = 3. Let now (g, A, θ) be a triple on M satisfying (4.6 4.6) and (4.7 4.7) with A being a cubic differential. The connection form of the Levi-Civita connection on the tangent bundle T M is 0 −ψ ψ 0 .
We define a 1-form on SM with values in gl(2, R) It is a consequence of the equivariance properties that the 1-form Υ is the connection 1-form of a unique (torsion-free) connection ∇ on the tangent bundle T M . Moreover, since the interior product i F Υ 2 1 vanishes identically for λ = a − V θ, it follows that the geodesics of the connection ∇ can be reparametrised to agree with the projections to M of the orbits of the thermostat flow defined by λ, see [32 32, Lemma 3.1] for details. Moreover, if θ is closed the connection ∇ admits an interpretation as a Lagrangian minimal surface, see [31 31]. If A is holomorphic so that θ vanishes identically, then the connection ∇ defines a properly convex projective structure on M , see the work of Labourie [25 25] and [30 30, 31 31]. This means that the universal cover Ω of M is a properly convex open subset of the real projective plane RP 2 for which there exists a discrete group Γ of projective transformations which acts cocompactly on Ω and so that M = Ω/Γ. Thus, (Ω, Γ) is a divisible convex set. Moreover, the segments of the projective lines RP 1 contained in Ω project to M to agree with the (unparametrised) geodesics of ∇. The universal cover Ω being a convex set, it is equipped with the Hilbert metric. The geodesic flow of the Hilbert metric descends to SM and by a result of Benoist [3 3], is Anosov if and only if Ω is strictly convex. In [3 3], it is also shown that a divisible convex set is strictly convex if and only if the group dividing it is word-hyperbolic. Since the fundamental group of a closed surface of negative Euler characteristic is word-hyperbolic, it thus follows from known results that the thermostat flow associated to a holomorphic cubic differential is a reparametrisation of an Anosov flow. However, since the Anosov property is invariant under reparametrisation of the flow, we conclude that the thermostat flow associated to a holomorphic cubic differential is Anosov, which is the statement of our Theorem 5.1 5.1 for the special case m = 3.

Regularity of weak foliations
As we previously mentioned, the case of m = 2 has the distinctive feature of having weak bundles of class C ∞ . It is natural to ask what happens for m ≥ 3. One approach to this question would be to compute the Godbillon-Vey invariant following [33 33]. Unfortunately for m ≥ 3 this calculation does not yield information conducive to an answer. However, for the case m odd, we can use reversibility of the flow combined with Theorem 5.5 5.5 to derive: Theorem 7.1. Suppose an Anosov thermostat given by the coupled vortex equations has a weak foliation of class C 2 and m is odd. Then A vanishes identically.
Proof. When m is odd there is an important additional symmetry in the flow: the flip σ given by (x, v) → (x, −v). We note that this map is isotopic to the identity. If φ denotes the thermostat flow then, σ • φ t = φ −t • σ. This relation easily implies that σ maps the weak stable foliation to the unstable one. Hence, if one of them is of class C 2 , the other one is also of class C 2 .
As we have already mentioned, Theorem 4.6 in [ 15 15] asserts that a smooth Anosov flow on a closed 3-manifold with weak stable and unstable foliations of class C 2 , is smoothly orbit equivalent to a quasi-fuchsian flow ψ that depends on a pair of points ([g 1 ], [g 2 ]) in Teichmüller space. The flow ψ has smooth weak stable foliation C ∞ -conjugate to the weak stable foliation of the constant curvature metric g 1 and smooth weak unstable foliation C ∞ -conjugate to the weak unstable foliation of the constant curvature metric g 2 . But since σ is isotopic to the identity we must have [g 1 ] = [g 2 ] and ψ is an ordinary geodesic flow preserving a volume form. Thus our thermostat flow preserves a volume form and by Theorem 5.5 5.5 we must have A = 0.
Remark 7.2. It is instructive to discuss Theorem 7.1 7.1 in the light of the remarks in Section 6 6 for m = 3. As pointed out, in this case, the thermostat flow is a C ∞ parametrisation of the geodesic foliation of a Hilbert metric. Benoist observes in [3 3] that the regularity of the weak foliations of the Hilbert geodesic flow coincides with the regularity of the boundary. Hence if the boundary of the strictly convex domain defining the Hilbert metric is C 2 , then the associated thermostat flow also has C 2 weak foliations and therefore A = 0. This implies that the convex domain is an ellipsoid, thus recovering a result of Benzécri [5 5] for the case of 2-dimensional domains (note however, that the proof in [5 5] is more direct and straightforward).

The path geometry defined by a thermostat
A thermostat naturally defines a path geometry and in this final section we show that the path geometry associated to the thermostat coming from a holomorphic differential A of degree m 2 is flat if and only if A vanishes identically or m = 3. The former case corresponds to the paths being the geodesics of a hyperbolic metric and the latter case to the paths being the geodesics of a convex projective structure. We first recall some elementary facts about path geometries while referring the reader to [7 7] for further details.
An (oriented) path geometry on an oriented surface M is given by an oriented line bundle L on the projective circle bundle SM := (T M \ {0}) /R + having the property that L together with the vertical bundle of the projection map ν : SM → M spans the contact distribution of SM . The paths of L are the projections of its integral curves to M . Note that the orientation of L naturally equips its paths with an orientation.
Example 8.1. Taking M to be the oriented 2-sphere S 2 , we obtain a canonical path geometry L 0 whose paths are the great circles. In this case Let now F = X + λV be a thermostat on the unit tangent bundle SM of a oriented Riemannian 2-manifold (M, g). We henceforth identify SM SM in the obvious way. In doing so, we obtain a path geometry by defining L := RF and by declaring vectors in L to be positive if they are positive multiples of F .
Clearly, if a path geometry is flat, then it must have the property that its paths agree with the geodesics of some projective structure. In [32 32,Proposition 3.4] it is shown that the path geometry defined by a thermostat X + λV shares its paths with the geodesics of some projective structure if and only if Using this fact we immediately obtain: Consequently, a and hence A must vanish identically or m = 3. Conversely, assume A is a cubic differential satisfying ∂A = 0 and K g = −1 + 2|A| 2 g . The path geometry associated to (g, A) defines a properly convex projective structure on the oriented surface M . An oriented properly convex projective surface is an example of a surface carrying a (G, X)-structure where X = S 2 is the oriented projective 2-sphere and G = SL(3, R) its group of projective transformations, cf. [21 21]. In particular, it follows that the path geometry associated to (g, A) is flat.