## Relative Helicity and Tiling Twist

Boris Khesin

Nicolau C. Saldanha

## Abstract

We consider domino tilings of 3D cubiculated regions. The tilings have two invariants, flux and twist, often integer-valued, which are given in purely combinatorial terms. These invariants allow one to classify the tilings with respect to certain elementary moves, flips and trits. In this paper we present a construction associating a divergence-free vector field ξ t to any domino tiling t , such that the flux of the tiling t can be interpreted as the (relative) rotation class of the field ξ t , while the twist of t is proved to be the relative helicity of the field ξ t .

## 1 Introduction

Domino tilings of 3D regions traditionally have two invariants associated to them, flux and twist. Flux is understood as the homology class associated with a certain cycle constructed for the tiling. For a cubiculated region (a topological 3-manifold with boundary) the flux Flux( ) of a tiling t t is an element of the first homology group of the region, defined up to an additive constant. The ambiguity can be removed by considering the relative flux RFlux( ) as discussed below. t For tilings with zero relative flux there is an integer invariant, the twist , associated with the tiling and measuring the 'mutual linking of tiles' around each other. It is invariant with respect to flips , local moves which consist of removing two neighboring parallel dominoes and placing them back after a rotation. The twist changes under another move, a trit , which replaces a frame-like triple of tiles to the one pointing in the opposite way (see Section 2 and a detailed discussion in [8]).

While there have been pointed out similarities of the twist invariant with the Hopf invariant, the correspondence remained at either an intuitive level or in the continuous limit for turning tiles into vector fields via a broadly understood tiling's refinement.

In this paper we present a construction of a smooth divergence-free vector field associated to an arbitrary tiling ('5-pipe construction', see Section 5.1) so that the twist invariant becomes, up to a factor, the relative helicity of that vector field. Furthermore, we extend the notion of relative helicity [5] to divergence-free vector fields on arbitrary three-manifolds and not necessarily tangent to their boundaries. This allows one to compare relative helicity with twists of tilings in nonsimply-connected regions (see Section 4.3). The toolbox includes an introduction of 'an isolating shell' for a cubiculated region, the use of refinements and appropriate connectivity of the spaces of tilings. Finally, we relate the flux invariant of a tiling to the rotation class of the associated vector field.

In a nutshell, the results of the paper are as follows. Recall that for a divergence-free vector field ξ in a simply-connected domain M ⊂ R 3 and tangent to its boundary ∂M , its helicity is the quantity Hel( ξ ) = ∫ M ( ξ, curl -1 ξ ) d x 3 . For a field whose support consists of several linked tubes the helicity reduces to the mutual linking of those tubes, the self-linking and fluxes in the tubes (see [12] and Section 4.2 and Appendix).

0 2010 Mathematics Subject Classification . Primary 05B45; Secondary 57K12, 52C22, 05C70. Keywords and phrases: Three-dimensional tilings, dominoes, dimers, flux, helicity, rotation class, relative helicity.

The helicity notion can be extended to null-homologous vector fields tangent to the boundary in arbitrary three-manifolds M equipped with a volume form. If the field ξ is not tangent to the boundary ∂M , one can define only its relative helicity, i.e., the difference of helicities of two vector fields ξ and η with the same behaviour at the boundary ∂M . Namely, one can extend ξ and η by the same field to vector fields in a bigger manifold ˜ M ⊇ M , where their helicities are well-defined, but depend on the extension. However, the difference of helicities will not depend on the extension, hence the name relative helicity , see Section 4.3.

On the other hand, given a 3D domino tile, consisting of two cubes (say, white and black), consider five smooth curves joining the symmetric faces at their centers and approaching those faces orthogonally, see Figure 4. Tubular neighborhoods of those curves will become supports of the five smooth divergence-free vector fields directed from the black to the white unit cube and each having flux φ . By performing this construction in each domino, to any tiling t of a cubiculated region R one associates a smooth divergence-free vector field ξ t in the whole region. By using an 'isolating shell' for the region, one can define the relative helicity of the vector field ξ t associated with the tiling. The construction of such an isolating shell turns out to be possible whenever the relative flux vanishes, RFlux( t ) = 0 ∈ H 1 ( R , ∂ R ). Our main result is the following:

Theorem 1.1. For two tilings t 0 and t 1 of the same flux and of zero relative flux in a 3D cubiculated region R , the difference of their twists is proportional to the relative helicity of the associated vector fields ξ t 0 and ξ t 1 :

$$\text{Hel} ( \xi _ { t _ { 1 } } ) - \text{Hel} ( \xi _ { t _ { 0 } } ) = 3 6 \varphi ^ { 2 } \left ( \text{Tw} ( \mathfrak { t } _ { 1 } ) - \text{Tw} ( \mathfrak { t } _ { 0 } ) \right ) \,,$$

where φ ∈ R is the flux in a single pipe of the vector fields ξ t i .

In order to clarify the condition of vanishing relative flux, recall that given a vector field ξ in a manifold M with a volume form µ and tangent to the boundary ∂M , the 2-form ω ξ := i ξ µ is closed iff ξ is divergence-free with respect to µ . If ξ is tangent to ∂M , the restriction of the 2-form ω ξ vanishes on the boundary ∂M . The cohomology class [ ω ξ ] ∈ H 2 ( M,∂M ) can be identified with a class in H M 1 ( ) (via the Poincar´ e isomorphism) and is called the rotation class of the dynamical system ξ , see [1]. Similarly, if a divergence-free vector field ξ is not tangent to ∂M , then [ ω ξ ] ∈ H 2 ( M ) ≃ H M,∂M 1 ( ). Thus such a vector field ξ defines a relative rotation class as an element of H M,∂M 1 ( ).

In Section 5.2 we prove the relation of flux Flux( t ) and relative flux RFlux( t ) for a tiling t to the (relative) rotation class of an appropriate vector field ξ t :

Theorem 1.2. Given a tiling t of a 3D cubiculated region R its relative flux RFlux( ) t ∈ H 1 ( R , ∂ R ) coincides modulo a factor with the relative rotation class [ ξ t ] of the vector field ξ t obtained via the 5-pipe construction:

$$[ \xi _ { t } ] = 6 \varphi \, R F l u x ( t ) \,.$$

In particular, the relative rotation class of the field ξ t vanishes iff the relative flux RFlux( t ) vanishes. By normalizing the flux as φ = 1 6, consistent with the 6-pipe setting of Section 5.2, / the two main theorems provide the exact matching: for any tiling its relative flux coincides with the relative rotation class of the associated vector field, while the relative helicity of the latter coincides with the tiling's twist.

Remark 1.3. There are numerous papers studying tilings in the context of quantum dimer models. In the papers [7, 4] the authors emphasized the existence of topologically non-trivial configurations called 'Hopfions' in those systems with the understanding that in a suitable continuous limit of the dimer model, those excitations would correspond to the Hopf map S 3 → S 2 . In [2] (see also [3]), Arnold proved that an asymptotic version of the Hopf invariant for a vector field is equal to

the field's helicity. Theorem 1.1 shows that even without taking the continuous limit, a tiling's twist itself is already equal (modulo the flux factor) to the corresponding relative helicity.

In a sense, the present paper shows that one can regard tilings as occupying an intermediate place between divergence-free vector fields on the one hand and links with framings on the other: by associating special divergence-free vector fields to tilings, one can compute their helicities via consideration of linking and self-linking of framed cycles. ⋄

The paper is organized as follows. In Sections 2 and 3 we recall the main moves and invariants for domino tilings in purely combinatorial terms. In Section 4 we introduce relative rotation class and relative helicity of vector fields not tangent to boundary. In Section 5 we present the 5-pipe construction associating a divergence-free vector field to a domino tiling, introduce an isolating shell, and prove Theorem 1.2 on the relative rotation class. Section 6 presents two key examples which are used in the proof of the main Theorem 1.1 in Section 7. One more example manifesting the main theorem is described in Section 8. Appendix 9 derives the formula for helicity via linking and self-linking numbers of pipes.

Acknowledgements. We are grateful to the anonymous referee for useful remarks. B.K. is indebted to PUC-Rio and IMPA in Rio de Janeiro and IHES in Bures-sur-Yvette for their kind hospitality. He was partially supported by an NSERC Discovery Grant. N.S. is thankful for the generous support of CNPq, CAPES and FAPERJ (Brazil).

## 2 Domino tilings, combinatorial flux and local moves

