## Abstract

The cooling of boiling water all the way down to freezing, by thermally connecting it to a thermal bath held at ambient temperature without external intervention, would be quite unexpected. We describe the equivalent of a “thermal inductor,” composed of a Peltier element and an electric inductance, which can drive the temperature difference between two bodies to change sign by imposing inertia on the heat flowing between them, and enable continuing heat transfer from the chilling body to its warmer counterpart without the need of an external driving force. We demonstrate its operation in an experiment and show that the process can pass through a series of quasi-equilibrium states while fully complying with the second law of thermodynamics. This thermal inductor extends the analogy between electrical and thermal circuits and could serve, with further progress in thermoelectric materials, to cool hot materials well below ambient temperature without external energy supplies or moving parts.

## INTRODUCTION

According to the rules of classical thermodynamics, the flow of heat between two thermally connected objects with different temperatures is determined by Fourier’s law of heat conduction, stating that the rate of heat flow between these objects increases with growing temperature difference, and, more generally, by the second law of thermodynamics, which requires that heat can flow by itself only from a warmer to a colder body. The majority of measures for energy-efficient thermal management in everyday infrastructure are based on these laws, be it for the thermal insulation of buildings and heat accumulators, or for harvesting a maximum of mechanical work from a heat engine operating between two different temperatures. While the validity of these fundamental laws is undisputed, there have been exciting technical developments during the past few years, which appear to be at odds with popular interpretations of these laws. In the simplest version of Fourier’s law, for example, the rate of heat (*T*_{b} that is connected to an object with a different temperature *T*_{r} is given by *k* of the connection is expected to be independent of the sign of the temperature difference *T*_{b} − *T*_{r}. Nevertheless, a number of recent experiments demonstrated that a thermal rectification is possible to some extent, thereby opening up the way for a customized thermal management beyond the simple form of Fourier’s law (*1*). The operation of such a “thermal diode” (i.e., an analog to a diode rectifying electric current) has been demonstrated, for example, by the use of phononic devices (*2*–*6*) and phase change materials (*7*–*10*) as well as the application of quantum dots (*11*) and engineered metallic hybrid devices at low temperatures (*12*), paving the way for even more sophisticated thermal circuits such as thermal transistors and logic gates (*6*, *13*, *14*).

The second law of thermodynamics, on the other hand, imposes strict limits on the efficiency of heat engines, cooling devices, and heat pumps (*15*) and suggests a preferred direction of the flow of heat to reach a thermodynamic equilibrium (*16*). According to the latter interpretation, a hot body with temperature *T*_{b} that is connected to a colder object at temperature *T*_{r} approaches thermodynamic equilibrium with strictly *T*_{b} > *T*_{r}, and *T*_{b} is expected to monotonically fall as a function of time *t* because heat is not supposed to flow by itself from a cold to a warmer body (Fig. 1, A and C). Any undershooting or oscillatory behavior of *T*_{b} with respect to *T*_{r} (Fig. 1, B and D), with a reversing direction of the heat flow and transferring heat from cold to hot, would normally be ascribed to an active intervention to remove heat with external work to be done, or to a violation of the second law of thermodynamics.

Here, we describe a simple thermal connection containing a novel thermal element, namely, the equivalent of a “thermal inductor” that we had originally designed for precise heat capacity measurements in an actively driven thermal circuit (*17*). We show that it can also act in a passive way without any external or internally hidden source of power, and is able to drive the temperature difference between two massive bodies to change sign by imposing a certain thermal inertia on the flow of heat. We demonstrate in an experiment that such a process can occur through a series of quasi-equilibrium states, and show that it still fully complies with the second law of thermodynamics in the sense that the entropy of the whole system monotonically increases over time, albeit heat is temporarily flowing from cold to hot.

## RESULTS AND DISCUSSION

### Modeling of the oscillating thermal circuit

The considered thermal connection consists of an ideal electrical inductor with inductance *L* and a Peltier element with Peltier coefficient Π = α*T*, forming a closed electrical circuit (Fig. 2A) (*17*). Here, α stands for the combined Seebeck coefficient of the used thermoelectric materials and *T* is the absolute temperature of the considered junction between these materials. When an electric current *I* is flowing through, heat *Q* is generated or absorbed at a rate *C* and a thermal reservoir (or two finite bodies) are each in thermal contact with the opposite sides of the Peltier element, providing a thermal link by its thermal conductance *k*. The process is described by Kirchhoff’s voltage law in Eq. 1A containing the generated thermoelectric voltage α(*T*_{b} − *T*_{r}) and by the thermal balance equations (Eqs. 1B and 1C) for the rates of heat removed from (or added to) one body *R* is the internal resistance of the Peltier element and α is taken as a constant for simplicity. We also temporarily ignore parasitic effects due to other sources of electrical resistance or thermal transport through leads or convection, and the heat capacity contribution of the Peltier element is thought to be absorbed in *C*. We count *I* and α in Eq. 1 turns out to be unimportant. The dissipated joule heating power *RI*^{2} is regarded to be equally distributed to both sides of the Peltier element. The individual contributions to the flow of heat in Eqs. 1B and 1C are visualized in Fig. 2B. The set of equations (Eq. 1), but without the inductive term *18*, *19*).

