Soliton Fermionic number from the heat kernel expansion

TL;DR

This paper develops systematic heat kernel expansion methods to accurately compute the fractional fermion number of solitons, addressing convergence issues and enabling automated calculations in complex models.

Contribution

It introduces a refined heat kernel expansion approach and a formula for the localized eta function, improving the calculation of fermion numbers in soliton models.

Findings

Derived a systematic eta function formula for spectral asymmetry

Proposed an improved heat kernel expansion addressing convergence

Enabled automated computation of fermion numbers in multiflavour models

Abstract

We consider different methods of calculating the (fractional) fermion number of solitons based on the heat kernel expansion. We derive a formula for the localized eta function a more systematic version of the derivative expansion for spectral assymmetry and that provides a more systematic version of the derivative expansion for spectral asymmetry and compute the fermion number in a multiflavour extension of the Goldstone-Wilczek model.We also propose an improved expansionof the heat kernelthat allows the tackling ofthe convergence issues and permits an automated computation of the coefficients

Equations

$$N=-\frac{1}{2}\eta (0,H),$$

$$\eta (s,H)=\sum _{\lambda >0}{\lambda }^{-s}-\sum _{\lambda <0}(-\lambda {)}^{-s}.$$

$$\eta (s,H)=\mathrm{Tr}\left((H^2{)}^{-s/2}H/|H|\right)=\frac{1}{\Gamma \left(\frac{s+1}{2}\right)}{\int }_0^{\infty }d\tau {\tau }^{\frac{s-1}{2}}\mathrm{Tr}\left(H{e}^{-\tau H^2}\right).$$

$$\eta (s,H;\rho )=\mathrm{Tr}\left(\rho \cdot (H^2{)}^{-s/2}H/|H|\right).$$

$$H_{\rho }=H+\epsilon \rho .$$

$$\eta (s,H;\rho )=-\frac{1}{2\Gamma \left(\frac{s+1}{2}\right)}{\int }_0^{\infty }d\tau {\tau }^{\frac{s-3}{2}}{\frac{\mathrm{d}}{\mathrm{d}\epsilon }|}_{\epsilon =0}\mathrm{Tr}\left(e^{-\tau H_{\rho }^2}\right).$$

$$\eta (0,H;\rho )=\int d^{d}x\rho (x)\eta (0,H,x),\eta (0,H)=\int d^{d}x\eta (0,H;x).$$

$${\tilde{H}}_{\rho }^2:=H_{\rho }^2-M^2.$$

$$\mathrm{Tr}\,\left(Qe^{-\tau L}\right)=\sum_{k}\tau^{\frac{k-d}{2}}a_{k}\bigl{(}L;Q).$$

$$\left.\frac{\mathrm{d}}{\mathrm{d}\varepsilon}\right|_{\varepsilon=0}\,a_{k}\bigl{(}\tilde{H}^{2}_{\rho}\bigr{)}=0\quad\mbox{for}\quad k\leq d+1\,.$$

$$\eta(0,H;\rho)=-\frac{1}{2\sqrt{\pi}}\sum_{k}\Gamma\left(\frac{k-1-d}{2}\right)|M|^{1+d-k}\left.\frac{\mathrm{d}}{\mathrm{d}\varepsilon}\right|_{\varepsilon=0}\,a_{k}\bigl{(}\tilde{H}^{2}_{\rho}\bigr{)}.$$

$$L=-({\nabla }^2+E),{\nabla }_{\mu }={\partial }_{\mu }+{\omega }_{\mu }$$

$$a_{2p}(L;Q)=\frac{1}{(4\pi)^{d/2}p!}\int d^{d}x\,\mathrm{tr}\,\bigl{(}QE^{p}\bigr{)}.$$

$$\delta \eta (0,H)=-\frac{2}{\sqrt{\pi }}{a}_{d-1}(H^2,\delta H).$$

$$S=\int dtdx\bar{\psi }/D\psi ,/D=\mathrm{i}{\gamma }^{\mu }{\partial }_{\mu }-{\phi }_1-\mathrm{i}{\gamma }_{\ast }{\phi }_2,$$

$$H=-\mathrm{i}{\gamma }_{\ast }{\partial }_x+{\phi }_1{\gamma }^0+\mathrm{i}{\phi }_2{\gamma }^1.$$

$$H_{\rho}^{2}=-\partial_{x}^{2}+G^{2}+\mathrm{i}\gamma^{1}\varphi_{1}^{\prime}-\gamma^{0}\varphi_{2}^{\prime}+\varepsilon\bigl{(}2\rho G-\mathrm{i}\gamma_{*}(2\rho\partial_{1}+\rho^{\prime})\bigr{)},$$

$$G={\phi }_1{\gamma }^0+\mathrm{i}{\phi }_2.$$