A cubiculated region is a connected and oriented manifold R (usually with boundary) decomposed into finitely many unit cubes. A simple example of cubiculated region is a box : R = [0 , L ] × [0 , M ] × [0 , N ] ⊂ R 3 , LMN even; the unit cubes are [ a, a +1] × [ b, b +1] × [ c, c +1], ( a, b, c ) ∈ Z 3 . We assume that the interior of a cubiculated region is as in this example, so that interior edges (resp. vertices) are surrounded by four (resp. eight) unit cubes. We also assume that unit cubes are painted black and white, with adjacent cubes of opposite colors and the same number of cubes of each color. A (3D) domino is the union of two unit cubes with a common face, thus a 2 × 1 × 1 rectangular cuboid. A (domino) tiling of R is a set of dominoes with disjoint interiors whose union is equal to R .

Figure 1: A tiling of the box [0 , 4] × [0 , 4] × [0 , 2]. The orientation of R 3 is important: the z axis points upward away from the paper. Examples of dominoes in this tiling are [0 1] , × [0 , 2] × [0 , 1], [0 , 2] × [0 , 1] × [1 , 2] and [1 , 2] × [1 , 2] × [0 , 2]. This tiling admits no flips.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-8e09f0aa-ba1c-4283-a6f0-a4e2c4d1a334.jpg)

We follow [10], [14] and [8] in drawing tilings of cubiculated regions by floors, as in Figure 1. Vertical dominoes (i.e., dominoes in the z direction) appear as two squares, one in each of two adjacent floors; the top square of a vertical domino (which, in the figure, appears at the right) is left unfilled for visual facility.

Definition 2.1. Given a region R , let T ( R ) denote the set of domino tilings of R . A flip is a local move in T ( R ): two parallel and adjacent (3D) dominoes are removed and placed back in a different position. A trit is the only local move in T ( R ) involving three dominoes which does not reduce to flips. The three dominoes involved are in three different directions and fill a 2 × 2 × 2 box minus two opposite unit cubes. ⋄

Figure 2 shows in the 3 × 3 × 2 box a trit followed by a sequence of flips (the first tiling admits no flips). We shall come back to this example in Section 8.

Figure 2: Five tilings of the box [0 , 3] × [0 , 3] × [0 , 2]. The first move is a trit, the two other moves are flips. Five more flips take us from the fourth to the fifth tiling. The first tiling has twist -1, the others have twist 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-170a5f98-87e6-406b-9b3f-36e153960fb9.jpg)

To apply homology theory we will use a finer complex structure R ♯ . The vertices of R ♯ are the original vertices of R , together with the centers of edges, faces and unit cubes of R . Each edge of R is decomposed into two edges of R ♯ . Similarly, each square (resp. cube) of R is decomposed into four squares (resp. eight cubes) of R ♯ . Notice that the boundary ∂ R also acquires the structure of a polyhedral complex ∂ R ♯ .

A domino is identified with the sum of two oriented edges in R ♯ , from the center of the black cube to the center of the common face and from there to the center of the white cube. A domino tiling t of R is therefore identified as t ∈ C 1 ( R ♯ ). The boundary ∂ t ∈ C 0 ( R ♯ ) is the sum of the centers of the white cubes (of R ) minus the sum of the centers of the black cubes. In particular, for any two tilings t 0 , t 1 we have ∂ t 0 = ∂ t 1 or, equivalently, t 1 -t 0 ∈ Z 1 ( R ♯ ). We then write [ t 1 -t 0 ] ∈ H 1 ( R ) = H 1 ( R ; R ); in the present paper we almost always use coefficients in R .

Given an initial tiling t 0 , let T 0 ( R ) ⊆ T ( R ) be the set of tilings t belonging to the same 'homology class', i.e. such that [ t -t 0 ] = 0 ∈ H 1 ( R ). We will say that tilings t have the same flux as t 0 , Flux( ) = Flux( t t 0 ).

In the present paper we define a somewhat different object, the relative flux RFlux( ) t ∈ H 1 ( R , ∂ R ). It takes values in relative homology and does not require the choice of a base tiling. We work in the complex R ♯ and, given a tiling t , define t ∈ C 1 ( R ♯ ) as above. We construct q 1 ∈ C 1 ( R ♯ ) as follows: it is the sum over all black unit cubes of all edges from the cube centers to the centers of their faces minus the same sum (of all edges from cube centers to their faces) over all white cubes. Thus, for every tiling t we have that q 0 = ∂ (6 t -q 1 ) ∈ C 0 ( ∂ R ♯ ) is the sum of +1 (resp. -1) times the centers of the white (resp. black) squares in ∂ R . Thus, 6 t -q 1 ∈ Z 1 ( R ♯ , ∂ R ♯ ); we define

$$R \text{Flux} ( \mathfrak { t } ) = \left [ \mathfrak { t } - \frac { \mathfrak { q } _ { 1 } } { 6 } \right ] \in H _ { 1 } ( \mathcal { R }, \partial \mathcal { R } ).$$

Inclusion induces a homomorphism i : H 1 ( R ) → H 1 ( R , ∂ R ), while for any tiling t one has i (Flux( )) = RFlux( ) t t -RFlux( t 0 ).

Example 2.2. Already for cubiculated regions R ⊂ R 3 , there are examples with nontrivial homology or homotopy. The solid torus

$$\mathcal { R } _ { 0 } = ( [ 0, 1 2 ] \times [ 0, 1 2 ] \times [ 0, 4 ] ) \smallsetminus ( ( 4, 8 ) \times ( 4, 8 ) \times [ 0, 4 ] )$$

is not simply connected, with π 1 ( R 0 ) = Z , H 1 ( R 0 ; R ) = R and H 1 ( R 0 , ∂ R 0 ; R ) = 0. There exist tilings of R 0 with different values of flux, but the relative flux is always 0.

The cube with a hole

$$\mathcal { R } _ { 1, 0 } = ( [ 0, 1 2 ] \times [ 0, 1 2 ] \times [ 0, 1 2 ] ) \smallsetminus ( ( 4, 8 ) \times ( 4, 8 ) \times ( 4, 8 ) )$$

is simply connected but has nontrivial relative homology: H 1 ( R 1 0 , , ∂ R 1 0 , ) = R . For any tiling t ∈ T ( R 1 0 , ) we have RFlux( ) = 0. For a different cube with hole t

$$\mathcal { R } _ { 1, 1 } = ( [ 0, 1 3 ] \times [ 0, 1 3 ] \times [ 0, 1 3 ] ) \smallsetminus ( ( 4, 9 ) \times ( 4, 9 ) \times ( 4, 9 ) )$$

̸

we have RFlux( ) = 0 for all t t ∈ T ( R 1 1 , ), since the numbers of white and black cubes in the hole surrounded by R 1 1 , are different, see [8].

Quotients and similar identifications give us other interesting examples of regions. The quotient

$$\mathcal { R } _ { 2 } = \mathbb { R } ^ { 3 } / ( 6 \mathbb { Z } ) ^ { 3 }$$

is a cubiculated region homeomorphic to a 3-torus T 3 = ( S 1 ) 3 (see Figures 2, 7 and 8 in [8] for examples of tilings of R 2 ); since ∂ R 2 = , we have ∅ H 1 ( R 2 ) = H 1 ( R 2 , ∂ R 2 ).

The quotient

$$\mathcal { R } _ { 3 } = ( \mathbb { R } ^ { 2 } / ( 8 \mathbb { Z } ) ^ { 2 } ) \times [ 0, 4 ]$$

is a cubiculated region homeomorphic to T 2 × [0 , 1], so that H 1 ( R 3 ) = R 2 . It admits tilings with different values of the flux. We have H 1 ( R 3 , ∂ R 3 ) = R , while the map i : H 1 ( R 3 ) → H 1 ( R 3 , ∂ R 3 ) is the zero homomorphism, i = 0 and RFlux( ) = 0 for all t t ∈ T ( R 3 ).

We may also remove a few cubes from a larger region; for instance, R 4 = R 2 ∖ (2 4) , 3 is a 3-torus with a hole. ⋄

̸

Remark 2.3. The examples above demonstrate independence of the assumptions of zero relative flux and two tilings having the same flux. Indeed, for R 2 the relative and absolute homology groups coincide and one can have two tilings with the same non-zero flux, [ t 1 ] = [ t 0 ] = 0 ∈ H 1 ( R 2 ) = H 1 ( R 2 ; ∂ R 2 ). On the other hand, all tilings of R 3 have zero relative flux, RFlux( t ) = 0, while realizing various fluxes in H 1 ( R 3 ) = R 2 . ⋄

## 3 Twist of a tiling

Recall that given an initial tiling t 0 we denote by T 0 ( R ) the set of tilings t with Flux( ) = [ t t -t 0 ] = 0 ∈ H 1 ( R ). In general, the twist of a tiling is a map Tw : T 0 ( R → ) Z / m ( Z ) where m ∈ Z , m ≥ 0. In the present paper we are concerned with the case m = 0, where the twist map is a function Tw : T 0 ( R → ) Z well-defined up to an additive constant. In this case, one can define the twist of a tiling recursively and, if needed, using refinements as follows (we refer to [8] for a general combinatorial definition of the twist and more details).