To consider the situation for the actual experiment to be presented further below, where a finite heat capacity *C* is connected to an infinite thermal reservoir as shown in Fig. 1 (A and B), we combine *T*_{r} = const. with Eqs. 1A and 1B and obtain a nonlinear differential equation for *I*(*t*), namely

This equation can be easily numerically solved with high accuracy. The time evolution *T*_{b}(*t*) could then be obtained by inserting the corresponding solution *I*(*t*) into Eq. 1A. For a systematic analysis, we restrict ourselves to Δ_{0} = *T*_{b}(0) − *T*_{r} < < *T*_{r} with temperature-independent α, *k*, *R*, and *C*, valid within a sufficiently narrow temperature interval ± Δ_{0} around *T*_{r}. Then, the last two terms of the differential equation (Eq. 2) are negligible because with Eq. 1A,

### Discussion of the oscillatory solutions

It is very instructive to discuss at first the analytical solutions of this simplified equation. The initial conditions for *t* = 0, i.e., the time when the virgin thermal connection is inserted, are *I*(0) = 0 and *T*_{b}(0) − *T*_{r} = Δ_{0}, thereby fixing *T*_{b}(*t*) below *T*_{r}, we focus at first on oscillatory solutions of *I*(*t*). The condition for their occurrence, 4α^{2}*T*_{r} > *LC*(*k*/*C* − *R*/*L*)^{2}, can always be fulfilled for any value of α, if *L* is chosen as *L** = *RC*/*k*. While *R* and *k* are given by the characteristics of the Peltier element, *C* is limited only by the heat capacity of the Peltier element but can otherwise be chosen at will. The solution of interest is *I*(*t*) = *I*_{0} exp(− *t*/τ)sin(ω*t*), with τ and ω matched to fulfill the differential equation, and *I*_{0} = αΔ_{0}/(*L*ω). The corresponding phase-shifted solution for the temperature is *T*_{b}(*t*) − *T*_{r} = Δ_{0} exp(− *t*/τ)cos(ω*t* − δ)/cosδ with tanδ = (*R* − *L*/τ)/*L*ω. We now seek the particular solution realizing the weakest possible damping of *I*(*t*) within an oscillation cycle. This occurs for a maximum value of ωτ, where *L* = *L**, *C*/*k*. Introducing the standard definition of the dimensionless figure of merit for a Peltier element at *T* = *T*_{r} (*20*), *ZT*_{r} = α^{2}*T*_{r}/*kR* = (ω*τ*)^{2} (hereafter abbreviated as *Z*_{T}), the first minimum of *T*_{b}(*t*) is attained with these parameters at *t*_{min} = βτ* with

If expressed by the dimensionless quantities (*T*_{b}(*t*) − *T*_{r})/Δ_{0} and *t*/τ*, the time evolution and (*T*_{b,min} − *T*_{r})/Δ_{0}, a measure for the maximum obtained cooling effect, only depend on *Z*_{T} but are independent of the thermal load *C* and the other parameters of the system. In Fig. 3A, we show a series of corresponding *T*_{b}(*t*) curves that we numerically obtained with Eq. 3 for different values of *Z*_{T}, expressed in these dimensionless units, with *T*_{b}(0) − *T*_{r} = Δ_{0} = 80 K and *T*_{r} = 293 K to mimic a realistic scenario. However, despite the fairly large ratio Δ_{0}/*T*_{r} ≈ 0.27, the difference between the thus obtained values and the results calculated using the explicit solution of Eq. 3 with

According to Eq. 4, the minimum temperature decreases with increasing *Z*_{T} but is limited to *T*_{b,min} > *T*_{r} − Δ_{0} for a finite value of *Z*_{T} so that no catastrophic oscillation can occur. If the thermal connection were removed after the body has reached its minimum temperature, then *T*_{b} would stay at *T*_{b,min} < *T*_{r} under perfectly isolated conditions, as sketched in Fig. 1B. Removing it in a state where *I* = 0 even leaves the thermal connection in its original virgin state at *T*_{b} < *T*_{r} (inset of Fig. 3A). Any external work associated with the act of removing the thermal connection could be made infinitesimally small, for example, by opening a nanometer-sized gap between the body and the thermal connection. If the connection is not removed at all, *T*_{b}(*t*) exhibits a damped oscillation around *T*_{r}, approaching thermal equilibrium with eventually *T*_{b} = *T*_{r}. We note that the maximum possible cooling effect is not reached for the above parameters, but for *L*_{opt} = λ *L** with λ(*Z*_{T}) > 1 (see section S1). The corresponding solutions for *I*(*t*) are overdamped for *Z*_{T} < 1/3, but the temperature of the body still undershoots *T*_{r} for all values of *Z*_{T} > 0.

In a closely related scenario, two finite bodies with identical heat capacities 2*C* and different initial temperatures *T*_{b}(0) and *T*_{r}(0) are thermally connected in the same way and observed under completely isolated conditions (Fig. 1D). In the limit *I*(*t*) but with *T*_{r} replaced by *T*_{b}(*t*) and *T*_{r}(*t*) for *Z*_{T} = 5, together with the average temperature

### Oscillatory flow of heat

Each time when *T*_{b} − *T*_{r} changes its sign (this occurs for the first time as soon as *T*_{b} drops below *T*_{r}, for *L* = *L** after *t*_{0} = π/2ω*; Figs. 1 and 3), heat is still continuously flowing from the cold to the warmer object (dark/purple arrows in Figs. 1 and 2B) until |*T*_{b} − *T*_{r}| reaches its maximum, where the direction of the heat flow is reversed. The moving charge carriers drive virtually all of this heat directly away from the cold to the warm end without being temporarily stored as energy of the magnetic field residing in the inductor. The maximum amount of magnetic energy, _{0}/*T*_{r} < < 1 of the excess heat ~ *C*Δ_{0} that has been initially stored in the warmer body. In this sense, the electrical inductor acts, in interplay with the Peltier element, only as the driver of the temperature oscillation by imposing a certain thermal inertia that temporarily counteracts the flow of heat dictated by Fourier’s law. Thus, we can interpret the role of the circuit as that of a thermal inductor. In analogy to the self-inductance *L* of an electrical inductor generating a voltage difference Δ*V* according to *L*_{th} = *L*/(α^{2}*T*_{r}) (*17*), obeying

