---
abstract: 'The spectral eta function for certain families of Dirac operators on noncommutative $3$-torus is considered and the regularity at zero is proved. By using variational techniques, we show that $\eta_{D}(0)$ is a conformal invariant. By studying the Laurent expansion at zero of $\text{TR} (|D|^{-z})$, the conformal invariance of $\zeta''_{|D|}(0)$ for noncommutative $3$-torus is proved. Finally, for the coupled Dirac operator, a local formula for the variation $\partial_A\eta_{D+A}(0)$ is derived which is the analogue of the so called induced Chern-Simons term in quantum field theory literature.'
author:
- |
    $ $\
    Ali Fathi, Masoud Khalkhali
title: '**On Certain Spectral Invariants of Dirac Operators on Noncommutative Tori**'
---

Department of Mathematics, The University of Western Ontario\
London, ON, Canada[^1]\

Introduction
============

In this paper we study the variations of spectral zeta and eta functions $\zeta_{|D|} (z)= \TR(|D|^{-z})$ and $\eta_D (z)=\text{TR} (D|D|^{-z-1}) $ associated to certain families of Dirac operators on noncommutative $3$-torus. In the classical case the canonical trace TR [@Kontsevich-Vishik1995] provides a unified method of studying various spectral functions of elliptic operators and their variations. Connes’ pseudodifferential calculus for noncommutative tori makes it possible to define a suitable notion of noncommutative canonical trace [@det-work], and translate some of the properties of the canonical trace on manifolds to noncommutative settings. Among these, a fundamental result is the explicit description of the Laurent expansion at zero of the function $\text{TR} (AQ^{-z})$ where $A$ and $Q$ are classical elliptic operators [@Paycha-Rosenberg2006]. This result enables us to prove the regularity of $\zeta_{|D|}(z)$ and $\eta_D(z)$ at $z=0$, and also gives a local description for variations of $\eta_D(0)$ and $\zeta'_{|D|}(0)$. In particular, we show that $\eta_D(0)$ is constant over the family $\{e^{th}De^{th}\}$ and hence is a conformal invariant of noncommutative $3$-torus. Also, using the local description for conformal variation of $\zeta'_{|D|}(0)$, we prove that this quantity is a conformal invariant of noncommutative $3$-torus. This paper is organized as follows. In Section 2 we recall the definition of a spectral triple which is the basic ingredient in the definition of a noncommutative Riemannian space [@Connesbook1994]. Our main example is the spin spectral triple for noncommutative tori and its conformal perturbation first proposed in [@Cohen-Connes1992; @Connes-Tred2009]. In Section 3 we give a brief review of Connes’ pseudodifferential calculus for noncommutative tori from [@Connes1980; @Connes-Tred2009], and recall the extension of the Kontsevich-Vishik canonical trace to the setting of noncommutative tori from [@det-work]. It should be mentioned that this is also done in [@Paycha-Levy2014] where one works with toroidal symbols instead of Connes’ symbols.

In Section 4 we study the eta function associated to the Dirac operators of the conformally perturbed spectral triples $(\T^3,\mathcal{H},e^{th}De^{th}),$ and also to the coupled Dirac operator of the spectral triple $(\T^3,\mathcal{H},D+A)$. By exploiting the developed canonical trace, the regularity of the eta function at zero in above cases will be proved. Next, by using variational techniques we show that the value of the eta function at zero is constant over a conformally perturbed family. Also, by considering the spectral triple $(\T^3,\mathcal{H},D$) and the family $D_t=D+tu^*[D,u]$ for a unitary element $u\in\T^3$, we relate the difference $\eta_{D_1}(0)-\eta_{D_0}(0)$ to the spectral flow of the family $D_t$ and give a local formula for index of the operator $PuP$. This is the analogue of the result of Getzler [@Getzler-flow], in the case of noncommutative $3$-torus.

In Section 5 we consider the spectral zeta function $\zeta_{|D|}(z)=\TR(|D|^{-z})$ and study the conformal variation of the spectral value $\zeta'_{|D|}(0)$ within the framework of the canonical trace. We show that for the noncommutative $3$-torus this quantity is a conformal invariant. In even dimensions though, the conformal variation is not zero and hence conformal anomaly exists. Following [@Connes-Moscovici2014] we give a local formula for the conformal variation of $\zeta_{\Delta}'(0)$ in the case of noncommutative two torus.

Finally, in section 6 we consider the coupled Dirac operator $D+A$ and study the value $\zeta'_{D}(0)$ where $\zeta_D(z)=\TR(D^{-z})$. Since the spectrum of $D$ is extended along the real line, there is an ambiguity in the definition of the complex power $D^{-z}$ and hence in the value $\zeta'_D(0)$. In odd dimensions, this ambiguity can be expressed in terms of $\eta_{D+A}(0)$ and hence has a dependence on the coupled gauge field $A$. This dependence can be computed by a local formula and in physics literature it is usually referred to as the induced Chern-Simons term generated by the coupling of a massless fermion to a classical gauge field (cf. e.g. [@QFT-analytic]). We give an analogue of this computation and the local formula in the case of noncommutative $3$- torus.

Conformal and complex geometry of noncommutative two tori were first studied in the seminal work of Connes and Tretkoff [@Connes-Tred2009] where a Gauss-Bonnet theorem was proved for a conformally perturbed metric (cf. [@Cohen-Connes1992] for a preliminary version). This result was extended in [@Masoud-Farzad2012] where the Gauss-Bonnet theorem was proved for metrics in all translation invariant conformal structures. The problem of computing the scalar curvature of the curved noncommutative two torus was fully settled in [@Connes-Moscovici2014], and in [@Farzad-Masoud2013], and in [@Farzad-Masoud2014; @FF1] in the four dimensional case. Other works on the spectral geometry and heat kernel invariants of curved noncommutative tori include [@Marcolli-Tanvir2012; @Ali-Masoud2014; @dabsit1; @dabsit2; @FF1]. The computation of the curvature of the determinant line bundle in the sense of Quillen for certain families of Dirac operators on noncommutative tori was carried out in [@det-work].

We would like to thank Asghar Ghorbanpour for useful discussions on the subject of this paper. MK would like to thank the Hausdorff Institute in Bonn for its hospitality and support while this work was being completed.

Noncommutative geometry framework 
==================================

In this section we recall the basic ingredients in the definition of a noncommutative geometry. Our main example is the spin spectral triple for noncommutative tori on which all the material in this paper is based. The data of a noncommutative Riemannian geometry is encoded in a spectral triple [@Connes1980].

A spectral triple $(\mathcal{A},\mathcal{H},D)$ is given by an involutive unital (possibly noncommutative) algebra $\mathcal{A}$, a representation $\pi:\mathcal{A}\longrightarrow\mathcal{B}(\mathcal{H})$ on a Hilbert space $\mathcal{H}$, and a self-adjoint densely defined operator $D: Dom (D) \subset \mathcal{H} \longrightarrow\mathcal{H}$ with compact resolvent and the property that $[D, \pi(a)]$ is bounded for any $a\in\mathcal{A}$.

A spectral triple $(\mathcal{A},\mathcal{H},D)$ is called [*even*]{} if there exist a self-adjoint unitary operator $\Gamma:\mathcal{H}\to\mathcal{H}$, such that $a\Gamma=\Gamma a,$ for $a\in\mathcal{A}$ and $D\Gamma=-\Gamma D.$ The operator $\Gamma$ induces a grading $\mathcal{H}=\mathcal{H}^+\oplus\mathcal{H}^-$ with respect to which the Dirac operator is odd, $$D=
\begin{bmatrix}
0&D^-\\
D^+&0
\end{bmatrix}.$$