Given a cubiculated region R , the refinement R ′ of R is obtained by decomposing each unit cube in R into 5 × 5 × 5 smaller cubes. Given a tiling t of R , its refinement t ′ is a tiling of R ′ : each domino d in t is decomposed into 5 × 5 × 5 smaller dominoes, all parallel to d . We write R (0) = R and R ( k +1) = ( R ( k ) ) ′ and define t ( k ) ∈ T ( R ( k ) ) in a similar manner.

Theorem 3.1. ([8]) Consider a cubiculated region R and two tilings t 0 , t 1 of R . If Flux( t 1 ) = Flux( t 0 ) then there exist k such that the tilings t ( k ) 0 and t ( k ) 1 can be joined by a finite sequence of flips and trits.

Define the set of all refinements of T 0 ( R ) by

$$\mathcal { T } _ { 0 } ( \mathcal { R } ^ { ( * ) } ) = \bigsqcup _ { k \in \mathbb { N } } \mathcal { T } _ { 0 } ( \mathcal { R } ^ { ( k ) } )$$

where T 0 ( R ( k ) ) ⊆ T ( R ( k ) ) is the subset of tilings t with [ t -t ( k ) 0 ] = 0 ∈ H 1 ( R ). One can interpret T 0 ( R ( ∗ ) ) as the set of vertices of an infinite graph: two vertices t 0 and t 1 are connected by an edge if the tilings t 0 and t 1 differ by a flip, a trit or by one step of a refinement. The graph T 0 ( R ( ∗ ) ) is connected, as follows directly from Theorem 3.1.

The following result gives us a recursive definition of twist in T 0 ( R ). It is proved in [8] in different terms and (in certain cases before that) in [11]. Below we provide a different proof.

Theorem-Definition 3.2. (cf. [8]) Consider a cubiculated region R , with an initial tiling t 0 . If the relative flux of t 0 is zero then there exists a function Tw : T 0 ( R ( ∗ ) ) → Z (unique up to an additive constant) with the following properties:

- 1. If two tilings t a , t b ∈ T 0 ( R ( ∗ ) ) are joined by a flip then Tw( t b ) = Tw( t a ) .
- 2. If t a and t b are joined by a trit then Tw( t b ) = Tw( t a ) ± 1 ; sign is given by the orientation of the trit.
- 3. If t ′ is the refinement of t ∈ T 0 ( R ( ∗ ) ) then Tw( t ′ ) = Tw( ) t .

The function Tw is called the twist function .

We say that a trit from t a to t b is positive (respectively, negative ) if Tw( t b ) -Tw( t a ) = +1 (respectively, -1).

̸

Remark 3.3. If RFlux( t 0 ) = 0, there is m &gt; 0 such that the twist function becomes a map Tw : T 0 ( R → ) Z / m ( Z ), which, in addition to the properties listed above satisfies the following. There exist a finite sequence of tilings t 0 , t 1 , . . . , t ℓ = t 0 in T 0 ( R ( k ) ) (for some k ≥ 0) and a sequence of integers t 0 , t 1 , . . . , t ℓ such that for every i , t i and t i +1 are joined by a flip or trit, Tw( t i ) = t i mod m , and t ℓ = t 0 + m .

One should stress that all known definitions of the twist of a tiling involve a significant combinatorial part and are relatively complicated, cf. [8]; one of the main aims of the present paper is to give an alternative geometric description. ⋄

## 4 Rotation class and helicity of vector fields

## 4.1 Relative rotation class

Let M be a manifold with boundary ∂M and volume form µ . Consider a divergence-free vector field ξ on M and not necessarily tangent to its boundary.

Definition 4.1. ([1]) Consider the 2-form ω ξ := i ξ µ . It is closed, since ξ is divergence-free, and hence it defines the cohomology class [ ω ξ ] ∈ H 2 ( M ). By the Lefschetz isomorphism (the Poincar´ e isomorphism for manifolds with boundary) H 2 ( M ) ≃ H M,∂M 1 ( ). The relative homology class defined by [ ω ξ ] in H M,∂M 1 ( ) is called a relative rotation class of the field ξ .

If ∂M = ∅ or ξ is tangent to ∂M , then ω ξ | ∂M = 0 and hence the corresponding rotation class of the vector field ξ is an element in H M 1 ( ) ≃ H 2 ( M,∂M ). ⋄

This definition holds for M of any dimension, although our main application is in 3D.

## 4.2 Helicity of fields and pipes

Now we assume that M ⊂ R 3 is a three-dimensional domain and the field ξ is divergence-free, tangent to ∂M , and has zero rotation class, i.e. it is exact or null-homologous ).

Definition 4.2. ([12]) The helicity of a null-homologous field ξ in a domain M ⊂ R 3 is the number

$$\text{Hel} ( \xi ) \coloneqq \int _ { M } ( \xi, \text{curl} ^ { - 1 } \xi ) \ d ^ { 3 } x,$$

where the vector field curl -1 ξ is a divergence-free vector potential of the field ξ (which exists since ξ is exact), i.e., ∇× (curl -1 ξ ) = ξ and div(curl -1 ξ ) = 0. ⋄

Let a divergence-free field ξ be confined to two narrow linked flux tubes. Its helicity can be found explicitly as follows. Suppose that the tube cores are closed curves C 1 and C 2 , the fluxes of the field in the tubes are Flux 1 and Flux 2 (see Figure 3). The curves C i are oriented so that Flux i = Q i ≥ 0. Assume also that there is no net twist within each tube or, more precisely, that the field trajectories foliate each of the tubes into pairwise unlinked circles and the periods of those trajectories are equal. One can show that the helicity invariant of such a field is given by

$$\text{Hel} ( \xi ) = 2 \, l k ( C _ { 1 }, C _ { 2 } ) \cdot \text{Flux} _ { 1 } \cdot \text{Flux} _ { 2 },$$

where lk ( C , C 1 2 ) is the (Gauss) linking number of C 1 and C 2 , which explains the term 'helicity' coined in [12], as the measure of coiling one curve about the other. Recall, that the number lk ( C , C 1 2 ) for two oriented closed curves is the signed number of the intersection points of one curve with an arbitrary oriented surface spanning the other curve.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-28b0a04e-d7fd-4949-912d-c4769c6a6861.jpg)

/C9 /BD

Figure 3: Two linked tubes: the solid tori are thin tubular neighborhoods of the curves C i . A few cross-sections are shown, while the vector field is transversal to these sections.

Suppose now that a divergence-free field ξ is confined to several narrow linked oriented tubes with cores given by closed curves C i and fluxes Flux i ≥ 0. We still assume that the field trajectories foliate each of the tubes into topological circles. However we do not assume any longer that these field trajectories (circles) inside the same tubes are pairwise unlinked, while their periods can now be arbitrary. The corresponding linking of circles inside the i th tube (given by the linking number of the core curve C i with any satellite curve) is called the self-linking slk ( C i ) of the curve C i in the corresponding pipe.

## Proposition 4.3. The helicity of such a field ξ is

$$\text{Hel} ( \xi ) = 2 \sum _ { i < j } l k ( C _ { i }, C _ { j } ) \cdot \text{Flux} _ { i } \cdot \text{Flux} _ { j } + \sum _ { i } s l k ( C _ { i } ) ( \text{Flux} _ { i } ) ^ { 2 }.$$

We recall the proof of this folklore formula in Appendix, and will use it for computations related to tilings. Note that self-linking cannot be defined for an isolated closed curve (since there is no canonical choice of a satellite curve), but only for a curve with a framing, which delivers such a choice. For a field trajectory in a pipe, neighbouring trajectories provide the framing.

Remark 4.4. While helicity was defined above by using the Riemannian metric on M , it is actually a topological characteristic of a divergence-free vector field, depending only on the choice of a volume form on the manifold. Namely, consider a manifold M (possibly with boundary) with a volume form µ , and let ξ be a null-homologous vector field on M (tangent to the boundary). The divergence-free condition means that the Lie derivative of µ along ξ vanishes: L µ ξ = 0, or, which is the same, the substitution i ξ µ =: ω ξ of the field ξ into the 3-form µ is a closed 2-form: dω ξ = 0. If ξ is moreover null-homologous, then ω ξ is actually an exact 2-form: ω ξ = dα for some 1-form α , called a potential. (On a simply connected manifold M any divergence-free vector field is null-homologous.) ⋄

