arxiv: v2 [cond-mat.stat-mech] 19 Aug 2015

Transcription

1 Energy repartition for a harmonic chain with local reservoirs arxiv: v [cond-mat.stat-mech] 9 Aug 05 Gianmaria Falasco Marco Baiesi 3 Leo Molinaro Livia Conti 3 and Fulvio Baldovin 34 Institut für Theoretische Physi Universität Leipzig Postfach D Leipzig Germany Department of Physics and Astronomy University of Padova Via Marzolo 8 I-353 Padova Italy 3 IF Sezione di Padova Via Marzolo 8 I-353 Padova Italy 4 Sezione CISM Università di Padova Via Marzolo 8 I-353 Padova Italy (Dated: August 9 08) We exactly analyze the vibrational properties of a chain of harmonic oscillators in contact with local Langevin heat baths. onequilibrium steady-state fluctuations are found to be described by a set of mode-temperatures independent of the strengths of both the harmonic interaction and the viscous damping. Energy is equally distributed between the conjugate variables of a given mode but differently among different modes in a manner which depends exclusively on the bath temperatures and on the boundary conditions. We outline how bath-temperature profiles can be designed to enhance or reduce fluctuations at specific frequencies in the power spectrum of the chain length. PACS numbers: Ln e a I. ITRODUCTIO The enhancement of nonequilibrium fluctuations at low wavenumbers is a ey feature of systems driven by thermodynamic gradients (see [] for a review). For temperature gradients it has been thoroughly studied both theoretically [] and experimentally in systems ranging from simple fluids [3] to polymer solutions [4] and fluid layers also under the influence of gravity [5]. More recently fluctuations in nonisothermal solids have been the subject of experimental investigation fostered by the possibility of technological applications in fields as diverse as microcantilever-based sensors [6] and gravitational wave detectors [7]. For example the low frequency vibrations of a metal bar whose ends are set at different temperatures were found to be larger than those predicted by the equipartition theorem at the local temperature [8] thus corroborating the generality of the results obtained for nonequilibrium fluids [0]. (See also experiments with cantilevers [9]). Theoretical studies of nonequilibrium solids focused more on thermal conduction in low dimensions where crystals are usually modeled as Fermi-Pasta-Ulam oscillator chains coupled at the boundaries with heat baths at different temperatures [ 3]. Thans to their simplicity integrable and quasi-integrable models may be taen as a paradigm to describe more comprehensively the energetics of normal solids under nonisothermal conditions. For instance anomalous features are nown to disappear when in place of non-homogeneous boundary conditions at the borders a temperature gradient is generated by stochastic heat baths displaced along the system [3]. Specifically it has been shown that self-consistent heat baths such that no energy flows on average into or out of the reservoirs are sufficient to recover the Fourier s falasco@itp.uni-leipzig.de baldovin@pd.infn.it law of heat conduction in a harmonic chain [4 6]. Lifting the self-consistency condition one obtains a simple yet general model which describes a solid immersed in a locally equilibrated medium [7 9]. This can find application in all cases where the study of fluctuations is applied to an extended system with a complex thermal balance. As an example we may cite cryogenic gravitational wave detectors where thermal fluctuations of the systems composed by the test masses and their multistage suspension chains are of central importance. The latter are effectively coupled to different heat baths and flows [0]. Here we analyze the energy repartition among the elastic modes of a harmonic chain held in temperature gradient as setched in Fig. (a). In a coarse-grained picture the oscillator displacements can be thought of as the local strain of a (one-dimensional) elastic dispersive body such that the model describes the damped propagation of thermal phonons (Fig. (b)). Our approach is fully analytic and provides an explicit expression for the energy repartition among the modes in terms of their effective temperatures T. Exemplifying our results for temperature profiles with a defined concavity we show that T s depend only on this concavity and on the boundary conditions of the system. A naive expectation could be that deviations from energy equipartition are to be antic- FIG.. (a) Setch of the linear chain of harmonic oscillators held in a temperature gradient: Each oscillator is coupled to an independent heat bath at temperature T n n = 0... (in the picture =3). (b) Schematic interpretation in terms of sound propagation in a medium.

