arxiv: v1 [cond-mat.stat-mech] 1 May 2015

Transcription

1 Quantum Isothermal Reversible Process of Particles in a Box with a Delta Potential arxiv:155.68v1 [cond-mat.stat-mech] 1 May 15 Minho Park, Su Do Yi, and Seung Ki Baek Deartment of Physics, Pukyong National University, Busan , Korea (Received 5 November 14) Abstract For an understanding of a heat engine working in the microscoic scale, it is often necessary to estimate the amount of reversible work extracted by isothermal exansion of the quantum gas used as its working substance. We consider an engine with a movable wall, modeled as an infinite square well with a delta eak inside. By solving the resulting one-dimensional Schrödinger equation, we obtain the energy levels and the thermodynamic otentials. Our result shows how quantum tunneling degrades the engine by decreasing the amount of reversible work during the isothermal exansion. PACS numbers: 3.65.Ge,5.7.Ce,7..Pe Keywords: Schrödinger equation, Quantum thermodynamics, Isothermal exansion Electronic address: seungki@knu.ac.kr; Fax:

2 I. INTRODUCTION A fundamental idea behind classical statistical mechanics is that we are living in a macroscoic world, which means that we are unable to deal with microstates on a molecular level. This forces us to distinguish heat from work and to accet the increase of entroy on robabilistic ground, when we formulate the laws of thermodynamics. However, the border between our macroscoic world and the microscoic world of molecules is becoming vague due to technological develoments over the last century. This has led to a natural question; i.e., what would it mean to thermodynamics if it became ossible to access and maniulate microstates just as we do macrostates? The founders of statistical mechanics were already aware of this roblem: For examle, Maxwell imagined an intelligent being s intervention on a molecular level in his famous thought exeriment, now known as Maxwell s demon [1], and exressed dee concern about the foundation of the second law of thermodynamics. The Szilard engine has been regarded as the simlest imlementation of Maxwell s demon with a single-article gas []. Its cycle consists of four rocesses; inserting a wall into the center of the box, measuring the osition of the article, exanding the quantum gas isothermally, and finally removing the wall. Considering its microscoic nature, it is temting to translate the cycle into quantum-mechanical language as has been done in Refs. 3 and 4 for other engines. In fact, there are more reasons than that: Zurek [5], when he discusses Jauch and Baron s objection that the gas is comressed without the exenditure of energy uon inserting the wall, argues that one has to resort to quantum mechanics to understand the Szilard engine [6]. According to Ref. 5, this aarent inconsistency is removed by quantum-mechanical considerations. Therefore, in a sense, it is a matter of theoretical consistency. Kim and coworkers have resented a detailed account for its quantum-mechanical cycle in Refs. 7 9, with emhasis on the isothermal rocess. There seem to remain some subtle issues to settle, however, as shown in the debate on how to deal with the quantum tunneling effect through the wall in a three-boson case [1]. In this work, we describe the isothermal exansion of the quantum Szilard engine by considering a quantum gas confined in a one-dimensional cylinder [11, 1]. Although some revious studies, such as Refs , assume that the wall is imenetrable, we will relax that assumtion because it is, strictly seaking, exerimentally infeasible. Starting from the Schrödinger equation, we calculate the energy levels and thereby obtain how much work can 1