Definition 4.5. ([2]) The helicity Hel( ξ ) of a null-homologous field ξ on a three-dimensional manifold M equipped with a volume element µ is the integral of the wedge product of the form ω ξ := i ξ µ and its potential:

$$\text{Hel} ( \xi ) = \int _ { M } d \alpha \wedge \alpha, \text{ where } d \alpha = \omega _ { \xi }.$$

As discussed, this generalizes Definition 4.2.

⋄

An immediate consequence of this purely topological (i.e. metric-free) definition is the following Arnold's theorem: The helicity Hel( ξ ) is preserved under the action on ξ of a volume-preserving diffeomorphism of M . In this sense Hel( ξ ) is a topological invariant: it can be defined without coordinates or a choice of metric, and hence every volume-preserving diffeomorphism carries a field ξ into a field with the same helicity. The physical significance of helicity is due to the fact that it appears as a conservation law in both fluid mechanics and magnetohydrodynamics: Kelvin's law implies the invariance of helicity of the vorticity field for an ideal fluid motion.

## 4.3 Relative helicity for vector fields

First we recall the definition of relative helicity for vector fields (elaborating the definitions in [5] and [3]). Suppose that a domain in the space R 3 (or a closed oriented manifold M 3 ) is split into two regions A and B separated by a boundary surface S . Assume further that two divergence-free vector fields ξ and η in A coincide on the boundary S and have the same extension ζ into the region B . Call the extended fields in M respectively ˜ and ˜. ξ η Abusing notation we will denote them as the sums ˜ = ξ ξ + ζ and ˜ = η η + , where ζ ξ , η and ζ are regarded as the (discontinuous) vector fields in the entire manifold M with supports supp ξ, supp η ⊆ A and supp ζ ⊆ B . (Alternatively one can modify the fields ξ , η , ζ in a narrow neighborhood of the boundary S to avoid discontinuity; this will not affect the argument below.)

Assume that both the extended fields ˜ and ˜ are null-homologous in ξ η M and tangent to ∂M , while we do not impose any restrictions on the topology of the sets A B , , or M .

Definition 4.6. The difference ∆ Hel( ξ, η ) := Hel( ˜ ) ξ -Hel(˜) measures the η relative helicity of the fields ξ and η in A . ⋄

Theorem 4.7. The relative helicity ∆Hel( ξ, η ) = Hel( ˜ ) ξ -Hel(˜) η is independent of their common extension ζ in the region B .

Proof. Define the (closed) two-forms α, β, and ω by substituting the vector fields ξ, η, and ζ into the volume form µ on M i µ : ξ = α i µ , η = β , and i ζ µ = ω . Then one has to show that the difference

$$\Delta \, \text{Hel} = \text{Hel} ( \tilde { \xi } ) - \text{Hel} ( \tilde { \eta } ) = \int _ { M } ( \alpha + \omega ) \wedge d ^ { - 1 } ( \alpha + \omega ) - \int _ { M } ( \beta + \omega ) \wedge d ^ { - 1 } ( \beta + \omega )$$

does not depend on ω . Note that the 2-forms α + ω and β + ω are exact in M (because of the assumption that the fields ˜ and ˜ are null-homologous) and taking ξ η d -1 of them makes sense.

Assume first that the regions A and B are both simply-connected. One readily obtains

$$\Delta \text{Hel} = \int _ { M } \alpha \wedge d ^ { - 1 } \alpha - \int _ { M } \beta \wedge d ^ { - 1 } \beta + \int _ { M } ( \alpha - \beta ) \wedge d ^ { - 1 } \omega + \int _ { M } \omega \wedge d ^ { - 1 } ( \alpha - \beta ).$$

Here d -1 applied to a discontinuous 2-form is a continuous 1-form (the 'form-potential'). The terms in ∆ Hel containing ω are ∫ M ( α -β ) ∧ d -1 ω + ∫ M ω ∧ d -1 ( α -β ), and we want to show that their contribution vanishes.

Integrating by parts one of the terms, we come to 2 ∫ M ( α -β ) ∧ d -1 ω , which, in turn, is equal to 2 ∫ A ( α -β ) ∧ d -1 ω , since supp( α -β ) ⊆ A .

On the other hand, in the domain A the 1-form d -1 ω is the differential of a function, d -1 ω = dh . Indeed, it is closed (the differential d d ( -1 ω ) = ω vanishes in A due to the condition on supp ζ = supp ω ⊆ B ), and hence it is exact in the simply connected region A . Hence,

$$2 \int _ { A } ( \alpha - \beta ) \wedge d ^ { - 1 } \omega = 2 \int _ { A } ( \alpha - \beta ) \wedge d h = 2 \int _ { S } h ( \alpha - \beta ) = 0,$$

where the last equality is due to the assumption on the identity of the fields ξ and η on the boundary S . This proves that ∆ Hel is unaffected by the choice of the extension ζ .

Now consider the case of arbitrary domains A and B . In the latter case, the 2-form ω is not exact, but only closed (and such that α + ω and β + ω are exact, i.e. [ α + ] = [ ω β + ] = 0 ω ∈ H 2 ( M )). To emphasize the ambiguity in the choice of ω we represent is as ω = ω 0 + γ , where we fix some 'reference' closed 2-form ω 0 , while an arbitrary 2-form γ is exact, and both ω 0 and γ have their support in B . Now the difference we are studying is

$$\text{Hel} ( \tilde { \xi } ) - \text{Hel} ( \tilde { \eta } ) = \int _ { M } ( \alpha + \omega _ { 0 } + \gamma ) \wedge d ^ { - 1 } ( \alpha + \omega _ { 0 } + \gamma ) - \int _ { M } ( \beta + \omega _ { 0 } + \gamma ) \wedge d ^ { - 1 } ( \beta + \omega _ { 0 } + \gamma ) \,,$$

and we would like to show that it is independent of the choice of an exact 2-form γ . We claim that this is evident once we introduce new forms ¯ = α α + ω 0 and ¯ = β β + ω 0 and rewrite the difference as

$$\stackrel { \sim } { \text{Hel} } ( \tilde { \xi } ) - \text{Hel} ( \tilde { \eta } ) = \int _ { M } ( \bar { \alpha } + \gamma ) \wedge d ^ { - 1 } ( \bar { \alpha } + \gamma ) - \int _ { M } ( \bar { \beta } + \gamma ) \wedge d ^ { - 1 } ( \bar { \beta } + \gamma ) \,,$$

thus mimicking the expression above for simply-connected regions. Indeed, now it boils down to the same computation leading to the form

$$2 \int _ { A } ( \bar { \alpha } - \bar { \beta } ) \wedge d ^ { - 1 } \gamma = 2 \int _ { A } ( \alpha - \beta ) \wedge d ^ { - 1 } \gamma = 0 \,.$$

The latter expression vanishes for the same reason as above: α -β has the support in A and it is exact (since both α + ω and β + ω are exact), while d -1 γ is closed in A . This concludes the proof that relative helicity is well-defined in the general case.

## 5 Pipes, fluxes, and shells

In this section we present the constructions taking us from the context of cubiculated regions and domino tilings to that of vector fields and helicities.

## 5.1 The five pipes construction

In the first construction, for a domino tile d ⊂ R 3 , we construct a specific divergence-free vector field ξ d in d . The vector field ξ d is confined to five narrow flux tubes. The tubes are the boundaries of tubular neighborhoods of smooth curves approximating the polygonal lines. This is therefore the situation discussed in Section 4.2. If the domino is d = [0 2] , × [0 , 1] × [0 , 1], the tubes are thin neighborhoods of the following five oriented polygonal lines, also shown in Figure 4:

$$\L o n o w m g n v e o n n u c o p o y s o n a m n o, a o v o o n w m g n g \\ ( 0, \frac { 1 } { 2 }, \frac { 1 } { 2 } ) \longrightarrow ( 2, \frac { 1 } { 2 }, \frac { 1 } { 2 } ), \\ ( \frac { 1 } { 2 }, 0, \frac { 1 } { 2 } ) \rightarrow ( \frac { 1 } { 2 }, \frac { 1 } { 4 }, \frac { 1 } { 2 } ) \rightarrow ( \frac { 3 } { 2 }, 0, \frac { 1 } { 2 } ), \\ ( \frac { 1 } { 2 }, \frac { 1 } { 2 }, 0 ) \rightarrow ( \frac { 1 } { 2 }, \frac { 1 } { 2 }, \frac { 1 } { 4 } ) \rightarrow ( \frac { 3 } { 2 }, \frac { 1 } { 2 }, \frac { 1 } { 4 } ) \rightarrow ( \frac { 3 } { 2 }, \frac { 1 } { 2 }, 0 ), \\ ( \frac { 1 } { 2 }, 1, \frac { 1 } { 2 } ) \rightarrow ( \frac { 1 } { 2 }, \frac { 3 } { 4 }, \frac { 1 } { 2 } ) \rightarrow ( \frac { 3 } { 2 }, \frac { 3 } { 4 }, \frac { 1 } { 2 } ) \rightarrow ( \frac { 3 } { 2 }, 1, \frac { 1 } { 2 } ), \\ ( \frac { 1 } { 2 }, \frac { 1 } { 2 }, 1 ) \rightarrow ( \frac { 1 } { 2 }, \frac { 1 } { 2 }, \frac { 3 } { 4 } ) \rightarrow ( \frac { 3 } { 2 }, \frac { 1 } { 2 }, \frac { 3 } { 4 } ) \rightarrow ( \frac { 3 } { 2 }, \frac { 1 } { 2 }, 1 ). \\ \text{re is commetric under rotations} \, h w \, \frac { \pi } { 2 } \, \text{around the control} \, \, \text{ }$$