$$\delta\eta(0,H)=-\frac{2}{\sqrt{\pi}}a_{0}(H^{2},\delta H)=-\frac{1}{\pi}\int dx\,\mathrm{tr}\bigl{(}\gamma^{0}\delta\varphi_{1}+\mathrm{i}\gamma^{1}\delta\varphi_{2}\bigr{)}=0,$$

$$E=(M^2-G^2)-\mathrm{i}{\gamma }^1{\phi }_1^{\prime }+{\gamma }^0{\phi }_2^{\prime }-2\rho \epsilon G,{\omega }_1=\mathrm{i}{\gamma }_{\ast }\epsilon \rho$$

$$\frac{1}{\sqrt{4\pi}l!}\mathrm{tr}\,\bigl{[}(M^{2}-G^{2})^{l}(-2\varepsilon\rho G)(-\mathrm{i}\gamma^{1}\varphi^{\prime}_{1}+\gamma^{0}\varphi^{\prime}_{2})\bigr{]}$$

$$G^2={\phi }_1^2+{\phi }_2^2+\mathrm{i}{\gamma }_{\ast }[{\phi }_1,{\phi }_2].$$

$$\frac{2\rho\varepsilon}{\sqrt{4\pi}l!}\mathrm{tr}\,\bigl{[}(z_{+}^{l}+z_{-}^{l})(\varphi_{2}\varphi^{\prime}_{1}-\varphi_{1}\varphi^{\prime}_{2})+\mathrm{i}(z_{+}^{l}-z_{-}^{l})(\varphi_{1}\varphi^{\prime}_{1}+\varphi_{2}\varphi^{\prime}_{2})\bigr{]}$$

$$z_{\pm }=M^2-y_{\pm },y_{\pm }={\phi }_1^2+{\phi }_2^2\pm \mathrm{i}[{\phi }_1,{\phi }_2].$$

$$\sum _{l=0}^{\infty }\frac{z_{\pm }^l}{M^{2l}}={\left(1-\frac{z_{\pm }}{M^2}\right)}^{-1}.$$

$$y_{\pm }=({\phi }_1\mp \mathrm{i}{\phi }_2)({\phi }_1\mp {\phi }_2{)}^{†}$$

$$\eta(0,H)=-\frac{1}{2\pi}\int dx\,\mathrm{tr}\,\bigl{[}(y_{+}^{-1}+y_{-}^{-1})(\varphi_{2}\varphi^{\prime}_{1}-\varphi_{1}\varphi^{\prime}_{2})+\mathrm{i}(y_{+}^{-1}-y_{-}^{-1})(\varphi_{1}\varphi^{\prime}_{1}+\varphi_{2}\varphi^{\prime}_{2})\bigr{]}.$$

$$\eta(0,H)=-\frac{\mathrm{i}}{2\pi}\mathrm{tr}\,\left.\bigl{[}\ln(\varphi_{1}+\mathrm{i}\varphi_{2})-\ln(\varphi_{1}-\mathrm{i}\varphi_{2})\bigr{]}\right|_{-\infty}^{+\infty}$$

$$\eta (0,H)=\frac{1}{\pi }{\left(\mathrm{arctg}\frac{{\phi }_2}{{\phi }_1}\right)|}_{-\infty }^{+\infty },$$

$$H=H_{0}+\mu\sigma_{3}=\left(\begin{array}[]{cc}\mu&-\frac{d}{dx}+\Phi(x)\\ \frac{d}{dx}+\Phi(x)&-\mu\end{array}\right)=\left(\begin{array}[]{cc}\mu&D\\ D^{\dagger}&-\mu\end{array}\right)$$

Full text

1 Introduction

More than 40 years ago Jackiw and Rebbi [1] discovered that the Fermi number of solitonic ground states may be half-integer. Later, it was demonstrated that fractional and even irrational values of the fermion number are also allowed, see [2] and the review paper [3], where the whole structure of fermion fractionization by solitons was exhibited for several systems with solitons of different dimensions and properties, namely, 1D kinks, planar vortices, and 3D magnetic monopoles. A crucial step in understanding the fermion fractionization in solitonic backgrounds was given by Niemi and Semenoff in Reference [4] where it was reinterpreted in the context of QFT anomalies.

On the theoretical side this unexpected phenomenon showed that very subtle mathematical concepts arising in topological index theory played an important r$\hat{\rm o}$le in its understanding. Together with other phenomena, anomalies, instantons, etcetera, fermionic fractionization triggered by solitons laid the foundations for the extraordinary cross fertilization between Mathematics and Quantum Physics within the QFT framework successfully developed during the last forty years.