3 be extracted by an isothermal reversible exansion. We find the following: The amount of reversible work is significantly reduced comared to the claims in Refs. 5 and 7 if we assume full thermal equilibrium inside the box, as is usual in describing isothermal exansion. This conclusion should not be affected by either the wall s insertion or the measurement before the isothermal exansion because thermal equilibrium does not remember its history. Moreover, if we ut two or three fermions in the box, the free-energy landscae is not monotonic with resect to the wall s osition, so that we should sometimes erform work to exand the gas. These effects have not been exlicitly discussed in revious studies. This work is organized as follows: An exlanation of our basic setting in terms of the Schrödinger equation is given in Section II. A numeric calculation of the energy levels and the free energy and a comarison of a single-article case to the two- and three-article cases is resented in in Section III. This is followed by discussion of results and conclusions. II. PARTICLE IN A BOX Consider the following one-dimensional otential landscae with size L: aδ(x L) if < x < L V(x) = elsewhere, (1) where a is the strength of the delta otential and (,1) secifies its location [16]. We exect the Schrödinger equation ] [ d mdx +V(x) ψ(x) = Eψ(x) () to have a solution with a sinusoidal form. The boundary conditions at x = and x = L are satisfied if we set Asinkx if < x < L ψ(x) = Bsink(L x) if L < x < L, where the coefficients A and B, as well as the wavenumber k, are assumed to be nonzero. The eigenfunction ψ(x) has the following roerties: First, it is continuous over [, L]; i.e., (3) AsinkL = BsinqkL (4)

4 with q 1. Second, the change in the derivative of ψ(x) with resect to x aroundx = L is obtained by using { L+ǫ [ d lim ǫ L ǫ dx ψ(x) from which it follows that ( dψ lim ǫ dx dψ L+ǫ dx ] dx+ m L ǫ ) [ L+ǫ L ǫ Eψ(x)dx L+ǫ L ǫ ]} aψ(x)δ(x L)dx =, (5) = ma ψ(l). (6) Plugging Eq. (3) here and multilying both sides by sinqkl to use Eq. (4), we obtain ksinkl = ma sinklsinqkl. (7) This formula has also been derived in Ref. 1 in a different way. Equation (7) can be numerically solved by using the Newton-Rahson method to yield the allowed values of k. The method works as follows: First, let us define f(k) ksinkl + ma sinklsinqkl, and look for its zeros to solve Eq. (7). Starting from k =, we check whether the sign of f(k) changes with increasing k by a sufficiently small amount, say, dk 1. Every time the sign changes between k = k > and k +dk, we run the following iteration starting from k () = k : k (j+1) = k (j) f[k(j) ] f [k (j) ], (8) where j =,1,,... When k (j) has converged to a stationary value k ( ), i.e., if k (j) = k (j+1) within numerical accuracy, we take it as a solution of Eq. (7) and reeat the rocedure by checking the sign of f(k) from k = k ( ) + dk. Each of these solutions is identified with the wavenumber k n of the nth eigenmode in an ascending order. For the alication of the Newton-Rahson method numerically, it is convenient to exress the quantities in dimensionless units, so we divide Eq. () by E () g π /(ml ), the ground state energy for a =, to obtain [ ] d dη ψ(η)+φ(η) ψ(η) = ǫψ(η), (9) where η πx/l. Then, the otential inside the box is exressed as Φ(η) = ml aδ(x L) = αδ(η π), (1) π 3

5 κ ψ(η) 6 4 (a) (c) =.5 ε ε 1 ε ε ψ(η) (b) (d) =.4 ε ε 1 ε η π η π FIG. 1: (Color online) (a) Lowest wavenumbers in units of π/l as a function of, and (b) the corresonding energy levels in units of E () g = π /(ml ). The other two anels show the robability densities for the three lowest eigenmodes (c) when the wall is in the middle, =.5, and (d) when it has moved to =.4. The strength of the otential is set as α = 1. with α aml/( π) = a(l/π)/e () g. We can also define a dimensionless wavenumber κ Lk/π to rewrite Eqs. (7) as κsinκπ+αsinκπsinqκπ =. A little algebra shows that the energy level is given by ǫ = κ + α [ ( csc κπ 1 ) ( π κπ sinκπ +csc qκπ q 1 )] 1 κπ sinqκπ, (11) where the second term reresents the otential energy coming from the overla between the wavefunction and the delta eak otential. The numerical solutions for these κ and ǫ are deicted in Fig. 1, together with the robability density lots. III. REVERSIBLE WORK A. Single Particle We are interested in how much work can be extracted by erforming a reversible isothermal rocess with this system. Let ǫ n denote the nth energy level obtained by using the Newton-Rahson method with n = 1,...,M so that the ground-state energy is denoted by ǫ 1. For a given value of temerature T = (k B β) 1, where k B is the Boltzmann constant, we 4

6 F()-F(.5) T=1-3 α= FIG. : (Color online) Change in the free energy F in units of E () g = π /(ml ) when the wall, starting from the center of the box, moves to a new osition. The height of the wall is fixed at α = 1, and each curve reresents a different T in units of E () g /k B. The number of energy levels considered is ket as M = 1 3 throughout the free-energy calculations. run the calculation to get βǫ M 1. The artition function for a single quantum article can then be obtained as Z = lim M M e βǫn. n=1 (1) The free energy is F = k B T lnz, accomanied by a differential form df = SdT PdV, where S, P, and V are the entroy, ressure, and volume, resectively. The classical concets, such as ressure and work, need careful consideration in quantum mechanics [17]. To calculate work, we begin with the occuation robability of the nth level given by ρ n = e βǫn = k BT ( Z/ ǫ n ) n Z e βǫn = ǫ n ( k B T lnz) = F ǫ n. (13) The change in the internal energy U = n ρ nǫ n is then exressed as du = n ǫ ndρ n + ρ n dǫ n = dq+ dw, where dq and dw are differentials of heat and work, resectively. The quantum thermodynamic work is identified with dw = n ρ ndǫ n because heat is usually involved with changes in occuation robabilities for given energy levels whereas work is erformed when the energy levels themselves change [7]. Consequently, the amount of reversible work during this rocess is W = n with the classical case. F ǫ n dǫ n = df = F, which is consistent Figure shows F resulting from moving the wall, which was initially located at the center, =.5. Suose that T is so low that T E () g /k B. When the wall is at =, only the ground state contributes to the summation in Eq. (1), which means that Z( = ) e βe() g. When the wall is at =.5, the ground-state energy equals 4E () g, and the 5

7 first excited state also contributes to the artition sum in Eq. (1) because it lies close to the ground state [see Fig. 1(b)]. The artition function is, therefore, aroximated as Z( =.5) e 4βE() g. The free-energy difference in this low-t region is, thus, given as F() F(.5) 3E () g + k B T ln. When T = E () g /k B, for examle, the free-energy difference roughly amounts to.3e () g, which exlains the size of the free-energy dro in Fig.. As T increases, a lateau develos near the center of the box because F does not resond much to the wall s osition. However, the dro in F at the end of the rocess increases with increasing T if T O(1) in units of E () g /k B. Our question is how it grows with T. Let us take α and γ β π /(ml ) 1. Then, it is the artition function for the initial and final wall ositions can be aroximated. The latter case of the final wall osition is estimated as Z( = ) = = e γn (14) n=1 [ e γn e γ(n+1)] (15) ex( γn )dn 1 n= π 4γ 1, (16) where the summation in Eq. (15) is a first-order correction to the integral aroximation. Likewise, the former case of the initial wall osition gives Z( =.5) = n=1 [ e 4γn ex( 4γn )dn 1 ] π = 1, (17) 4γ where the factor of in front of the summation is due to the fact that every air of adjacent energy levels becomes degenerate in the limit of α. This effect cancels the factor of 4 inside the exonential arising from the reduced volume L/. Although F does not comletely vanish due to the correction, the imortant oint is that Z(.5)/Z() for γ 1 is far less than, in contrast with the claim that F = k B T ln in Refs. 5 and 7. Because α must be finite in any exerimental situation, the article can, in full thermal equilibrium, be observed on either side of the box, which reduces the amount of extracted work. On the other hand, the revious studies on the quantum Szilard engine [5, 7] have assumed that the rocess is erformed within a shorter time scale than required for tunneling through the wall and, therefore, concluded that Z( =.5) ex( 4γn )dn without the factor of in 6

8 F()-F(.5) F()-F(.5) (a) F()-F(.5) T= Two Bosons (c) T=1 3 1 Three Bosons F()-F(.5) -5 (b) T=1 3 1 Two Fermions (d) T=1 3 1 Three Fermions FIG. 3: (Color online) Free-energy changes in the resence of multile non-interacting quantum articles: (a) bosonic and (b) fermionic cases with two articles. With three articles, (c) the bosonic case remains qualitatively similar, but (d) the fermions behave differently due to the Pauli exclusion rincile. front. Strictly seaking, their engine undergoes an isothermal rocess only when artially equilibrated to maintain the confinement and extract W = k B T ln. B. Two- and Three-article Cases If two quantum articles are in the box, we should consider their symmetry, i.e., whether they are fermions or bosons. Without any consideration of the sin degree of freedom, the elements of the density matrix for the two articles are exressed in the bracket notation as follows: ǫ 1ǫ e βĥ ǫ 1 ǫ = e ( ) β(ǫ 1+ǫ ) 1 ǫ 1 ǫ ± 1 ǫ ǫ 1, where δ is the Kronecker delta symbol, and we have a lus sign for bosons and a minus sign for fermions. The two-article artition function is then obtained by taking the trace oeration: ( ) Z = Tr e βĥ = 1 e β(ǫ1+ǫ) (1±δ ǫ1 ǫ! ). (18) ǫ 1,ǫ The same alies to the three-article case. The density-matrix elements are written as ǫ 1ǫ ǫ 3 e βĥ ǫ 1 ǫ ǫ 3 = e ( β(ǫ 1+ǫ +ǫ 3 ) 1 ǫ 1 ǫ 3 ǫ 3 ±... ), and the artition function reads ( ) Z = Tr e βĥ = 1 e β(ǫ1+ǫ+ǫ3) (1±δ ǫ1 ǫ 3! ±δ ǫ ǫ 3 ±δ ǫ3 ǫ 1 +δ ǫ1 ǫ δ ǫ ǫ 3 ). (19) ǫ 1,ǫ,ǫ 3 7

9 Figure 3 shows the free-energy differences for the two- and the three-article cases. The bosonic cases are qualitatively similar to the single-article case. However, we see very different behavior for the fermionic cases where the Pauli exclusion rincile is in action: At low T, we need to erform ositive work to move the wall from the center to a certain osition [Figs. 3(b) and (d)]. Obviously, the reason is that it costs free energy to lace two fermions close to each other. IV. SUMMARY AND OUTLOOK In summary, we have considered an isothermal exansion rocess of a quantum gas, taking the tunneling effect into consideration. We have found that the amount of reversible work is smaller than k B T ln in the high-t region if the system is fully equilibrated at every moment of the isothermal exansion. The difference from Refs. 5 and 7 arises because they have searated the time scale of tunneling from that of the artial equilibration that occurs on only one side of the wall. Although it was not exlicitly stated in revious studies, this searation may be a lausible assumtion for the following reason: As we increase the otential height α, we may well exect the time scale for tunneling to grow whereas the artial equilibration before tunneling is achieved within a finite amount of time. The quantum Szilard engine will show its exected erformance only between these two time scales. The question is, then, how large a value of α one should have to ensure the searation of the time scales, which will be ursued in our future studies. Acknowledgments This work was suorted by a research grant from Pukyong National University (14). [1] K. Maruyama, F. Nori, and V. Vedral, Rev. Mod. Phys. 81, 1 (9). [] T. Sagawa and M. Ueda, in Nonequilibrium Statistical Physics of Small Systems: Fluctuation Relations and Beyond, edited by R. Klages, W. Just, and C. Jarzynski (Wiley, Weinheim, 13), [3] C. M. Bender, D. C. Brody, and B. K. Meister, J. Phys. A 33, 447 (). 8

10 [4] T. D. Kieu, Phys. Rev. Lett. 93, 1443 (4). [5] W. Zurek, in Frontiers of Nonequilibrium Statistical Physics, NATO ASI Series, Vol. 135, edited by G. T. Moore and M. O. Scully (Plenum, New York, 1986), [6] J. M. Jauch and J. G. Baron, Helv. Phys. Acta 45, (197). [7] S. W. Kim, T. Sagawa, S. D. Liberto, and M. Ueda, Phys. Rev. Lett. 16, 741 (11). [8] K.-H. Kim and S. W. Kim, Phys. Rev. E 84, 111 (11). [9] K.-H. Kim and S. W. Kim, J. Korean Phys. Soc. 61, 1187 (1). [1] M. Plesch, O. Dahlsten, J. Goold, and V. Vedral, Phys. Rev. Lett. 111, (13). [11] G. A. Vugalter, A. K. Das, and V. A. Sorokin, Phys. Rev. A 66, 114 (). [1] P. Pedram and M. Vahabi, Am. J. Phys. 78, 839 (1). [13] C. M. Bender, D. C. Brody, and B. K. Meister, P. Roy. Soc. A 461, 733 (5). [14] H. Dong, D. Z. Xu, C. Y. Cai, and C. P. Sun, Phys. Rev. E 83, 6118 (11). [15] H. Li, J. Zou, J.-G. Li, B. Shao, and L.-A. Wu, Ann. Phys. 37, 955 (1). [16] D. J. Griffiths, Introduction to Quantum Mechanics (Prentice Hall, Uer Side River, NJ, 1995), [17] P. Borowski, J. Gemmer, and G. Mahler, Eurohys. Lett. 6, 69 (3). 9