Notice that the figure is symmetric under rotations by π 2 around the central line, and also under reflection on the plane y = 1 2 . The narrow tubes are assumed to be likewise symmetric.

Figure 4: A domino, five polygonal lines inside it and the respective smooth tubes approximating the lines. Tubes are only very partially drawn in order to keep the figure simple. Lines are always oriented from black to white. At the right, the central square separating the two unit cubes of the domino, crossed by five lines and five tubes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-1165e44a-ad6a-4b74-ae57-7d1699cc287a.jpg)

For other dominoes, the figure is appropriately rotated and translated, while the arrows always point from a black cube to a white one. Thus, pipes enter the domino through a neighborhood of the center of a face of the black cube and exit the domino near the center of a face of the white cube.

If we have a tiling t of a region R , the same construction is performed in each domino. This produces a vector field ξ t in R . Assume that the flux through each pipe is equal to the same φ &gt; 0. We may assume that ξ t is smooth in the interior of R with the support in the union of the pipes. Notice that the restriction of ξ t to a thin neighborhood of the boundary ∂ R does not depend on the choice of the tiling t .

## 5.2 Proof of theorem on relative flux

Here we prove Theorem 1.2 that the relative flux RFlux( ) t ∈ H 1 ( R , ∂ R ) of a tiling t coincides modulo a factor with the relative rotation class [ ξ t ] the vector field ξ t obtained via the 5-pipe construction: [ ξ t ] = 6 φ RFlux( ). t

Proof of Theorem 1.2. Adjust the field ξ t as follows. Add to it one more pipe in the shape of a long torus, see Figure 5, which will follow the black-to-white direction as on Figure 4 and then returns after a U-turn, without self-linking or linking with any other pipe.

With this addition the new field ¯ ξ t in the tile: (1) is still divergence-free without singular points, (2) has 6 pipes in the black-to-white direction and 1 tube in the opposite direction, (3) has

Figure 5: An additional narrow toroidal pipe with a field flux φ and directed along the axis of the tile.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-3586b1b0-b4b3-4a92-8b89-06c577f2fb43.jpg)

one pipe from each face towards the vicinity of the center of the black cube and one pipe from (the vicinity of) the center of the white cube toward its faces, (4) the flux in each narrow pipe is φ .

With this construction, mimicking the definition of the chain 6 t -q 1 , one achieves that the intersection number of the latter (multiplied by φ ) with any closed 2D surface in R transversal to this chain is equal to the flux of the corresponding vector field ¯ ξ t through such a surface. Hence [ ¯ ] realizes the same homology class in ξ t H 1 ( R , ∂ R ) as φ [6 t -q 1 ].

Remark 5.1. While this adjusted '6-pipe' field ¯ ξ t more adequately represents the 6 t -q 1 , the contribution of the extra 'long solid torus' pipe vanishes to both the flux and the helicity, discussed below, so one can confine to the 5-pipe construction. However, it explains the origin of the factor 6 in the formulas. ⋄

Remark 5.2. The chain q 1 , as well as the part of the field ¯ ξ t defined by it, is determined solely by the coloring of the region R . It defines a net of pipes with sources at each white center and sinks at each black center, 'following' the chain q . Then the tiling t can be regarded as an additional system of pipes, 6 times stronger, joining the neighboring centers according to the tiling geometry and satisfying the Kirchhoff junction rule.

Note that the same claim about proportionality of the relative flux RFlux( t ) ∈ H 1 ( R , ∂ R ) and the relative rotation class [ ξ t ] of the corresponding field holds mutatis mutandis for tilings in any dimension. ⋄

## 5.3 An isolating shell

Consider first the case where R ⊂ R 3 is a well behaved region, perhaps a box. We apply the construction discussed in Section 4.3. Let M = R 3 , A = R (possibly rounded at the corners) and B = M ∖ A . Choose a divergence-free vector field ζ in B coinciding with ξ t near ∂ R . For simplicity, we may choose ζ to be also confined to a few tubes. We call ζ or, more precisely, ( M,A,B,ζ ), an isolating shell . Given a tiling t , let ˜ ξ t be the corresponding smooth, divergence-free vector field in M . Since M is contractible, the extended vector field ˜ ξ t is trivially null-homologous.

The more general case is similar, with a few adjustments. If ∂ R = ∅ , nothing needs to be done. Otherwise, recall the assumption that ∂ R is a topological manifold. We can therefore construct an open manifold M 0 with A = R ⊂ M 0 by taking M 0 to be the union of A with a tubular neighborhood of ∂ R . A minor difficulty is that M 0 ∖ A may well be disconnected. If this happens, we construct M ⊇ M 0 by adding tubes (disjoint from R ) connecting the components of M 0 ∖ A . We thus have an open manifold M A , = R⊂ M and a connected subset B = M ∖ A . As above, the vector fields ξ t for all tilings t coincide in a neighborhood of ∂ R ⊂ M . Again, choose a divergence-free vector field ζ in B , also confined to a few tubes and coinciding with ξ t near ∂ R . The desired isolating shell is ( M,A,B,ζ ). We need to check whether the extended vector field ˜ ξ t = ξ t + ζ is null-homologous in M .

Lemma 5.3. Consider a cubiculated region R and a domino tiling t . For M ⊃ R as above, the vector field ζ can be chosen so that ˜ ξ t = ξ t + ζ is null-homologous in M if and only if RFlux( ) = 0 t .

Proof. In this proof, we interpret the vector fields as thin tubular neighborhoods of weighted oriented curves. Thus, we can obtain a class in H M 1 ( ) (with real coefficients) from a vector field.

Assume first that ζ has been chosen so that ˜ ξ t is null-homologous in M . Inclusion defines a map i : H M 1 ( ) → H M,B 1 ( ) = H 1 ( R , ∂ R ). Since [ ˜ ] = 0 we have ξ t i ([ ˜ ξ t ]) = 0 ∈ H 1 ( R , ∂ R ). But i ([ ˜ ξ t ]) = [ ξ t ] = 6 φ RFlux( ) (from Theorem 1.2), proving one implication. t

Conversely, assume that RFlux( ) = 0. Construct a weighted oriented surface t S in R such that ∂S coincides with t -q 1 6 in the interior of R . Use the part of ∂S on ∂ R as guides to construct the required pipes in the intersection of B with a thin tubular neighborhood of ∂ R , i.e., to construct ζ . The surface S is a witness to the fact that ˜ ξ t is null-homologous in M .

## 6 Key examples

In this section we present the first examples of isolating shells and of the five pipes construction. The results of these examples will be used in later proofs.

Example 6.1. Let R = [0 2] , × [0 , 2] × [0 , 1], the 2 × 2 × 1 box: R admits exactly two tilings, t 0 and t 1 , shown in Figure 6.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-92fded31-2147-4d27-9495-2f9f141dd803.jpg)

Figure 6: The only two tilings of the 2 × 2 × 1 box are joined by a flip.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-d3ee6b20-6309-4f97-861f-cb42984131f9.jpg)

Figure 7: A valid isolating shell for the 2 × 2 × 1 box. We show the planes z = -1 4 , 1 4 , 1 2 , 3 4 , 5 4 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-ad25a65b-fb28-4734-bf79-7cbf1c511262.jpg)

Figure 7 shows a valid isolating shell for R . Figure 8 shows the two vector fields ˜ ξ t 0 and ˜ ξ t 1 : both have zero helicity. Indeed, a curve C 0 is trivial (for a family of curves ( C i ) including C 0 ) if there exists a contractible open neighborhood A ⊃ C 0 intersecting no other curve C i and slk ( C 0 ) = 0. In this case, C 0 can be discarded from the family without changing the helicity. In the present example, every curve is trivial. ⋄

Example 6.2. Let R be the region shown in Figure 9, i.e.,