On the experimental side Fermion fractionization meant electric charge fractionization violating the Millikan principle observed in Nature, regarding the indivisibility of the electron charge. This possibility is susceptible of being experimentally tested and a positive result on the existence of electric charges not being an integer multiple of the electron charge was found in exotic macromolecules, e.g. polyacetylene, see [5, 6]. The link between the existence of the fermion fractionization phenomenon in condensed matter physics and relativistic field theory was established in Reference [7] giving rise to one more connection between these two apparently distant areas of Physics. Also, much of the renewed interest in this problem was caused by various applications to Condensed Matter Physics, see e.g. the papers [8, 9] that treat the charge fractionization in graphene-related models. Various condensed matter applications bring up more and more complex systems [10, 11, 12] that require more and more sophisticated methods of analysis.

Most reliable methods use anomalies to compute the charge fractionization. Such methods however cannot be applied to all models. One thus needs certain kinds of the derivative expansion [2, 9]. This method is however not free of problems. On solitonic backgrounds the fields are usually of the same order as their derivatives. Thus application of the derivative expansion to Feynman diagrams is tricky, besides the combinatorial complexity if matrix-valued fields are involved. On top of this, convergence of the expansion is far from being obvious. 555We also like to mention at this point that the related trace identities method was successfully applied in [20] to the computation of fermi fractionization for magnetic monopoles.

The purpose of this work is to reconsider and improve the calculation methods of the fractional fermion number of solitons by systematically using heat kernel techniques (see, e.g., [13] for a review). We shall concentrate on combinatorial aspect, derive a new expression for the fractional fermion number in a model with matrix-valued fields, and analyze the convergence issues in an improved version of the heat kernel expansion.

Note, that for computations of the parity anomaly [4, 14], which is a close cousin of the problem considered here, the heat kernel methods proved to be very efficient [15, 16, 17, 18, 19]. The same can be said about the mass shift of solitons due to quantum corrections. The usual heat kernel expansion alone delivers the answer in the supersymmetric case [21, 22]. In non-supersymmetric case, however, certain modifications of the expansion are needed [23].

Our strategy is as follows. In the next section we derive localized expressions for the $\eta$ function and for the spectral asymmetry. This allows us to use the standard heat kernel expansion and derive an expression for the fermion number in a non-abelian (multiflavor) generalization of the Goldstone-Wilczek model. Then to address certain drawbacks of the standard heat kernel expansion we develop an improved version of the heat kernel. To set up the stage, we re-derive the Niemi-Semenoff formula for fermion number. In Section 4.3 we analyze convergence of the expansions.

2 Fermion number through the heat kernel expansion

Let us consider a theory of Dirac fermions in $d$ spatial dimensions described by a Hamiltonian $H$. The fermion number $N$ is given by the formula [3]

[TABLE]

where $\eta(0,H)$ is the value at [math] of the spectral $\eta$ function that measures the spectral asymmetry of the Dirac Hamiltonian $H$,

[TABLE]

Here $\lambda$’s are the eigenvalues of $H$. Through a Mellin transform, this formula can be converted to

[TABLE]

Note, that zero modes are excluded from the summation in (2). However, these modes can be included in the trace in the last expression in (3) since the operator $He^{-\tau H^{2}}$ vanishes on the null space of $H$.

We shall need a localized version of the formula (3). To this end we introduce a smooth function $\rho$ which plays the role of a chemical potential and define

[TABLE]

Let us define a shifted Hamiltonian

[TABLE]

so that the localized $\eta(s,H;\rho)$ function is obtained by differentiating $\eta(s,H_{\rho})$ with respect to $\varepsilon$,

[TABLE]

Smeared and unsmeared $\eta$ functions are related through the equations

[TABLE]

Let us assume that there is a natural mass gap parameter $M^{2}$ and define

[TABLE]

For any Laplace type operator $L$ and any smooth matrix-valued function $Q$ there is an asymptotic expansion of the corresponding heat trace

[TABLE]

For $Q\equiv 1$ we shall abbreviate $a_{k}(L;1):=a_{k}(L)$. If there are no boundaries, only even-numbered heat kernel coefficients $a_{k}(L)$ with $k=2p$ are non-zero. Clearly, the integral in (6) converges at $s=0$ if

[TABLE]

Let us assume that this condition is satisfied. Then, we can substitute (9) in (6), put there $s=0$ and integrate over $t$ to obtain

[TABLE]

General operator of Laplace type can be written as

[TABLE]

with a suitable choice of the connection one-form $\omega$ and a matrix-valued potential $E$. In a flat space, the heat kernel coefficients are integrals of polynomials constructed from $E$, $\Omega_{\mu\nu}=[\nabla_{\mu},\nabla_{\nu}]$ and their covariant derivatives. The general form of all heat kernel coefficients is not known. However, on solitonic backgrounds just a few structures containing a limited number of derivatives contribute to the spectral asymmetry. Such terms in $a_{k}(L)$ can be found for all $k$. For example, if one is interested in terms containing only $E$ and no derivatives the expression for heat kernel coefficients looks particularly simple,

