# Frequency thermal response and cooling performance in a microscopic system with a time-dependent perturbation.

###### Abstract

Following the nonequilibrium Green’s function formalism we study the thermal transport in a composite chain subject to a time-dependent perturbation. The system is formed by two finite linear asymmetric harmonic chains subject to an on-site potential connected together by a time-modulated coupling. The ends of the chains are coupled to two phononic reservoirs at different temperatures. We present the relevant equations used to calculate the heat current along each segment. We find that the system presents different transport regimes according the driving frequency and temperature gradients. One of the regimes corresponds to a heat pump against thermal gradient, thus a characterization of the cooling performance of the device is presented.

Natalia Beraha, Alejandro Soba ,
M. Florencia Carusela

Instituto de Ciencias, Universidad Nacional de Gral. Sarmiento, Los Polvorines, Argentina

Consejo Nacional de Investigaciones Científicas y Técnicas, Buenos Aires, Argentina

Email: nberaha@ungs.edu.ar

## 1 Introduction

Nowadays, the technological implementations in the meso and nanoscale require the management of substantial energies that can be generated, becoming harmful for a device. So significant improvements are needed in the direction of controlling energy, and thus heat, to avoid structural damages. However, many studies show that the management of heat can also display intriguing features which allow the design of devices with novel operating regimes [1, 2].

Spontaneously heat flows from objects at high temperatures to objects at lower temperatures. However heat pumps enables heat to flow against temperature gradient by means of an applied external work. At the meso and nanoscale several interesting applications have been developed in molecular electronics, thermometry and thermal machinary [3, 4, 5]. Many models of heat pumps have been proposed based on different mechanisms such as heat ratchets that periodically adjusts two baths’ temperatures while the average remains equal, brownian heat motors to shuttle heat across the system [6], heat pumps which directs heat against thermal bias in nanomechanical systems [7]. At molecular levels phonon pumps can also be induced by an external force or by mechanical switch on-off of the coupling between different parts of the system [8, 9, 10, 11, 12, 13, 14]. Experimentally the last can be done in molecular junctions or in molecular systems, for example, varying the distance among them or applying stretchings and compressions.

Cooling is another relevant feature because of applications in the quantum realm. Experimental implementations of quantum refrigerators have been developed based on electronic devices in the presence of ac driving fields [15, 16, 17, 18], absortion of phonons (heat) with electrons using nanomechanical devices or active feedback for cooling nanomechanical cantilevered beams [19, 20, 21, 22]. In turn models based on pumping phenomena have been proposed as moving barriers in a cavity to pump phonons from a cold reservoir to a hotter one, or driven two level systems or molecular junctions in contact with phononic baths [23, 24]. In the present work, the mechanism underlying cooling is mainly the phonon manipulation that can work both in insulators and electrical conductors, unlike the electronic cooling where charge, spin and coherence are important.

We present a microscopic model of a phonon pump based on a time dependent modulation of the contact between two one dimensional components. We show that this system can work not only as a heat engine, but it can also operate as a phononic refrigerator.

The paper is organized as follows. In section we present the microscopic model. In section we describe the theoretical framework and methodology used to solve numerically the problem. In section we present the different regimes for heat transport, discussing the role of external frequency and the size of system. In section we analyze the cooling regime comparing with the optimal case. The last section is devoted to summarize the results.

## 2 The model

We propose a microscopic model for energy transport assisted by acoustic phonons. We consider a one dimensional array of atoms, harmonically and bidirectionally coupled, sketched in Fig.1. The chain is divided in two segments (left) and (right) with different coupling intensities and between identical atoms or molecules (referred as ”masses”) and coupled together also harmonically with a coupling constant . The system is subject to a local harmonic pinning potential and connected to (left end) and (right end) semi-infinite chains of masses and respectively, coupled harmonically by spring constants and . This chains play the role of phonon reservoirs with temperatures and . We assume that the system can only vibrate longitudinally that is, we are modelling a heat pump assisted only by longitudinal acoustic modes.

The Hamiltonian of the system can be written as: where are the Hamiltonians of the left (L)/right (R) reservoirs respectively, describes the two central chains (I and II), represents the interaction between them and represents the contact between the two central chains and the reservoirs.

(1) |

with the total number of particles.

The two segments of the central chain have equal length . are the elastics constants in each segment and is the spring constant of the local pinning potential. They are harmonically coupled with a time dependent strength . Thus the interaction Hamiltonian can be written as:

(2) |