$$\mathcal { R } = [ 0, 2 ] ^ { 3 } \smallsetminus ( ( ( 1, 2 ] \times ( 1, 2 ] \times [ 0, 1 ) ) \cup ( [ 0, 1 ) \times [ 0, 1 ) \times ( 1, 2 ] ) ).$$

It is a 2 × 2 × 2 cube without two small, 1 × 1 × 1 cubes on a diagonal. The region R admits precisely two tilings t 0 and t 1 , also shown in Figure 9. Notice that the two tilings differ by a positive trit from t 0 to t 1 . They also differ by reflection on the plane x = y . In this case there is no natural definition of Tw( t ) as an integer, that is, we need to choose a base tiling in a more or less arbitrary manner. If we choose t 0 as a base tiling we then have Tw( t 0 ) = 0 and Tw( t 1 ) = 1.

Figure 8: The vector fields ˜ ξ t i for the only two tilings of the 2 × 2 × 1 box (with the same planes as Figure 7).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-9c3bfc00-d0e6-4945-a7be-49af8805c936.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-55f7f11e-52c5-4744-896d-6c64b620997d.jpg)

Figure 9: The region R and its two tilings t 0 and t 1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-68ae8543-3bb2-4a40-b4c7-031fb4eb1896.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-c9b2f87b-0925-47a1-a9f0-37889d42abb8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-0e92fc72-17e8-4dda-bc1b-38ec45470e1b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-3fea077d-e571-4bd5-a4e2-4fbe4f3abb1a.jpg)

Figure 10 shows a valid isolating shell for R and the two vector fields ˜ ξ t 0 and ˜ ξ t 1 . Notice that ˜ ξ t 1 is the mirror image of ˜ ξ t 0 under reflection on the plane x = y . This implies Hel( ˜ ξ t 1 ) = -Hel( ˜ ξ t 0 ).

̸

Figure 11 shows the several curves describing ˜ ξ t 0 . There are three nontrivial curves C i , 1 ≤ i ≤ 3 (one per row in Figure 11) and several trivial curves. An explicit computation shows that lk ( C , C i j ) = -2 for all i = j and slk ( C i ) = -2 for all i . Figure 12 shows a sample computation that lk ( C , C 2 3 ) = -2. Indeed, we draw projections of these two curves onto (a small perturbation of) the horizontal xy plane (red for C 2 , blue for C 3 ). This perturbed plane is chosen so that the curves' projections on it intersect transversally. There are eight intersection points between the projections of different curves (indicated in green). Computing the signs, we verify that six of them are negative and two are positive, yielding lk ( C , C 2 3 ) = -2. Similarly, Figure 13 shows that slk ( C 2 ) = -2. Here we project C 2 (red) and a satellite curve C sat 2 (blue). There are six intersection points between different curves (indicated in green, with larger disks indication pairs of twin intersection points). There are five negative and one positive intersection and therefore slk ( C 2 ) = lk ( C , C 2 sat 2 ) = -2. Other cases are treated similarly.

̸

It is convenient to keep this data in the 3 × 3 tabulation matrix L whose entries are l ij = lk ( C , C i j ) for i = j and l ij = slk ( C i ) for i = : j

$$L = \begin{pmatrix} - 2 & - 2 & - 2 \\ - 2 & - 2 & - 2 \\ - 2 & - 2 & - 2 \end{pmatrix}.$$

Proposition 4.3 implies Hel( ˜ ξ t 0 ) = -18 φ 2 and Hel( ˜ ξ t 1 ) = 18 φ 2 . We then have (for all t ∈ T ( R ))

$$\text{Hel} ( \tilde { \xi } _ { t } ) = 3 6 \varphi ^ { 2 } \, \text{Tw} ( t ) + C, \quad C = - 1 8 \varphi ^ { 2 }.$$

The value of the constant C depends on two arbitrary choices made above: the choice of base tiling and the choice of isolating shell. ⋄

Remark 6.3. Strictly speaking, the region R in Example 6.2 must be massaged to make its boundary ∂ R a manifold, since the point (1 1 1) , , is singular. This can be easily achieved by considering the boundary ∂ R ϵ of a small ϵ -neighborhood R ϵ of R , as this will not affect the computations of Hel and Tw. Another solution is to include in R a domino covering either of the missing cubes, such as [1 , 3] × [1 , 2] × [0 , 1]. The isolating shell needs to be slightly modified but we arrive at the same conclusions. ⋄

Figure 10: A valid isolating shell and the vector fields ˜ ξ t 0 and ˜ ξ t 1 . We show the following planes: -1 4 , 1 4 , 1 2 , 3 4 , 5 4 , 3 2 , 7 4 , 9 4 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-d984212d-9c87-45af-94a1-a1b2648c3514.jpg)

## 7 Proof of the main theorem on twist and helicity

Our aim is to prove Theorem 1.1, i.e., to relate Hel( ˜ ) with the twist Tw( ξ t t ) for a tiling of zero relative flux and satisfying the properties of Theorem 3.2. We prove two lemmas in this direction.

Lemma 7.1. Let R be a cubiculated region with a fixed isolating shell. Let t 0 , t 1 be domino tilings of R of zero relative flux and let ξ t i be the corresponding divergence-free vector fields.

- 1. If t 0 and t 1 differ by a flip then Hel( ˜ ξ t 1 ) = Hel( ˜ ξ t 0 ) .
- 2. If there is a positive trit from t 0 to t 1 then Hel( ˜ ξ t 1 ) = Hel( ˜ ξ t 0 ) + 36 φ 2 .

Proof. To see the invariance under the flip (item 1) we start with A being the union of the four unit cubes involved in the flip (as in Example 6.1) and B = M ∖ A its any isolating shell. The fields ξ t i are defined and different in A , but continued into B as the same field ζ . One can see that the relative helicities of ˜ ξ t i = ξ t i + ζ coincide, Hel( ˜ ξ t 1 ) = Hel( ˜ ξ t 0 ). (For instance, for the choice of continuation ζ in Figure 7 both helicities vanish.) For a general region R and a fixed isolating shell B with M := R∪ B we still define A as the four cubes 'participating in the flip' from ξ t 0 to ξ t 1 , while regard its complement B ′ := M ∖ A as a new isolating shell. The same example shows the invariance of relative helicity, Hel( ˜ ξ t 1 ) = Hel( ˜ ξ t 0 ).

The change of relative helicity under the trit (item 2) we define A to be the union of the six unit cubes involved in the trit. The rest of the proof is similar, except that we now use Example 6.2.

Our next lemma considers refinements. We must first define how to refine an isolating shell. Consider a cubiculated region R , a tiling t and a valid isolating shell ( M,A,B,ζ ). The boundary ∂ R is quadriculated: in the middle of each square, ζ draws a pipe with flux ϕ , pointing in or out according to color. That pipe is connected to somewhere in ∂ R to another square of the opposite color.

When we refine R to obtain R ′ , each old large square of ∂ R is decomposed into 25 new small squares of ∂ R ′ . Assume without loss of generality that the old pipe in ζ is now in the central small square. Construct 12 new short pipes matching the other 24 small squares as in Figure 14. (Note that different choices of the horizontal direction on a square face of ∂ R define equivalent isolating shells.) This defines ζ ′ and the desired isolating shell for R ′ .

Figure 11: Curves for the vector fields ˜ ξ t 0 from Example 6.2. The first three rows show three nontrivial curves; the last row shows several trivial curves.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-b6d3532d-6473-4149-b2bc-fd7f10fee50e.jpg)

Figure 12: Example of computation of the linking number, between the curves in the second and third rows of Figure 11; see Example 6.2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-97a0de9e-f073-41cf-b51c-fc384a03bf63.jpg)

Lemma 7.2. Let R be a cubiculated region with a fixed isolating shell. Refine R to obtain R ′ and a corresponding isolating shell. Let t be a domino tiling of R and t ′ be its refinement, a tiling of R ′ . Let ξ t and ξ t ′ be the corresponding divergence-free vector fields. We have Hel( ˜ ) = Hel( ˜ ξ t ξ t ′ ) .

Proof. We first perform a few flips on t ′ to obtain a more drawable tiling t ⋆ of R ′ : by Lemma 7.1, Hel( ξ t ⋆ ) = Hel( ξ t ′ ). The tiling t ⋆ is drawn in Figure 15. Notice that each domino of t is tiled by 125 = 5 3 small dominoes of t ⋆ .

We now construct the vector field ξ t ⋆ , that is, the pipe system corresponding to the tiling t ⋆ as in Figure 16.

The old pipes (that is, the pipes which were already in ξ t ) remain essentially as they were (in ξ t ), with no change of linking or self-linking numbers. The new pipes define only trivial links. This implies Hel( ˜ ) = Hel( ˜ ξ t ξ t ⋆ ), completing the proof.