It can be shown that the triple $(C^{\infty}(M), L^2(M,S), \slashed{D})$, consisting of $\mathcal{A}=C^{\infty}(M)$, the algebra of smooth functions on a closed Riemannian spin manifold $(M,g)$, and the spin Dirac operator $\slashed{D}$ on the Hilbert space $\mathcal{H}= L^2(M,S)$ of $ L^2$-spinors satisfies the requirements of a spectral triple (cf. e.g. [@Green-book']). In even dimensions, the spinor bundle admits a grading and we have $S=S^{+}\oplus S^{-}$ with respect to which the Dirac operator is odd. Therefore in even dimensions the spectral triple $(C^{\infty}(M), L^2(M,S), \slashed{D})$ is even.

Geometry of noncommutative tori
-------------------------------

Let $\Theta\in M_n(\mathbb{R})$ be a skew symmetric matrix. The noncommutative $n$-torus $C(\mathbb{T}^n _{\Theta})$ is defined to be the universal $C^{*}$-algebra generated by the unitaries $U_k$ for $k\in\mathbb{Z}^n$ and relations, $$U_kU_l=e^{\pi i\Theta(k,l)}U_{k+l},\quad k,l\in\mathbb{Z}^n.$$ Consider the standard basis $\left\{e_i\right\}$ for $\mathbb{R}^n$ and let $u_i=U_{e_i}$. Then it follows that $$u_ku_l=e^{2\pi i\theta_{kl}}u_lu_k,$$ where $\theta_{kl}=\Theta(e_k,e_l)$. The smooth noncommutative $n$-torus $C^{\infty}(\mathbb{T}^n_{\Theta})$ is defined to be the Fréchet $*-$subalgebra of elements with Schwartz coefficients in the Fourier expansion, that is all the $a\in C^{\infty}(\mathbb{T}^n_{\theta})$ that can be written as $$a=\sum_{k\in\mathbb{Z}^n}a_k U_k,$$ where $\{a_k\}\in\mathcal{S}(\mathbb{Z}^n)$. In fact, $C^{\infty}(\mathbb{T}^n_{\theta})$ is a deformation of $C^{\infty}(\mathbb{T}^n)$ and consists of the smooth vectors under the periodic action of $\mathbb{R}^n$ on $C(\mathbb{T}^n _{\Theta})$ given by $$\alpha_s(U_k)=e^{is.k}U_k,~~s\in\mathbb{R}^n, k\in\mathbb{Z}^n.$$ The algebra $C^{\infty}(\mathbb{T}^n_{\theta})$ is equipped with a tracial state given by $$\tau(\sum_{p\in\mathbb{Z}^n}a_pU_p)=a_0.$$ We also denote by $\delta_{\mu}$ the analogues of the partial derivatives $ \frac{1}{i}\frac{\partial}{\partial x^{\mu}}$ on $C^{\infty}(\mathbb{T}^n)$ which are derivations on the algebra $C^{\infty}(\mathbb{T}^n_{\theta})$ defined by $$\delta_{\mu}(U_k)=k_{\mu} U_k.$$ These derivations have the following property $$\delta_{\mu}(a^*)=-(\delta_{\mu}a)^*,$$ and also satisfy the integration by parts formula $$\tau(a\delta_{\mu}b)=-\tau((\delta_{\mu}a)b),\quad a,b\in C^{\infty}(\mathbb{T}^n_{\theta}).$$

By GNS construction one gets the Hilbert space $\mathcal{H}_{\tau}$ on which $C^{\infty}(\mathbb{T}^n_{\theta})$ is represented by left multiplication which will be denoted by $\pi(a)$ for $a\in C^{\infty}(\mathbb{T}^n_{\theta})$. The spectral triple describing the noncommutative geometry of noncommutative $n$-torus consists of the algebra $C^{\infty}(\mathbb{T}^n_{\theta})$ , the Hilbert space $\mathcal{H}=\mathcal{H}_{\tau}\otimes\mathbb{C}^{N}$, where $N=2^{[n/2]}$ with the inner product on $\mathcal{H}_{\tau}$ given by $$\left<a,b\right>_{\tau}=\tau(b^*a),$$ and the representation of $C^{\infty}(\mathbb{T}^n_{\theta})$ given by $\pi\otimes 1$.

The Dirac operator is $$D=\slashed{\partial}:=\partial_{\mu}\otimes\gamma^{\mu},$$ where $\partial_{\mu}=\delta_{\mu}$, seen as an unbounded self-adjoint operator on $\mathcal{H}_{\tau}$ and $\gamma^{\mu}$s are Clifford (Gamma) matrices in $M_N(\mathbb{C})$ satisfying the relation $$\gamma^i\gamma^j+\gamma^j\gamma^i=2\delta^{ij} I_N.$$ In $3-$dimension the Clifford matrices are given by the Pauli spin matrices, $$\gamma^{1}=
\begin{bmatrix}
0&1\\
1&0
\end{bmatrix} ,
\gamma^{2}=
\begin{bmatrix}
0&-i\\
i&0
\end{bmatrix},
\gamma^{3}=
\begin{bmatrix}
1&0\\
0&-1
\end{bmatrix}.$$

Consider the chirality matrix $$\gamma=(-i)^m\gamma^1\cdots\gamma^n,$$ where $n=2m$ or $2m+1$. It is seen that $\gamma\cdot\gamma=1$, and $\gamma$ anti-commutes with every $\gamma^{\mu}$ for $n$ even. The operator $\Gamma=1\otimes\gamma$ on $\mathcal{H}=\mathcal{H}_{\tau}\otimes\mathbb{C}^{2^{[n/2]}}$ defines a grading on $\mathcal{H}$ and for even $n$ one has $$\Gamma D=-D\Gamma.$$ Therefore $(C^{\infty}(\mathbb{T}^n_{\theta}),\mathcal{H},D)$ is an even spectral triple for $n$ even.

Real structure 
---------------

A real structure on a spectral triple $(\mathcal{A},\mathcal{H},D)$ is an anti-unitary operator $J:\mathcal{H}\to\mathcal{H}$ such $ J (Dom D) \subset Dom D$, and $[ a,  Jb^*J^{-1}]=0$ for all $a, b\in\mathcal{A}$. We say $(\mathcal{A},\mathcal{H},D)$ is of KO-dimension $n$ if the operator $J$ satisfies the following commutation relations $$J^2=\epsilon 1,~JD=\epsilon' DJ,~J\Gamma=\epsilon''\Gamma J,$$ where the $\epsilon,\epsilon',\epsilon''$ depend on $n\in\mathbb{Z}_8$ according to the following table: \[real-table\] $$

  n              0   1    2    3    4    5    6    7
  -------------- --- ---- ---- ---- ---- ---- ---- ---
  $\epsilon$     1   1    -1   -1   -1   -1   1    1
  $\epsilon'$    1   -1   1    1    1    -1   1    1
  $\epsilon''$   1        -1        1         -1   

$$

Note that the real structure operator switches left action to right action. More precisely, Since $\mathcal{A}$ commutes with $J\mathcal{A}J^*$ the real structure gives a representation of the opposite algebra $\mathcal{A}^{op}$ by $b^{\circ}=Jb^*J^*$ and turns the Hilbert space $\mathcal{H}$ into an $\mathcal{A}-$bimodule by $$a\psi b=aJb^*J^*(\psi)~~a,b\in\mathcal{A},~\psi\in\mathcal{H}.$$

\[real-tori\] The spectral triples $(C^{\infty}(\mathbb{T}^n_{\theta}),\mathcal{H}=H_{\tau}\otimes\mathbb{C}^N,D)$ for $N=2^{[n/2]}$ are all equipped with real structure. The real structure operator $J:H_{\tau}\otimes\mathbb{C}^N\rightarrow H_{\tau}\otimes\mathbb{C}^N$ is given by $$J=J_0\otimes C_0,$$ where $J_0$ is the Tomita conjugation map associated to the GNS Hilbert space $H_{\tau}$ given by $$J_0(a)=a^*,~a\in C^{\infty}(\mathbb{T}^n_{\theta})$$ and $C_0$ is the charge conjugation operator on $\mathbb{C}^N$ (cf. e.g. [@Green-book']). For $N=2$ we have $$J=
\begin{bmatrix}
0&-J_0\\
J_0&0
\end{bmatrix}.$$ $J_0^2=1$ for all $n$-mod 8, so by extending the Hilbert space $H_{\tau}$ and coupling $J_0$ with $C_0$ one can check that the resulting operator $J$ satisfies the requirements for a real structure (see [@Green-book']).

The coupled Dirac operator
--------------------------

In this section We consider the spectral triples obtained by coupling (twisting) the Dirac operator by a gauge potential. We begin by recalling the construction in commutative case.

Consider a compact Riemannian spin manifold $(M,g)$ with the spin Dirac operator $D:L^2(S)\rightarrow L^2(S)$ on spinors. Let $V\rightarrow M$ be a Hermitian vector bundle equipped with a compatible connection $\nabla^V$. Then one constructs the coupled (twisted) spinor bundle $S\otimes V\rightarrow M$ with the extended Clifford action $$\pi(\omega).(\psi\otimes f)= (\omega.\psi)\otimes f,$$ where $\pi(\omega)$ is a section of the Clifford bundle over $M$ represented on the space of spinors and $\psi$ and $f$ are sections of spinor bundle and the vector bundle $V$ respectively. Also, the vector bundle $S\otimes V$ is equipped with the coupled (twisted) spin connection $$\nabla^{S\otimes V}=\nabla^S\otimes1+1\otimes\nabla^V,$$ where $\nabla^S$ is the spin connection on $S$. This connection gives rise to the coupled Dirac operator $D_{\nabla^V}$ acting on the sections of $S\otimes V$. When $V$ is a trivial vector bundle, the connection $\nabla^V$ can be globally written as $\nabla^V=d+A$ where $A$ is a matrix of one forms (vector potential) and an easy computation shows that $$D_A=D+\pi(A).$$

The above construction can be generalized to the setting of spectral triples. Starting with a spectral triple $(\mathcal{A},\mathcal{H},D)$, one can construct a new spectral triple $(\mathcal{A},\mathcal{H},D+A)$ by adding a gauge potential to the Dirac operator, this corresponds to picking a trivial projective module over the algebra $\mathcal{A}$.

More precisely, $A$ is a self-adjoint element of $$\Omega^1_{D}=\left\{\sum_j a_j^0[D,a_j^1],~a_j^i\in\mathcal{A}\right\}.$$ In fact, $\Omega^1 _{D}$ is the image inside $\mathcal{B}(\mathcal{H})$ of the noncommutative 1-forms on $\mathcal{A}$ under the induced map $$\pi: a^0da^1\longrightarrow a^0[D,a^1].$$

Below, we explicitly write down the coupled Dirac operator for the spectral triple $(C^{\infty}(\mathbb{T}^n_{\theta}),\mathcal{H},D+A)$. First note that for any element $a=\sum_{k\in\mathbb{Z}^n}a_kU_k$ in $C^{\infty}(\mathbb{T}^n_{\theta})$ we have, $$[D,a]=\partial_{\mu}(a)\otimes\gamma^{\mu}.$$ Any $A\in\Omega^1(C^{\infty}(\mathbb{T}^n_{\theta}))$ is of the form $A=\sum_ia_idb_i$ where $a_i,b_i$ are in $C^{\infty}(\mathbb{T}^n_{\theta})$ and we have $$\pi(A)=\sum_ia_i\partial_{\mu}(b_i)\otimes\gamma^{\mu}.$$ We denote the elements $ \sum_ia_i\partial_{\mu}b_i$ by $A_{\mu}$ and hence, $$\pi(A)=A_{\mu}\otimes\gamma^{\mu},$$ also self adjointness of $A$ gives $A_{\mu}^*=A_{\mu}$. Therefore, the coupled Dirac operator is given by $$D+A=\slashed{\partial}+\slashed{A},$$ where again by Feynman slash notation $\slashed{A}=A_{\mu}\otimes\gamma^{\mu}$.

Elliptic theory on noncommutative tori
======================================

Our aim in this section is to recall the extension of the Kontsevich-Vishik canonical trace to the setting of noncommutative tori from [@det-work]. Alternatively, this is also done in [@Paycha-Levy2014] where they work with toroidal symbols instead of Connes’ symbols. We begin by a brief review of the basics of Connes’ pseudodifferential calculus for noncommutative tori from [@Connes1980; @Connes-Tred2009].

Matrix pseudodifferential calculus on $C^{\infty}(\mathbb{T}^n_{\theta})$
-------------------------------------------------------------------------

We shall use the multi-index notation $\alpha=(\alpha_1,..,\alpha_n)$, $\alpha_i\geq0$ , $|\alpha|=\alpha_1+...+\alpha_n$, $\alpha!=\alpha_1!...\alpha_n!$, $\delta^{\alpha}=\delta_1^{\alpha_1}...\delta_n^{\alpha_n}$ and $\partial_{\xi}^{\beta}=
\partial_{\xi_1}^{\beta_1}...\partial_{\xi_n}^{\beta_n}$.

A matrix valued symbol of order $m$ on noncommutative $n-$torus is a smooth map $$\sigma:\mathbb{R}^n \longrightarrow C^{\infty}(\mathbb{T}^n_{\theta})\otimes M_N(\mathbb{C}),$$ such that $$||\delta^{\alpha}\partial_{\xi}^{\beta}\sigma(\xi)||\leq C_{\alpha, \beta} (1+|\xi|)^{m-|\beta|},$$ and there exists a smooth map $k:\mathbb{R}^n \setminus \{0\}\longrightarrow C^{\infty}(\mathbb{T}^n_{\theta})\otimes M_N(\mathbb{C})$ such that $$\lim_{\lambda\rightarrow\infty}\lambda^{-m}\sigma(\lambda\xi_1,\lambda\xi_2,...,\lambda\xi_n)=k(\xi_1,\xi_2,...,\xi_n).$$ We denote the symbols of order $m$ by $\mathcal{S}^m(C^{\infty}(\mathbb{T}^n_{\theta}))$.

A matrix pseudodifferential operator associated with $\sigma\in\mathcal{S}^m(C^{\infty}(\mathbb{T}^n_{\theta}))$ is the operator $A_{\sigma}:C^{\infty}(\mathbb{T}^n_{\theta})\otimes\mathbb{C}^N
\longrightarrow C^{\infty}(\mathbb{T}^n_{\theta})\otimes\mathbb{C}^N$ defined by $$A_{\sigma}(a)=\int_{\mathbb{R}^n}\int_{\mathbb{R}^n} e^{-is.\xi}\sigma(\xi)\alpha_s(a)dsd\xi,$$ where $\alpha_s$ is the extended action of $\mathbb{R}^n$ on $C^{\infty}(\mathbb{T}^n_{\theta})\otimes\mathbb{C}^N$.

Two symbols $\sigma$, $\sigma'\in {\mathcal S}^m(C^{\infty}(\mathbb{T}^n_{\theta}))$ are considered equivalent if $\sigma-\sigma'\in {\mathcal S}^m(C^{\infty}(\mathbb{T}^n_{\theta}))$ for all $m$. The equivalence of the symbols will be denoted by $\sigma \sim \sigma'$. We denote the collection of pseudodifferential operators by $\Psi^*(C^{\infty}(\mathbb{T}^n_{\theta}))$. The order gives a natural filtration on $\Psi^*(C^{\infty}(\mathbb{T}^n_{\theta}))$ and the following proposition [@Connes-Tred2009] gives an explicit formula for the symbol of the product of pseudodifferential operators as operators on $\mathcal{H}$ modulo the above equivalence relation.

Let $P$ and $Q$ be pseudodifferential operators with the symbols $\sigma$ and $\sigma'$ respectively. The product $PQ$ is a pseudodifferential operator with the following symbol, $$\sigma(PQ)\sim\sum_{\alpha}\frac{1}{\alpha!}\partial_{\xi}^{\alpha}(\sigma(\xi))\delta^{\alpha}(\sigma'(\xi)).$$

A symbol $\sigma\in\mathcal{S}^m(C^{\infty}(\mathbb{T}^n_{\theta}))$ is called elliptic if $\sigma(\xi)$ is invertible for $\xi\neq0,$ and for some $c$ $$||\sigma(\xi)^{-1}||\leq c(1+|\xi|)^{-m},$$ for large enough $|\xi|$.

Consider the Dirac operator of the spectral triple $\left(\T^3,\mathcal{H}_{\tau}\otimes\mathbb{C}^2,D=\slashed{\partial}\right)$, $$\slashed{\partial}=\partial_{\mu}\otimes\gamma^{\mu}:\mathcal{H}_{\tau}\oplus\mathcal{H}_{\tau}\longrightarrow\mathcal{H}_{\tau}\oplus\mathcal{H}_{\tau}.$$ The symbol reads $$\sigma(D)(\xi)=\xi_{\mu}\otimes\gamma^{\mu}=
\begin{bmatrix}
\xi_3&\xi_1-i\xi_2\\
\xi_1+i\xi_2&-\xi_3
\end{bmatrix},$$ which is clearly elliptic.

A smooth map $\sigma:\mathbb{R}^n\to C^{\infty}(\mathbb{T}^n_{\theta})\otimes M_N(\mathbb{C})$ is called a classical symbol of order $\alpha \in \mathbb{C}$ if for any $L$ and each $0\leq j\leq L$ there exist $\sigma_{\alpha-j}:\mathbb{R}^n\backslash\{0\}\to C^{\infty}(\mathbb{T}^n_{\theta})\otimes M_N(\mathbb{C})$ positive homogeneous of degree $\alpha-j$, and a symbol $\sigma^L\in\mathcal{S}^{\Re(\alpha)-L-1}(C^{\infty}(\mathbb{T}^n_{\theta}))$, such that $$\label{sigma}
\sigma (\xi)=\sum_{j=0}^{L}\chi(\xi)\sigma_{\alpha-j}(\xi)+\sigma^L(\xi)\quad\xi\in\mathbb{R}^n.$$ Here $\chi$ is a smooth cut off function on $\mathbb{R}^n$ which is equal to zero on a small ball around the origin, and is equal to one outside the unit ball. The homogeneous terms in the expansion are uniquely determined by $\sigma$. The set of classical symbols of order $\alpha$ on noncommutative n-torus will be denoted by $\mathcal{S}^{\alpha}_{cl}(C^{\infty}(\mathbb{T}^n_{\theta}))$.

The analogue of the Wodzicki residue for classical pseudodifferential operators on the noncommutative $n$-torus is defined in [@Farzad-Wong2011].

The noncommutative residue of a classical pseudodifferential operator $A_{\sigma}$ is defined as $$\Wres(A_{\sigma})=\tau \left(\res (A_\sigma)\right),$$ where $\res(A_\sigma):=\int_{|\xi|=1} \tr\, \sigma_{-n}(\xi)d\xi$.

It is evident from the definition that noncommutative residue vanishes on differential operators, operators of order $<-n$ as well as non-integer order operators

The canonical trace
-------------------

In what follows, we recall the analogue of Kontsevich-Vishik canonical trace [@Kontsevich-Vishik1995] on non-integer order pseudodifferential operators on the noncommutative tori from [@det-work]. For an alternative approach based on toroidal noncommutative symbols see [@Paycha-Levy2014]. For a thorough review of the theory in the classical case we refer to [@Paycha2012; @Paycha-Scott2007].

The existence of the so called cut-off integral for classical noncommutative symbols is established in [@det-work].

\[cut-off\] Let $\sigma\in\mathcal{S}_{cl}^{\alpha}(C^{\infty}(\mathbb{T}^n_{\theta}))$ and $B(R)$ be the ball of radius $R$ around the origin. One has the following asymptotic expansion $$\int_{B(R)}\sigma(\xi)d\xi\sim_{R\rightarrow\infty}\sum_{j=0,\alpha-j+n\neq0}^{\infty}\alpha_j(\sigma)R^{\alpha-j+n}+\beta(\sigma)\log R+c(\sigma),$$ where $\beta(\sigma)=\int_{|\xi|=1}\sigma_{-n}(\xi)d\xi$ and the constant term in the expansion, $c(\sigma)$, is given by $$\int_{\mathbb{R}^n}\sigma^L+\sum_{j=0}^{L}\int_{B(1)}\chi(\xi)\sigma_{\alpha-j}(\xi)d\xi-\sum_{j=0,\alpha-j+n\neq0}^L \frac{1}{\alpha-j+n}\int_{|\xi|=1}\sigma_{\alpha-j}(\omega)d\omega.$$ Here we have used the notation of (\[sigma\]).

The cut-off integral of a symbol $\sigma\in\mathcal{S}_{cl}^{\alpha}(C^{\infty}(\mathbb{T}^n_{\theta}))$ is defined to be the constant term in the above asymptotic expansion, and we denote it by $\Int \sigma(\xi)d\xi$.

Now the canonical trace of a classical pseudodifferential operator of non-integer order on $C^{\infty}(\mathbb{T}^n_{\theta})$ is defined as follows [@det-work]:

The canonical trace of a classical pseudodifferential operator $A $ of non-integral order $\alpha$ is defined as $$\TR(A) =\tau \left(\Int \tr \, \sigma_A(\xi)d\xi\right).$$

The relation between the TR-functional and the usual trace on trace-class pseudodifferential operators is established in [@det-work]. Note that any pseudodifferential operator $A$ of order less that $-n$ is a trace class operator and its trace is given by $$\Tr(A)=\tau \left(\int_{\mathbb{R}^n}\tr\sigma_{P}(\xi)d\xi\right).$$ The symbol of such operators is integrable and we have $$\label{TRTr}
\Int\sigma_A(\xi)=\int_{\mathbb{R}^n}\sigma_A(\xi)d\xi.$$ Therefore, the $\TR$-functional and operator trace coincide on classical pseudodifferential operators of order less than $-n$.

The $\TR$-functional is in fact the analytic continuation of the operator trace and using this fact we can prove that it is actually a trace.

A family of symbols $\sigma(z)\in\mathcal{S}_{cl}^{\alpha(z)}(C^{\infty}(\mathbb{T}^n_{\theta}))$, parametrized by $z\in W\subset \mathbb{C}$, is called a holomorphic family if

-   The map $z\mapsto \alpha(z)$ is holomorphic.

-   The map $z\mapsto \sigma(z)\in \mathcal{S}_{cl}^{\alpha(z)}(\A)$ is a holomorphic map from $W$ to the Fréchet space $\mathcal{S}_{cl}(C^{\infty}(\mathbb{T}^n_{\theta})).$

-   The map $z\mapsto \sigma(z)_{\alpha(z)-j}$ is holomorphic for any $j$, where $$\label{sigma2}
     \sigma(z)(\xi)\sim\sum_j\chi(\xi)\sigma(z)_{\alpha(z)-j}(\xi)\in\mathcal{S}_{cl}^{\alpha(z)}(C^{\infty}(\mathbb{T}^n_{\theta})).$$

-   The bounds of the asymptotic expansion of $\sigma(z)$ are locally uniform with respect to $z$, i.e, for any $L\geq1$ and compact subset $K\subset W$, there exists a constant $C_{L,K,\alpha,\beta}$ such that for all multi-indices $\alpha,\beta$ we have $$\left|\left|\delta^{\alpha}\partial^{\beta}\left(\sigma(z)-\sum_{j<L}\chi\sigma(z)_{\alpha(z)-j}\right)(\xi)\right|\right|<C_{L,K,\alpha,\beta}|\xi|^{\Re(\alpha(z))-L-|\beta|}.$$

A family $\left\{A_z\right\}\in\Psi_{cl}(C^{\infty}(\mathbb{T}^n_{\theta}))$ is called holomorphic if $A_z=A_{\sigma(z)}$ for a holomorphic family of symbols $\left\{\sigma(z)\right\}$.

Complex powers of elliptic operators are an important class of holomorphic families. Let $Q\in\Psi^q_{cl}(C^{\infty}(\mathbb{T}^n_{\theta}))$ be a positive elliptic pseudodifferential operator of order $q>0$. The complex power of such an operator, $Q_\phi^z$, for $\Re(z)<0$ can be defined by the following Cauchy integral formula. $$\label{Qz}
Q_\phi^z=\frac{i}{2\pi}\int_{C_\phi}\lambda_\phi^z (Q-\lambda)^{-1} d\lambda.$$ Here $\lambda^z_\phi$ is the complex power with branch cut $L_\phi =\{re^{i\phi}, r\geq 0\}$ and $C_\phi$ is a contour around the spectrum of $Q$ such that $$C_\phi\cap {\rm spec}(Q)\backslash\{0\}=\emptyset,\qquad L_\phi\cap C_\phi=\emptyset,$$ $$C_\phi\cap\{ {\rm spec}(\sigma(Q)^L(\xi)),\,\xi\neq 0\}=\emptyset.$$

More generally, an operator for which one can find a ray $L_\phi$ with the above property, is called an admissible operator with the spectral cut $L_\phi$ and its complex power can be defined as above. Self-adjoint elliptic operators are admissible (see [@Paycha2012; @det-work]).

The following Proposition is the analogue of the result of Kontsevich and Vishik [@Kontsevich-Vishik1995], for pseudodifferential calculus on noncommutative tori.

\[thmLaurent\] Given a holomorphic family $\sigma(z)\in\mathcal{S}_{cl}^{\alpha(z)}(C^{\infty}(\mathbb{T}^n_{\theta}))$, $z\in W \subset\mathbb{C}$, the map $$z\mapsto \Int\sigma(z)(\xi)d\xi,$$ is meromorphic with at most simple poles located in $$P=\left\{z_0\in W;~\alpha(z_0)\in\mathbb{Z}\cap[-n,+\infty]\right\}.$$ The residues at poles are given by $$\Res_{z=z_0} \Int\sigma(z)(\xi)d\xi=-\frac{1}{\alpha'(z_0)}\int_{|\xi|=1}\sigma(z_0)_{-n}d\xi.$$

By definition, one can write $\sigma(z)=\sum_{j=0}^L\chi(\xi)\sigma(z)_{\alpha(z)-j}(\xi)+\sigma(z)^L(\xi)$, and by Proposition (\[cut-off\]) we have, $$\begin{aligned}
\Int\sigma(z)(\xi)d\xi
&=\int_{\mathbb{R}^n}\sigma(z)^L(\xi)d\xi
+\sum_{j=0}^{L}\int_{B(1)}\chi(\xi)\sigma(z)_{\alpha(z)-j}(\xi)\\
&-\sum_{j=0}^{L}\frac{1}{\alpha(z)+n-j}\int_{|\xi|=1}\sigma(z)_{\alpha(z)-j}(\xi)d\xi.\end{aligned}$$ Now suppose $\alpha(z_0)+n-j_0=0$. By holomorphicity of $\sigma(z)$, we have $\alpha(z)-\alpha(z_0)=\alpha'(z_0)(z-z_0)+o(z-z_0)$. Hence $$\Res_{z=z_0} \Int\sigma(z)=\frac{-1}{\alpha'(z_0)}\int_{|\xi|=1}\sigma(z_0)_{-n}(\xi)d\xi.$$

\[analyticcontin\] The functional $\TR$ is the analytic continuation of the ordinary trace on trace-class pseudodifferential operators.

First observe that, by the above result, for a non-integer order holomorphic family of symbols $\sigma(z)$, the map $z\mapsto\Int\sigma(z)(\xi)d\xi$ is holomorphic. Hence, the map $\sigma\mapsto\Int\sigma(\xi)d\xi$ is the unique analytic continuation of the map $\sigma\mapsto\int_{\mathbb{R}^n}\sigma(\xi)d\xi$ from $\mathcal{S}_{cl}^{<-n}(C^{\infty}(\mathbb{T}^n_{\theta}))$ to $\mathcal{S}^{\notin\mathbb{Z}}_{cl}(C^{\infty}(\mathbb{T}^n_{\theta}))$. By we have the result.

Let $A\in\Psi ^{\alpha}_{cl}(C^{\infty}(\mathbb{T}^n_{\theta}))$ be of order $\alpha\in\mathbb{Z}$ and let $Q$ be a positive elliptic classical pseudodifferential operator of positive order $q$. We have $$\Res_{z=0}\TR(AQ^{-z})= \frac{1}{q}\Wres(A).$$

For the holomorphic family $\sigma(z)=\sigma(AQ^{-z})$, $z=0$ is a pole for the map $z\mapsto\Int\sigma(z)(\xi)d\xi$ whose residue is given by $$\Res_{z=0}\left(z\mapsto \Int\sigma(z)(\xi)d\xi\right)=-\frac{1}{\alpha'(0)}\int_{|\xi|=1}\sigma_{-n}(0)d\xi=-\frac{1}{\alpha'(0)}\res(A).$$ Taking $\tau$-trace on both sides gives the result.

Now we can prove the trace property of $\TR$-functional.

We have $\TR(AB)=\TR(BA)$ for any $A,B\in\Psi^*_{cl}(C^{\infty}(\mathbb{T}^n_{\theta}))$, provided that $ord(A)+ord(B)\notin\mathbb{Z}$.

Consider the families $\{A_z\}$ and $\{B_z\}$ such that $A_0\sim A$, $B_0\sim B$ , $\ord(A_z)=\ord(A)+z$ and $\ord(B_z)=\ord(B)+z$. For $z\in W=-(\ord(A)+\ord(B))+\mathbb{Z}$ the families $\{A_zB_z\}$ and $\{B_zA_z\}$ have non-integer order . For $\Re(z)\ll 0$, the two familes are trace-class and $\Tr(A_zB_z)=\Tr(B_zA_z)$. now by analytic continuation we have $TR(A_zB_z)=TR(B_zA_z)$, for $z\in\mathbb{C}-W$. Putting $z=0$ gives $TR(AB)=TR(BA)$.

The above result provides another proof for the trace property of the non-commutative residue on $\Psi^{\mathbb{Z}}_{cl}(C^{\infty}(\mathbb{T}^2_{\theta}))$ given in [@Farzad-Wong2011], namely, for $A,B\in\Psi^{\mathbb{Z}}_{cl}(C^{\infty}(\mathbb{T}^2_{\theta}))$, $$\Wres([A,B])=0.$$ On can write, $$\Wres([A,B])=\Res_{z=0}\TR([A,B]Q^{-z})=\Res_{z=0}\TR(C_z)+\Res_{z=0}\TR([AQ^{-z},B]),$$ where $C_z=ABQ^{-z}-AQ^{-z}B$. For $Re(z)\gg0$, the operator $AQ^{-z}$ is trace-class and $\Tr([AQ^{-z},B])=0$, so by analytic continuation,$\TR([AQ^{-z},B])=0$ and therefore, $\Res_{z=0}\TR([A,B]Q^{-z})=\Res_{z=0}\TR(C_z)$. Finally, $C_0=ABQ^0-AQ^0B\in\Psi_{cl}^{-\infty}(C^{\infty}(\mathbb{T}^2_{\theta}))$, so $$\Wres([A,B])=\Res_{z=0}\TR(C_z)=\Wres(C_0)=0,$$ where in the last equality we used the fact that the noncommutative residue of a smoothing operator is zero.

Log-polyhomogeneous symbols
---------------------------

In general, $z$-derivatives of a classical holomorphic family of symbols are not classical anymore and therefore we introduce log-polyhomogeneous symbols which include the $z$-derivatives of the symbols of the holomorphic family $\sigma(AQ^{-z})$.

A symbol $\sigma$ is called a log-polyhomogeneous symbol if it has the following form $$\label{logpolysym}
\sigma(\xi) \sim \sum_{j\geq 0}\sum_{l=0}^\infty\sigma_{\alpha-j,l}(\xi)\log^l|\xi|\quad |\xi|>0,$$ with $ \sigma_{\alpha-j,l}$ positively homogeneous in $\xi$ of degree $\alpha - j$.

An important example of an operator with such a symbol is $\log Q$ where $Q\in\Psi^q_{cl}(C^{\infty}(\mathbb{T}^n_{\theta}))$ is a positive elliptic pseudodifferential operator of order $q>0$. The logarithm of $Q$ can be defined by $$\log Q
 =Q\left.\frac{d}{dz}\right|_{z=0}Q^{z-1}
 =Q\left.\frac{d}{dz}\right|_{z=0}\frac{i}{2\pi}\int_C\lambda^{z-1}(Q-\lambda)^{-1}d\lambda.$$ It is a pseudodifferential operator with symbol $$\label{symboflog}\sigma(\log Q)
\sim \sigma(Q)\star \sigma\Big(\left.\frac{d}{dz}\right|_{z=0} Q^{z-1}\Big),$$ where $\star$ denotes the product of symbols. One can show that is a log-polyhomogeneous symbol of the form $$\sigma(\log Q)(\xi)=q\log|\xi| I+ \sigma_{cl}(\log Q)(\xi),$$ where $\sigma_{cl}(\log Q)$ is a classical symbol of order zero (see [@Paycha2012]).

By adapting the proof of Theorem 1.13 in [@Paycha-Scott2007] to the noncommutative case, we have the following theorem for the family $\sigma(AQ^{-z})$.

\[Laurent exp\] Let $A\in\Psi_{cl}^\alpha(C^{\infty}(\mathbb{T}^n_{\theta}))$ and $Q$ be a positive (or more generally, an admissible) elliptic pseudodifferential operator of positive order $q$. If $\alpha\in P$ then $z=0$ is a possible simple pole for the function $z\mapsto \TR(AQ^{-z})$ with the following Laurent expansion around zero, $$\begin{aligned}
\TR(AQ^{-z})&=\frac{1}{q}\Wres(A)\frac{1}{z}\\
 &+\tau\left(\Int\sigma(A)- \frac{1}{q}\res(A\log Q)\right)-\Tr(A\Pi_Q)\\
& +\sum_{k=1}^K(-1)^k\frac{(z)^k}{k!} \\
&\times \left(\tau\left( \Int\sigma(A(\log Q)^k)d\xi-\frac{1}{q(k+1)}\res(A(\log Q)^{k+1})\right)-\Tr(A\log^k Q\Pi_Q)\right)\\
&+o(z^{K}).\end{aligned}$$ Where $\Pi_Q$ is the projection on the kernel of $Q$.

The term $\res(A\log Q)$ appearing in above Laurent expansion is an extension of Wodzicki residue density to operators with Log-polyhomogeneous symbols [@Lesch1999]. For an operator $P$ with log-polyhomogeneous symbol, by $\res(P)$ we mean, $$\res(P)=\int_{|\xi|=1}\sigma_{-n,0}(\xi)d\xi,$$ (see ).

The spectral eta function
=========================

In this section we study the eta function associated with the family of spectral triples $(\T^3,\mathcal{H},e^{th}De^{th})$ where $h\in\T^3$ is a self-adjoint element [@Connes-Moscovici-twisted], and also the coupled spectral triple $(\T^3,\mathcal{H},D+A)$. By exploiting the developed pseudodifferential calculus, the regularity of the eta function at zero in above cases will be proved.

Regularity at zero
------------------

The spectral eta function was first introduced in [@APS] where its value at zero appeared as a correction term in the Atiyah-Patodi-Singer index theorem for manifolds with boundary. It is defined as $$\eta_D(z)=\sum_{\lambda\in\spec(D),\lambda\neq0}\sgn(\lambda)|\lambda|^{-z}=\TR\left(D|D|^{-z-1}\right),$$ where $D$ is a self-adjoint elliptic pseudodifferential operator. Unlike the spectral zeta functions for positive elliptic operators, proving the regularity of eta function at zero is difficult. This was proved in [@APSIII] and [@Gilkey-eta] using K-theoretic arguments and in [@Bismut-Freed], Bismut and Freed gave an analytic proof of the regularity at zero of the eta function for a twisted Dirac operator on an odd dimensional spin manifold.

Note that for an even spectral triple $(\mathcal{A},H,D)$ we have $D\Gamma=-\Gamma D$, therefore the spectrum of the Dirac operator is symmetric and $\eta_D(z)$ is identically zero. Also the same vanishing happens if the spectral triple admits a real structure with KO-dimensions $1$ or $5$ (see definition \[real-table\]).

The meromorphic structure of eta function for Dirac operator can be studied by the pole structure of the canonical trace for holomorphic families.

\[eta-res\] Let $D$ be an elliptic self-adjoint first-order differential operator on $\T^3$. The poles of the eta function $\eta_D(z)$ are located among $\{3-i,~i\in\mathbb{N}\}$, and $$\Res_{z=0}\eta_D(z)=\Wres(D|D|^{-1}).$$

By using the result of Proposition \[thmLaurent\], the family $\left\{\sigma\left(D|D|^{-z-1}\right)\right\}$ has poles within the set $\left\{z;-z\in\mathbb{Z}\cap[-3,\infty]\right\}$ or $\left\{z=3-i,i\in\mathbb{N}\right\}$. Also, by Proposition \[Laurent exp\] we have $$\eta_D(z)=\TR\left(D|D|^{-z-1}\right)=\Wres(D|D|^{-1})\frac{1}{z}+a_0+a_1z+\cdots,$$ Hence, $$\Res_{z=0}\eta_D(z)=\Wres(D|D|^{-1}).$$

We now prove the regularity at $z=0$ of eta function for the 1-parameter family $\{e^{th}De^{th}\}$ for the spectral triple $(\T^3,H_{\tau}\otimes\mathbb{C}^2,D=\partial_{\mu}\otimes\gamma^{\mu})$ on noncommutative 3-torus.

\[eta-reg\] Consider the family of operators $\{e^{th}De^{th}\}$ on $\mathcal{H}_{\tau}\otimes\mathbb{C}^2$ where $h=h^*\in\T^3$, then $$\Res_{z=0} \eta_{e^{th}De^{th}}(z)=0.$$

By definition, $\eta_{D_t}(z)=\TR(D_t|D_t|^{-z-1})$. Using Proposition \[eta-res\] we have $$\Res_{z=0} \eta_{D_t}(z)=\Wres(D_t|D_t|^{-1}),$$ where the right hand side is the Wodzicki residue on noncommutative $3$-torus.

Now for each element of the family $D_t=e^{\frac{th}{2}}De^{\frac{th}{2}}$, $D_t^2=e^{\frac{th}{2}}De^{th}De^{\frac{th}{2}}$ and $|D_t|=\sqrt{D_t^2}$. By using the product formula for the symbols we have, $$\sigma(D_t)\sim\xi_{\mu}e^{th}\otimes\gamma^{\mu}+e^{\frac{th}{2}}\delta_{\mu}(e^{\frac{th}{2}})\otimes\gamma^{\mu},$$ and $$\begin{aligned}
\sigma(D_t^2)&=\sigma(e^{\frac{th}{2}}De^{th}De^{\frac{th}{2}})\sim\xi_{\lambda}\xi_{\mu}e^{2th}\otimes\gamma^{\lambda}\gamma^{\mu}\\
&+\xi_{\lambda}e^{\frac{3th}{2}}\delta_{\mu}(e^{\frac{th}{2}})\gamma^{\lambda}\gamma^{\mu}+\xi_{\mu}e^{\frac{th}{2}}\delta_{\lambda}(e^{\frac{3th}{2}})\otimes\gamma^{\lambda}\gamma^{\mu}\\
&+e^{\frac{th}{2}}\delta_{\lambda}(e^{th})\delta_{\mu}(e^{\frac{th}{2}})\otimes\gamma^{\lambda}\gamma^{\mu}+e^{\frac{3th}{2}}\delta_{\lambda}\delta_{\mu}(e^{\frac{th}{2}})\otimes\gamma^{\lambda}\gamma^{\mu}.\\\end{aligned}$$ To compute the homogeneous terms in the symbol of $|D_t|$, we observe that $|D_t|=\sqrt{D_t^2}$ and hence, $$\sigma(|D_t|)=\sqrt{\sigma(D_t^2)}\sim\sigma_1(\xi)+\sigma_0(\xi)+\sigma_{-1}(\xi)+\cdots.$$ We compute the first three terms, which we need for computing the symbol of $|D_t|^{-1}$. $$\begin{aligned}
\sigma_1(\xi)&=\lim_{k\to\infty}\frac{\sigma(|D_t|)(k\xi)}{k}=\sqrt{\left(\xi^2\right)}e^{th}\otimes I,\\
\sigma_0(\xi)&=\lim_{k\to\infty}\sigma(|D_t|)(k\xi)-k\sigma_1(\xi)\\
&=\left(\xi_{\lambda}e^{\frac{3th}{2}}\delta_{\mu}(e^{\frac{th}{2}})\right)(\frac{1}{2\sqrt{\xi^2}}e^{-th})\otimes\gamma^{\lambda}\gamma^{\mu}+\left(\xi_{\mu}e^{\frac{th}{2}}\delta_{\lambda}(e^{\frac{3th}{2}})\right)(\frac{1}{2\sqrt{\xi^2}}e^{-th})\otimes\gamma^{\lambda}\gamma^{\mu},\\
\sigma_{-1}(\xi)&=\lim_{k\to\infty}\frac{\left(\sigma(|D_t|)(k\xi)-\sigma_1(\xi)-\sigma_0(\xi)\right)}{k^{-1}}\\
&=\frac{1}{2\sqrt{\xi^2}}\left(e^{\frac{th}{2}}\delta_{\lambda}(e^{th})\delta_{\mu}(e^{\frac{th}{2}})\right)e^{-th}\otimes\gamma^{\lambda}\gamma^{\mu}+\frac{1}{2\sqrt{\xi^2}}\left(e^{\frac{3th}{2}}\delta_{\lambda}\delta_{\mu}(e^{\frac{th}{2}})\right)e^{-th}\otimes\gamma^{\lambda}\gamma^{\mu}\\
&-\frac{1}{8\xi^2\sqrt{\xi^2}}\left(\xi_{\lambda}e^{\frac{3th}{2}}\delta_{\mu}(e^{\frac{th}{2}})\right)e^{-th}\left(\xi_{\nu}e^{\frac{3th}{2}}\delta_{\rho}(e^{\frac{th}{2}})\right)e^{-2th}\otimes\gamma^{\lambda}\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\\
&-\frac{1}{8\xi^2\sqrt{\xi^2}}\left(\xi_{\lambda}e^{\frac{3th}{2}}\delta_{\mu}(e^{\frac{th}{2}})\right)e^{-th}\left(\xi_{\rho}e^{\frac{th}{2}}\delta_{\nu}(e^{\frac{3th}{2}})\right)e^{-2th}\otimes\gamma^{\lambda}\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\\
&-\frac{1}{8\xi^2\sqrt{\xi^2}}\left(\xi_{\mu}e^{\frac{th}{2}}\delta_{\lambda}(e^{\frac{3th}{2}})\right)e^{-th}\left(\xi_{\nu}e^{\frac{3th}{2}}\delta_{\rho}(e^{\frac{th}{2}})\right)e^{-2th}\otimes\gamma^{\lambda}\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\\
&-\frac{1}{8\xi^2\sqrt{\xi^2}}\left(\xi_{\mu}e^{\frac{th}{2}}\delta_{\lambda}(e^{\frac{3th}{2}})\right)e^{-th}\left(\xi_{\rho}e^{\frac{th}{2}}\delta_{\nu}(e^{\frac{3th}{2}})\right)e^{-2th}\otimes\gamma^{\lambda}\gamma^{\mu}\gamma^{\nu}\gamma^{\rho},\\\end{aligned}$$ where we have used the notation $\xi^2=\sum_{k=0}^3\xi_k^2$.

Now by using the relation $|D_t|^{-1}|D_t|=1$ and by recursive computation we find the homogeneous terms in the symbol of the inverse, $$\sigma(|D_t|^{-1})\sim\sigma_{-1}(|D_t|^{-1})(\xi)+\sigma_{-2}(|D_t|^{-1})(\xi)+\cdots,$$ where, $$\sigma_{-1}(|D_t|^{-1})(\xi)=\frac{1}{\sqrt{\xi^2}}e^{-th}\otimes I,$$ $$\begin{aligned}
\sigma_{-2}(|D_t|^{-1})(\xi)&=\\
&-\sigma_{-1}(|D_t|^{-1})\{ \sigma_0(|D_t|)\sigma_{-1}(|D_t|^{-1})\\
&+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D_t|)\delta^{\alpha_1}\delta^{\alpha_2}\delta^{\alpha_3}\sigma_{-1}(|D_t|^{-1})\},\\
\sigma_{-3}\{|D_t|^{-1})(\xi)&=\\
&-\sigma_{-1}(|D_t|^{-1})\{ (\sigma_{-1}(|D_t|)\sigma_{-1}(|D_t|^{-1})+\sigma_0(|D_t|)\sigma_{-2}(|D_t|^{-1})\\
&+\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D_t|)\delta^{\alpha_1}\delta^{\alpha_2}\delta^{\alpha_3}\sigma_{-1}(|D_t|^{-1})\\
&+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D_t|)\delta^{\alpha_1}\delta^{\alpha_2}\delta^{\alpha_3}\sigma_{-2}(|D_t|^{-1})\},\end{aligned}$$ $$\begin{aligned}
\sigma_{-4}(|D_t|^{-1})(\xi)&=\\
-\sigma_{-1}(|D_t|^{-1})&\{ \sigma_{-2}(|D_t|)\sigma_{-1}(|D_t|^{-1})
+\sigma_{-1}(|D_t|)\sigma_{-2}(|D_t|^{-1})
+\sigma_0(|D_t|)\sigma_{-3}(|D_t|^{-1})\\
&+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_0(|D_t|)\delta^{\alpha_1}\delta^{\alpha_2}\delta^{\alpha_3}\sigma_{-2}(|D_t|^{-1})\\
&+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{-1}(|D_t|)\delta^{\alpha_1}\delta^{\alpha_2}\delta^{\alpha_3}\sigma_{-1}(|D_t|^{-1})\\
&+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_1(|D_t|)\delta^{\alpha_1}\delta^{\alpha_2}\delta^{\alpha_3}\sigma_{-3}(|D_t|^{-1})\\
&+\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_0(|D_t|)\delta^{\alpha_1}\delta^{\alpha_2}\delta^{\alpha_3}\sigma_{-1}(|D_t|^{-1})\\
&+\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_1(|D_t|)\delta^{\alpha_1}\delta^{\alpha_2}\delta^{\alpha_3}\sigma_{-2}(|D_t|^{-1})\\
&+\sum_{\alpha_1+\alpha_2+\alpha_3=3}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_1(|D_t|)\delta^{\alpha_1}\delta^{\alpha_2}\delta^{\alpha_3}\sigma_{-1}(|D_t|^{-1})\}.\end{aligned}$$ Therefore, the symbol of $\sigma(D_t|D_t|^{-1})$ reads $$\begin{aligned}
       &\sigma(D_t|D_t|^{-1})\sim\left(\sigma_1(D_t)+\sigma_0(D_t)\right)\star\left(\sigma_{-1}(|D_t|^{-1})+\sigma_{-2}(|D_t|^{-1})+\sigma_{-3}(|D_t|^{-1})+\cdots\right)\\
  &\sim\left(\sigma_{1}(D_t)\star\sigma_{-1}(|D_t|^{-1})\right)+\left(\sigma_{1}(D_t)\star\sigma_{-2}(|D_t|^{-1})\right)\\
  &+\left(\sigma_{1}(D_t)\star\sigma_{-3}(|D_t|^{-1})\right)+\left(\sigma_1(D_t)\star\sigma_{-4}(|D_t|^{-1})\right)+\cdots\\
  &+\left(\sigma_{0}(D_t)\star\sigma_{-1}(|D_t|^{-1})\right)+\left(\sigma_{0}(D_t)\star\sigma_{-2}(|D_t|^{-1})\right)+\left(\sigma_0(D_t)\star\sigma_{-3}(|D_t|^{-1})\right)+\cdots.\\ 
   \end{aligned}$$ For $a=1,0$ and $b=-1,-2,-3.\cdots$, one has $$\sigma_a(D_t)\star\sigma_b(|D_t|^{-1})=\sum_{\alpha}\frac{1}{\alpha!}\partial_{\xi}^{\alpha}\sigma_a(D_t)\delta^{\alpha}\sigma_b(|D_t|^{-1}),$$ and each term is of order $a-|\alpha|+b$. By collecting the terms of order $-3$ we obtain, $$\begin{aligned}
   \sigma_{-3}(D_t|D_t|^{-1})&\sim\left(\sigma_1(D_t)\sigma_{-4}(|D_t|^{-1})\right)+\left(\sigma_0(D_t)\sigma_{-3}(|D_t|^{-1})\right)\\
   &+\left(\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(D_t)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-3}(|D_t|^{-1})\right)\\
   &+\left(\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(D_t)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-2}(|D_t|^{-1})\right)\\
   &+\left(\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{0}(D_t)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-2}(|D_t|^{-1})\right)\\
   &+\left(\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{0}(D_t)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-1}(|D_t|^{-1})\right).\\
   \end{aligned}$$ Now, one should notice that Wodzicki residue of a matrix pseudodifferential operator involves a trace taken over the matrix coefficients. By analyzing the terms involved in $\sigma_{-3}(D_t|D_t|^{-1})$ and using the trace identities for gamma matrices (cf. e.g. [@Folland-book]), we see that the only contribution is from the the terms whose matrix coefficient consist of either one $\gamma$ matrix or product of three. Next, we use the following identities, $$\begin{aligned}
