---
abstract: 'Measuring and quantifying non-equilibrium dynamics in active biological systems is a major challenge, because of their intrinsic stochastic nature and the limited number of variables accessible in any real experiment. We investigate what non-equilibrium information can be extracted from non-invasive measurements using a stochastic model of soft elastic networks with a heterogeneous distribution of activities, representing enzymatic force generation. In particular, we use this model to study how the non-equilibrium activity, detected by tracking two probes in the network, scales as a function of the distance between the probes. We quantify the non-equilibrium dynamics through the cycling frequencies, a simple measure of circulating currents in the phase space of the probes. We find that these cycling frequencies exhibit power-law scaling behavior with the distance between probes. In addition, we show that this scaling behavior governs the entropy production rate that can be recovered from the two traced probes. Our results provide insight in to how internal enzymatic driving generates non-equilibrium dynamics on different scales in soft biological assemblies.'
author:
- Federica Mura
- Grzegorz Gradziuk
- 'Chase P. Broedersz'
title: 'Non-equilibrium scaling behaviour in driven soft biological assemblies'
---

[^1]

[^2]

Cells and tissue constitute a class of non-equilibrium many-body systems [@Gnesottoreview; @Fodorreview; @Needleman2017; @Schmidt2010; @Jülicherreview]. Indeed, non-equilibrium activity has been observed in various biological systems, including membranes [@Betz2009; @Turlier2016], chromosomes [@Weber], and the cytoplasm [@Lau2003; @Guo2014; @Fakhri2014]. A distinguishing physical feature of such biological assemblies is that they are driven out of equilibrium collectively by internal enzymatic processes that break detailed balance at the molecular scale. The active nature of living matter on larger scales can be determined non-invasively by observing the steady-state stochastic dynamics of mescoscopic degrees of freedom using time-lapse microscopy experiments: The non-equilibrium dynamics of these systems can manifest as circulating probability currents in a phase space of mesoscopic coordinates [@Battle2016; @Gladrow2016; @Gladrow2017; @Gnesottoreview]. However, it remains unclear how such non-equilibrium measures depend on the spatial scale on which the measurement is performed. This issue is not only of practical relevance in an experiment, it is also of fundamental importance: a theoretical understanding of the spatial scaling behavior of broken detailed balance in internally driven systems may reveal how to extract quantitative information from measurable phase space currents about the active nature of the system.

![Schematic illustrating soft viscoelastic networks with heterogeneous driving for various types of cellular systems. A) chromosome B) red blood cell membrane C) cytoskeletal network with in D-F associated bead-spring models with heterogeneous active driving. The color of the bead indicates the intensity of activity, representing the variance (increasing from blue to red) of the associated active noise process. []{data-label="fig:bead_FDT"}](fig1.png){width="7cm"}

Here we consider a simple, yet general model for an internally driven elastic assembly to study non-equilibrium scaling behavior. This assembly is driven out of equilibrium by heterogeneously distributed stochastic forces, representing internal enzymatic activity (Fig. \[fig:bead\_FDT\]). We quantify the non-equilibrium dynamics of such an assembly by the cycling frequencies associated to steady-state circulating currents in phase space [@Gladrow2016; @Gladrow2017]. To study how broken detailed balance manifests on different scales in a given system, we investigate how the cycling frequency of a pair of tracer probes depends on the spatial distance between these probes. Interestingly, the cycling frequencies in our model exhibit a power-law scaling with the distance between probes with an exponent that depends on the dimensionality of the system. To provide a conceptual understanding of this scaling behavior, we develop an analytical calculation of these exponents. Furthermore, we show that the exponent associated to the power law of the cycling frequencies also underlies the scaling behavior of the entropy production rate that can be recovered from measured trajectories. Therefore, we provide a framework to study the spatial scaling behavior of non-equilibrium measures in soft elastic assemblies.

Our model consists of a $d$-dimensional elastic network of $N$ beads, immersed in a simple Newtonian liquid at temperature $T$ [@Yucht2013; @Broedersz2014; @Mao2018]. We assume a lattice structure where each bead is connected to its nearest neighbours by springs of elastic constant $k$, as illustrated in Fig. \[fig:bead\_FDT\]. For simplicity, we model internal enzymatic activity by a Gaussian white noise with variance $\alpha _i$ at bead $i$. By assuming white noise, we effectively consider the dynamics of biological systems on time scales much longer than the characteristic timescales of the active processes [@Gladrow2016; @MacKintosh2008; @Ruostekoski2003]. Importantly, these activity amplitudes, $\alpha _i\geq0$, are spatially heterogeneous, reflecting a spatial distribution of active processes in the system. These activity amplitudes are drawn independently from a distribution $p_\alpha$ with mean $\bar{\alpha}<\infty$ and standard deviation $\sigma_\alpha<\infty$ for each realization of the system. This description of a heterogeneously driven assembly is similar to bead-spring models in which the beads are coupled to distinct heat baths at different temperatures [@Falasco2015; @Bonetto2004; @Rieder1967].

The temporal evolution of the probability distribution, $p({\mathbf{x}},t)$, of the beads’ displacements ${\mathbf{x}}$, relative to their rest positions, is governed by a Fokker-Planck equation: $$\begin{aligned}
\begin{split}
\frac{\partial p({\mathbf{x}},t)}{\partial t} &= - \nabla \cdot [{\mathbf{A}}{\mathbf{x}}p({\mathbf{x}},t)]+\nabla \cdot {\mathbf{D}}\nabla p({\mathbf{x}},t),\\
&=- \nabla \cdot  {\mathbf{j}}({\mathbf{x}},t) 
\label{eq:Fokker}
\end{split}\end{aligned}$$ where ${\mathbf{j}}({\mathbf{x}},t) ={\mathbf{A}}{\mathbf{x}}p({\mathbf{x}},t) + {\mathbf{D}}\nabla p({\mathbf{x}},t) $ is the probability current. Here, ${\mathbf{A}}$ is the elastic interaction matrix, incorporating all nearest neighbor spring interactions between beads; the mobility matrix is assumed to be diagonal to exclude hydrodynamic interactions between the beads, and is absorbed in ${\mathbf{A}}$. The diffusion matrix, ${\mathbf{D}}$, is diagonal with elements $d_{ij}= \delta _{ij}\frac{k_B(T+ \alpha _i)}{\gamma}$, where $\gamma$ is the damping coefficient describing the viscous interaction between a bead and the immersing liquid. The steady-state dynamics of this active network is described by $$p({\mathbf{x}}) = \frac{1}{\sqrt{(2\pi)^{dN} \det {\mathbf{C}}}} e^{-\frac{1}{2}{\mathbf{x}}^T{\mathbf{C}}^{-1}{\mathbf{x}}},$$ where ${\mathbf{C}}= \langle {\mathbf{x}}\otimes {\mathbf{x}}\rangle $ is the covariance matrix, which can be obtained by solving the Lyapunov equation ${\mathbf{A}}{\mathbf{C}}~+~{\mathbf{C}}{\mathbf{A}}^T~ =~ -2{\mathbf{D}}$ [@Risken]. In the simplest limit, the activities are spatially homogeneous: $\alpha _i= \alpha $ $\forall \, i$, resulting in a dynamics that reduce to that of an effective equilibrium scenario with ${\mathbf{j}}=0$. By contrast, in heterogeneously driven systems with non-identical $\alpha _i$’s, we obtain Non-Equilibrium Steady-State (NESS) dynamics with ${\mathbf{j}}\neq 0$ [@Risken].

