# Polarized heat current generated by quantum pumping in two-dimensional topological insulators

## Abstract

We consider transport properties of a two dimensional topological insulator in a double quantum point contact geometry in presence of a time-dependent external field. In the proposed setup an external gate is placed above a single constriction and it couples only with electrons belonging to the top edge. This asymmetric configuration and the presence of an ac signal allow for a quantum pumping mechanism, which, in turn, can generate finite heat and charge current in an unbiased device configuration. A microscopic model for the coupling with the external time-dependent gate potential is developed and the induced finite heat and charge current are investigated. We demonstrate that in the non-interacting case, heat flow is associated with a single spin component, due to the helical nature of the edge states, and therefore a finite and polarized heat current is obtained in this configuration. The presence of e-e interchannel interactions strongly affects the current signal, lowering the degree of polarization of the system. Finally, we also show that separate heat and charge flows can be achieved, varying the amplitude of the external gate.

## I Introduction

The fast development in nanotechnologies has spurred the investigation of quantum effects in electronic devices at submicron scale (1); (2); (3); (4); (5). In past years, the attention has been devoted mostly to the study of charge transport and related phenomena in electronic nanocircuits (6); (7); (8); (9); (10); (11). Here spectacular effects have been predicted and also experimentally observed, including quantum interference pattern (1), single electron injection and control in quantum conductors (7); (8); (12); (13) and coherent behavior in hybrid devices (4); (5); (9); (14); (15); (16).

Conversely, the question how thermodynamic aspects such as heat transport, power conversion and energy exchange work at the nanoscale has posed new interesting challenges (17); (18); (14); (19); (20); (21); (22); (23). A precise control and manipulation of heat flows in quantum conductors can lead to interesting applications, i.e., new logic devices based purely on thermal transport. Towards this goal, pioneering experiments have achieved the phase control of coherent thermal transport in hybrid quantum systems (5); (14); (24); (25); (26); (27).

In this respect, the possibility to tune and control some physical parameters by means of external fields and time-dependent potential is of great interest (28); (29); (30); (31); (32); (33). Indeed, periodically driven system can be immediately associated to thermal machines or heat engines which are based on cyclic operations. Different papers focused on the effects induced by an external driving field both in closed and open quantum systems.
For example, the influence of an ac field on the efficiency of thermoelectric setups or heat pumping mechanism on quantum dot based devices have been investigated (28); (29); (31); (34); (35); (36).

Dynamical aspects of time-dependent energy and heat flows have been addressed in order to understand fundamental thermodynamical problems such as energy exchange, work distribution and entropy production (see reviews (20); (31) and references therein).

Interestingly, locally applied time-dependent voltages can lead to heat pumping mechanisms in quantum systems (9); (31); (37); (38); (39); (40); (41).
The peculiar characteristic of heat pumps is that a direct heat current is generated by a purely ac drive which acts against some present thermal gradients.
Moreover, it can also generate dc finite current signals even in a purely equilibrium situation, in absence of any thermal or voltage bias (31); (42); (43).

Recently, a lot of attention has also been put on the interplay between heat and spin transport properties, in the emerging research field of spin caloritronics (44); (45); (46). Thermally driven spin sources are very interesting for the improvement of spintronic devices, which are believed to be more efficient with respect to their electronic counterparts (14); (20); (47).
Striking experimental results in this field have extended thermoelectric concept to spin transport by the observation of spin-Seebeck (48) and spin-Peltier (49) effects in magnetic systems. Intriguingly, spin-dependent effective temperatures of electrons were observed in nanopillar spin valves powered by a heat source, thus showing that electron heat conduction can depend on the spin degree of freedom (50); (51).

Among all mesoscopic systems for the investigation of spin transport properties, a relevant role is certainly played by two-dimensional topological insulators (2dTIs)(52); (53); (54); (55). These system are constituted by a wide bulk gap and charge and energy transport are mediated by the presence of edge states.
These are constituted by two counter-propagating electronic states with opposite spin polarization (spin-momentum locking), and are also called helical edge states (52); (56). Moreover, as long as time-reversal symmetry is preserved, electrons flowing along these edge states are topologically protected against backscattering and transport occurs in the so-called ballistic regime (52); (53); (54).
Experimental evidences of 2dTIs were reported in CdTe/HgTe (57) and InAs/GaSb (58); (59) quantum wells and have been predicted for a large class of new materials which relies on the presence of strong spin-orbit coupling (60); (61); (62).
These systems constitute an interesting playground to study coherent heat transport in which the spin degree of freedom can play a non-trivial role (31); (63).
Moreover, the unavoidable presence of e-e interactions in helical edge states (52); (64); (65); (66); (67) can dramatically affect the dynamics and have important consequences on thermal transport properties (63).