We are ready to state and prove our main result, Theorem 1.1. We restate it in a slightly different way.

Figure 13: Example of computation of the self-linking number for the curve in the second row of Figure 11; see Example 6.2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-df796866-21c9-447d-b4fa-c64a24f15c1c.jpg)

Figure 14: A square in ∂ R is decomposed into 25 small squares in ∂ R ′ . We refine an isolating shell by adding short pipes joining the centers of the unmatched new squares.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-a8171e67-5733-4162-979b-9638d84248cd.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-8e3d1158-fa60-467e-8fce-7255a7443487.jpg)

Theorem 7.3. (=Theorem 1.1 ) ′ Let R be a cubiculated region. Let t 0 be a tiling of R such that RFlux( t 0 ) = 0 ∈ H 1 ( R , ∂ R ) . Construct an isolating shell for t 0 . There exists a contant C ∈ R such that

$$T w ( \mathbf t ) = \frac { 1 } { 3 6 \varphi ^ { 2 } } \, \text{Hel} ( \tilde { \xi } _ { \mathbf t } ) + C$$

for all t ∈ T 0 ( R ) .

Proof. Recall that given an initial tiling t 0 , the set of all refinements of tilings T 0 ( R ) with the same flux is denoted by T 0 ( R ( ∗ ) ). Define Tw : ˜ T 0 ( R ( ∗ ) ) → R by

$$\widetilde { T w } ( t ) = \frac { 1 } { 3 6 \varphi ^ { 2 } } \left ( \text{Hel} ( \tilde { \xi } _ { t } ) - \text{Hel} ( \tilde { \xi } _ { t _ { 0 } } ) \right ).$$

We claim that Tw assumes integer values and is the twist function, as in Theorem-Definition 3.2. ˜ Indeed, properties 1 and 2 follow from items 1 and 2 of Lemma 7.1, while its property 3 follows from Lemma 7.2. The fact that Tw assumes integer values now follows from the connectivity of ˜ graph T 0 ( R ( ∗ ) ) (Theorem 3.1). Uniqueness of Tw up to an additive constant (which also follows from Theorem 3.1) completes the proof.

̸

Remark 7.4. The above theorem also gives an alternative proof of Theorem-Definition 3.2. Indeed, we constructed a twist function Tw assuming values in Z . In [8] also the setting with m = 0 is discussed from a combinatorial viewpoint, where m is related to the relative flux of the initial tiling. One can prove that m = 0 if and only if RFlux( t 0 ) = 0, but we do not discuss this here as the proof uses different (combinatorial) tools. ⋄

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-8cb2a90f-7ae7-4589-9858-771053ccee5d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-a80e891e-54eb-4897-8235-9480ef09bb5c.jpg)

Figure 15: The tiling t ⋆ , a minor modification of the refinement t ′ . We show here two floors of a domino of t . The large domino is drawn horizontally. At the left, we draw the central floor. At the right, any of the other 4 floors.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-a09e00fb-4ac8-4820-918e-15641de82663.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-5ea452a3-d0b4-429b-a7b5-438b5c804cb7.jpg)

Figure 16: The vector field ξ t ⋆ . The position is the same as in Figure 15. The pipes which were already present in ξ t are drawn in red and remain essentially unchanged. The remaining pipes correspond to new curves, all trivial.

## 8 A 3D 'crosses-and-toes' example

In this section we consider yet one more, larger example of tiling, which is of interest by itself. One can use this example and a slightly different construction of the set A new for the proof of the second part of Lemma 7.1.

Example 8.1. Let R = [0 3] , × [0 , 3] × [0 , 2], the 3 × 3 × 2 box, the same shown in Figure 2. A valid isolating shell is shown in Figure 17. Figure 18 shows ˜ (i.e., the tubes) for the fifth tiling in ξ Figure 2, the one with nine vertical dominoes. Even though there are many curves, they are all trivial and helicity is again zero.

Figure 17: A valid isolating shell for the 3 × 3 × 2 box. We show the following planes: -1 4 , 1 4 , 1 2 , 3 4 , 5 4 , 3 2 , 7 4 and 9 4 . The red dots are vertical tubes. There is yet another central vertical tube, not entirely shown.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-b3bccbc1-9083-4122-b2fe-6a7abf4050b1.jpg)

Figure 18: The five-tube version of the vertical tiling of the 3 × 3 × 2 box (Example 8.1).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-57592719-3256-4df0-bdad-f10b1731d022.jpg)

Let t 0 be the first tiling in Figure 2, the one with nonzero twist, so that Tw( t 0 ) = -1. Figure 19

shows t 0 again; on the second row, the tubes for ξ t 0 . The other rows show one curve at a time: call them C 1 to C 6 , from top to bottom. There is a central vertical curve C 7 , not shown.

̸

Computing linking and self-linking numbers by using signed intersections of curve projections is not hard, but laborious, and we omit the details for now. The results are shown in a 7 × 7 tabulation matrix L with entries l ij = lk ( C , C i j ) for i = j and l ij = slk ( C i ) for i = : j

$$L = \begin{pmatrix} 0 & 0 & - 1 & - 1 & - 1 & - 1 & - 1 \\ 0 & 0 & 0 & 0 & 0 & 0 & - 1 \\ - 1 & 0 & - 1 & - 1 & - 1 & - 1 \\ - 1 & 0 & - 1 & - 1 & - 1 & - 1 \\ - 1 & 0 & - 1 & - 1 & - 1 & - 1 \\ - 1 & 0 & - 1 & - 1 & - 1 & - 1 \\ - 1 & - 1 & - 1 & - 1 & - 1 & 0 \end{pmatrix} ; \\ \vdots \ r r \, 1 / \vdots \ \partial _ { \rho } \, 2 \, \partial _ { \rho } \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rangle \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rror \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rro \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rq \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rp \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rgo \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rhop \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rphoto \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rho \, 2 \, \rno \, 2 \, \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \, \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rror \, \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 3 \rho \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \lho \rror \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \r 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 0 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \ rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho \, 2 \rho. 2 \$$

Equation (2) then implies Hel( ξ t 0 ) = -36 φ 2 = 36 φ 2 Tw( t 0 ).

Note that our box R has precisely 229 tilings. These are t 0 (with twist -1), the mirror image of t 0 (on the plane z = 1, for instance; this tiling has twist +1) and another 227 tilings of twist 0. In this example, any two tilings of twist 0 can be joined by a finite sequence of flips. It can be verified (and it follows from Theorem 1.1) that

$$\text{Hel} ( \xi _ { t } ) = 3 6 \varphi ^ { 2 } \, T w ( t )$$

for all t ∈ T ( R ). If t 1 is the last tiling in Figure 2 it is easy to verify that Hel( ξ t 1 ) = 0 (all curves are trivial; see Figure 18). ⋄

## 9 Appendix: Helicity via linking and self-linking

Here we recall the derivation of the explicit formula (2):

Proposition 9.1. For the helicity of a divergence-free field ξ confined to several narrow nonintersecting linked flux pipes ∪ T i in a simply-connected manifold M :

$$\text{Hel} ( \xi ) = 2 \sum _ { i < j } l k ( C _ { i }, C _ { j } ) \cdot \text{Flux} _ { i } \cdot \text{Flux} _ { j } + \sum _ { i } s l k ( C _ { i } ) \cdot ( \text{Flux} _ { i } ) ^ { 2 }.$$

Here closed curves C i are cores of the pipes T i and Flux i are the fluxes in those pipes, while lk ( C , C i j ) is the Gauss linking number of the curves. The self-linking number slk ( C i ) is the 'average self-linking' of trajectories of ξ in the given pipe. Namely, if all trajectories of the field ξ in a pipe are closed and have the same pairwise linking (for instance, the Poincar´ map in the pipe e cross-section is a solid rotation by a rational angle), then slk is equal to that pairwise linking number. For an arbitrary field, the self-linking number slk must be replaced by the integral of the asymptotic winding number of the field's trajectories introduced by Fathi [6] and discussed below.

This formula was a folklore statement since Moffatt's discovery in [12] of the formula

$$\text{Hel} ( \xi ) = 2 \text{ } l k ( C _ { 1 }, C _ { 2 } ) \cdot \text{Flux} _ { 1 } \cdot \text{Flux} _ { 2 }$$

for helicity Hel( ξ ) = ∫ M ( ξ, curl -1 ξ ) d x 3 of the field ξ supported in two tubes T 1 ∪ T 2 ⊂ R 3 with all orbits closed, with equal periods, and having no internal twist inside the tubes, and since Arnold's generalization of it to arbitrary divergence-free fields as the asymptotic linking number of its trajectories. It explicitly appeared e.g. in [16] with slightly stronger assumptions and a special interpretation of slk via the Calugareanu invariant, cf. [13].