If we were able to observe the stochastic motion of all beads in the network, we could directly measure the full probability current ${\mathbf{j}}({\mathbf{x}})$ and extract information about the complete non-equilibrium dynamics of the system. However, in an actual experiment typically only a small subset of the degrees of freedom can be tracked (Fig. \[fig:reduceddfig\]A). What information on the non-equilibrium dynamics of the system can be extracted from such limited observations? To address this question, we investigate a scenario where only a few degrees of freedom are accessible.

We start by reducing our description to the marginal distribution, $p_{\rm r}({\mathbf{x}_{\rm r}})~=~\int dx_{k\not\in [r]} p(x_1,x_2,..,x_{dN})$, of a subset $[\rm r]$ of $n$ tracked degrees of freedom ${\mathbf{x}_{\rm r}}$. By integrating out the subset $[{\rm l}]$ of $m$ unobserved degrees of freedom ${\mathbf{x}_{\rm l}}$ on both sides of Eq.  in the steady-state limit, we obtain (see supplementary material): $$0= - \nabla \cdot[{\mathbf{A}_\textup{eff}}{\mathbf{x}_{\rm r}}p_{\rm r}({\mathbf{x}_{\rm r}})]+\nabla \cdot {\mathbf{D}_{[{\rm r},{\rm r}]}}\nabla p_{\rm r} ({\mathbf{x}_{\rm r}}),
\label{eq:effectivefokker}$$ where the sub-index $[{\rm r},{\rm r}]$ of a matrix indicates the sub-matrix corresponding to the reduced set of observed variables. In addition, we introduce the effective linear interaction (Fig. \[fig:reduceddfig\]B), which can be written as ${\mathbf{A}_\textup{eff}}{\mathbf{x}_{\rm r}}$, with ${\mathbf{A}_\textup{eff}}= {\mathbf{A}_{[{\rm r},{\rm r}]}}+ {\mathbf{A}}_{[{\rm r},{\rm l}]} {\mathbf{C}}_{[{\rm l},{\rm r}]}{\mathbf{C}_{[{\rm r},{\rm r}]}}^{-1} $. Here, ${\mathbf{A}}_{[{\rm r},{\rm l}]}$ and ${\mathbf{C}}_{[{\rm l},{\rm r}]}$ are rectangular matrices of sizes $[n \times m]$ and $[m \times n]$, given by the elements of indices $\rm [r,l]$ of ${\mathbf{A}}$ and $\rm [l,r]$ of ${\mathbf{C}}$, respectively.

![Reduced system of tracked probed. A) Schematic of two fluorescently labelled probe beads in a larger system. B) Elastic force acting on bead $j$ obtained at different time steps of a simulation of the Langevin dynamics of the full system (blue points), and the effective linear force, ${\mathbf{A}_\textup{eff}}{\mathbf{x}_{\rm r}}$, from analytical calculations (light blue plane). C) Probability density (color map) and probability current (white arrows) calculated analytically from the effective 2D system, together with results from simulating the full system in the inset. D) The non-conservative part of the effective force field: $\frac{({\mathbf{A}_\textup{eff}}- {\mathbf{A}_\textup{eff}}^T)}{2}{\mathbf{x}_{\rm r}}$ (black arrows) can contribute to the rotation in phase space in non-equilibrium systems.[]{data-label="fig:reduceddfig"}](fig2.png){width="8"}

Thus, we obtain an effective stationary Fokker-Planck equation for the reduced system (Eq. ). By solving this equation, we obtain the exact steady-state reduced probability distribution $p_{\rm r}({\mathbf{x}_{\rm r}})$ and probability current density: $${\mathbf{j}_{\rm r}}({\mathbf{x}_{\rm r}})={\mathbf{A}_\textup{eff}}{\mathbf{x}_{\rm r}}p_{\rm r}({\mathbf{x}_{\rm r}}) +{\mathbf{D}_{[{\rm r},{\rm r}]}}{\mathbf{C}_{[{\rm r},{\rm r}]}}^{-1}{\mathbf{x}_{\rm r}}p_{\rm r}({\mathbf{x}_{\rm r}}), 
 \label{eq:effectivecurrent}$$ which can, in principle, be measured directly from the trajectories of the observed degrees of freedom at steady state (Fig. \[fig:reduceddfig\]C).

We can use this reduced description to investigate how broken detailed balance manifests at different scales in the network. In particular, we consider the simplest case of a reduced system of only two tracked beads in a larger system, as illustrated in Fig. \[fig:reduceddfig\]A. It is convenient to quantify the probability currents in the 2D phase space of these two tracer beads by a pseudoscalar quantity: the average cycling frequency around the origin [@Weiss2003; @Gladrow2016; @Gladrow2017]. For linear systems, we can express the reduced probability current as ${\mathbf{j}_{\rm r}}({\mathbf{x}_{\rm r}})={\mathbf{\Omega}_{\rm r}}{\mathbf{x}_{\rm r}}p_{\rm r}({\mathbf{x}_{\rm r}})$, where ${\mathbf{\Omega}_{\rm r}}$ is a 2D matrix with purely imaginary eigenvalues $\lambda= \pm i \omega$, with $\omega$ representing the cycling frequency.