In this work we focus exactly on these issues. We consider a 2dTI device in a double quantum point contact geometry (68); (69); (70); (71); (72); (73). In this setup tunneling of electrons is allowed by the presence of constrictions and can lead to quantum interference effects.
We also introduce a model to couple electrons of the top edge to a local time-dependent gate potential located above a single QPC. Our aim is to investigate the generation of a direct heat current by the time-dependent drive along the edge of a 2dTI. We demonstrate that the presence of an ac field can generate dc finite current signals in absence of thermal or voltage bias, relying on a quantum pumping mechanism.
We inspect the interplay with the peculiar spin properties of helical edge states, demonstrating that in the proposed setup a polarized heat current can be enstaured, carried solely by a single spin species.
Finally, the presence of e-e interactions can strongly modify the behavior of the heat pumped current in this device. Interestingly, it turns out that by looking at the resulting oscillating patterns one can identify the presence of e-e interactions in the system.

The paper is organized as follows. In Sec. II, we present the setup and describe the various contributions to the Hamiltonian of the system. In Sec. III, we define charge and heat currents and evaluate their average values. Sec. IV is devoted to the discussion of our main results. Finally, in Sec. V, we draw our conclusions.

## Ii Model and setup

We consider a 2dTI connected to two reservoirs at equilibrium kept at the same chemical potential and temperature .

Here, electrons can propagate along the edges and, in the absence of interactions, have well-defined chirality of right () and left () moving particles (52); (53). Moreover, the spin-momentum locking property of 2dTI constrains the direction of propagation and the spin projection of electrons: in the top edge the spin of -electrons is and the spin of -electrons is and viceversa for the bottom edge. Top and bottom edges are separated by macroscopic distances and inter edge tunneling is only possible in presence of some constrictions. Here, we consider a double quantum point contact (QPC) geometry (70); (71); (72), where the two constrictions are placed at , as schematically shown in Fig. 1. Electrons in the top edge are capacitively coupled to an external gate potential , confined in a region , see Fig. 1. The presence of a time-dependent gate potential, assumed to be periodic in time for simplicity, is crucial in order to generate dc finite current signals, without any external dc bias (31); (41). The total system can be described by the following Hamiltonian

(1) |

where describes the interacting edge states of the 2dTI. The free contribution associated to fermionic channels can be written as

(2) |

with (from now on )

(3) |

where and are the Fermi velocity and Fermi momentum respectively, , and the corresponding electron creation operator. In the following, electron interactions on each edge are treated assuming the breaking of galileian invariance, i.e., in the absence of intra channel interactions (3); (56); (4); (74); (75), and the corresponding term reads

(4) |

where is the coupling constant describing the interchannel interaction strength (for sake of convenience, we have adopted the same label as in the Luttinger liquid language(3)). In the above equation the density operator is denoted as . Following standard bosonization prescriptions(3); (4); (74) (see Appendix A for more details), can be written in a diagonal form in terms of chiral bosonic field operators as

(5) |

where represents the renormalized velocity, expressed in terms of the dimensionless interaction parameter (75) and identifies the top/bottom edge respectively. It is worth to note that in this work we restrict the analysis to the range of interaction strengths given by , for which it is possible to show that single electron tunneling is the dominant process (in the renormalization group sense), both in the single and double QPC geometry(56); (71); (72); (63).

Tunneling events between the two edges are allowed in the double QPC geometry under consideration. For sake of simplicity, here we focus on spin-preserving tunneling events described by (68); (69); (67); (70)

(6) |

with the space dependent amplitude describing the double constrictions depicted in Fig.1 (here represents the constant tunneling amplitude).An analogous analysis can be carried out considering spin-flipping tunneling at the QPCs: however, we expect this contribution to be smaller with respect to the spin-preserving one(68); (76).

Finally, we assume that the external gate potential is capacitively coupled to the electron density of the top edge via