&\tr(\gamma^{\lambda})=0,\\
&\tr(\gamma^{\lambda}\gamma^{\mu}\gamma^{\nu})=i\epsilon^{\lambda\mu\nu}\tr(I),\end{aligned}$$ where $\epsilon^{\lambda\mu\nu}$ is the Levi-Civita symbol. By observing that $$\sum_{\lambda,\mu,\nu}\epsilon^{\lambda,\mu,\nu}=0,$$ we conclude that $$\Wres\left(D_t|D_t|^{-1}\right)=\tau\left(\int_{|\xi|=1}\tr\left(\res(D_t|D_t|^{-1}\right)d\xi\right)=0.$$

In [@Bismut-Freed], Bismut-Freed showed that in fact, for a twisted Dirac operator on a Riemannian Spin manifold, the [*local eta-residue*]{}, $\tr \left(\res_x\left(D|D|^{-1}\right)\right)$ vanishes. The above and the following results confirm the same vanishing in the noncommutative case $\T^3$.

The following result proves the regularity at zero of eta function for the coupled Dirac operator $D+A$ on $\T^3$.

\[eta-coupled\] Consider the coupled Dirac operator $\slashed{\partial}+\slashed{A}$ on $\mathcal{H}_{\tau}\otimes\mathbb{C}^2$ over $\T^3$, then $$\Res_{z=0} \eta_{D+A}(z)=0.$$

The coupled Dirac operator $D$ and $D^2$ are given by $$\begin{aligned}
D&=\partial_{\mu}\otimes\gamma^{\mu}+A_{\mu}\otimes\gamma^{\mu}\\
D^2&=\partial_{\mu}\partial_{\lambda}\otimes\gamma^{\mu}\gamma^{\lambda}+\partial_{\mu}(A_{\lambda})\otimes\gamma^{\mu}\gamma^{\lambda}+
   A_{\mu}\partial_{\lambda}\otimes\gamma^{\mu}\gamma^{\lambda}+A_{\mu}A_{\lambda}\otimes\gamma^{\mu}\gamma^{\lambda}.\\
   \end{aligned}$$ Note that $\partial_{\mu}A_{\lambda}(a)=\partial_{\mu}(A_{\lambda})a+A_{\lambda}\partial_{\mu}(a)$ and also $A_{\mu}A_{\lambda}\neq A_{\lambda}A_{\mu}$. Hence, $$D^2=\sum_{\mu}\partial_{\mu}^2\otimes I+2A_{\mu}\partial_{\lambda}\otimes\gamma^{\mu}\gamma^{\lambda}+\partial_{\mu} (A_{\lambda})\otimes\gamma^{\mu}\gamma^{\lambda}+A_{\mu}A_{\lambda}\otimes\gamma^{\mu}\gamma^{\lambda},$$ where $\partial_{\mu} (A_{\lambda})$ in the third term above is a multiplication operator. The symbols are: $$\sigma(D)=\slashed{\xi}+\slashed{A},$$ and $$\sigma(D^2)=\xi^2+2A_{\mu}\xi_{\lambda}\otimes\gamma^{\mu}\gamma^{\lambda}+\partial_{\mu} (A_{\lambda})\otimes\gamma^{\mu}\gamma^{\lambda}+A_{\mu}A_{\lambda}\otimes\gamma^{\mu}\gamma^{\lambda},$$ where $\xi^2=\sum_{\mu}\xi_{\mu}^2\otimes I$. Now, $$\begin{aligned}
   \sigma(|D|)&=\sqrt{\xi^2+2A_{\mu}\xi_{\lambda}\otimes\gamma^{\mu}\gamma^{\lambda}+\partial_{\mu} (A_{\lambda})\otimes\gamma^{\mu}\gamma^{\lambda}+A_{\mu}A_{\lambda}\otimes\gamma^{\mu}\gamma^{\lambda}}\\
   &\sim \sigma_1(\xi)+\sigma_0(\xi)+\sigma_{-1}(\xi)+\cdots,
   \end{aligned}$$ where $$\begin{aligned}
\sigma_1(\xi)&=\sqrt{\xi^2}\otimes I,\\
\sigma_0(\xi)&=\frac{\slashed{A}\slashed{\xi}}{\sqrt{\xi^2}},\\
   \sigma_{-1}(\xi)&=\frac{1}{2\sqrt{\xi^2}}\left(\partial_{\mu}(A_{\lambda})\otimes\gamma^{\mu}\gamma^{\lambda}+A_{\mu}A_{\lambda}\otimes\gamma^{\mu}\gamma^{\lambda}\right)\\
&-\frac{1}{2\xi^2\sqrt{\xi^2}}\left(A_{\mu}\xi_{\lambda}A_{\nu}\xi_{\rho}\otimes\gamma^{\mu}\gamma^{\lambda}\gamma^{\nu}\gamma^{\rho}\right).\end{aligned}$$ By using $\sigma(|D|)\star\sigma(|D|^{-1})\sim1$ and by recursive computation, we get: $$\sigma_{-1}(|D|^{-1})=\frac{1}{\sqrt{\xi^2}}\otimes I.$$ $$\begin{aligned}
   \sigma_{-2}(|D|^{-1})&=\\
   -\frac{1}{\sqrt{\xi^2}}\otimes I&\left(\sigma_0(|D|)\sigma_{-1}(|D|^{-1})+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-1}(|D|^{-1})\right)\\
   &=-\frac{1}{\sqrt{\xi^2}}\otimes I\left(\frac{\slashed{A}\slashed{\xi}}{\sqrt{\xi^2}}.\frac{1}{\sqrt{\xi^2}}\otimes I\right)=-\frac{\slashed{A}\slashed{\xi}}{\xi^2\sqrt{\xi^2}}.\\
   \end{aligned}$$ $$\begin{aligned}
   \sigma_{-3}(|D|^{-1})&=-\frac{1}{\sqrt{\xi^2}}\otimes I\{ \sigma_{-1}(|D|)\sigma_{-1}(|D|^{-1})+\sigma_0(|D|)\sigma_{-2}(|D|^{-1})\\
   &+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{0}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-1}(|D|^{-1})\\
   &+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-2}(|D|^{-1})\\
  &+\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-1}(|D|^{-1})\}\\
   &=-\frac{1}{\sqrt{\xi^2}}\otimes I\left(\frac{\slashed{A}\slashed{\xi}}{\sqrt{\xi^2}}.\frac{-\slashed{A}\slashed{\xi}}{\xi^2\sqrt{\xi^2}}+\partial_{\xi}^1(\sqrt{\xi^2}\otimes I).\delta^1(\frac{-\slashed{A}\slashed{\xi}}{\xi^2\sqrt{\xi^2}})\right)\\
   &=\frac{\xi_{\mu}\xi_{\nu}A_{\rho}A_{\lambda}}{\xi^4\sqrt{\xi^2}}\otimes\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\gamma^{\lambda}+\frac{1}{\sqrt{\xi^2}}\partial_{\xi}^1(\sqrt{\xi^2})\delta^1(\frac{A_{\mu}\xi_{\nu}}{\xi^2\sqrt{\xi^2}})\otimes\gamma^{\mu}\gamma^{\nu}.\\
   \end{aligned}$$ $$\begin{aligned}
   \sigma_{-4}(|D|^{-1})&=-\frac{1}{\sqrt{\xi^2}}\otimes I\{\sigma_{-2}(|D|)\sigma_{-1}(|D|^{-1})+\sigma_{-1}(|D|)\sigma_{-2}(|D|^{-1})\\
   &+\sigma_0(|D|)\sigma_{-3}(|D|^{-1})\\
   &+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{-1}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-1}(|D|^{-1})\\
   &+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{0}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-2}(|D|^{-1})\\
   &+\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{0}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-1}(|D|^{-1})\\
  &+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-3}(|D|^{-1})\\
   &+\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-2}(|D|^{-1})\\
   &+\sum_{\alpha_1+\alpha_2+\alpha_3=3}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-1}(|D|^{-1})\}\\
   &=-\frac{1}{\sqrt{\xi^2}}\otimes I\{\sigma_0(|D|)\sigma_{-3}(|D|^{-1})+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{0}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-2}(|D|^{-1})\\
   &+\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-3}(|D|^{-1})\\
   &+\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(|D|)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-2}(|D|^{-1})\}.
   \end{aligned}$$ Therefore the symbol of $\sigma(D|D|^{-1})$ reads $$\begin{aligned}
       &\sigma(D|D|^{-1})\sim\left(\sigma_1(D)+\sigma_0(D)\right)\star\left(\sigma_{-1}(|D|^{-1})+\sigma_{-2}(|D|^{-1})+\sigma_{-3}(|D|^{-1})+\cdots\right)\\
  &\sim\left(\sigma_{1}(D)\star\sigma_{-1}(|D|^{-1})\right)+\left(\sigma_{1}(D)\star\sigma_{-2}(|D|^{-1})\right)\\
  &+\left(\sigma_{1}(D)\star\sigma_{-3}(|D|^{-1})\right)+\left(\sigma_1(D)\star\sigma_{-4}(|D|^{-1})\right)+\cdots\\
  &+\left(\sigma_{0}(D)\star\sigma_{-1}(|D|^{-1})\right)+\left(\sigma_{0}(D)\star\sigma_{-2}(|D|^{-1})\right)+\left(\sigma_0(D)\star\sigma_{-3}(|D|^{-1})\right)+\cdots.\\ 
   \end{aligned}$$ For $a=1,0$ and $b=-1,-2,-3.\cdots$, one has $$\sigma_a(D)\star\sigma_b(|D|^{-1})=\sum_{\alpha}\frac{1}{\alpha!}\partial_{\xi}^{\alpha}\sigma_a(D)\delta^{\alpha}\sigma_b(|D|^{-1}),$$ and each term is of order $a-|\alpha|+b$. By collecting the terms of order $-3$ we get $$\begin{aligned}
   \sigma_{-3}(D|D|^{-1})&\sim\left(\sigma_1(D)\sigma_{-4}(|D|^{-1})\right)+\left(\sigma_0(D)\sigma_{-3}(|D|^{-1})\right)\\
   &+\left(\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(D)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-3}(|D|^{-1})\right)\\
   &+\left(\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(D)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-2}(|D|^{-1})\right)\\
   &+\left(\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{0}(D)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-2}(|D|^{-1})\right)\\
   &+\left(\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{0}(D)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-1}(|D|^{-1})\right).\\
   \end{aligned}$$ Since $\sigma_0(D)$ has no $\xi$ dependence, the last two terms vanish and therefore we have: $$\begin{aligned}
   \sigma_{-3}(D|D|^{-1})&\sim\left(\sigma_1(D)\sigma_{-4}(|D|^{-1})\right)+\left(\sigma_0(D)\sigma_{-3}(|D|^{-1})\right)\\
   &+\left(\sum_{\alpha_1+\alpha_2+\alpha_3=1}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(D)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-3}(|D|^{-1})\right)\\
   &+\left(\sum_{\alpha_1+\alpha_2+\alpha_3=2}\frac{1}{\alpha_1!\alpha_2!\alpha_3!}\partial_{\xi_1}^{\alpha_1}\partial_{\xi_2}^{\alpha_2}\partial_{\xi_3}^{\alpha_3}\sigma_{1}(D)\delta_1^{\alpha_1}\delta_2^{\alpha_2}\delta_3^{\alpha_3}\sigma_{-2}(|D|^{-1})\right).\\
   \end{aligned}$$ Finally, similar to the proof of Proposition , one observes that the matrix coefficients of all terms in $\sigma_{-3}(D|D|^{-1})$ consist of either one $\gamma$ matrix or product of three $\gamma$ matrices. Again, by using the similar trace identities of $\gamma$ matrices identities we obtain, $$\Wres(\frac{D}{|D|})=\tau\left(\res(D|D|^{-1})\right)=\tau\left(\int_{|\xi|=1}\tr(\sigma_{-3})d\xi\right)=0.$$

