Bio-moleculear thermal oscillator and constant heat current source

The demand for materials and devices that are capable of controlling heat fl ux has attracted many interests due to desire to attain new sources of energy and on-chip cooling. Excellent properties of DNA make it as an interesting nanomaterial in future technologies. In this paper, we aim to investigate the thermal fl ow through two sequence combinations of DNA, e.g, (AT)4 (CG)4 (AT)4 (CG)4 and (CG)8 (AT)8. Two interesting phenomena have been observed respectively. In the fi rst confi guration, an oscillatory thermal fl ux is observed. In this way, an oscillating heat fl ux from a stationary spatial thermal gradient is provided by varying the gate temperature. In the second confi guration, the system behaves as a constant heat current source. The physical mechanism behind each phenomenon is identifi ed. In the fi rst case, it was shown that the transition between thermal positive conductance and negative diff erential conductance implies oscillatory heat current. In the latter, the discordance between the phonon bands of the two coupled sequences results in constant thermal fl ow despite of increasing in temperature gradient.


Introduction
The Recent decade has been faced with the great research interests in heat low management in nano scale devices [1,2]. The ability to thermal control may have interesting applications such as heat waste recovery, energy storage, and phononics. In the phononics (understanding and controlling the thermal low), it is desirable to develop the thermal circuit ingredients, such as thermal diodes [3][4][5], thermal transistors [6,7], thermal memory [8] and thermal logic gate [9], analogous to electrical circuit components. The irst report of heat manipulation was represented in 2002 [3]. It was an elementary presentation of nano thermal diode or recti ier. In 2006, Li, et al. demonstrated the irst model of thermal transistor [6]. It was a three-gate system in which by tuning the temperature of the third gate, one can manage the thermal current in the two gates.
Moreover, constant current sources and current oscillators, which have the signi icant roles in electronic circuits, have been founded for many years. Their thermal peer, namely, heat current oscillators and constant heat current sources are poorly realized. Thermal oscillator provides the oscillatory thermal current from a stationary temperature bias. While, a steady heat current source prepares a constant thermal flow contrary to a variation in the temperature bias. They are two main bases of phononics.
The purposes of the current paper are two folds. First, a thermal oscillator for creating an oscillating heat lux from a stationary spatial temperature gradient between a hot reservoir and a cold reservoir is suggested. Second, a new design of constant thermal current sources is presented.
In the irst part, the system con iguration includes two portion, e.g. "Source" (S) and "Drain" (D). At the right end, there is a hot reservoir and the cold reservoir is at the left end. The source and drains are attached to each other. The source is joined to the hot reservoir from right. Drain is connected to the cold reservoir from left. At the interface, the system is coupled to the gate reservoir. The gate reservoir play the role of control thermostat equilibrated at T G . The con iguration is shown in igure 1. Each portion includes an asymmetric sequence combination of DNA molecule. Superior attributes of DNA propose employing it as an outstanding nanomaterial in future technologies. DNA molecule bene its from high mechanical lexibility, so it will not be imposed any break down during an operation.

Materials and Method
Many efforts as much as mathematical models have been devoted to characterize DNA dynamic [10]. It is stated that Peyrard-Bishop-Daxois model (PBD) can successfully describe DNA dynamics [11][12][13]. In this model, the strong covalent bonds were utilized to describe the nucleotides which are aligned in the same strand. An anharmonic nearest-neighbor interactions demonstrates the covalent bonds. A nonlinear onsite potential is employed to characterize the coupled strands through the weak hydrogen bonds.
The Hamiltonian of the model writes  When a system is connected to the heat baths at different temperatures, one can consider the heat transport characteristics in a Hamiltonian system at a stationary state. Various models have been employed to control the temperature of the system as a thermostat, e.g, the Anderson thermostat [14], Berendsen thermostat [15], Nos´e-Hoover thermostat [16], and Langevin thermostat [17]. In the Langevin approach, the heat bath is modeled by using the stochastic process, so the reservoirs is simulated by adding a dissipative term and a stochastic Gaussian white noise term in the Newtonian equation of motion. Due to stochastic nature of Langevin heat bath, its utilization is accompanied by complication and consuming. To prevail the dif iculties originated from stochastic processes, a new approach is developed by Nose and Hoover with a deterministic nature. It includes only one imaginary variable. This variable simulates system's achieve to desired temperature. It is known that in the nearly linear chains, Nos´e-Hoover thermostat is not effective and will result to arti icial results [18,19]. In such cases, it is essential to use Langevin approach despite of its dif iculties. Fortunately, DNA has a strong nonlinear nature and one can use the Nose-Hoover approach with con ident. Now, we thermalize the r irst particles at T R and l last particles at T L and g intermediate particles at T G by using Nose-Hoover thermostat chains [16]. The dynamic of thermostats is described as: , ,  (1). ODE45 solver is used to integrate the differential equations of motion. To reach the steady state, long enough time is regarded. Local temperature at site n is computed by Virial theorem: So, the total heat lux is de ined as n n J j   . High enough time (200 ns) is regarded so that the system attains to the steady state. Figure 2 shows the steady state low of heat lux through a homogenous chain (AT) after 2 × 10 7 tu.