(7) |

with the time-dependent potential

(8) |

Here, are linked to the dimension of the gate, see Fig.1, and is a periodic time-dependent drive, with period . In particular, we will consider the gate located around , i.e., () and , as shown in Fig. 1. Finally we note that the chosen configuration, with the external gate located around a single QPC at , naturally guarantees the required asymmetry needed for a quantum pumping mechanism (9); (40); (41).

## Iii Average pumped currents

We are interested in the study of transport properties in the setup introduced above. We recall that the 2dTI is connected to two reservoirs, which we assume to be at equilibrium at the same chemical potential and same temperature . Nevertheless finite dc current signals can be generated due to the presence of the time-dependent gate potential. Indeed, as already mentioned, the ac field supplied by the external gate above one single QPC can lead to a quantum pumping mechanism, which, in turn, generates dc finite currents.
Here, we investigate this aspect, by studying the pumped currents which flow through the 2dTI.
In particular we will consider the generation of heat current and charge current , focusing on their dc components, extracted from the finite temperature current contributions by further averaging the signals over one period of the gate potential .

It is possible to show, see Appendix A, that in the dc limit the only possible contributions to pumped currents are due to backscattering. The latter are present whenever the edges are coupled and tunneling events occur, i.e., in the proximity of the two QPCs.
Therefore, in this section we focus on the evaluation of the backscattering contribution to the average currents (we thus introduce a related index ).
Backscattering charge current can be defined as

(9) |

where is the electron charge () and we have identified the particle current per spin component, which is given by

(10) |

with the particle number operator . Using Eq. (6) we can thus write

(11) |

Analogously, the backscattering heat current flowing between the two reservoirs can be written as the sum of the heat current per spin direction

(12) |

with

(13) |

where . The expression in Eq. (13), related to the heat current contribution per spin component, in terms of fermionic fields is

(14) |

It is worth to underline that Eq. (12) corresponds to the net amount of heat flow exchanged between -channel and -channel due to backscattering processes and it does not give any information on the heat which locally enters into each reservoir separately.

Using standard perturbative approach (4); (70); (63) in the tunneling Hamiltonian, the average pumped currents can be calculated.
At the lowest order in the tunneling one has the finite temperature currents (here for particle-number and heat contributions)

(15) |

where the average is taken with respect to .

In the above equation we have introduced the index to indicate the time evolution of operators evaluated in absence of tunneling, i.e. with respect to for the top edge and to for the bottom one.
We remind the reader that in our model the external gate potential couples only to top edge electrons, see Fig. 1.
In order to evaluate the expression in Eq. (15), it is sufficient to know the time evolution of fermionic operators, obtained from . In particular we are interested only in their time-evolution close to the QPC positions at , since we are assuming local tunneling with a - like spatial shape for the QPCs centered around .
The explicit derivation of the time-evolution of these operators is reported in Appendix B for the case of a sinusoidal voltage , and here we quote only the relevant results:

(16) | ||||

(17) |

in the region close to the QPC1 (right) at and

(18) | ||||

(19) |

in the region around the QPC2 (left) located at . In the above equations we have denoted with the time-evolution of operators with respect to the bare , without the external gate. The amplitudes and phases introduced just above are given by

(20) | ||||

(21) | ||||

(22) | ||||

(23) |

where the amplitudes have no spatial dependence, contrary to what happens for which is affected by the gate voltage . The coefficients and are connected to the geometry of the setup and are given by

(24a) | ||||

(24b) |

One can immediately note that in the non-interacting case, i.e. , the amplitude vanishes. This will have significant implications on the pumped charge and heat currents as we will discuss below.

We recall that the tunneling is local at the QPC positions, with amplitude , and we introduce the following series representation(77); (78)

(25) |

where is the Bessel function of -th order. Inserting this relation into Eq. (15) one can obtain the expressions for the average pumped currents. We stress that the only contributions come from positions located at the QPCs . We thus introduce the simplified notations

(26) |

and we obtain the dc components:

(27) |

for the average particle current per spin component, and

(28) |

for the contribution associated to the heat current per spin component. Notice that the contributions in the infinite sums cancel out and the series expressions present in Eqs.(27)-(28) start from . In the above equations we have introduced the usual short-distance cut-off of HLL (4), and , where is the bosonic correlation function (4); (10); (79); (80)