[TABLE]

Let us consider a smooth localized variation $\delta H$ of the Hamiltonian $H$. Suppose that $\delta H$ is zero order, i.e. it does not contain derivatives. Then, one can show [24] that

[TABLE]

Strictly speaking, this formula is valid only if under the variation $\delta H$ eigenvalues of the Hamiltonian do not cross the zero value. In physics’ literature Eq. (14) is usually used without any restrictions (and without running into troubles because of this), see e.g. [15]. Here, we assume that the background is ”generic”, i.e. the number of zero modes is stable under infinitesimal localized variations.

The main novelty of the approach suggested above is the use of a localizing function $\rho$. This enables us to utilize standard expressions for the heat kernel coefficients that were derived under an assumption that all integrations by parts were possible. Some other benefits of introducing $\rho$ will be discussed below, see Sec. 4.3.

3 A non-abelian extension of the Goldstone-Wilczek model

Let us consider a model in $1+1$ dimensions described by the action

[TABLE]

where $\gamma_{*}=\gamma^{0}\gamma^{1}$ is the chirality matrix in $1+1$ dimensions, while $\varphi_{1}$ and $\varphi_{2}$ are scalar fields that are hermitian matrices in the ”flavor” space. In the abelian one-flavor case the action (15) is the model considered by Goldstone and Wilczek [2].

The Hamiltonian for this problem reads

[TABLE]

This yields the following relation for shifted Hamiltonian (5)

[TABLE]

where prime denotes derivative with respect to $x$ and

[TABLE]

We omitted the terms with $\varepsilon^{2}$ since they do not contribute to the spectral asymmetry.

Take infinitesimal localized variations $\delta\varphi_{1}$ and $\delta\varphi_{2}$ so that $\delta H=\gamma^{0}\delta\varphi_{1}+\mathrm{i}\gamma^{1}\delta\varphi_{2}$. From (14) one computes the corresponding variation of spectral asymmetry,

[TABLE]

due to the trace of gamma matrices. We conclude that $\eta(0,H)$ is a topological invariant. I.e., it depends on the asymptotic values of $\varphi_{1}$ and $\varphi_{2}$ only.

Let us introduce a mass gap parameter $M^{2}$ whose value we do no specify now. The operator $\tilde{H}_{\rho}^{2}$ (see eq. (8)) can be brought to the canonical form (12) with

[TABLE]

We consider only solitonic backgrounds, meaning that $\varphi_{1}$ and $\varphi_{2}$ go fast to their asymptotic values as $x\to\pm\infty$. Thus $\eta(0,H)$ depends on these asymptotic values and not on the derivatives of $\varphi_{1,2}$. Consequently, in the expansion (11) we have to keep only the terms that contain at most a single first derivative of $\varphi_{1}$ or $\varphi_{2}$. This excludes the terms in the heat kernel expansion containing derivatives of $E$. Consequently, only the terms listed in (13) contribute. Besides, in the powers of $E$ one should select the terms that contain $\rho$ only once. Taken together, all relevant contributions to $a_{2(l+1)}$ read

[TABLE]

Note, that

[TABLE]

After taking the trace over $\gamma$-matrices in (21) we obtain

[TABLE]

with

[TABLE]

As a next step, one has to substitute the expression (23) for $a_{2(l+1)}$ in (11) and compute the sum with the help of

[TABLE]

This sum is convergent if $\|z_{\pm}\|=\|M^{2}-y_{\pm}\|<M^{2}$. The matrices

[TABLE]

are non-negative. Thus, the problem with convergence appears only if $\varphi_{1}$ and $\varphi_{2}$ vanish simultaneously in some direction of the flavor space either in asymptotics or at some points of the bulk (meaning just a breakdown of the derivative expansion). If the problem appears at the asymptotics, one cannot do much about it. We assume that this does not happen. However, since $\eta(0,H)$ is topological, we may use the freedom of smooth variations of background fields to shift the roots of $\varphi_{1}$ away from the roots of $\varphi_{2}$ in the bulk thus making the sum (24) convergent everywhere. Thus,

[TABLE]

As expected, this expression can be converted to a boundary term

[TABLE]

This result is new.

If $\varphi_{1}$ and $\varphi_{2}$ are one-component fields rather than matrices, the expression (27) can be simplified,

[TABLE]

which is the Goldstone and Wilczek result [2].

We see that the method proposed works well and relatively easily. However, there are some points that have to be clarified or improved.

The heat kernel expansion reflects the behavior of heat trace at small $\tau$. To put the integration over $\tau$ on firm grounds one needs to understand the behavior at large $\tau$ as well. 2. 2.

It is hard to control the rate of convergence of the expansion (26), especially since the individual terms in (26) do not give total derivatives after the variation with respect to $\rho$. 3. 3.