Conformal invariance of $\eta_D(0)$
-----------------------------------

In this section we study the variations of $\eta_{D}(0)$ for the spectral triple $(\T^3,\mathcal{H}, D=\slashed{\partial})$. We consider the conformal variation of the Dirac operator and show that $\eta_D(0)$ will remain unchanged.

Consider the commutative spectral triple $(C^{\infty}(M), L^2(M,S)_g,D_g)$ encoding the data of a closed $n$-dimensional spin Riemannian manifold with the spin Dirac operator on the space of spinors. By varying $g$ within its conformal class, we consider $\tilde{g}=k^{-2}g$ for some $k=e^{h}>0$ in $C^{\infty}(M)$. The volume form for the perturbed metric is given by $\dvol_{\tilde{g}}=k^{-n}\dvol_{g}$ and one has a unitary isomorphism $$U:L^2(M,S)_g\longrightarrow L^2(M,S)_{\tilde{g}}$$ by $$U(\psi)=k^{\frac{n}{2}}\psi,$$ It can be shown that (cf. [@Hitchin-harmonic]) $D_{\tilde{g}}=k^{\frac{n+1}{2}}D_g k^{\frac{-n+1}{2}}$ and hence

$$U^*D_{\tilde{g}}U=k^{\frac{-n}{2}}(k^{\frac{n+1}{2}}D_gk^{\frac{-n+1}{2}})k^{\frac{n}{2}}=\sqrt{k}D_g\sqrt{k}.$$