(29) |

Here, the energy cut-off, which represents the highest energy scale of the problem. Note that this quantity can be related to the energy gap between bulk conduction and valence bands of the 2dTI. The dimensionless parameter depends on the interaction strength and reads

(30) |

Performing the integrals present in Eqs. (27)-(28) (see Refs. (71) and Appendix C for details) the average pumped currents can be written as

(31) | ||||

(32) |

where and and are modulating functions (see Appendix C for the explicit expressions) and is the Fourier representation in energy domain of . In passing we note that and do not depend on the voltage amplitude , but this dependence is crucially present in the argument of the Bessel functions.
Before discussing the behavior of the average pumped currents some comments are in order.
First of all, contrary to previous works(70); (71), here time-dependent tunneling amplitudes and phase shifts of the two QPCs are not introduced a priori in the system. Indeed, the solution of the microscopic model of the coupling to the external gate have lead to the time-dependence of the tunneling amplitude and to the microscopic derivation of the phase shifts and .
One can note that dc finite current signals strongly rely on quantum interference effects (14); (63); (73), induced by the presence of the two constrictions in this geometry. Indeed, both currents in Eqs. (31)-(32) vanish in the limit .
From a physical point of view, each contributions in the sum over account for tunneling processes in which electrons absorb or emit photons. Therefore the tunneling amplitudes, for each QPCs, associated to these processes result weightened by the corresponding Bessel function present in Eq. (25).

We also note that quantum interference effects are suppressed if the tunneling amplitude through the left () QPC or through the right () QPC vanish , resulting in a vanishing pumped current signals. Finally we underline that the phases and acquired by electrons, have to be different in order to generate finite dc currents. This fact is consistent with the general prescription of the parametric pumping mechanism (9); (40); (41); (77); (78).
Once Eqs. (31) and (32) have been obtained, all stationary currents can be computed by using the relations in Eq.(9) and Eq.(12).

## Iv Results and discussion

In this section we will discuss our main results. We will mainly focus on the generation of finite heat current induced by the pumping mechanism in absence of external bias.
We investigate the net amount of heat current flowing through the 2dTI, studying the behavior of as a function of various external parameters, inspecting also the role played by e-e interactions.

In the following all energies are rescaled with respect to the chemical potential (dimensionless quantities are thus indicated with a bar , i.e., , , , ).
For sake of convenience, we also introduce

(33) |

a dimensionless parameter directly connected to the presence of quantum interference effects, as it appears in the expressions for the pumped current Eqs. (31)-(32) and is linearly proportional to the distance between the two QPCs.

We restrict the analysis to a physically reasonable range of parameters, remaining as close as possible to experimentally accessible values. We fix the chemical potential to meV and the temperature to mK ().
The energy cut-off, which represents the largest energy scale involved, is related to the typical energy bulk gap in 2dTI. Recent proposal have predicted very wide bulk gap for 2dTI in topological materials with strong spin-orbit coupling (60); (61); (62). Recent measurement on novel devices have reported bulk gap values up to meV (60). Here we fix the energy cut-off to .

Concerning the characteristic parameters of the external gate, here we fix which corresponds to a drive frequency of GHz. We vary the gate amplitude in the range . In this case, we can investigate a large range of parameters, including the non-linear regime where interesting features appear.

Let us start considering the non-interacting case (). Interestingly enough we get which, in turn, lead to the vanishing of the heat current component associated to electrons . This fact lead to the conclusion that, in the non-interacting case, a finite heat current is produced and is associated only to a single spin species

(34) |

Therefore this setup can allow for the generation of fully polarized heat current.
Note that the vanishing of the component associated to electrons with is intimately connected to the spin-momentum locking property of 2dTI.
In the end this is also connected to the chosen geometry. Indeed if the external gate is placed on top of the other QPC, this would imply with an opposite polarization of the pumped heat current.
Indeed, physically, the vanishing of spin component can be understood in terms of tunneling paths of electrons. When electrons reach the right QPC, they have not increased or reduced their energy yet by passing through the gate. Therefore, each contribution to transport that could arise from tunneling at is exactly compensated by contribution tunneling at the right QPC.
Only the tunneling path through the other QPC is due to electrons and interference effects would completely vanish.