In the example above we had to sum up the simplest contributions $E^{k}$ to the heat kernel coefficients. In more complicated cases, especially in higher dimensions, more complicated invariants may become relevant. Thus, it is useful to have a method to find all higher heat kernel coefficients for some types of the operators.

These points will be addressed below.

4 The kink fermion number through (generalized) heat kernel expansion

4.1 The spectrum of fluctuations

In this Section, we demonstrate how the heat kernel expansion can be improved by taking into account some global information on the background fields. Below, we describe the properties of Dirac spectrum on kink backgrounds (cf [3]). We take

[TABLE]

where $\mu$ is the mass, and $\Phi(x)$ refers to a yet unspecified topological kink with a smooth behavior. One can think of $\Phi=\nu\mathrm{tanh}(x)$, though we shall not restrict ourselves to this simplest example. The Hamiltonian (29) is a particular case of (16). The eigenspinors of (29) satisfy the Dirac equation:

[TABLE]

Finding the solution of the ODE system (30) if $|E|\neq|\mu|$ is equivalent to solving one spectral Schr$\ddot{\rm o}$dinger problem in two ways:

[TABLE]

where the Schr$\ddot{\rm o}$dinger operators are defined in terms of the topological kink configuration:

[TABLE]

In (32)-(33) we have taken into account that kink configurations have very definite asymptotic behavior:

[TABLE]

where $v$ is a constant. Thus, $V_{1}(x)=\Phi^{2}(x)-v^{2}-\frac{d\Phi}{dx}$ and $V_{2}(x)=\Phi^{2}(x)-v^{2}+\frac{d\Phi}{dx}$ also tend to zero at $x=\pm\infty$:

[TABLE]

We summarize the spectrum of the Dirac operator:

If $E=\mu$ then $D^{\dagger}u_{1}(x)=2\mu u_{2}(x)$ and $Du_{2}(x)=0$. Henceforth $DD^{\dagger}u_{1}(x)=0$, which means that $u_{1}$ is a zero mode of the operator $DD^{\dagger}$ if it is normalizable. 2. 2.

If $E=-\mu$ then $D^{\dagger}u_{1}(x)=0$ and $Du_{2}(x)=-2\mu u_{1}(x)$. Henceforth $D^{\dagger}Du_{2}(x)=0$, which means that $u_{2}(x)$ in this case is a zero mode of the operator $D^{\dagger}D$ if normalizable. We stress that only the spinor annihilated by either $D^{\dagger}$ or by $D$ is normalizable. Therefore, generically, for one-scalar field theory there will only be one normalizable zero mode of the Dirac operator governing the spinorial fluctuations around one kink. 3. 3.

If $E\neq\pm\mu$ the first component of the spinor $u_{1}(x)$ is obtained by solving the spectral problem

[TABLE]

whereas the second component $u_{2}$ is fixed in (31)-(left). Clearly this solution for $u_{2}$ is an eigenfunction of the spectral problem

[TABLE]

because the intertwining of the $DD^{\dagger}$ and $D^{\dagger}D$ operators which , except the zero mode, are isospectral:

[TABLE]

Denoting as $\widetilde{u}_{1}$ a normalized eigenfunction of the spectral problem (34) $DD^{\dagger}u_{1}(x)=\omega^{2}u_{1}(x)$ one easily checks that the spinor

[TABLE]

is a normalized eigenfunction of the Dirac Hamiltonian, since via a shrewd partial integration one finds that: $\int dx\,u^{\dagger}(x)u(x)=\int dx\,\widetilde{u}_{1}^{*}(x)\widetilde{u}_{1}(x)$.

The function

[TABLE]

is obviously an eigenfunction of the spectral problem (35) $D^{\dagger}Du_{2}(x)=(E^{2}-\mu^{2})u_{2}(x)$ because it is proportional to $u_{2}$ and, after another partial integration, one checks that its norm is one: $\int dx\,\widetilde{u}_{2}^{*}\widetilde{u}_{2}=\int dx\,\widetilde{u}_{1}^{*}\widetilde{u}_{1}$.

We thus may characterize the eigenspinors of the Dirac Hamiltonian as

[TABLE]