with the mass of the atom in the chain , denotes the displacement from the equilibrium position , where is the equilibrium distance between particles and is the momentum. The index labels the atom of a segement that is in contact with a reservoir

The reservoir Hamiltonians corresponding to the semi-infinity chains are

(3) |

with () number of particles.

It is convenient to express this Hamiltonian in terms of normal modes for open boundary conditions

(4) |

The corresponding transformations are:

(5) |

and

(6) |

with

(7) |

The contact is described by the Hamiltonian

(8) |

where elements and are the sites of the central chain in contact with the L/R bath.

In order to study the energy current that flows between the system and the reservoirs, we focus on the energy balance. Conservation of energy implies that the total power invested by all the external fields or mechanical agents must be dissipated into the reservoirs at a rate given by:

(9) |

where and are the rates of heat absorbed by reservoir and respectively (defined positive when heat flows into the reservoir). In the case of periodic fields or drivings ”” means time average over one period.

The energy current that flows to and out a given elementary volume should cancel each other when the system reaches the steady state and in the case there is not local power injection, since there are not temporal variations of the mean local density of energy. As our Hamiltonians involve only nearest-neighbors terms, the minimum volume that we consider is one that encloses nearest-neighbors sites. From the continuity equation and energy conservation, the local time-dependent heat current for the incoming energy current from connecting site of the central chain towards each reservoir gives:

(10) |

with .

We define the net heat current in each segment averaging over an integer number of periods after a transient time as:

(11) |

In the steady state the value of the current is independent of the site in each segment of the central part, therefore and .

## 3 Theoretical framework: Non-equilibrium Green’s function method and Dyson equation.

Following the procedure define in [25], we define the greater, lesser, advanced and retarded Green’s functions

(12) |

(13) |

(14) | ||||

(15) |

(16) | ||||

(17) |

The Dyson equation for the retarded function corresponding to the central chains (without the contacts) is:

(18) |

In turn, when writing the Dyson equation for the retarded Green’s function along the contacts, it is possible to integrate the degrees of freedom of the reservoirs, then:

(19) |

The force matrix includes the interparticle interactions, the pinning interaction and the system-reservoir interaction and is defined as : , where the matrix elements are:

(20) |

(21) |

(22) |

On the other hand, in Eq.(19) corresponds to the the self energy:

(23) |

with a positive infinitesimal and the spectral density of the reservoir :

(24) |

where .

When Langreth’s rules are applied recursively, we can obtain the Dyson’s equation for the lesser Green’s function corresponding to coordinates along the central chain:

(25) |

and to one coordinate on the chain and the other on the reservoir.

(26) |

The corresponding lesser Green’s function of the uncoupled reservoirs is:

(27) |

with a Fourier transform:

(28) |

with being the Bose-Eistein distribution, dependent on temperature of reservoir .

The main goal is to obtain the exact function , solving Eq.(19).

Following the strategy described in [24], we first perform the Fourier transform with respect the ”delayed time” (t’) in .

(29) |

Then, considering the Fourier expansion of the time dependent Hamiltonian, we substitute in the Dyson’s equation (Eq. (19)):

(30) |

where

(31) |

is the stationary component of the retarded Green’s function of the central chains connected to reservoirs, when there is no time-dependent perturbation. Due to the temporal periodicity of we can expand Eq.(30) in terms of the Floquet components:

(32) |

Generally, to obtain an exact solution of Eq.(32) it must be solved numerically. However under some approximations or conditions it can be solved more easily, as it is our case. As we are considering amplitudes of smaller than the energies of the time independent part of and a time-dependent perturbation that contains only one harmonic , we can evaluate the Green’s function up to first under as:

(33) | ||||

(34) |

Now we can we calculate and express the dc component of the heat current in terms of the obtained Green’s functions. Starting from Eq.(10), the time dependent heat current is:

(35) |

Using the function obtained in Eq.(26) and the expansion given in Eq.(32) we obtain the dc heat current flowing in or out of the reservoir :

(36) |

It is interesting to observe that the heat current may not be zero even in the absence of a temperature gradient, just mediated by the energy contribution due to the power injected into the system.

## 4 Results

In the linear chain, the energy flowing in/out of a lead can be calculated from Eq.(36). In a first step, we work with a modulated coupling with a strength smaller than the other coupling constants. Under this approximation the use of perturbative calculations would be sufficiently. However as we explore broad parameter ranges where these calculations are not expected to produce accurate results, we calculate currents exactly by numerically solving the Dyson equations.