In Fig. 2 we show the average pumped heat current in the non-interacting case, where an oscillating behavior is clearly visible as a function of the gate voltage .
The latter is due to the Bessel functions and present in Eq. (32).
In order to understand this oscillating behavior, it is instructive to look at the limit of small of (as it is in our case with the chosen parameters)

(35) |

proportional to .
One can see that for , as it is in our case, maxima and minima of , as a functions of the gate amplitude, are several in the considered range of . This basically would explain the oscillating pattern present in Fig. 2 for the polarized heat current in the non-interacting case. Indeed the other present a slower oscillating behavior, since is linked to always smaller than . We underline that the same qualitative features are present also for higher values of , where in general the full expression of governs the observed pattern. Finally, the magnitude of is related to the weigth function present in Eq. (32) and to the product of Bessel functions itself. The latter are responsible of the increasing magnitude for increasing gate voltage in the chosen range of parameters, where Bessel functions of higher order become relevant.

The presence of e-e interactions strongly modifies the above picture, as we now discuss.

The interacting case is reported in Fig. 3. The first important difference with respect to the non-interacting case is that now both spin components contribute to the average pumped heat current , and thus the heat current is not polarized anymore, see in particular Fig. 3b).

The behaviours of and can be ascribed qualitatively to the difference between the Bessel functions and . These functions are plotted in Fig. 3b) versus the gate voltage .
As one can see displays more oscillations than .
Indeed, looking at Eq. (22) one can note that , with , which is for all interaction strengths in the interval .
Therefore, the number of zeros of are reduced of a factor with respect to as a function of . We underline that this factor , related to interaction parameter , arises as a consequence of the particular form for the phases of the bosonized expression for and in Eqs. (53) and is a peculiar property of HLL.

The difference between and are also reflected in the different amplitudes of spin components of the heat current. Summing over higher -th order of the Bessel functions and considering the faster oscillation related to , the amplitude of decreases with respect to , see Fig. 3a).
This trend is reflected in the total heat current, i.e. the sum of the two spin contributions, which is thus dominated by the slower oscillation associated to , as one can see in Fig. 3c). This fact is clearer if one consider stronger interactions (see the green dashed line for ), where the period of oscillations is dominated by the slow spin component, presenting some beats due to the modulation with the fast oscillation associated to . Due to all these different features, the pumped heat current can be also used as a sensitive probe of the presence of e-e interactions. Indeed, as we have shown the heat current is no more fully polarized and presents a different pattern of oscillations, with characteristic modulations related to the interaction strength.

A similar behavior is expected also for the average charge current. Indeed the same pumping mechanism would lead to a finite charge current signal. This is reported in Fig. 4, where is plotted as a function of the gate voltage for different interaction strength. It turns out that in the non-interacting case also pumped charge current is fully polarized and carried solely by spin electrons. The presence of interactions induces the same qualitative changes discussed for the heat current and the associated oscillations can be explained in full analogy to what happens for .

It is interesting, however, to look at the sign of the product of the two pumped currents, i.e.

(36) |

as a function of the gate amplitude as shown in Fig. 5. Here we show the non-interacting case () in the upper panel and in the lower panel. In general, the fact that this quantity has not a definite sign indicates that by varying the gate voltage one can achieve regimes in which charge and heat currents flows in opposite directions. The difference in the sign of charge and heat current, as a function of the gate voltage, can be related to the modulating functions and , which weight the contribution of each photo-assisted tunneling occuring at energy . In the presence of interactions, when , the lower number of oscillations of with respect to allows to better distinguish the effect of processes of each order on the sign of charge and heat currents. On the other hand, for , only remains and brings into each -th process many oscillations, thus resulting in a strong reduction of the modulating function. In the interacting case (see the lower panel in Fig. 5) the regions with a definite sign of as a function of become wider with respect to the non-interacting one. Here charge and heat flow in the same - or opposite - directions for larger range of the external gate voltage .

## V Conclusions