instead of (36). To distinguish among the different eigenspinors we remark that $\widetilde{u}_{1}(x,E)$ and $\widetilde{u}_{2}(x,E)$ are, respectively, non-null normalized eigenfunctions of $DD^{\dagger}$ and $D^{\dagger}D$. They both are characterized by either a real vector number $k\in{\mathbb{R}}$ if one search for scattering states and/or purely imaginary momenta, $k=i\kappa_{j}$ with $\kappa_{j}>0$ and $j=1,2,\cdots,N_{B}$, if $N_{B}$ bound states exist. The corresponding energies are: $E^{(\pm)}=\pm\sqrt{k^{2}+v^{2}+\mu^{2}}$, and $E_{j}^{(\pm)}=\pm\sqrt{v^{2}+\mu^{2}-\kappa_{j}^{2}}$. Besides there is a zero mode $\kappa_{0}=v$ of either $DD^{\dagger}$ or $D^{\dagger}D$ which contributes to the spectral asymmetry selecting either $E_{\mu}=\mu$ or $E_{-\mu}=-\mu$. We denote by $B_{1}=\{\widetilde{u}_{1k}\}_{k\in I}$ the complete orthogonal set in $L^{2}({\mathbb{R}})$ formed by the eigenfunctions of $DD^{\dagger}$. $I$ is a set of indices that symbolically label the scattering states, the bound states, and, eventually, one zero mode. Also we dispose of $B_{2}=\{u_{2k}\}_{k\in I}$ as a second complete orthogonal basis in $L^{2}({\mathbb{R}})$ formed by the eigenfunctions of the operator $D^{\dagger}D$, except for the possible a priori presence of a zero mode.

Most of the properties established here are also valid in higher dimension [9].

4.2 Niemi-Semenoff formula and the spectral heat trace

The fermionic quantum field in the Schr$\ddot{\rm o}$dinger picture is expanded in terms of the Dirac eigenspinors:

[TABLE]

Here $[dk]$ denotes an integration measure on the spectrum of the Dirac operator. The fermionic creation and annihilation operators of particle and antiparticles $\hat{b}_{\pm}(k)$, $\hat{b}^{\dagger}_{\pm}(k)$ satisfies the anticommutation relations:

[TABLE]

All other anticommutators between these operators are zero. The one-particle states may be therefore only occupied by one Fermion or unoccupied and one state is distinguished as the ground state: \big{|}0\big{\rangle}, characterized as the vacuum state where all the one-particle states are unoccupied, \hat{b}_{+}(k)\big{|}0\big{\rangle}=\hat{b}_{-}(k)\big{|}0\big{\rangle}=0, $\forall k$. From the expectation value at the vacuum state of the normal ordered Fermi density, one has

[TABLE]

This calculation has been performed using the anticommutation relations. Further use of the anticommutation relations allows us to obtain the fermionic number of the kink ground state:

[TABLE]

Formula (38) permits to rewrite (42) in the form:

[TABLE]

Each “integral” in this formula for $N$ is ultraviolet divergent, $v^{2}$ and $\mu^{2}$, however, prevent infrared divergences. To control the ultraviolet divergences we shall use spectral zeta function/heat trace methods associated with the operators $DD^{\dagger}$ and $D^{\dagger}D$.

We put at work this strategy in a more general context. Translate the spectral zeta function corresponding to the second order differential operator666 Note that this differential operator is a 1D member of the family of Laplace type operators (12).

[TABLE]

to its spectral heat trace via the Mellin’s transform:

[TABLE]

To tame ultraviolet divergences we subtracted the contribution from ${\cal K}_{0}$, which is obtained from ${\cal K}$ by putting $V=0$. The operator ${\cal K}_{0}$ is the same for $DD^{\dagger}$ and $D^{\dagger}D$, so that the correction term is canceled in $N$. We denote $h_{\cal K}(\tau)={\rm Tr}\left[e^{-\tau\,{\cal K}}-e^{-\tau\,{\cal K}_{0}}\right]$. The kink fermion number as given in formula (43) may be written in terms of the spectral heat functions of the paired Schr$\ddot{\rm o}$dinger operators as follows:

[TABLE]

At this point we use the asymptotic expansion of the kink spectral trace, see [25],

[TABLE]

Here, $N_{\rm zm}^{\cal K}$ is the number zero modes in the spectrum of ${\cal K}$. The error function in equation (45) is introduced to guarantee the right behavior of the heat trace at low temperature, $\tau\to\infty$, when zero modes are present, whereas the high-temperature asymptotics is preserved. Any regular function with the same limits as the error function at $\tau\to\infty$ and $\tau\to 0$ could do the job but our choice optimizes the solution of the recurrence relations, see [25]. In addition, this procedure implies a natural relation between the Seeley coefficients of two intertwined Laplace type operators Darboux factorizable in terms of one ladder operator and its adjoint. The coefficients $c_{n}({\cal K})$ are obtained by solving the recurrence relations arising in the power series expansion solution of the ${\cal K}$-heat equation in the so-called modified Gilkey-de Witt approach, see [25]. The coefficients $c_{n}$ play a role similar to that of $a_{2n}$ is Sec. 2, apart from normalization and modification of the expansion.

Application to the kink fermion number formula gives

[TABLE]

or, after performing the asymptotic expansion and subsequently the Mellin transform we obtain

[TABLE]

Finally, taking the limit in the previous expression, one obtains the formula

[TABLE]

which is closely related to the Niemi-Semenoff formula.

