Dimensional transitions in thermodynamic properties of ideal Maxwell-Boltzmann gases
An ideal Maxwell-Boltzmann gas confined in various rectangular nano domains is considered under quantum size effects. Thermodynamic quantities are calculated from their relations with partition function which consists of triple infinite summations over momentum states in each direction. To get analytical expressions, summations are converted to integrals for macro systems by continuum approximation which fails at nanoscales. To avoid both from the numerical calculation of summations and the failure of their integral approximations at nanoscale, a method which gives an analytical expression for single particle partition function (SPPF) is proposed. It’s shown that dimensional transition in momentum space occurs at certain magnitude of confinement. Therefore, to represent SPPF by lower-dimensional analytical expressions becomes possible rather than numerical calculation of summations. Considering rectangular domains with different aspect ratios, comparison of the results of derived expressions with those of summation forms of SPPF is done. Itâs shown that analytical expressions for SPPF give very precise results with maximum relative errors of around , and at just the transition point for single, double and triple transitions respectively. Based on dimensional transitions, expressions for free energy, entropy, internal energy, chemical potential, heat capacity and pressure are given analytically valid for any scale.
Keywords: Quantum size effects; Dimensional transitions; Nano thermodynamics; Dimensional superpositions
Istanbul Technical University, Energy Institute, 34469, Istanbul, Turkey.
Advancements in nanoscience and nanotechnologies increased the importance of studies related to the nanoscale properties of matter, in recent years. On the presence of quantum confinement effects, thermodynamic properties of systems differ from those of macro scale. Due to wave character of particles, quantum size effects (QSE) appear at nanoscale systems as a consequence of quantum confinement.
Thermodynamics under QSE became an active research area and many novel aspects such as quantum boundary layer, anisotropic gas pressure, loss of additivity in extensive properties, quantum forces and gas diffusion due to size difference in both Maxwell-Boltzmann (MB) and quantum (Fermi-Dirac and Bose-Einstein) gases have been studied in literature[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17]. Even though size and shape dependencies in thermodynamics of nano systems has been recently studied, dimensional transitions due to QSE in thermodynamic quantities and transition conditions have not been examined in literature.
Thermodynamic state functions are represented by summations from their definitions in statistical mechanics[18, 19]. In thermodynamics of macro systems, it is common to replace summations by integrals to decrease the calculation burden and to obtain analytical expressions. Conversion of summations to integrals is done by using continuum approximation and density of states concept, under the assumption of discrete energy levels are so close to each other that they can be considered as continuous variables. Continuum approximation works extremely well in macro scale. However, in nanoscale, ascending divergence appears between summations and their integral approximations[7, 19].
At meso scales, instead of integral approximations, Weyl’s conjecture or first two terms of Poisson summation formula (PSF) are used in literature[2, 3, 4, 5, 7, 8, 9, 10, 11, 12, 19]. On the other hand, when characteristic sizes of the domain are smaller than half of the most probable de Broglie wavelengths of particles, even these approximations considered in literature become insufficient. In the case of such strong confinements, one should either consider the numerical calculation of exact summations in thermodynamic properties or look for a better solution.
In this article, a method based on dimensional transitions in momentum space due to quantum confinement is proposed for the calculation of single particle partition function (SPPF) of MB statistics. It is seen that the method leads to analytical expressions which give nearly the same results of exact summations. Then from the analytical expressions of SPPF, thermodynamic quantities such as free energy, entropy, internal energy, heat capacity and pressure are derived analytically. Consequently, the results successfully represent thermodynamic behavior of strongly confined systems without need to calculate relevant triple summations.
2 Determination of Dimensional Transition Point in Momentum Space
For a high temperature and/or low density ideal gas, SPPF is given in triple summation form by using MB statistics as follows
where are momentum state variables with subscript of indicating directions, are the energy eigenvalues from the solution of the Schrödinger equation for a rectangular domain, is Boltzmann’s constant and is temperature. By defining a confinement parameter that indicates the magnitude of confinement in ’th direction, Eq. (1) can also be written as follows
In the definition of , is the size of the domain in the ’th direction and is the half of the most probable de Broglie wavelength of particles where is Planck’s constant and is mass of the particle.
Although triple summations are the exact way to calculate SPPF, in macro scales where values are much smaller than unity (), summations are replaced by integrals making continuum approximation. On the other hand, when domain size in particular direction () is on the order of , like in nanoscale (), there is a considerable deviation between the results of summations and integrals. Therefore, a more precise evaluation method, such as Poisson summation formula (PSF), is necessary to calculate summations. PSF is written for even functions, , as
By evaluating each term of PSF, single particle partition function in 1D can be exactly represented as
The first term of PSF is the conventional integral term and it’s accuracy for macro domains is quite well. When domain size gets closer to , contribution of the second term (called zero correction term) of PSF becomes appreciable. First two terms of PSF can be called one-dimensional (1D) representation, since it represents the contributions of whole momentum state variables, . Contribution of the zero correction term is studied as QSE in literature[7, 8, 10, 11, 12, 15, 16] especially for weakly confined domains in which the contribution of the third term of PSF is negligible. To fully represent SPPF in any scale, even for strongly confined domains , whole terms of PSF have to be considered. However, since there is no known analytical solution for the third term which include another infinite summation, either we have to calculate infinite summations of SPPF numerically or find a method to give analytical expressions representing the exact behavior of SPPF successfully. For strongly confined cases , fortunately, contribution of just the first term () of SPPF’s summation becomes extremely dominant. Therefore, it can be called as zero-dimensional (0D) representation since it represents almost only the contribution of the ground state.
To see the success of 1D and 0D representations, variations of their analytical expressions with confinement parameter are compared with that of exact 1D SPPF based on numerical calculation of summation, in Figure 1. The solid curve shows the variation of exact form of 1D SPPF (Eq. 4a) with confinement parameter . Plus signs show the results of 1D analytical representation based on first two terms of PSF, while cross signs demonstrate the results of 0D representation based on just the first term () of SPPF’s summation.
It is understood that 1D analytical expression matches perfectly with the exact form of SPPF until the certain value of , after that it considerably deviates. In the same manner, 0D expression substantially deviates before , while it perfectly matches with the solid curve after that. This certain value of , which corresponds to the transition point from 1D to 0D representation, can be found as by equating 1D expression to 0D one. Therefore, an analytical expression representing SPPF for any values of can be given as follows
where the transition value is . This approximation gives error even at the transition point. The transition represents the dimensional change (contraction/expansion) of momentum space in particular direction. Henceforward, the domain in the ’th direction will be considered as free and confined for and respectively.
2.1 Dimensional Transitions and Superpositions in Single Particle Partition Function
By using the dimensional transition point in momentum space, we may now consider SPPFs for the domains which are confined at nanoscale in certain directions. 3D domain is free in all three directions, so that confinement values in all directions are smaller than , the transition value. For 2D domain, and for 1D domain, and are chosen as greater than respectively. 0D domain is considered as confined in all directions, which means all of their confinement parameters are larger than . In this sense, using Eq. (5), we can write analytical expressions for SPPF of an ideal MB gas for 3D, 2D, 1D and 0D domains as
Although dimensional transition point has a unique value, by engaging the sizes of the domain to each other in each direction, it is possible to obtain different transition points correspond to different dimensional transitions for anisometric domains. In order to engage sizes to each other, we may define constant aspect ratios between sizes in three directions. Therefore, we can change the size of the domain by changing only one of the sizes, which is chosen here as , while preserving its shape. Aspect ratio between sizes in first and second directions is denoted by and while the second and third directions by . In that case, confinement parameters of second and third directions can be written in terms of the confinement parameter in the first direction as and respectively. In Figure 2, relative errors () of analytical expressions given by Eq. (6) are given for different aspect ratios.
Peaks in Figure 2 correspond to the values , and for transitions from 3D-2D, 2D-1D and 1D-0D respectively. As it is seen from Figure 2a, each dimensional transition causes an error of even at the transition points. If the sizes in two different directions are the same, dimensional transitions in these directions superpose and it causes error as it is seen in Figures 2b and 2c. Similarly, for an isometric domain, superposition of three different dimensional transitions cause error, Figure 2d.
3 Dimensional Transitions in Thermodynamic Properties
Free energy of an ideal MB gas consisting of particles is simply given by [18, 19]. Dimensionless free energy can be defined as . Then, we can obtain analytical expressions for dimensionless free energies of a MB gas for 3D, 2D, 1D and 0D cases by considering dimensional transitions
respectively, where with as direction index. It is clear that is always smaller than since is necessary only for the case of . Note that, we can use MB statistics if and only if condition holds. This condition holds even if one of the confinement parameters is smaller enough than unity. Because in that condition the system contains enough particles to do statistics and interparticular distance is much greater than to use MB statistics. Since 0D case does not satisfy this condition, one could not use MB statistics for 0D conditions. Nevertheless, in this article, we gave also 0D expressions of thermodynamic properties, just to give a mathematical intuition. Therefore, it is worth to note that thermodynamic expressions given for 0D do not represent a meaningful physical situation since the condition of leads for 0D condition.
From free energies, we can derive all other thermodynamic quantities such as entropy, internal energy, chemical potential, heat capacity at constant sizes and pressure as
Thermodynamic properties can be written in their dimensionless forms as , , , and where . Even though pressure is a scalar quantity in macro scale, it becomes a tensorial quantity with zero off-diagonal elements for anisometric domains. In other words, for an anisometric rectangular domain where .
Dimensionless entropy, internal energy, chemical potential, heat capacity and pressure expressions for MB gases are derived by considering dimensional transitions and given in Table 1.
|Heat capacity at constant sizes|
It is easy to check that, all expressions derived in this article go to their classical expressions when confinement parameter ’s go to zero. So they are valid not only for nano systems but also for meso and macro systems. They are generic in that sense, within the error limit at dimensional transition points.
Note for heat capacity that, instead of referring to constant volume, we refer here constant sizes which means size of each direction , and are constant. Otherwise, combination of different values can also correspond to the same volume, although they correspond different heat capacity values, different from that in macro scale.
Unlike the scalar behavior of pressure in macroscopic scale, quantum confinement in one direction brings a positive contribution to the pressure in that direction which make the pressure a tensorial quantity in nanoscale. Besides, when confinement of a particular direction is larger than , pressure at that direction grows quadratically.
By introducing dimensional transitions of thermodynamic properties for MB gases, it’s shown that it is possible to obtain analytical expressions within a negligible error to represent the exact forms of thermodynamic properties based on triple infinite summations. In this regard, analytical results are obtained for SPPF, free energy, entropy, internal energy, chemical potential, heat capacity and pressure, for 3D, 2D, 1D and 0D cases. It is seen that thermodynamic quantities have strict size dependency which cause dimensional contraction in free momentum space for nano systems.
Critical value of confinement parameter is found as . When the confinement parameter in a particular direction becomes greater than critical value of transition, momentum component in that direction can only take its ground-state value. Then, degree of freedom of momentum space decreases one unit. Therefore, depending on the values of confinement parameter in each direction, degree of freedom takes the values of and . By considering this dimensional transition model, it is possible to express the thermodynamic properties analytically with error for any confined system.
As it is expected, classical expressions can be recovered in case of the absence of confinement, . Dimensional transitions in ideal Fermi and Bose gases are under consideration. Thermodynamic properties of strongly confined ideal Fermi and Bose gases in nano domains will then be analytically obtained for any scale.
- thanks: NOTICE: This is an author-created, un-copyedited version of an article accepted for publication in Physica Scripta. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at doi:10.1088/0031-8949/90/4/045208.
- Y. Volokitin, J. Sinzig, L. J. De Jongh, G. Schmid, M. N. Vargaftik, and I. I. Moiseevi. Quantum-size effects in the thermodynamic properties of metallic nanoparticles. Nature, 384:621–623, 1996.
- M. I. Molina. Ideal gas in a finite container. Am. J. Phys., 64:503–505, 1996.
- R. K. Pathria. An ideal quantum gas in a finite-sized container. Am. J. Phys., 66:1080–1085, 1998.
- W. S. Dai and M. Xie. Quantum statistics of ideal gases in confined space. Phys. Lett. A, 311:340–346, 2003.
- W. S. Dai and M. Xie. Geometry effects in confined space. Phys. Rev. E, 70:016103, 2004.
- H. Yu, J. Li, R. A. Loomis, L. W. Wang, and W. E. Buhro. Two- versus three-dimensional quantum confinement in indium phosphide wires and dots. Nature Materials, 2:517–520, 2003.
- A. Sisman and I. Müller. The casimir-like size effects in ideal gases. Phys. Lett. A, 320:360–366, 2004.
- A. Sisman. Surface dependency in thermodynamics of ideal gases. J. Phys. A- Math. Gen., 37:11353–11361, 2004.
- H. Pang, W. S. Dai, and M. Xie. The difference of boundary effects between bose and fermi systems. J. Phys. A- Math. Gen., 39:2563–2571, 2006.
- A. Sisman, Z. F. Ozturk, and C. Firat. Quantum boundary layer: a non-uniform density distribution of an ideal gas in thermodynamic equilibrium. Phys. Lett. A, 362:16–20, 2007.
- C. Firat and A. Sisman. Universality of quantum boundary layer for a maxwellian gas. Phys. Scripta, 79:065002, 2009.
- C. Firat, A. Sisman, and Z. F. Ozturk. Thermodynamics of gases in nano cavities. Energy, 35:814–819, 2010.
- E. O. Karabulut, M. Koyuncu, U. Atav, and M. Tomak. Thermodynamics of a two-dimensional interacting bose gas trapped in a quartic potential. Phys. Lett. A, 375:339–345, 2011.
- H. Pang, W. S. Dai, and M. Xie. Thermodynamics of a two-dimensional interacting bose gas trapped in a quartic potential. J. Phys. A-Math. Theo., 44:365001, 2011.
- A. Sisman and G. Babac. Quantum size effects on classical thermosize effects. Continuum Mech. Thermodyn., 24:339–346, 2012.
- C. Firat and A. Sisman. Quantum forces of a gas confined in nano structures. Phys. Scr., 87:045008, 2013.
- A. Aydin and A. Sisman. Discrete nature of thermodynamics in confined ideal fermi gases. Phys. Lett. A, 378:2001–2007, 2014.
- C. Kittel and H. Kroemer. Thermal Physics, 2nd Ed. W. H. Freeman, 1980.
- R. K. Pathria. Statistical Mechanics, 2nd Ed. Pergamon Press, 1996.