This cycling frequency can be measured for a pair of probe beads at a distance $r$, and this frequency will depend on the specific configuration of the activity amplitudes $\alpha_i$ at all beads in the system. We aim to compute how this cycling frequency depends on $r$ after averaging over all activity configurations. Since $\omega$ is expected to be distributed symmetrically around $0$, we calculate $\sqrt{\langle\omega^2(r)\rangle_\alpha}$ for pairs of beads separated by a distance $r$. Here, the average $\langle...\rangle_\alpha$ is taken over an ensemble of activities $\{\alpha_i\}$ drawn from the distribution $p_\alpha$. Intuitively, the magnitude of the circulation of currents in phase space typically decreases with the distance between the probes, as shown in Fig. \[fig:freqfig\]A. This reduction of the circulation is reflected by a decrease of the cycling frequency $\omega$ with distance. Remarkably, $\sqrt{\langle\omega^2(r)\rangle_\alpha}$ appears to depend on the distance between the tracer beads, $r$, as a power law, $\sqrt{\langle\omega^2(r)\rangle_\alpha}\propto r^{-\mu}$, with $\mu\approx 1.9$ for a 1D chain with a folded Gaussian or an exponential distribution of activities, as depicted in Fig. \[fig:freqfig\]B.

To investigate how the architecture of the system affects the scaling behavior of the cycling frequencies, we considered different network structures, including square, triangular, and cubic lattices. In particular, we calculated $\sqrt{\langle\omega^2(r)\rangle_\alpha}$, where we also averaged over different lattice directions. Interestingly, we find that the characteristic exponent $\mu$ appears to depend strongly on the dimensionality of the lattice, but not on its geometry, as shown in Fig. \[fig:freqfig\]B-C. These results suggest that the distance dependence of the cycling frequency is determined in part by the long wavelength elastic properties of the system. Importantly, however, the scaling of cycling frequency is sensitive to the spatial structure of the activities. For example, in the simple case of a delta-distributed (single-source) activity on a 1D chain, we find $\mu_{\rm single}\approx 2.4$ (Fig. \[fig:freqfig\]B) in contrast to the value $1.9$ obtained above for spatially distributed activities.

To obtain more insight into the scaling behavior of the cycling frequencies, we derive an analytical expression for the cycling frequency as a function of the distance between the observed beads, $\omega (r)$. In general, it can be shown that for a linear system described by a Fokker-Planck equation, the cycling frequencies are given by (see supplementary materials): $$\omega_{ij}=\frac{1}{2\gamma}\frac{\langle \tau_{ij} \rangle}{\sqrt{\det{\mathbf{C}_{[{\rm r},{\rm r}]}}}}
\label{eq:torque}$$ where $\tau_{ij}\coloneqq {\mathbf{x}}\times{\mathbf{f}_{\rm r}}({\mathbf{x}})= x_if_j({\mathbf{x}})-x_jf_i({\mathbf{x}})$ is a generalized phase space torque in the $x_i$-$x_j$ plane, with $f_i({\mathbf{x}})$ denoting the deterministic force acting on the $i^{\rm th}$ bead. This result is intuitive: for an overdamped system the mean angular velocity is proportional to the mean torque and the factor $1/\sqrt{\det{\mathbf{C}_{[{\rm r},{\rm r}]}}}$ ensures coordinate invariance. For the 1D chain of beads (Fig. \[fig:bead\_FDT\]D), Eq.  reduces to: $$\omega_{ij}=\frac{k}{\gamma}\frac{{\widetilde{\partial}_{2}^2}c_{ij}}{\sqrt{\det{\mathbf{C}_{[{\rm r},{\rm r}]}}}},
\label{eq:omchain}$$ where $c_{ij}$ is the $i,j^{\rm th}$ element of the covariance matrix ${\mathbf{C}}$, and with the discrete second derivative across rows denoted as: ${\widetilde{\partial}_{2}^2}c_{ij}=c_{i,j+1}-2c_{i,j}+c_{i,j-1}$. Thereby, we have reduced the problem of calculating $\omega(r)$ to finding the covariance matrix of the system.

![Spatial scaling behavior of cycling frequencies. A) Steady-state current cycles in phase space of two tracer beads for a nearby pair of probes (left) and distant pair of probes (right). B) Scaling behavior of the cycling frequencies, $\sqrt{\langle\omega^2(r)\rangle}$, of a pair of probes beads as a function of their spatial distances, obtained for a 1D chain and different activity distributions, as indicated in the legend. C) Scaling behavior of the cycling frequencies, $\sqrt{\langle\omega^2(r)\rangle}$, obtained for different lattices and a folded Gaussian activity distribution. Triangular and square markers represent triangular and square/cubic lattices, respectively. Light/dark blue triangles represent triangular networks with zero/finite rest length springs. []{data-label="fig:freqfig"}](fig3.png){width="8cm"}

The structure of ${\mathbf{D}}$ suggests a natural decomposition of the covariance matrix ${\mathbf{C}}$ into equilibrium (${\overline{\mathbf{C}}}$) and non-equilibrium (${\mathbf{C}^{*}}$) contributions: ${\mathbf{C}}=(k_BT/k){\overline{\mathbf{C}}}+(k_B\bar{\alpha}/k){\mathbf{C}^{*}}$, such that ${\overline{\mathbf{C}}}$ and ${\mathbf{C}^{*}}$ are dimensionless. Both ${\overline{\mathbf{C}}}$ and ${\mathbf{C}^{*}}$ can be found by solving the Lyapunov equation, which for the 1D chain is given by $$\begin{aligned}
{\widetilde{\partial}_{1}^2}{\overline{c}}_{ij}+{\widetilde{\partial}_{2}^2}{\overline{c}}_{ij}&=-2\delta_{ij} \label{eq:diffec} \\
{\widetilde{\partial}_{1}^2}{c^{*}}_{ij}+{\widetilde{\partial}_{2}^2}{c^{*}}_{ij}&=-2\delta_{ij}\frac{\alpha_i}{\bar{\alpha}} \label{eq:diffnc},\end{aligned}$$ where ${\widetilde{\partial}_{1}^2}$ indicates the discrete second derivative across columns. These equations represent discrete stationary diffusion equations, with sources of divergence given by $\delta_{ij}$ and $\delta_{ij}(\alpha_i/\bar{\alpha})$, respectively. This result prescribes how a spatial distribution of activities structures the covariance matrix.