The above property of the Dirac operator is usually referred to as being conformally covariant. It is a known fact that $\eta_D(0)$ for a Dirac operator on an odd dimensional manifold is a conformal invariant [@APS], i.e. it is invariant under the conformal changes of the metric. In a more general context, in [@Rosenberg-covariant; @Parker-Rosenberg] and [@Paycha-Rosenberg2006], by using variational techniques it was shown that for a coformally covariant self-adjoint differential operator $A$, $\eta_A(0)$ is a conformal invariant.

In the framework of noncommutative geometry, conformal perturbation of the metric is implemented by changing the volume form [@Cohen-Connes1992], namely, by fixing a positive element $k=e^h$ for $h^*=h$ in $C^{\infty}(\mathbb{T}^n_{\theta})$, one constructs the following positive functional $$\varphi_k(a)=\tau (ak^{-n}).$$ By normalizing the above functional one obtains a state which we denote by $\varphi$.\
The state $\varphi$ defines an inner product $$\left<a,b\right>_{\varphi}=\varphi(b^*a)\quad a,b\in C^{\infty}(\mathbb{T}^n_{\theta})$$ and hence one obtains a Hilbert space $\mathcal{H}_{\varphi}$ by GNS construction. The algebra $C^{\infty}(\mathbb{T}^n_{\theta})$ acts unitarilly on $\mathcal{H}_{\varphi}$ by left regular representation and the right multiplication operator $R_{k^{n/2}}$ extends to a unitary map $U_0:\mathcal{H}_{\tau}\longrightarrow\mathcal{H}_{\varphi}$. In fact one has, $$\left<U_0a,U_0b\right>_{\varphi}=\varphi(k^{n/2}b^*ak^{n/2})=\tau(k^{n/2}b^*ak^{n/2}k^{-n})=\tau(b^*a)=\left<a,b\right>_{\tau}.$$