From a refrigeration engineering point of view, we can divide a full period of an oscillation cycle of *T*_{b}(*t*) into four stages (see Fig. 2B). In the first quarter of a full period (i), the operation corresponds to that of a thermoelectric generator. In the second quarter (ii), the circuit acts as a thermoelectric cooler even for *T*_{b} < *T*_{r}, driven by the electric current that is still flowing in the original direction due to the action of the electric inductor. During the third quarter (iii), the electric current changes sign (thermoelectric generator), while heating persists in the fourth quarter (iv) even for *T*_{b} > *T*_{r}, and the device is then operating as a thermoelectric heater.

The order of magnitude of the rate of heat flowing between the objects can be chosen arbitrarily low by adjusting the thermal load *C*. In this way, the electronic oscillator circuit reaches the typically very large time scales encountered in thermal systems (i.e., seconds, minutes, or even longer). This guarantees that the processes, albeit irreversible, can run in a quasistatic way and pass through a series of quasi-equilibrium states with well-defined thermodynamic potentials and state variables of the bodies. This is in marked contrast to nonequilibrium oscillatory processes such as the Belousov-Zhabotinsky chemical reaction (*21*, *22*), other oscillations in complex systems far from thermodynamic equilibrium (*23*), to thermal inductor type of behaviors associated with the convection of heated fluids (*24*), or to transient switching operations in light-emitting diodes (*25*).

### Experiments