We can make further progress by noting that the principle of detailed balance imposes $\omega_{ij}=0$ at thermal equilibrium, which together with Eq.  implies ${\widetilde{\partial}_{2}^2}{\overline{c}}_{ij}=0$. We can, therefore, substitute ${\widetilde{\partial}_{2}^2}c_{ij}$ in Eq. by ${\widetilde{\partial}_{2}^2}{c^{*}}_{ij}$, and then expand this equation up to linear order in $\bar{\alpha}/T$ to obtain $$\omega_{ij}=\frac{k}{\gamma}\frac{\bar{\alpha}}{T}\frac{{\widetilde{\partial}_{2}^2}{c^{*}}_{ij}}{\sqrt{\det{\overline{\mathbf{C}}_{[{\rm r},{\rm r}]}}}}.
\label{eq:omegalin}$$ We proceed by calculating ${\mathbf{C}^{*}}$ for a given distribution of activities $\{\alpha_i\}$. Because of the linearity of Eq. , ${\mathbf{C}^{*}}$ is a superposition of steady-state solutions to single-source problem, i.e. a delta-distribution for which all but one of the activities would be set to zero. Denoting the element of ${\mathbf{C}^{*}}$ at a distance $r$ from the single activity source by ${c^{*}}(r)$, we obtain the “covariance current" $\partial_r {c^{*}}(r)\sim1/r$. Here we employed a continuous approximation of the discrete diffusion problem in Eqs.  and . Thus, ${c^{*}}(r)= -a\ln(r)+b$ for a single-source problem with integration constants $a$ and $b$. Using this expression for ${c^{*}}(r)$ together with Eq. , we obtain for the single source case: $\omega^2_{\rm single}(r)=\frac{k^2}{\gamma^2}\frac{\alpha^2}{T^2}\frac{a^2}{r^4}\frac{1}{\det{\overline{\mathbf{C}}_{[{\rm r},{\rm r}]}}(r)}$, where $\alpha$ is the source’s activity.

Next, we use a superposition of single source solutions for ${c^{*}}(r)$ to obtain the non-equilibrium contribution of the covariance matrix ${\mathbf{C}^{*}}$ for a specific configuration of many activity sources $\{\alpha_i\}$. Using this result in conjunction with Eq.  and performing an ensemble average over the distribution of activity realizations, we arrive at the central result $$\langle\omega^2(r)\rangle_\alpha=\frac{k^2}{\gamma^2}\frac{ \sigma ^2 _{\alpha}}{T^2}\frac{\pi a^2}{2r^3}\frac{1}{\det{\overline{\mathbf{C}}_{[{\rm r},{\rm r}]}}(r)}.
\label{eq:scaling}$$ Finally, we note that the elements of the equilibrium covariance matrix are given by ${\overline{c}}_{i,j}=\textup{min}(i,j) -ij/(N+1)$, and find that for $r\ll N$, $\det{\overline{\mathbf{C}}_{[{\rm r},{\rm r}]}}(r)$ exhibits a power law behavior, $\det{\overline{\mathbf{C}}_{[{\rm r},{\rm r}]}}(r)\sim r$. Therefore, from this analysis we find for a 1D chain with heterogenous activities $\mu=2$, independent of the activity distribution $p_\alpha$. Furthermore, we find $\mu_{\rm single}=2.5$ for a single-source activity, in accord with our numerical result (see Fig.\[fig:freqfig\]B). This calculation provides insight into how a combination of features of the equilibrium and non-equilibrium contributions to the covariance matrix determine the spatial scaling behavior of cycling frequencies.

Non-zero cycling frequencies directly reflect broken detailed balance, suggesting a connection between $\omega$ and measures of the internal driving, including the rate of entropy production. For a Markovian system described by a Fokker-Planck equation, the total entropy production rate under steady-state conditions is given by [@Seifertreview]: $$\Pi_\textup{tot}=k_{\rm B}\int d{\mathbf{x}}\frac{{\mathbf{j}}^T({\mathbf{x}}){\mathbf{D}}^{-1}{\mathbf{j}}({\mathbf{x}})}{p({\mathbf{x}})},
\label{eq:entropyFokker}$$ where $k_{\rm B}$ is Boltzmann’s constant. The validity of this result relies on the equivalence between the Fokker-Planck and Langevin descriptions. However, we have seen that the marginal probability density of the reduced system is only described by a Fokker-Planck equation (see Eq. ) at steady state, reflecting the loss of Markovianity after coarse-graining. Nonetheless, we can define an effective dynamics of the reduced set of variables through the Langevin equation $$\frac{d{\mathbf{x}_{\rm r}}(t)}{dt}=  {\mathbf{A}_\textup{eff}}{\mathbf{x}_{\rm r}}(t) + \sqrt{2{\mathbf{D}_{[{\rm r},{\rm r}]}}}\,\mathbf{\xi}_{\rm r}(t),
\label{eq:effectivelangevin}$$ with Gaussian white noise $\mathbf{\xi}_{\rm r}(t)$. This equation of motion results in the exact steady-state probability and current densities, but with an approximate stochastic dynamics. In particular, the effective interaction matrix ${\mathbf{A}_\textup{eff}}$ (see Eq. ) captures only the average interaction between the traced variables, as illustrated in Fig. \[fig:reduceddfig\]B. Furthermore, in contrast to the full deterministic forces (${\mathbf{A}}{\mathbf{x}}$), these effective interactions (Fig. \[fig:reduceddfig\]C) need not to derive from a potential and, thus, may contain a non-conservative component (Fig. \[fig:reduceddfig\]D).

The entropy production rate associated with the effective Markovian dynamics in Eq.(\[eq:effectivelangevin\]) is given by

![Spatial scaling behavior of the entropy production rate, $\Pi_{\rm r}^{\rm 2D}$, of a pair of probe beads as a function of their spatial distance $r$, obtained for different lattices and a folded Gaussian activity distribution. Note the entropy production rate of the reduced system is scaled by the total entropy production rate of the whole network, $\Pi_{\rm tot}$. Triangular and square markers represent triangular and square/cubic lattices, respectively. Light/dark blue triangles represent triangular networks with zero/finite rest length springs.[]{data-label="fig:entropyfig"}](fig4.png){width="8cm"}

$$\Pi_{{\rm r}}=k_B\int d{\mathbf{x}_{\rm r}}\frac{{\mathbf{j}_{\rm r}}^T({\mathbf{x}_{\rm r}}){\mathbf{D}_{[{\rm r},{\rm r}]}}^{-1}{\mathbf{j}_{\rm r}}({\mathbf{x}_{\rm r}})}{p_{\rm r}({\mathbf{x}_{\rm r}})}\leq \Pi_\textup{tot},
\label{eq:entropyred}$$