We use the following dimensionless parameters: spring constants in units of , moments in units , frequencies in units and temperatures in . For a typical atom and a typical situation these units corresponds to frequencies and temperatures . Thus the nondimensional temperatures to correspond to temperatures of the order to .

We explore several situations of operation for our device. First, we sketch the case when strength of interaction between segments takes a constant value. Transport is induced by the constant thermal bias between reservoirs, that in the absence of time-dependent perturbations the dc heat current is given by the well-known Landauer-Büttiker formula for phononic systems [24, 26]:

(37) |

with

(38) |

the transmission function between reservoirs and .

In order to show how the presence of the on-site potential and the coupling constant affects the transmission we plot in Fig.2, for a weak and strong coupling with and without pinning potential.

As our model only consider harmonic interactions, the on-site pinning potential plays a crucial role in establishing a steady-state current that can be estimated from Eq.(36). The dc heat current corresponds to the quasi-steady state current established after a transient that is established earlier as long the on-site pinning potential is stronger [27].

On the other hand the heat transmission is mainly mediated by low frequency modes, the main responsible for heat conduction. The on-site potential produce a shift in the minimum cutoff frequency of the phonon bands from zero to , so the transmission function takes appreciable values if [28]. When the strength of the pinning potential is increased, the device has a behavior closer to a thermal insulator. For higher values of the potential, large current oscillations can be suppressed, even for high temperatures. In Fig.2 it is exemplified how the transmission function depends on the pinning potential. The transmission for is independent of the value of the local coupling between segments, as it is expected for non localized modes. The system becomes a good thermal conductor. The effect of is to increase the contribution of the localized modes to the transmission along the interface. That is the contact acts as a scatter point for high frequency phonons but not for low frequency ones. These soft phonon modes (including zero-mode) play a fundamental role in the heat conduction.

In a similar way we can define from Eq.(36) the component of the transmission function for the perturbed system

(39) |

In Fig.3 we compare the normalized for two different frequencies. Changing the frequency of modulation the phonon spectra of the channel can be narrowed or widen, thus more o less higher frequency phonons (localized) can be activated affecting the heat conduction.

For time dependent coupling but in the limit of adiabatic driving or for very fast modulation so the system reduces to a static case. When is very low (adiabatic driving limit), the system reduces to two coupled segments with a coupling constant . In the very fast oscillating limit () the coupling convers to a time average constant value . Our interested is focus in the regime in between these two limits.

The transmission function shows that the system presents a multiresonant energy transport [8]. Therefore with a suitable modulation, the allowed phonon band can be changed or shifted and so the phonons involved in the transport. Consequently heat transfer along the device can be enhanced or restrict controlling and tuning dynamically the phononic thermal channels. And in turn the contact acts as a scatter point or interface where phonons can gain or lose a phonon of frequency .

From now on and without loss of generality, we take the reference temperature, being the variable one. We define and calculate and from Eq. 36, where the sign for rectified currents is positive when heat flows from the system to the reservoir. In Fig.4 we show versus for a resonant frequency.

We find four heat transport regimes:

Regime : , and This is a regime mainly ruled by the temperature gradient, with heat flowing from the hot to the cold reservoir. The difference correspond to the mean power injected according to Eq. 9. Although the currents have different absolute values, the power injected is low enough so the direction of the heat current is determined by . The system works as a heat conductor.

Regime : , This is a regime mainly governed by the injected power that is dissipated into both reservoirs. The onset of this regime is given by a current reversal in one of the segments. When or are equal to zero the device can operate as a local insulator, thus heat is inhibited to flow in one part of the system.

Regime : , and The system acts as a refrigerator pumping energy against temperature gradient. Heat flows from the cold reservoir to the hot one.

Regime : , and Heat flows from hot to cold but in the reverse direction than Regime . Transport is again mainly ruled by the temperature gradient.

Depending on the driving frequency the system may display all or some of these regimes. Thus, adjusting temperature gradients and ac-frequencies it is possible to tune changeovers between them.

The transmission function of the system is strongly dependent on the structural parameters such as , who plays an important role on the potential regimes and transitions between them. In Fig.5 we show the contour phase diagram versus and for and . In both cases the weak coupling condition is fulfilled and takes non zero values. Nevertheless we find that the system cannot pump energy against temperature bias in both cases. For the refrigeration regime extends over a considerable region of the parameter space. However for this regime is absent. As it was mentioned above, one possible reason is related to the differences between the transmission functions. In fig.3 we compare for both cases finding that high frequency channels associated to localized modes are not available for .