In a real experiment, measurements of sizeable temperature oscillations face certain challenges. At present, the most efficient Peltier elements have a maximum *Z*_{T} ≈ 2 (*26*). In a scenario of cooling an amount of boiling water from 100°C down to its freezing temperature at 0°C by connecting it to a heat sink at room temperature (20°C), for example, a *Z*_{T} ≈ 5 would be required (Fig. 3A), which is out of reach of the present technology. A further challenge is the choice of the inductance *L*. It should be large enough to allow the cooling of a substantial amount of material while keeping *k* as small as possible to maximize *Z*_{T}. With τ* = *C*/*k* of the order of several seconds and typical internal resistivities of Peltier elements *R* ≈ 0.1 ohm and higher, *L*_{opt} > *L** = *R*τ* must be of the order of 1 H or larger, although a useful cooling effect could still be achieved for *L* < *L**. Normal conducting inductors with these large inductances suffer from a considerable internal resistance *R*_{s}, however, thereby reducing the cooling performance well below the theoretical expectations. The incorporation of a corresponding finite resistivity *R*_{s} switched in series, in addition to the internal resistance *R* of the Peltier element as it is sketched in the inset of Fig. 4A, leads again to Eq. 3 but with *R* replaced by *R* + *R*_{s} (see section S3). This will reduce the effective dimensionless figure of merit *Z*_{T} of the Peltier element by a factor *R*/(*R* + *R*_{s}), which can be substantial as soon as *R*_{s} becomes of the order of *R* or even larger, with an associated decrease of the cooling performance, as illustrated in Fig. 3A and in section S1.

Therefore, we performed an experiment with two configurations of resistanceless superconducting coils (with *L* = 30 and 58.5 H, respectively) connected to a commercially available Peltier element. With a cube of ≈1 cm^{3} copper as the thermal load at temperature *T*_{b} and a massive copper base acting as the thermal reservoir held at *T*_{r} = 295 K (≈22°C), we are thereby implementing an entirely passive oscillating thermal circuit that should show the predicted temperature oscillations (see Materials and Methods). In Fig. 4A, we present the results of these measurements according to the experimental scheme shown in Fig. 1B and analyzed in Fig. 3A. Initially heated to 377 K (≈104°C), *T*_{b} dropped for *L* = 58.5 H by ≈1.7 K below the base temperature *T*_{r}, verifying that heat has been flowing from the chilling copper cube to the thermal reservoir as soon as *T*_{b} fell below 295 K around *t*_{0} ≈ 410 s. Using the known values *k*, *R*, and α of the Peltier element that we had previously obtained by the procedure described further below, we modeled these experimental *T*_{b}(*t*) data by taking into account the parasitic resistance *R*_{s} in series due to the long electric lines to the superconducting coils (inset of Fig. 4A), according to the corresponding analogon to Eq. 2 (see eq. S5 in section S3). The data are best reproduced with *R*_{s} = 21 and 43 milliohms for *L* = 30.0 and 58.5 H, respectively (solid lines in Fig. 4A), in very good agreement with our direct measurements for the total resistivity of the leads and connections, *R*_{s} ≈ 20 and 45 milliohms (see Materials and Methods).