In the case of kink type operators where $DD^{\dagger}$ and $D^{\dagger}D$ are intertwined operators it may be shown that $c_{n}(DD^{\dagger})=c_{n}(D^{\dagger}D)$, $\forall n$ when the modified Gilkey-de Witt approach is applied. We recall that this scheme is adapted to the presence of zero modes in the spectrum of a operator. We shall prove this fact for a pair of intertwined operators having an spectral structure identical to the spectral structures of $DD^{\dagger}$ and $D^{\dagger}D$ in the kink background. For the sake of simplicity it will be assumed that the operator $DD^{\dagger}$ encompasses a zero mode whereas $D^{\dagger}D$ has a positive spectrum. Otherwise, the following result remains true swapping the r$\hat{\rm o}$le of the first-order operators $D$ and $D^{\dagger}$.

Lemma: Let $D=-\frac{d}{dx}+\Phi(x)$ be a first-order differential operator such that $\lim_{x\rightarrow\pm\infty}\Phi(x)=\pm v$ and $\psi_{0}(x)=e^{-\int_{x}\Phi(x)dx}$ is a square-integrable function annihilated by $D^{\dagger}$. The Seeley coefficients of the operators $DD^{\dagger}$ and $D^{\dagger}D$ of the spectral $DD^{\dagger}$-and $D^{\dagger}D$-heat functions in the improved Gilkey-de Witt expansion coincide: $c_{n}(D^{\dagger}D)=c_{n}(DD^{\dagger})$ $\forall n\in\mathbb{N}$.

Proof: We shall use the notation given by (32) and (33). Let $\{f_{k}(x)\}$ and $\{g_{k}(x)\}$ be complete bases of orthonormal eigenfunctions of the operators $DD^{\dagger}$ and $D^{\dagger}D$, respectively. We summarize the spectra of these operators without taking into account yet that they are intertwined and, thus, except for the zero modes, they are isospectral

[TABLE]

In terms of the associated eigenfunctions the heat kernels of these operators are expressed as the (formal) series:

[TABLE]

where $f_{0}(x)$ is the zero mode of the operator $DD^{\dagger}$.

The action of the “ladder” operators $D^{\dagger}$ and $D$ respectively on the spectral equations for $DD^{\dagger}$ and $D^{\dagger}D$ easily show isospectrality of the positive spectra: $N_{1}=N_{2}$ and $\omega_{1i}^{2}=\omega_{2i}^{2}=\omega_{i}^{2}$ with $i\geq 1$. In fact $g_{k}(x)=\frac{1}{\omega(k)}D^{\dagger}f_{k}(x)$ are normalized eigenfunctions of the spectral the operator $D^{\dagger}D$. A caveat: since the mapping is via a differential operator the spectral densities of the scattering states are different. The heat equation kernel of the $D^{\dagger}D$ operators can be now written in terms of the eigenfunctions of the $DD^{\dagger}$ operator as

[TABLE]

The trace of this kernel $h_{D^{\dagger}D}(\tau)=\int dx\,K_{D^{\dagger}D}(x,x;\tau)$ may be written thus

[TABLE]

Moreover $\rho_{DD^{\dagger}}(k)=\rho_{D^{\dagger}D}(k)-\frac{1}{\pi}\frac{v}{k^{2}+v^{2}}$, see [26], which implies

[TABLE]

that is,

[TABLE]

Plugging the modified Gilkey-de Witt expansion of the heat function for these operators

[TABLE]

into the previous expression

[TABLE]

leads to the result

[TABLE]

as claimed $\Box$.

In sum, the Niemi-Semenoff formula for this type of Dirac operators hold:

[TABLE]

where $N_{zm}^{DD^{\dagger}}-N_{zm}^{D^{\dagger}D}=\pm 1$.

The fermion number computed at the kink ground state, with all the fermionic fluctuations around the kink being unoccupied, is an angle rather than an integer number- in fact zero- as it would be in the usual vacuum state. In the $\mu=0$ limit the Jackiw-Rebbi half-integer fermion number is recovered.

4.3 Fermion number from the spectral eta function

Here we evaluate the fermion number by using the expression (6) for spectral asymmetry. Let us write

[TABLE]

where $V(x)$ and $Q(x)$ are matrix-valued functions that vanish at $x\to\pm\infty$, $v^{2}$ is a diagonal matrix representing asymptotic values of the potential, $v^{2}={\rm diag}(v_{1}^{2},v_{2}^{2})$, and $\mathbb{I}$ is a unit matrix. All matrices are $2\times 2$. An extension to matrices of arbitrary size is straightforward.

Let us take an expansion

[TABLE]

which essentially coincides with the one that we used in Sec. 2 after obvious changes in the notations. In the expansion (4.3), however, the mass gap $v^{2}$ need not be proportional to a unit matrix. After plugging this expansion in (6), integrating over $\tau$ and passing to $s=0$, we obtain