The derivation of the formula above can be obtained as a following sequence of statements.

Figure 19: On the first row, the tiling t 0 from Example 8.1. On the second, the vector field ξ t 0 , i.e., the five-tube version of t 0 . Curves are shown separately on the other rows.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00522v2-figure-7518740c-279b-479d-96f6-b01fbf652f1d.jpg)

Proof. 1) According to Arnold's theorem [2] the total helicity integral Hel( ξ ) = ∫ M ( ξ, curl -1 ξ ) d x 3 is equal to the average linking number of trajectories,

$$\text{Hel} ( \xi ) = \iint _ { M \times M } l k ( x, y ) \mu _ { x } \mu _ { y },$$

where lk ( x, y ) is the asymptotic linking of the trajectories g s ( x ) and g t ( y ) of the field ξ starting at the points x, y ∈ M , and µ is the volume form on M . For a field supported on the union of pipes ∪ T i it splits into pairwise integrals.

2) For two tubes T i and T j carrying a field ξ , all whose orbits are closed with equal periods and with no internal twist, their cross-linking is given by Moffatt's formula. Note that this also holds in the more general setting of not necessarily closed trajectories and for arbitrary internal twist inside the tubes. Indeed, by taking the average linking of one closed curve with trajectories in the other tube (the case studied in [2]), the corresponding linking is proportional to the flux through a Seifert surface for that curve and it does not depend in the internal twist in the tube. Hence in the formula for cross-linking two tubes T i and T j , only the mutual linking of their cores C i with C j and the product of their fluxes Flux i · Flux j appears.

3) What remains to prove is that the helicity for a field inside the i th tube is given by slk ( C i ) · (Flux ) i 2 . This is essentially the statement from [9], see also [17].

In more detail, use some Riemannian metric on M to introduce trivialization along the tube, which now can be regarded as a solid torus D 2 × S 1 with globally defined cylindrical coordinates ( r, ϕ, z ), where r ≥ 0 and z mod1. (The invariants introduced will not depend on that choice.) Following [6] for any pair of points x, y ∈ D one considers the asymptotic winding number

$$W _ { \phi } ( x, y ) = \lim _ { T, S \rightarrow \infty } \frac { \phi ( g ^ { t } ( y ) ) - \phi ( g ^ { s } ( x ) ) } { 2 \pi \cdot T \cdot S }$$

for a pair of trajectories g t ( y , t ) ∈ [0 , T ] and g s ( x , s ) ∈ [0 , S ]. It is defined for almost all pairs x, y ∈ D and it is an integrable function on D × D . It was proved in [9, 17] that its integral over D × D is equal to the helicity invariant of the field. It is also the Calabi invariant of the corresponding Poincar´ map Ψ : e D → D for the flow of ξ in the tube:

$$\text{Hel} ( \xi ) = \text{Cal} ( \Psi ) = \iint _ { D \times D } W _ { \phi } ( x, y ) \, \omega _ { x } \wedge \omega _ { y }$$

where ω is the invariant area form induced by the divergence-free vector field ξ in the pipe crosssection, ω x = i ξ µ x . Note that for a field ξ in a tube along C with closed trajectories of period 1 the value of W x,y ϕ ( ) is constant and W x,y ϕ ( ) = slk ( C ). On the other hand, the integral of the induced area form on the cross-section is equal to Flux ξ of the corresponding vector field ξ through that cross-section.

Remark 9.2. Note that the above general formula is, in particular, valid for the case of closed curves with different periods . For instance, rescale a divergence-free field ξ as fξ for a function f on a tube T . The divergence-free constraint on fξ implies that f is a first integral of the field ξ , while the flux Flus fξ of the new field must be the same in all cross-sections. Then Hel( fξ ) = ( ∫ D f ( x ω ) x ) 2 · Hel( ξ ), see e.g. [9], while the latter expression transforms as follows:

$$\text{Hel} ( f \xi ) = ( \int _ { D } f ( x ) \, \omega _ { x } ) ^ { 2 } \cdot s l k ( C ) \cdot ( \text{Flux} _ { \xi } ) ^ { 2 } = s l k ( C ) \cdot ( \text{Flux} _ { f \xi } ) ^ { 2 } \,,$$

which emphasizes universality of the formula (2) via fluxes.

## References

- [1] V.I. Arnold: On the one-dimensional cohomology of the Lie algebra of divergence-free vector fields and rotation numbers of dynamical systems, Funct. Anal. Appl. , 3:4 (1969), 319-321.
- [2] V.I Arnold: The asymptotic Hopf invariant and its applications, Proc. Summer School in Diff. Equations at Dilizhan (Erevan) , 1973, English translation in: Selecta Math. Soviet. , 5:4 (1986), 327345.
- [3] V. Arnold and B. Khesin: Topological Methods in Hydrodynamics, Applied Mathematical Sciences, vol. 125, Springer-Verlag, New York 1998, pp. xv+374; second extended edition: Springer-Nature Switzerland 2021, xx+456pp.
- [4] G. Bednik: Hopfions in lattice dimer model, Phys. Rev. B
- 100, 024420 (2019), 9pp.; arXiv:1901.04527
- [5] M.A. Berger and G.B. Field: The topological properties of magnetic helicity, J. Fluid Mech. , 147 (1984), 133-148.
- [6] A. Fathi: Transformations et hom´omorphismes pr´ eservant la mesure, e Syst`mes dynamiques minimaux e , Th` ese Orsay (1980)
- [7] M. Freedman, M.B. Hastings, C. Nayak, and X.L. Qi: Weakly coupled non-Abelian anyons in three dimensions, Phys. Rev. B 84, 245119 (2011), 25pp.; arXiv:1107.2731
- [8] J. Freire, C. J. Klivans, P. H. Milet and N. C. Saldanha: On the connectivity of spaces of threedimensional tilings, Trans. Amer. Math. Soc. , 375 (2022), 1579-1605.
- [9] J.M. Gambaudo and E. Ghys: Enlacements asymptotiques, Topology , 36:6 (1997), 1355-1379.
- [10] P.H. Milet and N.C. Saldanha: Flip invariance for domino tilings of three-dimensional regions with two floors. Discrete &amp; Computational Geometry , 53:4 (2015), 914-940.
- [11] Pedro H Milet and Nicolau C Saldanha. Domino tilings of three-dimensional regions: flips and twists. arXiv:1410.7693 .
- [12] H.K. Moffatt: The degree of knottedness of tangled vortex lines. J. Fluid. Mech. , 35:1 (1969), 117-129.
- [13] H.K. Moffatt and R.L. Ricca: Helicity and the Calugareanu invariant, Proc. R. Soc. Lond. A , 439 (1992), 411-429.
- [14] N.C. Saldanha: Domino tilings of cylinders: the domino group and connected components under flips, Indiana Univ. Math. J. , 71:3 (2022), 965-1002.
- [15] N.C. Saldanha, C. Tomei, M.A. Casarin Jr, and D. Romualdo: Spaces of domino tilings, Discrete &amp; Computational Geometry , 14:1 (1995), 207-233.
- [16] M.W. Scheeler, W.M. van Rees, H. Kedia, D. Kleckner, W.T.M. Irvine: Complete measurement of helicity and its dynamics in vortex tubes, Science , 357 (2017), 487-491.
- [17] E. Shelukhin: 'Enlacements asymptotiques' revisited, Annales math´ ematiques du Qu´bec e , 39 (2015), 205-208.
- [18] W.P. Thurston: Conway's tiling groups, The American Mathematical Monthly , 97:8 (1990), 757-773.
- B.K.: Department of Mathematics, University of Toronto
- 40 St. George Street, Toronto, ON M5S 2E4, Canada

khesin@math.toronto.edu

N.S.: Departamento de Matem´ atica, PUC-Rio

Rua Marquˆ es de S˜o Vicente 225, Rio de Janeiro, RJ 22451-900, Brazil a saldanha@puc-rio.br
Boris Khesin Nicolau C. Saldanha
Relative Helicity and Tiling Twist
We consider domino tilings of 3D cubiculated regions. The tilings have two
invariants, flux and twist, often integer-valued, which are given in purely
combinatorial terms. These invariants allow one to classify the tilings with
respect to certain elementary moves, flips and trits. In this paper we present
a construction associating a divergence-free vector field $\xi_t$ to any domino
tiling $t$, such that the flux of the tiling $t$ can be interpreted as the
(relative) rotation class of the field $\xi_t$, while the twist of $t$ is
proved to be the relative helicity of the field $\xi_t$.