We put $\mathcal{\tilde{H}}=H_{\varphi}\otimes\mathbb{C}^N$, the action of $C^{\infty}(\mathbb{T}^n_{\theta})$ on $\mathcal{\tilde{H}}$ is given by $$a\longrightarrow a\otimes1$$ and the map $$U=U_0\otimes I:\mathcal{H}\longrightarrow\mathcal{\tilde{H}}$$ is a unitary equivalence between the two Hilbert spaces.\
Now we consider the operator $\tilde{D}=R_{k^{\frac{n+1}{2}}}DR_{k^{\frac{-n+1}{2}}}$.

$(C^{\infty}(\mathbb{T}^n_{\theta}),\mathcal{\tilde{H}},\tilde{D})$ and $(C^{\infty}(\mathbb{T}^n_{\theta}),\mathcal{H},R_{\sqrt{k}}DR_{\sqrt{k}})$ are spectral triples, and the map $U$ is a unitary equivalence between them.

Note that left multiplication by an element $a\in C^{\infty}(\mathbb{T}^n_{\theta}) $ commutes with right multiplication operators $R_{k^2}$, $R_{k^{-1}}$ and $R_{\sqrt{k}}$ . Also the two norms coming from $<,>_{\varphi}$ and $<,>_{\tau}$ are equivalent. So both commutators $[a,\tilde{D}]$ and $[a,\sqrt{k}D\sqrt{k}]$ are bounded. The unitary equivalence easily follows from definition of $U$ and $\tilde{D}$.

In next step, we convert the right multiplications in $R_{\sqrt{k}}DR_{\sqrt{k}}$ to left multiplication using the real structure on noncommutative tori (see example \[real-tori\]) . It is easily seen that $$J_0R_{\sqrt{k}}\partial_{\mu}R_{\sqrt{k}}J_0=-\sqrt{k}\partial_{\mu}\sqrt{k},$$ and $$JR_{\sqrt{k}}DR_{\sqrt{k}}J=\sqrt{k}D\sqrt{k}.$$

Since $\sqrt{k}D\sqrt{k}$ is iso-spectral to $R_{\sqrt{k}}DR_{\sqrt{k}}$ (being intertwined by $J$), the following definition is reasonable.

The conformal perturbation of the Dirac operator $D$ is the 1-parameter family $(C^{\infty}(\mathbb{T}^n_{\theta}),\mathcal{H},D_t)$, where $$D_t=e^{\frac{th}{2}}De^{\frac{th}{2}},~~~h=h^*\in C^{\infty}(\mathbb{T}^n_{\theta}).$$

We need the following formula for variation of eta function.

\[eta-var\] Let $\{D_t\}$ be a smooth 1-parameter family of invertible self-adjoint elliptic operators of order $d$, then $$\frac{d}{dt}\eta_{D_t}(z)=-z\TR(\dot{D_t}(D^2_t)^{-(\frac{z+1}{2})}).$$

For $k>0$ odd, $\eta_{D^k}(z)=\eta_D(kz)$, hence we can replace $D$ by $D^k$ for $k$ large enough.\
For $d$ large enough, $(D-\lambda)^{-1}$ is trace class, so one can write $$\eta_{D_t}(z)=\frac{1}{2\pi i}\int_{\Gamma_1}\lambda^{-z}\Tr(D_t-\lambda)^{-1}d\lambda-\int_{\Gamma_2}(-\lambda)^{-z}\Tr(D_t-\lambda)^{-1}d\lambda,$$ where $\Gamma_1$ and $\Gamma_2$ are appropriate contours around positive and negative eigenvalues of $D$.

Now, $$\frac{d}{dt}(D_t-\lambda)^{-1}=-(D_t-\lambda)^{-1}\dot{D_t}(D_t-\lambda)^{-1},$$ so $$\Tr(\frac{d}{dt}(D_t-\lambda)^{-1})=-\Tr(\dot{D_t}(D_t-\lambda)^{-2}),$$ therefore $$\frac{d}{dt}\eta_{D_t}(z)=\frac{1}{2\pi i}\Tr(\dot{D_t}(\int_{\Gamma_1}-\lambda^{-s}\Tr(D_t-\lambda)^{-2}d\lambda
 +\int_{\Gamma_2}(-\lambda)^{-z}\Tr(D_t-\lambda)^{-2}d\lambda)).$$ Now, integration by parts in both integrals and the formula $\frac{d}{d\lambda}(D_t-\lambda)^{-1}=(D_t-\lambda)^{-2}$ gives the result.

By applying the above lemma to the Dirac operator $D$ on noncommutative 3-torus we get, $$\begin{aligned}
   \partial\eta_D(0)=\left[\left.\frac{d}{dt}\right|_{t=0}\eta_{D_t}(z)\right]_{z=0}&=\left[-z\TR(\partial D|D|^{-1}(D ^2)^{\frac{-z}{2}})\right]_{z=0}\\
   =-\Wres(\partial D|D|^{-1})&.
  \end{aligned}$$

By the properties of Wodzicki residue, one sees that $\eta(0)$ is constant under smoothing perturbations, so if $\ker(D)\neq0$, on can replace $D$ by $D+\Pi$ where $\Pi$ is the projection on the finite dimensional kernel of $D$ and hence the result of above lemma still makes sense for an operator with nontrivial kernel.

Consider the spectral triple $(\T^3,\mathcal{H},D=\slashed{\partial})$, then $\eta_D(0)$ is invariant under the conformal perturbations.