[TABLE]

where we denoted

[TABLE]

To extract the Seeley coefficients $[c_{n}(H_{\rho}^{2})]_{ii}$ we first solve the recurrence relations777See Reference [27], Section §.3, to find the notational details and the developments of the main steps in this derivation. for the derivatives of the diagonal densities ${}^{(k)}C_{n}(x)=\lim_{y\to x}\frac{\partial^{k}}{\partial x^{k}}c_{n}(x,y)$

[TABLE]

starting with the initial conditions ${}^{(k)}C_{0}(x)=\delta_{0k}\,\mathbb{I}$. The Seeley coefficients are thus obtained from the diagonal densities via integration

[TABLE]

and the fermionic number is estimated as the $n_{\rm max}\to\infty$ limit of the partial sums:

[TABLE]

If the $v_{i}$ are all equal we have

[TABLE]

Recall that the square of the Niemi-Semenoff-Dirac Hamiltonian (29) is

[TABLE]

where $\Phi(x)$ is the solitonic background which complies with the asymptotic behavior

[TABLE]

Therefore, we have $v^{2}={\rm diag}(\nu^{2}+\mu^{2},\nu^{2}+\mu^{2})$,

[TABLE]

and

[TABLE]

Plugging these expressions into the recurrence relations the solutions for the lower Seeley coefficients provided by Mathematica read:

[TABLE]

For example,

[TABLE]

where we have used that $\frac{\partial}{\partial x}[\Phi(x)(\Phi^{\prime}(x))^{2}]=(\Phi^{\prime}(x))^{3}+2\Phi(x)\Phi^{\prime}(x)\Phi^{\prime\prime}(x)$ and therefore

[TABLE]

These results lead to the general formula and we obtain for the Fermi number:

[TABLE]

In these computations a very important, but expected point, must be stressed: The Seeley coefficients are independent of the specific dependence on $x$ of the Kink profile $\Phi$.

Denoting the partial sum of the series (47) as

[TABLE]

we find that this partial sum may be written in terms of special functions as:

[TABLE]

We then plot these partial sums

[TABLE]

as functions of $z=\frac{\mu}{\nu}=\frac{m}{ga}$, see Figure 1.

This Figure allows us to grasp an analytical understanding of the application of the Gilkey-de Witt asymptotic method to this kind of calculations. The behavior of $N^{(n_{\rm max})}$ is displayed in Figure 1. We observe that when $\mu\neq 0$ the partial sums converge to the exact response (shown as the lower blue solid line in the Figure). For small values of $m$ this convergence is very slow. For $\mu=0$ (which involves the presence of zero modes in the operator $H$) the response of the formula (47) is zero in contrast with the exact value $N=\pm\frac{1}{2}$. We can impose that when zero modes are present the exact response can be obtained by demanding the continuity of $N(H)$ as a function of $\mu$ at $\mu=0$.

Let us stress some features the method presented above:

This method singles out a value of the mass gap parameter that is consistent with global properties of the operator and that ensures convergence of the perturbation series. 2. 2.

The terms that improve the large $\tau$ behavior of the heat kernel, see (45) could have been included in (4.3). As one can easily see, such terms vanish upon variation with respect to $\rho$. This observation explains a posteriori why the expansion in Sections 2 and 3 worked well without any improvement. 3. 3.

This method is actually quite general. An extension to the abelian Goldstone-Wilczek model is quite straightforward. (We do not present it here). In other generalizations one depends just on the computer power.

5 Conclusions

In this paper, we suggested a method to compute the fermion number of solitons based on the heat kernel expansion. A crucial step was an introduction of a localizing function $\rho$ in the $\eta$ function that allowed us to use the standard heat kernel coefficients and simplified combinatorics of the problem. With the use of this method we obtain an expression for the fermion number in a multiflavor extension of the Goldstone-Wilczek model.

This method, being just a more systematic version of the usual derivative expansion, shares some drawbacks of the latter. It does not take into account the large distance and short proper time behavior of the heat kernel. Thus, it is not easy to explore the convergence. However, there is an improvement of the expansion [25] which allows to address this problem. With the help of this improved expansion we first reconfirmed the Niemi-Semenoff formula and then demonstrated the convergence of heat kernel series for fermion number.

We would like to stress, that the methods presented here are extendable to more complicated and higher-dimensional models. We are going to consider these models in a future work.

Acknowledgments

This work was supported in parts by the São Paulo Research Foundation (FAPESP), projects 2016/03319-6 and 2017/50294-1 (SPRINT), by the grants 303807/2016-4 and 428951/2018-0 of CNPq, by the RFBR project 18-02-00149-a and by the Tomsk State University Competitiveness Improvement Program. AAI and JMG also thank to the JCyL for partially supporting their research under Grant BU 229P18