In Fig. 5, we have also plotted the different fractions of the rates of energy flow in this experiment, as we have depicted them in Fig. 2B. The corresponding data were derived from the measured temperature difference Δ*T* = *T*_{b}(*t*) − *T*_{r}, the electric current *I*, and the known parameters of the Peltier element. Except for the comparably small electromagnetic power term (*R*_{s}, all the relevant exchange of energy occurs directly between the Peltier element, the thermal load and the thermal bath, with the thermoelectric contributions α*T*_{b}*I* and α*T*_{r}*I* dominating the other terms appearing in Eqs. 1B and 1C.

To be better able to quantitatively analyze our model and to precisely determine the parameters of the Peltier element, we had previously performed a series of additional experiments using an active gyrator-type substitute of a real inductor (*27*) that allows simulation of almost arbitrary values of *L* with negligible effective internal resistivity, *R*_{s} ≈ 0. Although the thermal connection can then no longer be considered as strictly operating “without external intervention,” no net external work is performed on the system. In Fig. 4B, we present the results of a series of corresponding measurements for four different values of the nominal inductance *L* of the gyrator (*L* = 33.2, 53.6, 90.9, and 150 H), with *T*_{b} reaching a temperature of ≈2.7 K below *T*_{r} for *L* = 90.9 H. In the same figure, we also show the results of a global fit to both the measured *T*_{b}(*t*) and *I*(t) data according to the relations given in this work, with only four free parameters: *C*, *R*, *k*, and *Z*_{T}. We obtain an excellent agreement for an effective *Z*_{T} = 0.432, *C* = 4.96 J/K (this value includes a heat capacity contribution of the Peltier element), *R* = 0.22 ohm, and *k* = 0.0318 W/K, resulting in a Seebeck coefficient ^{−3} V/K. The used values *L* = 33.2 and 90.9 H are very close to the resulting *L** = 34.3 H and the optimum inductance *L*_{opt} = 94.4 H, respectively (see section S1).

### Relation to the second law of thermodynamics

The fact that heat temporarily flows from cold to hot in the processes described above inevitably calls for a further discussion in view of the second law of thermodynamics. The proof that these processes do not violate this fundamental law is surprisingly simple. The total rate of entropy production is *k*(*T*_{b} − *T*_{r})/*T*_{b} + *k*(*T*_{b} − *T*_{r})/*T*_{r} ≥ 0, and we have *T*_{b}(*t*) and *T*_{r}(*t*) counter-oscillate in the second scenario (Figs. 1D and 3B) and repeatedly match around *T*_{b} = *T*_{r} must be larger than with *T*_{b} ≠ *T*_{r}. As *T*_{av} is not constant but lowest for *S*_{tot}(*t*) is a monotonously increasing function of time *t* (inset of Fig. 3B). For *t* → ∞ and in the limit *LCR* circuit.

Although no external intervention is generating the inversion of the temperature gradients and forcing heat to flow from cold to hot, the inductor is, as an integrated part of the thermal connection, perpetually changing its state due to the oscillatory electric current. Therefore, the described processes are also in full conformity with the original postulate of the second law of thermodynamics by Clausius, stating that a flow of heat from cold to hot must be associated with “some other change, connected therewith, occurring at the same time” (*16*).

We may finally analyze the considered thermal oscillator in view of the thermodynamic efficiency of heat engines. If we assume an ideal Carnot efficiency (*15*) for stage (i) in the configuration of Fig. 1B where the circuit is acting as a thermoelectric generator, a maximum amount of work d*W* ≤ *C*d*T*_{b}(1 − *T*_{r}/*T*_{b}) can be extracted from the heat capacity *C* by successively removing infinitesimally small amounts of heat *C*d*T*_{b} upon cooling it by d*T*_{b}, summing up to *W* ≤ *C*(Δ_{0} − *T*_{r}ln[*T*_{b}(0)/*T*_{r}]) for the full range of temperatures from *T*_{b}(0) = *T*_{r} + Δ_{0} down to *T*_{r}. By completely recycling this amount of work during stage (ii), where the thermal load is chilled from *T*_{b} = *T*_{r} further down with a corresponding Carnot engine acting as a cooler, we reach the lowest possible temperature *T*_{b,min} that can be achieved in this way, which is given by the implicit equation*Z*_{T} → ∞ and Δ_{0} < < *T*_{r}, where *T*_{r} − *T*_{b,min} → Δ_{0}. It is independent of any model assumptions and therefore represents an ultimate limit for the maximum possible cooling effect during any thermal oscillation cycle without external work done on a system. However, this thermodynamic constraint is crucial neither for the experiments reported here nor for practical applications using the present technology, with *Z*_{T} values still of the order of unity (*26*). In the temperature range used for the present experiments, for example, the lowest achievable temperature allowed by thermodynamics could still be as low as *T*_{b,min} ≈ 226 K (≈−47°C).

## CONCLUSIONS

We have shown both theoretically and experimentally that the use of a thermal connection containing a variant of a thermal inductor can drive the flow of heat from a cold to a hot object without external intervention. We have analyzed a corresponding oscillatory thermal process in detail, from both the electrical and thermodynamic points of view. While it fully complies with the second law of thermodynamics, we identify a general thermodynamic limit for the maximum possible cooling effect that can occur during such a thermal oscillation cycle. Despite the conceptual simplicity of the described experiment and based on the laws of classical physics, it has, to our knowledge, never been considered in the literature. With future progress in materials research, the technique may become technically useful and allow the cooling of hot materials well below the ambient temperature without the need for an external energy supply or any moving parts.

## MATERIALS AND METHODS

We used a commercial Peltier element (module TB-7-1.4-2.5, Kryotherm Inc.) for all the experiments. The temperature difference between the thermal load (a copper cube of dimensions ≈ 10 mm × 10 mm × 10 mm with a mass *m* ≈ 9 g) and the thermal bath (a block of 3.65 kg of copper) was measured with a copper-constantan differential thermocouple, thermally connected to the experiment with silver paint. The corresponding voltage was measured using a Keysight multimeter (model 34465A) and converted to a temperature difference according to a standard type T conversion table.

The superconducting coils used in the experiment were part of two physical property measurement system (PPMS) experimental platforms (Quantum Design Inc.) and were immersed in liquid helium at *T* = 4.2 K. The resistivities *R*_{s} of the electric lines and connections to these coils were measured with a PeakTech 2705 milliohm meter with an accuracy of ±2 milliohms. For one superconducting coil with inductance *L* = 30 H, we measured *R*_{s} = 20 milliohms, and for two coils in series with a total *L* = 58.5 H, we obtained *R*_{s} = 45 milliohms. The inductances were determined with a ≈2% accuracy by measuring the time constant of a decaying current in a closed loop with a known 20-ohm resistor. The resistivities *R*_{hs} ≈ 35 ohms of the heat switches (“persistent switches”) in parallel to the superconducting coils have been included in this analysis. However, for the operation in a thermal oscillator experiment, their presence is quantitatively negligible.

The gyrator was built on the basis of the scheme described in (*27*). The values for *L* can be adjusted by an appropriate choice of a resistor in the gyrator circuit, and they have been crosschecked by measuring the time constants of a respective *LR* test circuit. The current *I* was determined by monitoring the voltage across an internal shunt resistor inside the gyrator. All experiments were done in open air without any thermal shielding. In a global fit to all the available *T*_{b}(*t*) and *I*(*t*) data to Eq. 2, we obtained an excellent agreement for *C* = 4.96(1) J/K (this value includes a heat capacity contribution of the Peltier element), *R* = 0.220(1) ohm, and effective values for *k* = 0.0318(1) W/K and *Z*_{T} = 0.432(1), resulting in *L** = 34.3 H and *L*_{opt} = 94.4 H. We independently confirmed the ratio of the fitted values *C* and *k*, τ = *C*/*k* ≈ 156.1(1) s, in a direct measurement by monitoring the decay of the temperature of the copper cube without the inductor connected and the Peltier circuit open, with τ ≈ 162 s. The effective *Z*_{T} = 0.432 obtained from the global fitting procedure has to be interpreted as a constant *Z*_{T} = α^{2}*T*_{r}/*kR*, where *k* and *R* in Eq. 1 may include extrinsic contributions due to parasitic heat transport and internal resistance from the wiring. The true intrinsic *Z*_{T} value of the Peltier element may therefore be somewhat larger, and we estimated it to *Z*_{T} ≈ 0.5 near room temperature. The fitted value *R* = 0.22 ohm of the circuit including the internal wiring is also compatible with the resistance specifications of the Peltier element (*R* = 0.18 ohm ± 10%), which represents the main source of electrical resistivity. Pictures of the experimental setups are provided in section S5.

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/5/4/eaat9953/DC1

Section S1. Maximum possible cooling effect

Section S2. Two finite bodies with different temperatures

Section S3. Effect of a finite resistance *R*_{s} in series

Section S4. Analogy to a thermal inductor

Section S5. Pictures of the experimental setups

Fig. S1. Optimized values for a maximum temperature undershoot.

Fig. S2. Experimental setup of the oscillating thermal circuit.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is **not** for commercial advantage and provided the original work is properly cited.

## REFERENCES AND NOTES

**Acknowledgments:**We thank A. Vollhardt, D. Florin, and D. Wolf for their technical support.

**Funding:**The authors acknowledge the financial support by the UZH and by the Schweizerische Nationalfonds zur Förderung der Wissenschaftlichen Forschung (grant no. 20-131899).

**Author contributions:**O.B. introduced the concept of using a Peltier element in combination with an inductor for thermal experiments. X.Z. modified and characterized the setup for use with the superconducting coils. A.S. wrote the manuscript and performed all the measurements and the full analysis of the circuit with regard to its use as a passive thermal oscillator and its compliance with thermodynamics.

**Competing interests:**The authors declare that they have no competing interests.

**Data and materials availability:**All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

- Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).