2 ipated at long wavelengths only since local equilibrium conditions should hold at short scales. On the contrary we find that both long and short wavelength modes can either heat up or cool down well beyond the averagetemperature. We also study a reverse-engineering approach in which the heat bath temperatures are inferred starting from a desired energy repartition. T n n α=0. α= α=0 II. MODEL AD GEERAL RESULTS Consider a linear chain of + equal oscillators locatedatpositionsq n (n = 0...). Successivemasses are connected through a harmonic potential of equilibrium length l 0. Each of them is in contact with a specific Langevin bath at temperature T n [] providing viscous damping with coefficient γ and thermal noise ξ n. Setting masses to unity the equations of motion in the displacement coordinate R n q n nl 0 read R n = γṙn κ A nm R m +ξ n () m=0 where A nm is a tridiagonal matrix accounting for firstneighbors interactions via the potential κ (R m R m ). In Eq. () the standard Gaussian white noise ξ n has an amplitude given by the fluctuation-dissipation theorem at the local temperature (in units of B ): ξ m (t)ξ n (t ) = γt n δ(t t ) δ mn. () In the following we first consider the case of free boundary conditions (A 00 = A = ); fixed (A 00 = A = 0) and mixed (A 00 = A = 0) boundary conditions are discussed later in Appendix A. With free boundaries the matrix A nm is diagonalized by the linear transformation Φ AΦ with Φ n = ( π ( cos n+ ) ) (3) mapping the spatial coordinates R n into the coordinates of the normal modes X n Φ n R n for which Ẍ = γẋ ω X +η (4) ) where ω = 4κ sin( π (+) is the (squared) eigenfrequency of the -th mode. In this dynamics the only source of correlation between modes is contained in the transformed Gaussian white noises η n Φ n ξ n η (t)η (t ) = γt δ(t t )/() (5) These correlation include a temperature matrix T () Φ n Φ n T n (6) FIG.. results. Heat-bath profiles utilized for exemplifying our which is certainly diagonal only in the equilibrium case T n = T n whereenergyequipartitionisrecovered. In a nonequilibrium state generated by heterogeneous bathtemperatures the diagonal T still encodes information about how energy is distributed among the modes. onzero off-diagonal T emerge in connection with energy fluxes. To show this we consider the average inetic energy (K ) and potential energy (V ) of the -th mode K () Ẋ ( δ 0 /) V ()ω X where expectation values are taen over different realizations of the thermal noise ξ n. We get the variances X Ẋ from the solution of Eq. (4) (Appendix B) X (t) = α= t dt ( ) α λ e λα (t t ) η (t ) (7) λ where λ (α) with α = are the roots of the characteristic equation for the unforced harmonic oscillator; namely λ (α) = [γ + ( )α γ 4ω ]. For each mode both the average inetic and potential energy turn out to coincide with one half of the mode temperature (Appendix B): K X = V X = T / ( 0). (8) This relation establishes a form of energy equipartition between the conjugate variables of a single mode. Interestingly from (6) and (3) one sees that T does not depend on the details of both the harmonic interaction (κ) and the damping (γ). Therefore the amount of energy stored in the -th mode is directly determined by the choice of the bath temperature profile T n for given boundary conditions. Put in other words properly designing thermal profiles it is in principle possible to enhance or reduce the thermal vibrations of specific modes. All these findings are confirmed by numerical integration of ().

3 3 III. ROLE OF THE BOUDARY CODITIOS In the case of free boundaries Eq. (6) gives ) T = T + T ncos(n+ + π m=0 T ( 0) (9) m where T T n/() is the average imposed temperature. The center-of-mass inetic energy () X 0 / is equal to T00 = T. otice that Eq. (9) is valid in particular when T n corresponds to a selfconsistent profile [4 6]. In (9) the energy stored by the mode under stationary nonequilibrium conditions emerges lie a correction to the average temperature T which at most amounts to ±T. This correction can be viewed as a weighted average of a cosine function over the temperature profile: For parity it vanishes for all temperature profiles which are odd with respect to (/ T). The relevant physical consequence is that with free boundary conditions energy equipartition is extended to all nonequilibrium temperature profiles which are odd-symmetric with respect to (/ T) lie linear profiles. At variance if the temperature profile has a definite upwards (downwards) concavity in the interval [0] low- (high-) modes heat up and high- (low-) modes freeze down. We exemplify these findings assuming heat-bath temperatures T n = T 0 +(n/) α (T T 0 ) with T 0 = 0 T = and = 99 (see Fig. ): α = corresponds to a linear temperature profile whereas α < (α > ) corresponds to a profile with upwards (downwards) concavity. In Fig. 3(a) one finds the resulting T for free boundary conditions. Transport properties might depend crucially on the boundary conditions []. We show that the latter strongly influences also the repartition of energy among the normal modes. For fixed boundary conditions T becomes (Appendix A) [ ( )T n= T = T ncos ( n π) ] m= T (0) m (0 < < ) with T n= T n/( ). Fig. 3(b) shows the mode energy repartition for the same profiles T n used for open boundary conditions in Fig. 3(a). otably the low- behavior is inverted. For instance while free boundaries enhance the long-wavelength energy storage for concave-up T n fixed boundaries do the opposite. Hence if the aim were to store energy at low s a convenient strategy would be to heat up the boundaries and cool down the middle of a free chain and vice versa with a fixed chain. For mixed boundaries in which we leave free the mass atn = 0andfixthemassatn = wehave(appendixa) ) T = T + T ncos ( (n+)(+) + π m=0 T m () T / T α = 0. α = α = 0 (a) free (b) fixed (c) mixed free end hotter (d) mixed free end colder FIG. 3. Modes normalized temperatures T /T calculated through Eqs. (9) (0) () for the heat-bath profiles mentioned earlier. The four panels refer to different boundary conditions: (a) free (b) fixed (c) mixed with free end hotter than the fixed end and (d) vice versa. Insets enlarge the plots at low s. ( < ) with T T n/. Due tothe broensymmetry upon profile reflection with respect to the vertical axis passing through / in our exemplification we may distinguish two cases for each temperature profile: One in which the hotter temperatures are applied at the side ofthefreeendinn = 0(asinFig.)andoneinwhichhot temperatures are applied at the side of the fixed mass in n = (perform the transformation T n T n to the profiles in Fig. ). Results are respectively depicted in Fig. 3(c) and 3(d). In both cases even the linear temperature profile does not lead to equipartition. From the plots one notices that low- modes store more energy if the free end is hotter. This alludes to suggestive implications: The mixed boundary is the case considered in Ref. [8] where an experiment with a solid bar and a numerical study of an anharmonic chain displayed behaviors qualitatively consistent with that of Fig. 3(c). Our results thus suggest that the noise at lowest s would be lowered by letting the free end to float in a colder environment. This also points out a conceivable indication for reducing the measured thermal noise in experiments

4 4 T / T T n free n FIG. 4. Reconstruction of the temperature profile (lower panel) through Eq. () starting from T /T displayed in the upper panel and T = 5.5. passible to schematizations analogous to those in Fig.. IV. REVERSE EGIEERIG The expression T (T n ) may be inverted thus determining whichheat-bathtemperatureprofilest n maycorrespond to a given mode energy repartition. For definiteness let us focus on the case of free boundaries. Thans to simple identities(appendix C) the inversion of Eq.(9) gives T n +T n = = ( ) n+ cos π T +T. () otice that given T the temperature profile T n is not uniquely identified. In fact the relation T (T n ) is manyto-one for instance already on the basis of symmetry one can figure out that temperature profiles T n T n related by the transformation T n = T n produce the same energy repartition T. In the lower panel of Fig. 4 we display a profile reconstruction originated from the specific choice for T reported in the upper one. For simplicity we complemented Eq. () with the condition T n = T n ; this means in particular T 0 = T. Our example points out that in principle it is possible to design heat-bath temperature profiles so that the energy stored in the normal modes of the chain is arbitrarily distributed in the range [0T] consistently with the condition = T = T. V. POWER SPECTRUM To show that T also encodes the dynamics offluctuations we compute the power spectrum S(ω) of the chain length R R 0 + l 0 in the frequency domain a quantity typically monitored in experiments [8]. According to the Wiener-Khinchin theorem [3] under stationary conditions S(ω) is given by the Fourier transform of the chain length s autocorrelation function. Referring again to free boundary conditions in terms of normal modes we have R R 0 = = (Φ Φ 0 )X. Hence (Appendix D) S(ω) = γ 6γ = odd (Φ Φ 0 ) (Φ Φ 0 ) T (ω ω iγω) (ω ω +iγω) (3) ( ) cos (+) π T (ω ω ) +γ ω (4) (ω 0) where even modes do not contribute owing to the symmetry of the boundaries. Eq. (4) neglects the cross-correlations between modes at different. Such cross-correlation terms are instead responsible for the heat flux along the chain J n = l 0 κ ṘnR n []. In terms of normal modes we have (Appendix D) in fact J n = il 0κ (π) dω γ Φ n Φ n T ω (ω ω iγω) (ω ω +iγω) (5) (0 < n < ). We have checed that the contribution of terms with in Eq. (3) can only be appreciated in proximity of the negative peas of the power spectrum away from the resonances ω. VI. COCLUSIOS In summary our analytic study of energy repartition in a harmonic chain in contact with independent heat baths shows that both long and short wavelength modes may have energies which deviate significantly from the level expected if equipartition were to hold. This enhanced or reduced storage of energy depends critically on the shape of the temperature profile and on the boundary conditions. Other dynamical properties such as the damping or the elastic coupling are instead totally irrelevant. Thus for a generic harmonic chain information encoded in the temperature profile is mapped into a sequence of vibrational mode temperatures which in turn shape the power spectrum of the chain length. Investigations about the influence on the above picture of nonlinearities originating thermo-mechanical couplings is the next important step.

5 5 ACKOWLEDGMETS We would lie to than G. Benettin S. Lepri and R. Livi for useful discussions. We also than the Galileo Galilei Institute for Theoretical Physics for the hospitality and the IF for partial support during the completion of this wor. Appendix A: Boundary conditions. Free boundary conditions In the case of free boundary conditions with + oscillators the Laplacian matrix is A (A nm ) nm= (A) which is diagonalized by the linear transformation Φ AΦ with { = 0 Φ n = ( ) cos (n+) (+) π 0 (A) of oscillators becomes and the Laplacian matrix reads A (A nm ) nm= (A6) which now is diagonalized by (n) = (00) or (n) = () 0 (0 < n = 0) or Φ n = ( ) (0 < n = ) n sin π otherwise Φ = Φt. Also in this case it is straightforward to show that [ ( )T n= T = T ncos ( n π) ] m= T m (0 < < ) (A7) (A8) (A9) with T n= T n/( ). We have n= cos( n π) = so that at equilibrium we again recover T = T. Φ = Φt. (A3) 3. Mixed boundary conditions It is straightforward to see that the definition T () Φ n Φ n T n (A4) leads to ) T = T + T ncos(n+ + π m=0 T (0 < ) m (A5) and T 00 = T T n/( + ). otice that ) (n+ cos + π = 0 so that at equilibrium T n = T n we recover T = T.. Fixed boundary conditions In the case of fixed boundary we can stic to our notationsby fixingthe twomassesat theborder. The number In the case of mixed boundary we fix only the mass at n =. The number of oscillators becomes and the Laplacian matrix is A (A nm ) nm= (A0) which is diagonalized by the linear transformation Φ AΦ with Φ n = + (n) = () 0 (n = 0 < ) or (0 n < = ) cos ( (n+) ( +) () ) π otherwise (A)

6 6 Φ = Φt. As for the previous cases it is easy to prove that ( T = T + T (n+)(+) ncos m=0 T m + π (A) ) (A3) (0 < ) with T T n/. In this case ( ) cos (n+)(+) + π = / and again one recovers T = T at equilibrium. For the average inetic and potential energy per mode K () Ẋ ( δ 0 /) V ()ω X we thus obtain the basic result K X = V X = T / ( 0). (B) (B) (B3) Appendix C: Reconstructing the temperature profile Appendix B: Energy repartition among the modes Equation Ẍ = γẋ ω X +F (B) is a first-order linear differential equation in the vector Y (X Ẋ ). Its stationary solution is formally given by Y (t) = t with the definitions Λ = dt exp[(t t )Λ ] F (t ) ( ) 0 ω γ F = (B) ( 0. (B3) F) The matrix exponential in Eq. (B) is computed by diagonalizing Λ. Its eigenvalues λ are the two solutions of the characteristic equation for the unforced harmonic oscillator namely λ α = ) ( γ +( ) α γ 4ω α =. (B4) Therefore from the solutions X (t) = t dt A α exp( λα (t t ))F (t ) (B5) Ẋ (t) = with α= t α= dt B α exp( λα (t t ))F (t ) (B6) A = λ λ B = λ B = λ A λ = A (B7) we can evaluate the stationary equal-time correlations X X = ζt ẊẊ = ζt αβ= αβ= A α Aβ λ α +λβ B α Bβ λ α +λβ. (B8) (B9) (B0) Expressing the cosine in complex notation it is easy to prove the following identities: = ( ) ( ) (m+) (n+) cos π cos π = = = Hence from Eq. (6) we obtain or = (δ mn +δ m m ) (C) ( ) (m+) cos π = 0. (C) ( ) (m+) cos π T = (C3) ( ) (m+) = cos π T = ) + T ncos(n+ + π m=0 T (C4) m T n +T n = = ( ) n+ cos π T +T. (C5) Appendix D: Spectral density According to the Wiener-Khinchin theorem[3] under stationary conditions the spectral density S(ω) of the chain length R R 0 + l 0 is given by

7 S(ω) = dτ e iωt [R (t 0 ) R 0 (t 0 )] [R (t 0 +τ) R 0 (t 0 +τ)] +π l 0 δ(ω) 7 (D) We have R (t) R 0 (t) = (Φ Φ 0 ) X (t) (D) = so that S(ω) = = (Φ Φ 0 ) (Φ Φ 0 ) dτ e iωt X (t 0 ) X (t 0 +τ) +π l 0 δ(ω) (D3) We indicate the Fourier transform of a generic function h(t)asĥ(ω) dte iωt h(t) anddenoteitscomplexconjugate as ĥ (ω). The Fourier transform of Eq. (4) gives ω X (ω) = iωγ X (ω) ω X (ω)+ η (ω). (D4) Solving for X (ω) and using η (t)η (t ) = γt δ(t t )/() (D5) we obtain for ω 0 S(ω) = γ = (Φ Φ 0 ) (Φ Φ 0 ) T (ω ω iγω) (ω ω +iγω). (D6) The local heat flux J n along the chain [] is given by J n = l 0 κ Ṙn R n (0 < n < ). (D7) In terms of normal modes the local heat flux becomes J n = l 0 κ Φ n Φ n Ẋ X. (D8) Indeed stationarity implies Ẋ X = 0 for equal-time averages. We then have Ẋ X = (π) = (π) dω dω = i γ T (π) dω ( iω) e iωt e iω t X (ω) X (ω ) dω ( iω) e iωt e iω t γ dω (D9) T (ω ω iγω) (ω ω +iγω ) δ(ω ω ) (D0) ω (ω ω iγω) (ω ω +iγω). (D) Putting things together we obtain J n = il 0κ (π) γ Φ n Φ n T dω ω (ω ω iγω) (ω ω +iγω) (0 < n < ). (D) [] J. M. Ortiz de Zárate and J. V. Sengers Hydrodynamic Fluctuations in Fluids and Fluid Mixtures Elsevier ew Yor (006). [] R. Schmitz Fluctuations in onequilibrium Fluids Physics Reports 7 (988). [3] P.. Segrè R.W. Gammon J. V. Sengers and B. M. Law Rayleigh scattering in a liquid far from thermal equilibrium Phys. Rev. A (99). [4] W. B. Li K. J. Zhang J. V. Sengers R. W. Gammon and J. M. Ortiz de Zárate Concentration Fluctuations in

8 8 a Polymer Solution under a Temperature Gradient Phys. Rev. Lett (998). [5] C. J. Taacs A. Vailati R. Cerbino S. Mazzoni M. Giglio and D. S. Cannell Thermal Fluctuations in a Layer of Liquid CS Subjected to Temperature Gradients with and without the Influence of Gravity Phys. Rev. Lett (0). [6] A. C. Bleszynsi-Jayich W. E. Shans and J. G. E. Harris oise thermometry and electron thermometry of a sample-on-cantilever system below Kelvin Applied Physics Letters (008). [7] L. Conti M. Bonaldi L. Rondoni Rareoise: nonequilibrium effects in detectors of gravitational waves Cl. Quant. Grav (00). [8] L. Conti P. De Gregorio G. Karapetyan C. Lazzaro M. Pegoraro M. Bonaldi L. Rondoni Effects of breaing vibrational energy equipartition on measurements of temperature in macroscopic oscillators subject to heat flux J. Stat. Mech. P003 (03). [9] F. Aguilar Sandoval M. Geitner E. Bertin L. Bellon Resonance frequency shift of strongly heated microcantilevers J. Appl. Phys (05) [0] In order to identify low frequency fluctuations with low wavenumber fluctuations one has to assume a monotonic dispersion relation which is the case for acoustic phonons. [] S. Lepri R. Livi and A. Politi Heat Conduction in Chains of onlinear Oscillators Phys. Rev. Lett (997) [] S. Lepri R. Livi and A. Politi Thermal conduction in classical low-dimensional lattices Phys. Rep. 377 (003). [3] A. Dhar Heat transport in low-dimensional systems Adv. in Phys (008). [4] M. Bolsterli M. Rich W. M. Visscher Simulation of onharmonic Interactions in a Cristal by Self-Consistent Reservoirs Phys. Rev. A 086 (970). [5] F. Bonetto J. L. Lebowitz J. Luarinen Fourier s Law for a Harmonic Crystal with Self-Consistent Stochastic Reservoirs J. Stat. Phys (004). [6] D. Dhar and D. Roy Heat Transport in Harmonic Lattices J. Stat. Phys (006). [7] E. Pereira and R. Falcao ormal Heat Conduction in a Chain with a Wea Interparticle Anharmonic Potential Phys. Rev. Lett (006). [8] H.C. Fogedby and A. Imparato Heat fluctuations and fluctuation theorems in the case of multiple reservoirs J. Stat. Mech. P0 (04). [9] F. icacio A. Ferraro A. Imparato M. Paternostro and F. L. Semião Thermal transport in out-of-equilibrium quantum harmonic chains Phys. Rev. E (05). [0] K. Somiya Detector configuration of KAGRA - the Japanese cryogenic gravitational-wave detector Class. Quant. Grav (0). [] This can be viewed as a simplified description of a continuous bacground temperature field as long as the characteristic length scale of the temperature variations are large with respect to the typical fluctuations of the oscillator positions i.e. T/ T n R n. [] A. Kundu A. Chaudhuri D. Roy A. Dhar J. L. Lebowitz and H. Spohn Heat conduction and phonon localization in disordered harmonic crystals Europhys. Lett (00). [3] See e.g. M. Toda. Saito and R. Kubo Statistical Physics II: onequilibrium Statistical Mechanics (Springer Series in Solid-State Sciences 99).