In this work we have investigated a double quantum point contact geometry in a 2dTI. Based on this, we have proposed a setup in which a pumping mechanism can generate finite heat and charge currents without external dc bias or thermal gradient. The pumping mechanism relies on the presence of an external gate potential placed on top of one of the two constrictions which acts on the system with an external ac field. We have developed a microscopic model of the coupling between the external gate and the electrons of the edge states of the 2dTI. We have therefore evaluated the average heat and charge currents pumped into the system as a function of various parameters, taking also into account the presence of e-e interactions. In the non-interacting case, this setup generates finite heat current which is carried solely by a single spin species, resulting into a fully polarized heat flow. The presence of e-e interactions strongly modifies the behavior of pumped currents, which are no more fully polarized since both spin components now contribute to the net heat flow. Moreover, looking at the behavior of the heat current as a function of the gate voltage, one can distinguish between different pattern of oscillations. The latter are connected to the interaction strength, thus allowing to use this setup as a tool to identify the presence of e-e interactions. Finally we have studied the sign of the product of heat and charge pumped currents, showing that one can reach situation in which the two quantities flow in the same or in the opposite directions, by varying the external gate voltage.

###### Acknowledgements.

F.R., M.C., and M.S. acknowledge the support of the MIUR-FIRB2013 – Project Coca (Grant No. RBFR1379UX) and the COST Action MP1209. M.C. wish to thank A. Braggio, F. Taddei, and L. Arrachea for fruitful discussion. People from the Marseille group would like to acknowledge the support of Grant No. ANR-2010-BLANC-0412 (“1 shot”) and of ANR-2014-BLANC “one shot reloaded”. This work was carried out in the framework of Labex ARCHIMEDE Grant No. ANR-11-LABX-0033 and of A*MIDEX project Grant No. ANR-11-IDEX-0001-02, funded by the “investissements d’avenir” French Government program managed by the French National Research Agency (ANR).## Appendix A Zero order currents

This Appendix is devoted to the evaluation of the average currents in the absence of tunneling, taking into account the presence of the time-dependent external gate . In particular we are interested in average quantities in the dc limit or stationary regime. Here, we demonstrate that, in absence of tunneling events (i.e. without QPCs), all average currents vanish.

Since the gate potential couples only to electrons belonging to the top edge, we do not consider here operators related to the bottom edge.
Charge and energy densities associated to the top edge (denoted with the index ) can be written as

(37a) | ||||

(37b) |

The corresponding currents can be defined by using the generalized continuity equation ()

(38) |

Density operators, expressed in terms of bosonic fields, are given by

(39) | ||||

(40) |

We now use the results found in Appendix B, where the time-evolution of operators in presence of the external gate has been calculated. Given Eqs. (49) and (50), we can calculate average currents in absence of tunneling events, in the dc limit. After performing a thermal average with respect to , charge and heat currents are evaluated in the long-time limit (dc regime), using (38)

(41) |

for charge contribution,

(42) |

for heat current.

The expression in Eq. (41) is a periodic function of time and it is linear in the time-dependent gate potential. Since the average over one period of is zero, it is easy to realize that charge current has no dc contribution in the absence of tunneling.

On the other hand, heat current is quadratic in the gate potential, with a non-trivial spatial dependence. In order to calculate its contribution, we thus consider two different regions with and .
In the first region (), the heat current, averaged over one period, reads

(43) |

while in the other region () one has

(44) |

It is worth to note that corresponds to the dc heat current flowing into the right reservoir, while represents the dc heat current flowing into the left reservoir. Since these two contributions are equal, the net heat current globally flowing along the 2dTI is zero. This means that the heat current introduced by the pumping mechanism, without tunneling events and backscattering, is partitioned on two equal parts one to the left and one to the right, resulting into a zero net contribution to the pumped heat flow.

## Appendix B Presence of an external gate

This Appendix is devoted to the derivation of the time-evolution of operators in presence of a time-dependent external gate . As stated in the main text, we focus on a periodic time-dependent gate potential, whose explicit time-dependence is

(45) |

We start by evaluating the time-evolution of the operators in presence of the external gate. The Hamiltonian associated to the gate potential in Eq. (7) can be rewritten in terms of chiral bosonic fields as

(46) |

The equation of motion for the bosonic fields associated to the top edge (in the presence of ) is

(47) |

with the corresponding solution:

(48) |

Here, we have denoted with the solution in the absence of the external gate. Recalling the explicit form of the gate potential in Eq. (8), we have

(49) | ||||

(50) |

where . Recalling the explicit time-dependence in Eq. (45), the expression in Eqs. (49)-(50) become

(51) | |||

(52) |

### b.1 Time-evolution of fermionic operators

We now evaluate the time-evolution of fermionic operators