We also find another interesting feature that is related to the structure. For a reference value , there is a temperature for which cooling regime is achived with a temperature difference between reservoirs is . The last is strongly dependent on the size of the system, as it is shown in Fig.6. When is changed the spectral density also does, then the possible driving frequencies that will enable cooling in resonant condition in the parameter space. From an experimental point of view it can be crucial to decide the way to increase the temperature intervals for which cooling phenomena occurs. One possibility is related to structural features or ”phononic engineering”, varying the size of the system as in our work or adjusting characteristic parameters associated to interactions between atoms (or in other words the adequate selection of materials). Other possibility is related to operational features of a device, for example tuning the strength or frequencies of the time-dependent perturbations or adjusting the differences of temperature along the device. Figs. 4, 5 and 6 reveal the different possibilities.

## 5 Performance of the phonon refrigerator

The microscopic realization of a phonon refrigeration mechanism we are proposing is based on a dynamical cycle given by the periodic modulation of that produce the activation (or deactivation) of phonon modes or thermal channels.

Since the structure of the leads is asymmetric they produce different phonon speed in each segment. Additionally as they are finite, the system has only limited time before the reflected or transmitted phonons from the reservoirs begin to arrive and interfere with other phonons moving through the contact or propagate along it. Therefore for certain driving frequencies and intervals, the system can block the usual heat flowing in the direction of the temperature gradient, enabling the pumping of phonons against it.

Unlike other works based on a time-dependent on-site potential [7], our mechanism is based on an periodical contact interaction. Therefore our model can describe composite systems on substrates that can be realized by linear lattices with static pinning potentials. Examples of such systems are polymeric polar chains where the opening or locking of phonon channels can be adjusted under the action of electric fields and composed short alkane chains o molecular junctions that can be stretched or pulled under some defined protocol[29, 30, 31].

As the underlying mechanism is the shift and changes in the width of the phonon bands, the system size is one of the structural parameters that plays a relevant role affecting the number of available phonon modes, and therefore the cooling ability of the device. We calculate the ratio of cooling to energy consumption or Cooling Coefficient of Performance as:

(40) |

where is the heat extracted from the cold reservoir and is the heat injected on the hot one. The energy consumption or power injected is thus given by .

In Fig.7 we show the performance of the device plotting as a function of the heat pumped from the cooler reservoir. In our system is an operational control parameter for the efficiency when the device is subject to a temperature difference. The curves show the expected behaviour: as long the heat extracted is bigger so the efficiency of the device. However, the efficiency presents an upper and lower bound for the heat that can be extracted. The maximun and minimun heat extracted for an correspond to a going from zero to the maximun value respectively obtained from the phase diagram. In addition, the presents a monotonic behaviour, decreasing faster as long the heat is lower.

The multiresonant feature of the transport is also reflected in the cooling performance of the system as it is shown in Fig.8 where is plotted versus the size of the system. Therefor a similar non monotonic response can be also obtained by tuning adequately other structural features related to material properties (interaction constants) or operational factors.

## 6 Conclusions

We analyzed the heat current in a microscopic system consisting in a one-dimensional finite composed chain of atoms. It is built from two segments interacting with a substrate. They interact periodically by one of their ends, meanwhile the others are connected to phononic reservoirs. Therefore, our model considered heat transfer assisted by ’longitudinal’ vibrational modes. We have calculated the current using a nonequilibrium Green’s function approach, treatment capable of handling time-dependent perturbations and temperature gradients at the same time.

We found that the system displays different heat transfer regimes according to the temperature gradients, structural and operational factors. In one of them the device can act as a refrigerator, pumping energy against temperature gradient assisted by the power injected into the system. The underlying mechanism is based in the periodic modulation of the contact between different structures of the composite that opens or closes local phonon channels, enabling or suppressing the propagation of phonons along the interface. We showed two ways to increase or reduce the performance a system to operate as a heat pump under a temperature difference. One possibility is to design a device from a given material with the adequate size when the external action is fixed. Other possibility is to exert a controlled external action in order to tune the available thermal channels.

Nowadays, the energies involved in the micro and nanolevel applications as molecular electronics or nanomechanical devices can be important even producing damages or reduction of the operation performance of the object. We explore and proposed a physical mechanism for temperature reduction and control of the local heat flow along a device, topic of increasing and actual technological relevance.

## Acknowledgments

M.F.C., N.B., A.S are supported by PIO Conicet. We thank to cluster TUPAC of the Computational Simulation Center - CONICET.