where ${\mathbf{j}_{\rm r}}({\mathbf{x}_{\rm r}})$ is defined in Eq. . Note, estimating $\Pi_{\rm r}$ by using the Markovian formalism allows us to set a lower bound for the total entropy production rate $\Pi_{\rm tot}$ (see supplementary materials), similar to what already shown for discrete systems [@Bisker2017]. In the 2D case with two traced degrees of freedom that we consider here, Eq.  reduces to (see supplementary materials) $$\Pi_{\rm r}^{\rm 2D}=k_B\omega ^2\operatorname{Tr}{({\mathbf{C}_{[{\rm r},{\rm r}]}}{\mathbf{D}_{[{\rm r},{\rm r}]}}^{-1})}.
\label{eq:entropy2d}$$ This result provides an explicit relation between the partial entropy production rate and the cycling frequency $\omega$. Note, all quantities in the expression for $\Pi_{\rm r}^{\rm 2D}$ can be observed in an experiment, providing a direct way to non-invasively determine the reduced rate of entropy production for a set of traced degrees of freedom. Since $\operatorname{Tr}{({\mathbf{C}_{[{\rm r},{\rm r}]}}{\mathbf{D}_{[{\rm r},{\rm r}]}}^{-1})}$ depends only weakly on $r$, as long as $1\ll r\ll N$, we expect a scaling behavior $\langle \Pi_{r}^{\rm 2D}\rangle \sim r^{-2 \mu}$. This result shows that the spatial scaling behavior of the cycling frequencies directly determines the spatial scaling behavior of the entropy production rate.

In summary, we here demonstrate theoretically how experimental measures of non-equilibrium activity in internally driven linear networks are affected by the lengthscale at which the system is observed. Specifically, we developed a general framework to predict the scaling behavior of cycling frequencies and the entropy production rate that can be inferred by tracing pairs of degrees of freedom. We showed the exponent $\mu$ that governs this behavior for a system with heterogeneous random activities, is insensitive to the details of distribution of activities. However, this exponent depends sensitively on the dimensionality of the system. The predicted scaling behaviour can be tested in biological [@Betz2009; @Lau2003; @Guo2014; @Mizuno2007; @Fakhri2014; @Lieleg2010; @Koenderink2009; @Schaller2010; @Jensen2015] and artifical [@Palacci2013; @Buttinoni2013] systems under non-equilibrium steady-state conditions.

We thank E. Frey, J. Gladrow, F. Gnesotto, P. Ronceray, and C. Schmidt for many stimulating discussions. This work was supported by the German Excellence Initiative via the program NanoSystems Initiative Munich (NIM), by a DFG Fellowship through the Graduate School of Quantitative Biosciences Munich (QBM). Part of this work was performed at the Aspen Center for Physics, which is supported by National Science Foundation grant PHY-1607611.

[50]{}

Fodor E., Marchetti M. C. ArXiv [**1708**]{} 08652, (2017).

F. Gnesotto, F. Mura, J. Gladrow, C. P. Broedersz Rep. Prog. Phys. (in press) F. Jülicher, K. Kruse, J. Prost, J.-F. Joanny Physics Reports [**449**]{} 1, (2007).

D. Needleman, Z. Dogic Nat. Rev. Mater [**2**]{}, (2017).

C. F. Schmidt, and F. C. MacKintosh Curr. Opin. Cell Biol. 22, 29 (2010).

T. Betz, M. Lenz, J. F. Joanny, C. Sykes, Proc. Natl. Acad. Sci. [**36**]{}, 106 (2009).

H. Turlier, D. A. Fedosov, B. Audoly, T. Auth, N. S. Gov, J. -F. Joanny, G. Gompper, T. Betz Nature Physics [**12**]{}, 513-519 (2016).

S. C. Weber, A. J. Spakowitz, J. A. Theriot, Proc. Natl. Acad. Sci. [**19**]{}, 109 (2010).

A. W. C. Lau, B. D. Hoffman, A. Davies, J. C. Crocker, T. C. Lubensky Phys. Rev. Lett. [**19**]{}, 91 (2003).

M. Guo, A. J. Ehrlicher, M. H. Jensen, M. Renz, J. R. Moore, R. D. Goldman, J. L. Schwartz, F. C. Mackintosh, D. A. Weitz Cell [**4**]{}, 158 (2014).

N. Fakhri, A. D. Wessel, C. Willms, M. Pasquali, D. R. Klopfenstein, F. C. Mackintosh, C. F. Schmidt Science [**6187**]{}, 344 (2014).

C. Battle, C. P. Broedersz, N. Fakhri, J. Howard, C. F. Schmidt, F. C. MacKintosh Science [**6285**]{}, 352 (2016).

J. B. Weiss, Tellus A [**3**]{}, 55 (2003).

J. Gladrow, N. Fakhri, F. C. MacKintosh, C. F. Schmidt, C. P. Broedersz, Phys. Rev. Lett. [**24**]{}, 116 (2016).

J. Gladrow, C. P. Broedersz, C. F. Schmidt Phys. Rev. E [**2**]{}, 96 (2017).

M. G. Yucht, M. Sheinman, C. P. Broedersz, Soft Matter, [**9**]{} (29) (2013).

C. P. Broedersz, F. C. MacKintosh Rev. Mod. Phys, [**86**]{} (3) (2014).

X. Mao, T. C. Lubensky Annu. Rev. Condens. Matter Phys, [**9**]{} (0) (2018).

F. C. MacKintosh, A. J. Levine Phys. Rev. Lett, [**100**]{} (1) 018104 (2008).

J. Ruostekoski, J. R. Anglin Phys. Rev. Lett, [**91**]{} (19) 190402 (2003).

Z. Rieder, J. L. Lebowitz, E.Lieb J. Math. Phys [**8**]{} 1073 (1967).

F. Bonetto, J. L. Lebowitz, J. Lukkarinen J. Stat. Phys [**116**]{} (2004).

G. Falasco, M. Baiesi, L. Molinaro, L. Conti, F. Baldovin Phys. Rev. E [**92**]{} 022129 (2015).

H. Risken, Springer, Berlin, Heidelberg, 1996.

U. Seifert, Reports Prog. Phys. [**12**]{}, 75 (2012).

G. Bisker, M. Polettini, T. R. Gingrich, J. M. Horowitz J. Stat. Mech. Theory Exp [**2017**]{} 093210 (2017).

D. Mizuno, C. Tardin, C. F. Schmidt, F. C. MacKintosh Science [**5810**]{}, 315 (2007).