Thermal oscillator
An electronic relaxation oscillator creates an oscillating electronic current from a stationary spatial voltage. Thermal analogue of electronic oscillator is desirable for technological application like energy harvesting devices, switching devices or clocking devices and sensing devices. A few reports of thermal oscillator have been presented based on different mechanisms (Ref. 20 and references therein). Here, it was shown that the phase transition between the positive differential thermal conductivity and negative differential thermal conductivity is responsible to thermal oscillation mechanism. The system is con igured as follows: two coupled chain are the same, each one with the sequence (AT) 4 (CG) 4 . The system at the left end is attached to the cold thermostat at T C = 3 K. At the middle part, the chain is coupled to the gate thermostat as a control thermostat. The hot thermostat is connected to right end which is equilibrated at T H = 30 K. Figure 3 shows the oscillating heat current through two S and D segments. To consider the physical nature of phenomenon, we calculate the thermal. Figure 2 shows local heat flux j n through the system after 2 × 10 7 tu conductivity through the Fourier's law. The Fourier's law is expressed as where the heat flux j is the amount of transported heat through the unit surface per unit time. T is ×% the temperature gradient.
It was shown that the t r ansition between positive differential thermal conductivity and negative differential thermal conductivity implies oscillatory heat current. It is proved that negative differential thermal conductivity is possible due to nonlinear nature of system. Nonlinearity varies the phonon band's position with temperature. It makes phonon bands have match and mismatch dependent on temperature [6,21]. The result is shown in igure 4.

Constant heat current source
Constant current sources have been employed in electronic circuits for many years. A constant (electric) current source supplies a ixed charge flow, in spite of a variation in the power supply. In Ref. 22, by attaining complete NDTR, it has been presented the irst report of constant thermal current sources. Here we present new design of constant thermal current sources. The chain comprises of 32 base pairs. Sequence combination (CG) 8 (AT) 8 is used. The left end is coupled to the cold thermostat at 2 K. The right end is connected to the hot thermostat equilibrated at 20 K. It was shown that by varying right temperature, the thermal current in the sequence is nearly constant. The result is shown in igure 5. Under presented situation, the device can suffer the changes of the temperatures of heat baths and yields a fully stable value of heat current.
To investigate the mechanism which is resulted in occurrence of constant current source, we consider the power spectrum of AT and CG segments. The temperature of CG sequence (blue) is ixed at 2 K. The temperature of AT sequence is increasing T + = T -(1 + ). For AT sequence, power spectrums have been illustrated as a pink (18 K) and red (10 K) style in the igure 6a,b. For CG sequence which is equilibrated at 2 K, power spectrum is shown as a blue style.   Clearly it was seen from igure 6 that the discordance between the phonon power spectrums of the two attached sequences results in constant thermal low despite of increasing in temperature gradient. Moreover, it is obvious that by increasing the temperature of AT sequence e.g., by increasing the temperature difference, the overlap between phonon bands do not vary.
So, in the region determined at the igure 5, the thermal current through the system is constant. The temperature pro ile is shown in igure 7. It is clear the heat lux could not low from irst half part of the sequence to the left part. So, one can use this property to design the constant thermal source.

Conclusion
Thermal management in nanodevices may have a   fundamental role in the development of future nanoelectronic circuits due to the importance of the heat among nanoelements. Self-assembly capabilities and elastic properties of DNA make it as foremost nanowire materials. Ability of DNA for coding, prediction and its intrinsic diversity persuade one to construct the electronic and heattronic devices based on DNA. We proposed a model of the thermal oscillator and constant heat current source based on DNA nano wire. The underlying physical mechanisms of such behaviors are determined. It was shown that the change in material state between conducting and negative conducting phases is responsible for oscillation of heat current. Also, it was shown that the inconformity between the power spectrum of the two interface particles in some temperature range make constant heat current source and possible.