Consider the 1-parameter family of operators $$D_t=e^{\frac{th}{2}}De^{\frac{th}{2}},~~~h=h^*\in\T^3.$$ One computes $$\dot{D_t}=\frac{1}{2}(hD_t+D_th),$$ and $$\partial D=\frac{1}{2}(hD+Dh).$$ Therefore by trace property of noncommutative residue we get $$\partial\eta(D,0)=-\Wres(hD|D|^{-1}),$$ and it is seen that $$\Wres(hD|D|^{-1})=\tau\left( \tr( \res(hD|D|^{-1})\right)=\tau\left(\tr(h~\res(D|D|^{-1})\right).$$ Now by looking at the symbols $$\begin{aligned}
 \sigma(D)&\sim\xi_{\mu}\otimes\gamma^{\mu},\\
 \sigma(|D|)&\sim |\xi|^2\otimes I,
 \end{aligned}$$ we get, $$\sigma(D|D|^{-1})\sim\frac{\xi_{\mu}}{|\xi|}\otimes\gamma^{\mu}.$$ Therefore, there is no homogeneous term of negative order in the symbol of $\sigma(D|D|^{-1})$ and hence, $$\Wres(hD|D|^{-1})=0.$$

Spectral flow and odd local index formula
-----------------------------------------

The value of eta function at zero is intimately related to another spectral quantity called the spectral flow. To motivate this relation, consider a family $A_t$ of $N\times N$ Hermitian matrices. The spectral flow of the family $A_t$ is defined as the net number of the eigenvalues of $A_t$ passing zero. It is easily seen that the difference of the signatures of the end points is related to the spectral flow of the family by, $$\eta_{A_1}(0)-\eta_{A_0}(0)=2\SF(A_t).$$ The spectral flow of a family of self adjoint elliptic operators $A_t$ on manifolds can also be defined ([@APS]) but the above equality holds along with a correction term (see [@Paycha-Cardona] for a proof).

\[eta-flow\] Let $A_t$ be a self-adjoint family of elliptic operators of order $a$ with $A_0$ and $A_1$ invertible, then $$\eta_{A_1}(0)-\eta_{A_0}(0)=2\SF(A_t)-\frac{1}{a}\int_0^1\frac{d}{dt}\eta_{A_t}(0)dt.$$

It’s a remarkable fact due to Getzler ([@Getzler-flow]) that the spectral flow of the family of operators $D_t=D+tg^{-1}[D,g]$ interpolating $D$ and $g^{-1}Dg$ for the Dirac operator $D$ on an odd dimensional spin manifold and $g\in C^{\infty}(M,GL(N))$ in fact gives the index of $PgP$ where $P=\frac{1+F}{2}$ and $F=D|D|^{-1}$ is the sign operator. Therefore one obtains a pairing between the odd K-theory and K-homology: $$\label{flow-index}
\SF(D_t)=\index(PgP).$$

The above equality has been generalized to the framework of non commutative index theory [@Carey-Phillips98; @Carey-Phillips]. Consider a spectral triple $(\mathcal{A},\mathcal{H},D)$ with sufficiently nice regularity properties and let $u\in\mathcal{A}$ be a unitary element and $[u]\in K^1(A)$ the corresponding class in K-theory. The sign operator $F=\frac{D}{|D|}$ gives the projection $P=\frac{1+F}{2}$ and $PuP$ is a Fredholm operator. For the family $D_t=D+tu^*[D,u]$ interpolating between $D$ and $u^*Du$ one can show that(see [@Carey-Phillips98] and [@Carey-Phillips]): $$\SF(D_t)=\index(PuP).$$ The odd local index formula of Connes-Moscovici in turn expresses the right side of the above equality as a pairing between the periodic cyclic cohomolgy and K-theory given by a local formula in the sense of non commutative geometry [@CM-local-index]. Note that for the spectral triple $(\T^3,\mathcal{H},D=\slashed{\partial}$), by Proposition (\[eta-coupled\]), the eta function $\eta_{D_t}(0)$ over the family of Dirac operators $D_t=D+tu^*[D,u]$ makes sense and with a minor modification of the proof of Proposition (\[eta-flow\]) in the commutative case we get, $$\eta_{u^*Du}(0)-\eta_{D}(0)=2\SF(D_t)-\int_0^1\frac{d}{dt}\eta_{D_t}(0)dt.$$ Since $u^*Du$ and $D$ have the same spectrum, we have $\eta_{u^*Du}(0)=\eta_{D}(0)$ and therefore we obtain obtain the following integrated local formula for the index: $$\index(PuP)=\frac{1}{2}\int_0^1\frac{d}{dt}\eta_{D_t}(0)dt.$$

Bosonic functional determinant and conformal anomaly
====================================================

Consider a positive elliptic differential operator $A$ on a closed manifold $M$, the corresponding bosonic action is defined by $$\begin{aligned}
        S_{\bos}(\phi)=(\phi,A\psi),\end{aligned}$$ where $\phi:M\rightarrow\mathbb{C}^d$ is a bosonic field living on $M$. Upon quantization, one constructs the partition function $$\begin{aligned}
\mathcal{Z}_{\bos}=\int e^{\frac{-1}{2}S_{bos}(\phi)}[\mathcal{D}\phi],\end{aligned}$$ where $[\mathcal{D}\phi]$ is a formal measure on the configuration of bosonic fields. The corresponding effective action is given by $$\begin{aligned}
\mathcal{W}=\log\mathcal{Z}.\end{aligned}$$ One should think of $\mathcal{Z}_{\bos}$ formally as $(\det A)^{-\frac{1}{2}}$. Through zeta regularization scheme (see e.g. [@QFT-analytic]) one defines, $$\mathcal{Z}:=e^{\frac{1}{2}\zeta_A'(0)}$$ and hence, $$\mathcal{W}=\frac{1}{2}\zeta_A'(0),$$ where the spectral zeta function is defined by $$\zeta_A(z)=\TR(A^{-z}),$$ and we have dropped the dependence on the spectral cut $L_{\theta}$ in definition of the complex power. Note that since $A$ is a differential operator, by Proposition (\[Laurent exp\]), $\zeta_A(z)$ is regular at $z=0$ and $\zeta_A'(0)$ makes sense.

Let $D$ be the Dirac operator on a closed odd dimensional spin manifold, it’s a known fact that $\zeta_{|D|}'(0)$ is a conformal invariant (see [@Rosenberg-covariant] and [@Parker-Rosenberg]). We prove the analogue of this result for the Dirac operator on noncommutative n-torus $(C^{\infty}(\mathbb{T}^n_{\theta}),\mathcal{H},D=\slashed{\partial})$ in odd dimensions.

In even dimensions the conformal variation is not zero and hence conformal quantum anomaly exists. In dimension two, the Polyakov anomaly formula [@Polyakov] gives a local expression for conformal anomaly of the Laplacian on a Riemann surface. Following [@Connes-Moscovici2014] we derive an analogue of this formula for noncommutative 2-torus using the formalism of canonical trace.

For proving the conformal invariance of $\zeta_{|D|}'(0)$ on $C^{\infty}(\mathbb{T}^{2p+1}_{\theta})$ we will need to carefully analyze the coefficients in short time asymptotic of the trace of heat operator $e^{-tD^2}$.

\[heatasymp\] Consider the spectral triple $(C^{\infty}(\mathbb{T}^n_{\theta}),\mathcal{H},D)$. One has the following asymptotic expansion $$\Tr(e^{-tD^2})\sim\sum_{i=0}^{\infty}a_i(D^2)t^{\frac{i-n}{2}}\quad t\rightarrow0,$$ where coefficients $a_i$ are computed by integrating local terms.

We have, $$\Tr(e^{-tD^2})=\tau(\int_{\mathbb{R}^n}\tr\sigma(e^{-tD^2})(\xi)d\xi).$$ The operator $e^{-tD^2}$ is again a pseudodifferential operator of order $-\infty$ whose symbol can be computed using the symbol product rules. For any $\lambda$ in the resolvent of $D^2$, $(\lambda-D^2)^{-1}$ is a pseudodifferential operator of order $-2$ with symbol $r(\xi,\lambda)$ which can be written as its homogeneous parts $r=r_0+r_1+\cdots$, with $r_k(\xi,\lambda)$ homogeneous of order $-k-2$ in $(\xi,\lambda)$ i.e. $$r_k(t\xi,t^2\lambda)=t^{-2-k}r_k(\xi,\lambda)\quad \forall t>0$$ The operator $e^{-tD^2}$ is also a pseudodifferential operator with the symbol $e_0+e_1+e_2+\cdots$ where $e_n\in S^{-\infty}$ and are defined by $$e_n(t,\xi)=\frac{1}{2\pi i}\int_\gamma e^{-t\lambda}r_n(\xi, \lambda)d\lambda\quad t>0$$

$$\begin{aligned}
 \int_{\mathbb{R}^n} \sigma(e^{-tD^2})(t,\xi)d\xi\\
&=\sum_{i} \int  e_i(t,\xi)d\xi\\\end{aligned}$$

The integrands $\int_{\mathbb{R}^n}  e_i(t,\xi)d\xi$ are homogeneous in $\xi$. $$\begin{aligned}
\int_{\mathbb{R}^n} e_n(t,\xi) d\xi&= \frac{1}{2\pi i}\int_{\mathbb{R}^n}  \int_\gamma e^{-t\lambda}r_n(x,\xi, \lambda)d\lambda d\xi\\
&=\frac{1}{2\pi i}\int_{\mathbb{R}^n}  \int_\gamma e^{-\lambda'}r_n(t^{-1/2}\xi', t^{-1}\lambda')\frac{d\lambda'}{t} \frac{d\xi'}{t^{n/2}}\\
&=t^{\left(\frac{i-n}{2}\right)}\frac{1}{2\pi i}\int_{\mathbb{R}^n} \int_\gamma e^{-\lambda}r_n(\xi, \lambda)d\lambda d\xi\\\end{aligned}$$ In identity second identity we have used the change of variable $t\lambda=\lambda'$ and $t^{1/2}\xi=\xi'$.

Hence we have $$\int_{\mathbb{R}^n}\sigma(e^{-tD^2})(\xi)d\xi\sim\sum_{i=0}\beta_i t^{\frac{i-n}{2}},$$ Where $$\beta_i(x)=\frac{1}{2\pi i}\int_{\mathbb{R}^n} \int_\gamma e^{-\lambda}r_i(\xi, \lambda)d\lambda d\xi$$ and hence $$\Tr(e^{-tD^2})\sim\sum_{i=0}^{\infty}t^{\frac{i-n}{2}}a_i(D^2),~~t\rightarrow0,$$ Where $a_i(D^2)=\tau(\tr(\beta_i))$

The above asymptotic expansion holds for any self-adjoint elliptic differnetial operator $A\in\Psi_{cl}^{*}(C^{\infty}(\mathbb{T}^n_{\theta}))$. In fact, one has $$\Tr(e^{-tA^2})\sim\sum_{i=0}^{\infty}a_i(A^2)t^{\frac{i-n}{\alpha}}~~~t\rightarrow0,$$ where $\alpha=ord(A)$.

One also has the following asymptotic expansion $$\Tr(ae^{-tD^2})\sim\sum_{i=0}^{\infty}a_i(D^2,a)t^{\frac{i-n}{2}}~~~t\rightarrow0,$$ Where $a\in C^{\infty}(\mathbb{T}^n_{\theta})$ is considered as a multiplication operator and $$a_i(D^2,a)=\tau(\tr(a\beta_i))$$ in above computations.

Next, we have the following variational result.

\[zetavar\] Let $A_t$ be a 1-parameter family of positive elliptic differential operators of fixed order $\alpha$ on noncommutaive n-torus, then $$\frac{d}{dt}\zeta_{A_t}(z)=-z\TR(\dot{A_t}(A_t)^{-z-1}).$$

By using the contour integral formula for the complex powers we have, $$\begin{aligned}
\frac{d}{dt}\TR(A_t^{-z})=\frac{1}{2\pi i}\int_{\Gamma}\lambda^{-z}\frac{d}{dt}\TR(\lambda-A_t)^{-1}d\lambda\\=\frac{1}{2\pi i}\int_{\Gamma}\lambda^{-z}\TR(\dot{A_t(}\lambda-A_t)^{-2})d\lambda\\
=\TR\left(\dot{A_t}\frac{1}{2\pi i}\int_{\Gamma}\lambda^{-z}(\lambda-A_t)^{-2}d\lambda\right).\\\end{aligned}$$ Now using integration by parts formula gives the result.

Note that by above lemma we have: $$\frac{d}{dt}\zeta_{A_t}(0)=\left[-z\TR(\dot{A_t}A^{-1}A^{-z})\right]_{z=0}=-\Wres(\dot{A_t}A^{-1}),$$ and therefore, the value of zeta at zero is constant under smoothing perturbations and hence the above result makes sense for non invertible operators.

We also need the following result relating the constant term in Laurent expansion of $\TR(AQ^{-z})$ around $z=0$ to the constant term in the asymptotic expansion $\Tr(Ae^{-tQ})$ at $t=0$. Let $\f.p.\left.\TR(AQ^{-z})\right|_{z=0}$ denote the constant term in the Laurent expansion and also $\f.p.\left.\Tr(Ae^{-tQ})\right|_{t=0}$ denote the constant term in the heat trace asymptotic near zero. We refer the reader to [@Berlin-Getzler] for a proof of the following result in a more general form in commutative case.

\[BGVlemma\] Consider the differential operators $A,Q\in\Psi_{cl}^{*}(C^{\infty}(\mathbb{T}^n_{\theta}))$, where $Q$ is positive with positive order $q$. Then, $$\f.p.\left.\TR(AQ^{-z})\right|_{z=0}=\f.p.\left.\Tr(Ae^{-tQ})\right|_{t=0}.$$

Consider the spectral triple $(C^{\infty}(\mathbb{T}^{2p+1}_{\theta}),\mathcal{H},D=\slashed{\partial})$. The value $\zeta'_{|D|}(0)$ is a conformal invariant.

It is enough to show that $\zeta'_{\Delta}(0)$ for $\Delta=D^2$ is conformally invariant. By the Lemma (\[zetavar\]) we have, $$\partial\zeta_{\Delta}(z)=-z\TR(\partial \Delta\cdot \Delta^{-1}\Delta^{-z}),$$ so we have the following Laurent expansion around $z=0$: $$\partial\zeta_{\Delta}(z)=-z\left(a_{-1}\frac{1}{z}+a_0+a_1z+\cdots\right).$$ Therefore $$\label{zetavar2}
\partial\zeta_{\Delta}'(0)=\frac{d}{dz}\left.\partial\zeta_{\Delta}(z)\right|_{z=0}=-a_0=-\tau\left(\int\!\!\!\!\!\!-\sigma(\partial\Delta\Delta^{-1})-\frac{1}{2}res(\partial\Delta\Delta^{-1}\log\Delta)\right).$$

Consider the conformal perturbation of the Dirac operator, $$D_t=e^{\frac{th}{2}}De^{\frac{th}{2}}~~~h=h^*\in\mathcal{A}_{\theta}.$$ An easy computation gives $$\partial\Delta=\left.\frac{d}{dt}D_t^2\right|_{t=0}=\frac{h}{2}D^2+DhD+\frac{1}{2}D^2h,$$ and therefore by inserting the above identity into (\[zetavar2\]) we get $$\partial\zeta_{\Delta}'(0)=-2\tau\left(\int\!\!\!\!\!\!-\sigma(h)-\frac{1}{2}\res(h\log\Delta)\right)=-2\f.p.\left.\TR(h\Delta^{-z})\right|_{z=0}.$$

Now using the Lemma (\[BGVlemma\]) we have $$\label{zetavar3}
\partial\zeta_{\Delta}'(0)=-2\f.p.\left.\Tr(h e^{-t\Delta})\right|_{t=0}=-2a_{2p+1}(h,\Delta)=-2\tau(\tr(h\beta_{2p+1})).$$ By examining the proof of Lemma (\[heatasymp\]), we see that $a_n=0$ for odd $n$ since the integrand involved is an odd function, hence the result is obtained.

In the following we give an analogue of Polyakov anomaly formula [@Polyakov] for the spactral triple $(C^{\infty}(\mathbb{T}^2_{\theta}),\mathcal{H},\slashed{\partial})$. Note that although $\log\det(\Delta)=-\zeta_{\Delta}'(0)$ is not local[^2] (see Proposition (\[Laurent exp\])), the difference between $\log\det_{\Delta_h}$ of the conformally perturbed Laplacian and $\log\det_{\Delta}$ can be given by a local formula. This is of course an example of the local nature of anomalies in quantum field theory. Here we only express this difference as an integrated anomaly and refer the reader to [@Connes-Moscovici2014] for further computations and interpretation of the formula.

[**( A conformal anomaly formula)**]{} Consider the spectral triple $(C^{\infty}(\mathbb{T}^2_{\theta}),\mathcal{H},D=\slashed{\partial})$, $\Delta=D^2$ and $\Delta_h=D_h^2$ where $D=e^{\frac{h}{2}}De^{\frac{h}{2}}$. The difference between $\log\det_{\Delta_h}$ of the conformally perturbed Laplacian and $\log\det_{\Delta}$ can be given by a local formula: $$\label{logdetpolyakov}
\log\det(\Delta_h)-\log\det(\Delta)=-\left(\zeta_{\Delta_h}'(0)-\zeta_{\Delta}'(0)\right)=-\int_0^1\frac{d}{dt}\zeta_{\Delta_t}'(0)dt,$$ where $\Delta_t=D_t^2$, $D=e^{\frac{th}{2}}De^{\frac{th}{2}}$.

The proof follows from the Lemma (\[zetavar\]) and fundamental theorem of calculus along the family $\Delta_t$.

Note that, by using the similar argument used in derivation of the equations and we see that the integrand in equation is given by $$\begin{aligned}
  &\frac{d}{dt}\zeta_{\Delta_t}'(0)=-\tau\left(\int\!\!\!\!\!\!-\sigma(\frac{d}{dt}\Delta_t\Delta_{t}^{-1})-\frac{1}{2}\res(\frac{d}{dt}\Delta_t\Delta_{t}^{-1}\log\Delta_{t})\right)=\\
 & -2a_{2}(h,\Delta_t)=-2\tau\left(\tr(h\beta_{2}(\Delta_t))\right).
 \end{aligned}$$ of course, computing the density $\tr(\beta_{2}(\Delta_t))$ requires considerable amount of computations (see [@Connes-Moscovici2014] and [@Farzad-Masoud2013]).

Fermionic functional determinant and the induced Chern-Simons term
==================================================================

Consider a closed manifold $M$ and the classical fermionic action defined by $$S_{\fer}(\bar{\psi},\psi)=(\bar{\psi},D\psi),$$ where $\bar{\psi}$ and $\psi $ are fermion fields on $M$ and $D=i\gamma^{\mu}\nabla_{\mu}$ is the Dirac operator . The fermionic partition function is given by $$\mathcal{Z}_{\fer}=\int e^{-S_{\fer}(\bar{\psi},\psi)}[\mathcal{D}\bar{\psi}][\mathcal{D}\psi],$$ where $[\mathcal{D}\bar{\psi}]$ and $[\mathcal{D}\psi]$ are formal measures on the space of fermions. The partition function $\mathcal{Z}$ should be thought as the $\det(D)$. The one-loop fermionic effective action is defined as $$\mathcal{W}:=\log(\mathcal{Z})=\log\det(D).$$ Similar to bosonic case, after choosing a spectral cut $L_{\phi}$ [^3], by definition we have (see ): $$\zeta_D(z)=TR(D_{\phi}^{-z}),$$ and hence one can define $$\begin{aligned}
      \mathcal{W}=\log\det(D):=-\zeta_D'(0).
  \end{aligned}$$ Note that due to choice of a spectral cut, there is an ambiguity involved in above definition of one-loop effective action. Let $\zeta_D ^{\uparrow}(z)$ and $\zeta_D ^{\downarrow}(z)$ be the spectral zeta functions corresponding to a choice of $\phi$ in upper and lower half plane respectively. a quick computation shows that (cf. [@Ponge-assym]) $$\begin{aligned}
   \zeta_D ^{\uparrow}(z)-\zeta_D ^{\downarrow}(z)=(1-e^{-i\pi z})\zeta_D ^{\uparrow}(z)-(1-e^{-i\pi z})\eta_{D}(z),\end{aligned}$$ and the measure of ambiguity in the effective action is given by $$\begin{aligned}
\zeta_D ^{\uparrow}(0) '-\zeta_D ^{\downarrow}(0)'=i\pi\zeta_D ^{\uparrow}(0)+i\pi\eta_D(0).\end{aligned}$$ In $2p+1$ dimensions $\zeta_D ^{\uparrow}(0)=\zeta_D ^{\downarrow}(0)=0$ and hence $$\begin{aligned}
\zeta_D ^{\uparrow}(0) '-\zeta_D ^{\downarrow}(0)'=i\pi\eta_D(0).\end{aligned}$$ Therefore the ambiguity is given by the non local quantity $\eta_D(0)$ which also depends on the gauge field coupled to the Dirac operator. The measure of this dependence can be given by a local formula which in physics literature is referred to as the induced Chern-Simons term generated by the coupling of a massless fermion to a classical gauge field (cf. e.g. [@QFT-analytic]).

Here, we give an analogue of this local term for the coupled Dirac operator on noncommutative 3-torus. We consider the operator $D=\slashed{\partial}+\slashed{A}$ on $\T^3$ and compute the variation of the eta invariant $\eta_{D}(0)$ with respect to the vector potential.

First we state the following lemma. The proof in commutative case also works in noncommutative setting with minor changes and we will not reproduce it here (see [@Paycha-Rosenberg2006]).

\[res-asym\] In the asymptotic expansion of the heat kernel for a positive elliptic differential operator $A\in\Psi_{cl}^{*}(C^{\infty}(\mathbb{T}^n_{\theta}))$ of $ord(A)=\alpha$, $$\int_{\mathbb{R}^n}\sigma(e^{-tA})(\xi)d\xi\sim\sum_{i=0}^{\infty}\beta_i t^{\frac{i-n}{\alpha}},$$ one has $$\res(A^{-k})=\frac{\alpha}{(k-1)!}\beta_{n-\alpha k},$$ where $k\in\mathbb{Z}^{+}$ and $\beta_{n-\alpha k}=0$ if $\alpha k\notin\mathbb{Z}$.

Consider the family $D_t=\slashed{\partial}+\slashed{A}_t$ on $\T^3$. By performing the variation with respect to $A_{\mu}$ we have, $$\partial D=\gamma^{\mu}\partial A_{\mu}.$$

For the family of Dirac operators $\{D_t\}$, a proof similar to the proof of Proposition (\[eta-coupled\]) shows that $\eta_{D_t}(z)$ is regular at $z=0$ and therefore along this family, $\eta_{D_t}(0)$ makes sense. By using Lemma (\[eta-var\]) we obtain: $$\partial\eta(0)=\Wres\left(\partial D|D|^{-1}\right)=\Wres\left(\gamma^{\mu}\partial A_{\mu}|D|^{-1}\right).$$

Now by using Lemma (\[res-asym\]) and the fact that $\gamma^{\mu}\partial A_{\mu}$ is a zero order operator it follows that the variation of the eta invariant with respect to the gauge field $A_{\mu}$ is given by the following local formula, $$\partial\eta_{D}(0)=\Wres\left(\gamma^{\mu}\partial A_{\mu}|D|^{-1}\right)=\tau \left(\tr\left(\gamma^{\mu}\partial A_{\mu}\beta_{2}\right)\right).$$

\#1[0=]{}

[50]{}

M.F. Atiyah, V.K. Patodi, I.M. Singer. Spectral asymmetry and [R]{}iemannian geometry.[I]{}. , Volume77:43–69, 1975.

M.F. Atiyah, V.K. Patodi, I.M. Singer. , Volume79:71–99, 1976.

N. Berline, E. Getzler, M. Vergne. . , 1992.

T. A. Bhuyain and M. Marcolli. The [R]{}icci flow on noncommutative two-tori. , 101(2):173–194, 2012.

J.M. Bismut, D.  Freed. , 107: 103-163, 1986.

A.A. Bytsenko, G. Cognola, E. Elizalde, Moretti. . , 2003.

A. Cardona, C. Ducourtioux, S. Paycha . . , 242: 31–65, 2003.

A. Carey, J. Phillips. . , 50(4):673–718, 1998.

A. Carey, J. Phillips. . , 31(2):135–194, 2004.

P. Cohen, A. Connes. Conformal geometry of the irrational rotation algebra. Preprint of the MPI, Bonn. no. MPI/92-23, 1992.

A. Connes. algèbres et géométrie différentielle. , 290(13):A599–A604, 1980.

A. Connes. . Academic Press, Inc., San Diego, CA, 1994.

A. Connes, H. Moscovici. The local index formula in noncommutative geometry. , 5(2):174–243, 1995.

A. Connes, H. Moscovici. In . 57-71, Aspects Math., E38, Vieweg, 2008.

A. Connes, H. Moscovici. Modular curvature for noncommutative two-tori. , 27(3):639–684, 2014.

A. Connes, P. Tretkoff. The [G]{}auss-[B]{}onnet theorem for the noncommutative two torus. In [*Noncommutative geometry, arithmetic, and related topics*]{}, pages 141–158. Johns Hopkins Univ. Press, Baltimore, MD, 2011.

L. Dabrowski, A. Sitarz, [*Curved noncommutative torus and Gauss–Bonnet,*]{} J. Math. Phys. 54, 013518 (2013).

L. Dabrowski, A. Sitarz. Asymmetric noncommutative torus. arXiv:1406.4645, 2014.

A. Fathi, A. Ghorbanpour, M. Khalkhali. . , Arxiv:1410.0475, 2014.

F. Fathizadeh. On the Scalar Curvature for the Noncommutative Four Torus. ArXiv:1410.8705, 2014.

F. Fathizadeh, M. Khalkhali. The [G]{}auss-[B]{}onnet theorem for noncommutative two tori with a general conformal structure. , 6(3):457–480, 2012.

F. Fathizadeh, M. Khalkhali. Scalar curvature for the noncommutative two torus. , 7(4):1145–1183, 2013.

F. Fathizadeh, M. Khalkhali. Scalar curvature for noncommutative four-tori. , To appear.

F. Fathizadeh, M. Wong. Noncommutative residues for pseudo-differential operators on the noncommutative two-torus. , 2(3):289–302, 2011.

G.B. Folland. . AMS, Mathematical Surveys and Monographs, 2008.

E. Getzler. , 32(3):489–507, 1993.

P. B. Gilkey. , 240(2):183–189, 1979.

P. B. Gilkey. Volume 11 of [*Mathematics Lecture Series*]{}. Publish or Perish Inc., Wilmington, DE, 1984.

J.M. Gracia-Bond[í]{}a, J. V[á]{}rilly, H. Figueroa. . Birkhäuser Advanced Texts: Basler Lehrbücher, 2001.

N. Hitchin. . , 14(1):1 - 55, 1974.

M. Kontsevich, S. Vishik. Geometry of determinants of elliptic operators. In [*Functional analysis on the eve of the 21st century, [V]{}ol. 1 ([N]{}ew [B]{}runswick, [NJ]{}, 1993)*]{}, volume 131 of [*Progr. Math.*]{}, pages 173–197. Birkhäuser Boston, Boston, MA, 1995.

M. Khalkhali and A. Moatadelro. A Riemann-Roch theorem for the noncommutative two torus. , 86(0):19 – 30, 2014.

M. Lesch. On the noncommutative residue for pseudodifferential operators with log-polyhomogeneous symbols. , 17(2):151–187, 1999.

C. Levy, C. N. Jiménez, S. Paycha. . ArXiv:1303.0241 \[math.AP\], 2013.

T. Parker, S. Rosenberg. . , 25:199–222, 1987.

S. Paycha. . Volume 59, [ *University Lecture Series*]{}. American Mathematical Society, Providence, RI, 2012.

S. Paycha, S. Rosenberg. Conformal anomalies via canonical traces. In [*Analysis, geometry and topology of elliptic operators*]{}, pages 263–294. World Sci. Publ., Hackensack, NJ, 2006.

S. Paycha, S. Scott. A [L]{}aurent expansion for regularized integrals of holomorphic symbols. , 17(2):491–536, 2007.

S. Polyakov. Quantum geometry of bosonic strings. , 103:207–210, 1981.

R. Ponge. Spectral asymmetry, zeta functions, and the noncommutative residue. , 17(9): 1065–1090, 2006.

S. Rosenberg. . ,Volumes [2-36]{}(3):553-568, 1987.

S. Scott. . Oxford Mathematical Monographs. Oxford University Press, Oxford, 2010. M. Wodzicki. Noncommutative residue. I. Fundamentals. , Lecture Notes in Math., 1289, Springer, Berlin, 1987, p. 320-399.

[^1]: [*E-mail addresses*]{}: afathiba@uwo.ca, masoud@uwo.ca

[^2]: Roughly, it means that it can not be written as an integral of finite number of homogenous terms of the symbol of $\Delta$.

[^3]: There are two choices of spectral cut $L_{\phi}$ for the Dirac operator: in upper or lower half plane