## References

## References

- [1] J. Ren and B. Li, Phys. Rev. E 81, 021111 (2010)
- [2] N Li, J Ren, L Wang, G Zhang, p Hnggi, B Li, Rev. Mod. Phys. 84 , 1045-1066 (2012)
- [3] A. Nitzan, Science 317, 759 (2007)
- [4] L. Wang and B. Li, Phys. Rev. Lett. 99, 177208 (2007).
- [5] F. Giazotto, T. T. Heikkila, A. Luukanen, A. M. Savin, and J. P. Pekola, Rev. Mod. Phys. 78, 217-274 (2006)
- [6] N. Li, F. Zhan, p. Hnggi, and B. Li, Europhys. Lett 80, 011125 (2009).
- [7] L. Arrachea, E. R. Mucciolo, C.Chamon, and R.B. Capaz, Phys. Rev. B 86, 125424 (2012);
- [8] N. Beraha, A. Soba, R. Barreto, M.F. Carusela, Physica A, 433, 9-16 (2015)
- [9] L. S. Cao, R. W. peng and Mu Wang, Appl. phys. Lett. 84, 011908 (2008)
- [10] B.K. Agarwalla, J.S.Wang, and B.Li, phys. Rev. E 84, 041115(2011)
- [11] S. Zhang, J. Ren and B. Li, phys. Rev. E. 84, 031122 (2011); J. Ren and B. Li, phys. Rev. E 81, 021111 (2010)
- [12] E. Cuansing E., J.S. Wang, phys. Rev. B 81, 052302 (2010); E.Cuansing, J.S.Wang, phys. Rev. E 82, 021116 (2010);
- [13] N. Li, p. Hnggi, and B. Li, Europhys. Lett 84, 40009 (2008)
- [14] B.Q. Ai, D. He, and B. Hu, phys. Rev. E 81, 031124 (2010)
- [15] M. Switkes, C. M. Marcus, K. Campman, A. C. Gossard, Science 19, 283, 5409, pp. 1905-1908 (1999)
- [16] Susan K. Watson, R. M. Potok, C. M. Marcus, and V. Umansky, Phys. Rev. Lett. 91, 258301 (2003)
- [17] S.P. Giblin, M. Kataoka, J.D. Fletcher, P. See, T.J.B.M. Janssen, J.P. Griffiths, G.A.C. Jones, I. Farrer and D.A. Ritchie, Nature Communications 3, 930 (2012)
- [18] G. Granger, J. P. Eisenstein, and J. L. Reno, Phys. Rev. Lett. 102, 086803 (2009)
- [19] M. Montinaro, A. Mehlin, H. S. Solanki, P. Peddibhotla, S. Mack, D. D. Awschalom and M. Poggio, Appl. Phys. Lett. 101, 133104 (2012)
- [20] M. Poot, H. S.J. van der Zant, Physics Reports 511, Issue 5, Pages 273â335 (2012)
- [21] A. Naik, O. Buu, M. D. LaHaye, A. D. Armour, A. A. Clerk, M. P. Blencowe and K. C. Schwab, Nature 443, 193-196 (2006)
- [22] J. T. Muhonen, M. Meschke and J.P. Pekola, Reports on Progress in Physics 75, Number 4 (2012)
- [23] O. Entin-Wohlman, Y. Imry, and A. Aharony, Phys. Rev. B 91, 054302 (2015)
- [24] L.Arrachea and B. Rizzo, Journal of Physics: Conference Series , 012012 (2013)
- [25] L. Arrachea, Phys. Rev. B. 72, 125349 (2005)
- [26] J.S. Wang, J. Wang and N. Zeng, Phys. Rev. B 74, 033408 (2006)
- [27] E. C. Cuansing, H. Li and J.-S.Wang, Phys. Rev. E 86, 031132 (2012)
- [28] B. Li, J.H. Lan, L. Wang, Phys. Rev. Lett. 95 (2005) 104302.
- [29] D.X. Qi, R.W. Peng, L.S. Cao, Q. Hu, R. L. Zhang, X.R. Huang and M. Wang, Phys. Rev. B 85, 214123 (2012)
- [30] H. Wang and Y. Leng, J. Phys. Chem. C, 119 (27), 15216â15223 (2015)
- [31] Q. Li, I. Duchemin, S. Xiong, G.C. Solomon and D. Donadio, J. Phys. Chem. C, 119 (43), 24636â24642 (2015)