O. Lieleg, M. M. A. E. Claessens, A. R. Bauch Soft Matter, [**6**]{} (2) (2010).

M. H. Jensen, E. J. Morris, D. A. Weitz Biochim. Biophys, [**1853**]{} (11 0 0 ) (2015).

G. H. Koenderink, Z. Dogic, F. Nakamura, P. M. Bendix, F. C. MacKintosh, j. H. Hartwig, T. P. Stossel, and D. A. Weitz, Proc. Natl. Acad. Sci., [**106**]{}, 36. (2009)

V. Schaller, C. Weber, C. Semmrich, E. Frey, and A. R. Bausch, Nature, 467, 7311 (2010)

J. Palacci, S. Sacanna, A. P. Steinberg, D. J. Pine, P. M. Chaikin, Science, 1230020 (2013).

I. Buttinoni, J. Bialké, F. Kümmel, H. Löwen, C. Bechinger, T. Speck, Phys. Rev. Lett., 110(23), p.238301 (2013).

Y. L. Tong, Springer Series in Statistics, Springer-Verlag, 1990.

Derivation of Eq. (\[eq:effectivefokker\]) {#effectiveequation}
==========================================

Here, we derive Eq. , which describes the steady state distribution of traced variables. Integrating out the unobserved degrees of freedom on both sides of the Fokker-Plank equation (Eq. ), and using the Einstein notation for summing over repeated indexes, we obtain: $$\label{eq:FP3parts}
\overbracket{\int d{\mathbf{x}_{\rm l}}\partial_t p({\mathbf{x}})}^{(I)}=-\overbracket{\int d{\mathbf{x}_{\rm l}}\partial_i [a_{ij}x_j p({\mathbf{x}},t)]}^{(II)}+\overbracket{\int d{\mathbf{x}_{\rm l}}d_{ij}\partial_i\partial_j p({\mathbf{x}},t)}^{(III)}$$ where $a_{ij}$ and $d_{ij}$ are the elements of the interaction matrix ${\mathbf{A}}$ and the diffusion matrix ${\mathbf{D}}$, respectively. Rewriting the probability as $p({\mathbf{x}},t)=p({\mathbf{x}_{\rm l}}|{\mathbf{x}_{\rm r}},t)p_{\rm r}({\mathbf{x}_{\rm r}},t)$, we can separately calculate each term in Eq.$(\ref{eq:FP3parts})$. The first term $(I)$ gives: $$\begin{aligned}
\begin{split}
\int d{\mathbf{x}_{\rm l}}\partial_t p_{\rm r}({\mathbf{x}_{\rm r}},t)p({\mathbf{x}_{\rm l}}|{\mathbf{x}_{\rm r}},t)&=\partial_t p_{\rm r}({\mathbf{x}_{\rm r}},t)\int d{\mathbf{x}_{\rm l}}p({\mathbf{x}_{\rm l}}|{\mathbf{x}_{\rm r}},t)=\partial_t p_{\rm r}({\mathbf{x}_{\rm r}},t) 
\end{split}\end{aligned}$$ For the second term (II), we obtain $$\begin{aligned}
\begin{split}
\label{eq.partII}
\int d{\mathbf{x}_{\rm l}}\partial_i [p_{\rm r}({\mathbf{x}_{\rm r}},t)p({\mathbf{x}_{\rm l}}|{\mathbf{x}_{\rm r}},t)a_{ij}x_j]&=  \delta_{i,[{\rm r}]} \partial_i [p_{\rm r}({\mathbf{x}_{\rm r}},t)\int d{\mathbf{x}_{\rm l}}p({\mathbf{x}_{\rm l}}|{\mathbf{x}_{\rm r}},t)a_{ij}x_j] \\
 &=\delta_{i,[{\rm r}]} \partial_i [ p_{\rm r}({\mathbf{x}_{\rm r}},t)a_{ij} {\left\langle x_j|{\mathbf{x}_{\rm r}},t \right\rangle} ]
\end{split}\end{aligned}$$ where $\delta_{i,[{\rm r}]}=1$ if $x_i$ is one of the observed coordinates and zero otherwise. In the first line we use that the probability density vanishes at infinity faster than $1/x$. Similarly, the third term (III) can be written as $$\begin{aligned}
\begin{split}
\label{eq.partIII}
\int d{\mathbf{x}_{\rm l}}d_{ij}\partial_i\partial_j [p_{\rm r}({\mathbf{x}_{\rm r}},t)p({\mathbf{x}_{\rm l}}|{\mathbf{x}_{\rm r}},t)] &=\delta_{i,[{\rm r}]}\delta_{j,[{\rm r}]} d_{ij}\partial_i\partial_j [p_{\rm r}({\mathbf{x}_{\rm r}},t) \int d{\mathbf{x}_{\rm l}}p({\mathbf{x}_{\rm l}}|{\mathbf{x}_{\rm r}},t)]\\
&=\delta_{i,[{\rm r}]}\delta_{j,[{\rm r}]} d_{ij}\partial_i\partial_j p_{\rm r}({\mathbf{x}_{\rm r}},t) 
\end{split}\end{aligned}$$ An explicit calculation of the conditional averages appearing in Eq.(\[eq.partII\]) yields $ {\left\langle {\mathbf{x}_{\rm l}}|{\mathbf{x}_{\rm r}}\right\rangle} = {\mathbf{C}}_{[{\rm l},{\rm r}]}{\mathbf{C}}_{[{\rm r},{\rm r}]}^{-1} {\mathbf{x}_{\rm r}}$ [@Tong]. We can substitute contributions (I), (II) and (III) in Eq. (\[eq:FP3parts\]) under steady state conditions to obtain Eq. (\[eq:effectivefokker\]).

Derivation of Eq. (\[eq:torque\])
=================================

Here we derive the expression in Eq.  for the cycling frequencies. To this end, we first show that the right hand side of this equation is invariant under orientation preserving linear transformations restricted to the 2-dimensional reduced subspace. Let us consider such a transformation: ${\mathbf{x}_{\rm r}}'={\mathbf{B}}{\mathbf{x}_{\rm r}}$, ${\mathbf{f}_{\rm r}}'={\mathbf{B}}{\mathbf{f}_{\rm r}}$, and denote by ${\mathbf{C}_{[{\rm r},{\rm r}]}}'$ the reduced covariance matrix in the transformed coordinates. $$\label{eq:detB}
{\mathbf{B}}{\mathbf{C}_{[{\rm r},{\rm r}]}}{\mathbf{B}}^T={\mathbf{C}_{[{\rm r},{\rm r}]}}'\quad\Longrightarrow\quad \det{\mathbf{B}}=\sqrt{\frac{\det{\mathbf{C}_{[{\rm r},{\rm r}]}}'}{\det{\mathbf{C}_{[{\rm r},{\rm r}]}}}}$$ Using this result together with the transformation properties of the vector product, we obtain $$\frac{{\left\langle \tau_{ij} \right\rangle}}{\sqrt{\det{\mathbf{C}_{[{\rm r},{\rm r}]}}}}=\frac{{\left\langle {\mathbf{x}_{\rm r}}\times{\mathbf{f}_{\rm r}}({\mathbf{x}}) \right\rangle}}{\sqrt{\det{\mathbf{C}_{[{\rm r},{\rm r}]}}}}=\frac{{\left\langle {\mathbf{x}_{\rm r}}'\times{\mathbf{f}_{\rm r}}'({\mathbf{x}}') \right\rangle}}{\sqrt{\det{\mathbf{C}_{[{\rm r},{\rm r}]}}}}\frac{1}{\det {\mathbf{B}}}=\frac{{\left\langle \tau'_{ij} \right\rangle}}{\sqrt{\det{\mathbf{C}_{[{\rm r},{\rm r}]}}'}}.$$ The coordinate invariance of this term allows us to specifically consider the convenient coordinates in which ${\mathbf{C}_{[{\rm r},{\rm r}]}}={\mathbf{I}}$: $$\label{eq:taucalc}
\frac{1}{\gamma}{\left\langle \tau_{ij} \right\rangle}=\frac{1}{\gamma}{\left\langle {\mathbf{x}_{\rm r}}\times{\mathbf{f}_{\rm r}}({\mathbf{x}})  \right\rangle}=\frac{1}{\gamma}\int d{\mathbf{x}_{\rm r}}{\left\langle {\mathbf{x}_{\rm r}}\times{\mathbf{f}_{\rm r}}({\mathbf{x}})|{\mathbf{x}_{\rm r}}\right\rangle}p_{\rm r}({\mathbf{x}_{\rm r}})=\frac{1}{\gamma}\int d{\mathbf{x}_{\rm r}}\ {\mathbf{x}_{\rm r}}\times{\left\langle {\mathbf{f}_{\rm r}}({\mathbf{x}})|{\mathbf{x}_{\rm r}}\right\rangle}p_{\rm r}({\mathbf{x}_{\rm r}})$$ We can further expand this expression by using ${\mathbf{\Omega}_{\rm r}}={\mathbf{A}_\textup{eff}}+{\mathbf{D}_{[{\rm r},{\rm r}]}}{\mathbf{C}_{[{\rm r},{\rm r}]}}^{-1}$. (The expression for ${\mathbf{\Omega}_{\rm r}}$ follows immediately from Eq. , since we require ${\mathbf{j}_{\rm r}}({\mathbf{x}_{\rm r}})={\mathbf{\Omega}_{\rm r}}{\mathbf{x}_{\rm r}}p_{\rm r}({\mathbf{x}_{\rm r}})$.) $$\label{eq:fcond}
\frac{1}{\gamma}{\left\langle {\mathbf{f}_{\rm r}}({\mathbf{x}})|{\mathbf{x}_{\rm r}}\right\rangle}={\mathbf{A}_\textup{eff}}{\mathbf{x}_{\rm r}}={\mathbf{\Omega}_{\rm r}}{\mathbf{x}_{\rm r}}-{\mathbf{D}_{[{\rm r},{\rm r}]}}{\mathbf{C}_{[{\rm r},{\rm r}]}}^{-1}{\mathbf{x}_{\rm r}}.$$ Combining this result with Eq. , we arrive at $$\label{eq:meantau}
\frac{1}{\gamma}{\left\langle \tau_{ij} \right\rangle}=\int d{\mathbf{x}_{\rm r}}\ {\mathbf{x}_{\rm r}}\times ({\mathbf{\Omega}_{\rm r}}{\mathbf{x}_{\rm r}})p_{\rm r}({\mathbf{x}_{\rm r}})-\int d{\mathbf{x}_{\rm r}}\ {\mathbf{x}_{\rm r}}\times({\mathbf{D}_{[{\rm r},{\rm r}]}}{\mathbf{C}_{[{\rm r},{\rm r}]}}^{-1}{\mathbf{x}_{\rm r}})p_{\rm r}({\mathbf{x}_{\rm r}}).$$ Using the explicit form of ${\mathbf{\Omega}_{\rm r}}$ (see Eq. ), we evaluate the first term in this expression, $$\int d{\mathbf{x}_{\rm r}}\ {\mathbf{x}_{\rm r}}\times ({\mathbf{\Omega}_{\rm r}}{\mathbf{x}_{\rm r}})p_{\rm r}({\mathbf{x}_{\rm r}})=\int d{\mathbf{x}_{\rm r}}\ \omega_{ij}(x_i^2+x_j^2)p_{\rm r}({\mathbf{x}_{\rm r}})=\omega_{ij}(c_{ii}+c_{jj})=2\omega_{ij}.$$ In addition, we confirm by direct calculation, that, as expected, the second term in Eq.  vanishes: $$\begin{aligned}
-\int d{\mathbf{x}_{\rm r}}\ {\mathbf{x}_{\rm r}}\times({\mathbf{D}_{[{\rm r},{\rm r}]}}{\mathbf{x}_{\rm r}})p_{\rm r}({\mathbf{x}_{\rm r}})&=\int d{\mathbf{x}_{\rm r}}\ (-x_j,\ x_i)
\left(\begin{array}{cc}
d_{ii} & d_{ij} \\
d_{ij} & d_{jj} \\
\end{array}\right){\genfrac(){0pt}{0}{x_i}{x_j}}p_{\rm r}({\mathbf{x}_{\rm r}})=\\
&=\int d{\mathbf{x}_{\rm r}}\ [-d_{ii}x_ix_j-d_{ij}x_j^2+d_{ij}x_i^2+d_{jj}x_ix_j]p_{\rm r}({\mathbf{x}_{\rm r}})=\\
&=\underbrace{c_{ij}}_0 (d_{jj}-d_{ii})+d_{ij}\underbrace{(c_{ii}-c_{jj})}_0=0\end{aligned}$$ Altogether, this gives us the desired result: $$\frac{1}{2\gamma}\frac{{\left\langle \tau_{ij} \right\rangle}}{\sqrt{\det{\mathbf{C}_{[{\rm r},{\rm r}]}}}}=\omega_{ij}$$

Derivation of Eq. (\[eq:entropyred\])
=====================================

Here we show that $\Pi _\textup{tot} \ge \Pi_{\rm rr} $. $$\begin{aligned}
\begin{split}
\frac{\Pi _\textup{tot}-\Pi_{\rm r}}{k_{\rm B}}
&=\int d{\mathbf{x}}\frac{{\mathbf{j}}^T({\mathbf{x}}){\mathbf{D}}^{-1}{\mathbf{j}}({\mathbf{x}})}{p({\mathbf{x}})}-\int d{\mathbf{x}_{\rm r}}\frac{{\mathbf{j}_{\rm r}}^T ({\mathbf{x}_{\rm r}}){\mathbf{D}_{[{\rm r},{\rm r}]}}^{-1}{\mathbf{j}_{\rm r}}({\mathbf{x}_{\rm r}})}{p({\mathbf{x}_{\rm r}})}\\ 
&=\frac{\gamma}{k_{\rm B}}\sum_{j\in [l]}\int d{\mathbf{x}}\frac{v_j^2({\mathbf{x}})}{(T+\alpha_j)}p({\mathbf{x}}) + \frac{\gamma}{k_{\rm B}}\sum_{i\in [r]}\ \left[\left(\int d{\mathbf{x}}\frac{v_i^2({\mathbf{x}})}{(T+\alpha_i)}p({\mathbf{x}})\right)- \int d{\mathbf{x}_{\rm r}}\frac{\langle v_i({\mathbf{x}})|{\mathbf{x}_{\rm r}}\rangle^2}{(T+\alpha_i)}p({\mathbf{x}_{\rm r}})\right]   \\ 
&=\frac{\gamma}{k_{\rm B}}\left[\sum_{j\in [l]}\int d{\mathbf{x}}\frac{v_j^2({\mathbf{x}})}{(T+\alpha_j)}p({\mathbf{x}}) +\sum_{i\in [r]}\ \int d{\mathbf{x}_{\rm r}}\left[\left(\int d{\mathbf{x}_{\rm l}}\frac{v_i^2({\mathbf{x}})}{(T+\alpha_i)}p({\mathbf{x}_{\rm l}}|{\mathbf{x}_{\rm r}})p({\mathbf{x}_{\rm r}})\right)-\frac{\langle v_i({\mathbf{x}})|{\mathbf{x}_{\rm r}}\rangle ^2}{(T+\alpha_i)} p({\mathbf{x}_{\rm r}}) \right]\right]\\
&=\frac{\gamma}{k_{\rm B}}\left[\sum_{j\in [l]}\frac{\langle v_j^2({\mathbf{x}})\rangle}{(T+\alpha_j)}+ \sum_{i\in [r]}\ \int d{\mathbf{x}_{\rm r}}\underbrace{\left(\langle v_i^2({\mathbf{x}})|{\mathbf{x}_{\rm r}}\rangle - \langle v_i({\mathbf{x}})|{\mathbf{x}_{\rm r}}\rangle ^2 \right)}_{\geq 0} \frac{p({\mathbf{x}_{\rm r}})}{(T+\alpha_i)}\right]\geq 0
\label{eq.entineq}
\end{split}\end{aligned}$$ where in the second line we use that ${\mathbf{D}}$ is diagonal, $\mathbf{v}({\mathbf{x}})={\mathbf{j}}({\mathbf{x}})/p({\mathbf{x}})$, and ${\mathbf{j}_{\rm r}}({\mathbf{x}_{\rm r}})=p({\mathbf{x}_{\rm r}}) \int d{\mathbf{x}_{\rm l}}\ \mathbf{v}_{\rm r}({\mathbf{x}}) p({\mathbf{x}_{\rm l}}|{\mathbf{x}_{\rm r}})= p({\mathbf{x}_{\rm r}}) {\left\langle \mathbf{v}_{\rm r}({\mathbf{x}})|{\mathbf{x}_{\rm r}}\right\rangle}$, which follows from the derivation of Eq. .

Derivation of Eq. (\[eq:entropy2d\])
====================================

Here we derive the expression for the partial entropy production rate in terms of the cycling frequencies (see Eq.). It is convenient to substitute the current field $\mathbf{j}=\mathbf{\Omega}{\mathbf{x}}p({\mathbf{x}})$ in Eq. , which gives $$\begin{aligned}
\begin{split}
\Pi&=k_{\rm B}\int d{\mathbf{x}}({\mathbf{\Omega}}{\mathbf{x}})^T{\mathbf{D}}^{-1}({\mathbf{\Omega}}{\mathbf{x}})p({\mathbf{x}})= k_{\rm B} \int d{\mathbf{x}}\, x_i \Omega_{ij}^T ({\mathbf{D}}^{-1})_{jl} \Omega_{lm} x_m \\
&=k_{\rm B}\Omega_{ij}^T({\mathbf{D}}^{-1})_{jl}\Omega_{lm} c_{mi}=k_{\rm B} \operatorname{Tr}{({\mathbf{\Omega}}^T {\mathbf{D}}^{-1}{\mathbf{\Omega}}{\mathbf{C}})}.
\end{split}\end{aligned}$$ Since the entropy production is invariant under coordinate transformations, we can use a more suitable coordinate system. In particular, we choose a set of coordinates such that ${\mathbf{C}}=\mathbb{1}$. In this set of coordinates, the entries of the matrix $\Omega _{ij}$ correspond to the cycling frequencies in the coordinates space of the $i^{th}$ and $j^{th}$ coordinates [@Weiss2003]. Thus, in the 2D case ${\mathbf{\Omega}_{\rm r}}$ is given by $$\label{eq:omr}
{\mathbf{\Omega}_{\rm r}}=\begin{pmatrix}
0 &\omega \\
-\omega & 0
\end{pmatrix}$$ Furthermore, in this coordinate system ${\mathbf{C}_{[{\rm r},{\rm r}]}}$ and ${\mathbf{\Omega}_{\rm r}}$ commute, yielding $$\Pi_{\rm r}^{\rm 2D}=k_B\omega ^2\operatorname{Tr}{({\mathbf{C}_{[{\rm r},{\rm r}]}}{\mathbf{D}_{[{\rm r},{\rm r}]}}^{-1})}$$ Note, this expression is invariant under coordinate transformations.

[^1]: These authors contributed equally

[^2]: These authors contributed equally
