PACS: 05.20; 05.70; 65.80
Introduction
Due to the current boom of physics at a nanoscale, in this work we study a simple nanocrystal model inspired in the well-known Einstein model of solids 1-4. The use of highly idealized models is a common practice in physics because, even though they do not allow for precise quantitative predictions, their simplicity enables a clear qualitative understanding of a problem. The nanocrystal model proposed can be analyzed exactly, thus showing clearly which factors contribute to its differing behavior from that of the macroscopic solid. In a previous paper 5 we analyzed the Einstein solid and the two-level system with few elements (10-100); here, we want to analyze the effect of a surface. Therefore, we first need to compare the few-particle effect and the surface effect. Also, certain conceptual concerns should be taken into account, such as the additivity of the entropy and the equivalence between results obtained with the microcanonical and canonical ensembles. It is known that both ensembles lead to the same results in the thermodynamic limit, but not when dealing with finite systems. An introduction to this issue can be found in 6 and suitable teaching examples in Ref. 5.
This article continues in Sec. 2, which introduces the model and the equations that will be used throughout this paper. Section 3 presents the results: first we compare the few-particle and surface effects; then, we show how the specific heat converges to the value of the thermodynamic limit when the nanocrystal size is increased; the additivity of the entropy is dealt with next; and finally, we show how the results of the microcanonical ensemble differ from those of the canonical ensemble. Last, Sec. 4 summarizes the results.
1. Model
The Einstein solid, a well-known example for any advanced undergraduate student 1-4, is the starting point of our model. We intend to model a nanocrystal, i.e. a crystal formed by a few tens or hundreds of atoms. Therefore, we assume it is a simple cubic lattice with nL atoms per side, spaced by a distance a 0. Each atom is linked to its nearest neighbors through a quantum “spring” (oscillator) -see Fig. 1. It is assumed that all the oscillators have the same characteristic frequency ν 0. Then, the total number of atoms NA , the volume V of the nanocrystal and the number of oscillators N are:
It should be noted that the inner and surface atoms have a different number of neighbors, so that the number of quantum oscillators is not simply 3NA . A little thought shows that Eq. (1c) gives the right number of oscillators. When dealing with the Einstein solid, textbooks are interested in the thermodynamic limit and, therefore, the surface is not considered and the number of oscillators is simply 3NA .
The calculations will follow the microcanonical ensemble. The reason for this is that the nanocrystal is considered to be isolated and its energy constant. Using the canonical ensemble, the system would be connected to a reservoir at a constant temperature. However, it should be noted that, given the size of the nanocrystal, such interaction would get the system continuously out of equilibrium. Let us imagine the nanocrystal is in contact with a gas at room temperature: each impact of a gas molecule with the few-hundred-atom nanocrystal is a major event that would keep the system permanently out of equilibrium. This is why the proper way to study the nanocrystal is within the microcanonical ensemble.
The frequency of each oscillator is v 0, thus, the entropy S of a cubic crystal with N oscillators and a total energy E energy is 1,3,4:
Notice that (E/hν 0) is the total number of energy quanta and it is an integer M. A detailed derivation of Eq. (2) is a standard topic in any statistical mechanics course. The usual trick is to evaluate how many different ways there are to put M balls in N boxes and the problem can be reduced to calculate the number of permutations of (N - 1 + M) objects (M balls and N - 1 rods) taking into account that balls (rods) are indistinguishable. The usual treatment of the problem in textbooks considers N = 3NA , uses Stirling’s approximation with factorials, and then applies the usual thermodynamic equations for temperature (1/T = ∂S/∂E) and specific heat (Cυ = ∂E/∂T). However, a previous paper5 showed that, when working with few particles, factorials must be treated exactly and derivatives replaced by finite differences to account for the discrete nature of energy. On one side, Stirling’s approximation is valid for N large but it is not obvious that it is also right for small systems. For this reason it is necessary to calculate the factorial exactly. On the other side, a central result of quantum physics is usually forgotten and the energy of a system of oscillators is taken as a continuous variable. Such assumption is incorrect since it is known since the early twentieth century that the energy of an oscillator is quantized, that is, there are M, M + 1,…etc. energy quanta M being an integer given by M = E/hν 0. Strictly speaking one should deal with finite differences, then ∂S → ΔS = S(N, M + 1) - S(N,M) y ∂E → ΔE = (M + 1) - M = 1. Consequently, the temperature T and the specific heat Cυ need to be computed as follows:
Note that (3a) and (3b) are not approximations, but represent reality; only when working with large numbers for M, the discrete nature of energy can be disregarded, taking energy as a continuous variable and using derivatives.
One improvement of the model would be to consider different frequencies for the oscillators on the surface and for those inside the crystal. This assumption certainly makes sense since it remarks the difference between the surface and the volume inside the crytal; however this modified model cannot be studied in the microcanonical ensemble. Suppose that the surface oscillator frequency is ν'. Therefore two kinds of energy quanta should be considered: the M′ quanta with energy hν' and the M ones with energy hν 0. In the microcanonical formalism the system is isolated and M and M' are constant, consequently the volume and surface temperatures would be different but the system cannot reach thermal equilibrium since the energy quanta of the surface oscillators are not equivalent to those inside the crystal. If one is interested in analyzing this modified model, the canonical formalism has to be used: the nanocrystal is in contact with a heat reservoir that provides energy quanta of both frequencies and the surface and the inside volume can reach the same temperature.
2. Results
2.1. Few particle (FP) and finite size (FS) effects
We must distinguish now between two types of effects occurring when doing statistical mechanics with few elements. We analyzed a few-particle model in Ref. 5 without taking into account the surface, so that the number of oscillators was considered to be simply three times the number of atoms (N = 3NA ). The difference between these results and those obtained in the thermodynamic limit are due exclusively to the small number of particles in the system, so accordingly we will call this the few-particle (FP) effect. Now that we are considering the surface, the number of oscillators must be carefully recorded and, for this, the value for N is that given by Eq. (1c). Thus, we can speak of the finite-size (FS) effect, which includes the few-particle the surface effects. The first question that arises is, then, which effect is more relevant: that due to the few particles or the one due to the presence of a surface?
To answer this question, we evaluated numerically but exactly, i.e. without any approximations, the entropy given by Eq. (2). In one case, we consider N = 3NA and obtain an entropy SFP that only includes the few-particle effect. In the other case, the number of oscillators N is taken from Eq. (1c) and we obtain an entropy SFS that includes both the few particles and the surface. The difference between both entropies gives the surface effect. The reference value is the entropy Sth of the classical Einstein solid in the thermodynamic limit 1,4:
Once the entropies SFP
and SFS
are known, it is possible to find the temperatures TFP
and TFS
, and the specific heats
Now, we can define the relative effect for each of the thermodynamic quantities. If X represents the quantity of interest (S, T, Cυ ), the relative effect associated with the few particles (ΔXFP ) and that associated with the finite size (ΔXFS ) can be defined as follows:
Thus, it is clear which is the impact due to the small number of particles and that due to the surface. The results can be seen in Fig. 2.
Figure 2a shows the difference in entropy for two sizes of the few-particle system (nL = 3 and nL = 10). Consistently with the findings of 5, there is no discernible difference between the exact value and that of the thermodynamic limit for the system with 1000 atoms. On the contrary, Fig. 2b considers the surface and its effect is considerable in the small system (ΔSFS ∼ 22%) as well as in the large system (ΔSFS ∼ 5%). This demonstrates the impact of the presence of the surface in the thermodynamics of the nanocrystal.
The question arises of how this effect is shown in experimentally observable quantities. The answer is in Figs. 2c and2d. Figure 2c shows the impact on the specific heat due to the few particles effect. For the system with 1000 atoms, the specific heat is not different from the thermodynamic value. On the contrary, when considering the surface -Fig. 2d-, the specific heat is lower (
In conclusion, the surface plays a much more important role that the finite number of particles in the thermophysical properties of the nanocrystal.
One may ask whether it makes sense to evaluate the properties of a system using the exact approach used in this work. The answer is: it depends on the size of the system under study. For atomic clusters with tens of atoms, the calculation should be made exactly with the microcanonical formalism as there is a noticeable difference with the thermodynamic values. Conversely, if one has a system with thousands of atoms, the conventional results obtained with the canonical formalism can be used.
2.2. Specific heat scaling
Figure 3a illustrates the atomic specific heat of the nanocrystal as a function of temperature for the different sizes, and it also shows the curve of the usual Einstein solid. It is clear that the nanocrystal behavior is qualitatively similar to that of the macroscopic solid: with high temperatures, the specific heat reaches a constant value approaching the Dulong-Petit law. It should be highlighted that, in the units used in this article, Cυ /NA goes to 3 for the solid in the thermodynamic limit. Quantitatively, the specific heat of the nanocrystal is lower than that of the solid. It is interesting to illustrate the asymptotic value reached by the Cυ of the nanocrystal as a function of its size. As seen in Fig. 3b, it can be inferred the following scaling relationship:
This scaling is caused by the presence of the surface in the crystal. In fact, Eqs. (1a) and (1c) give that the number of oscillators can be written as
2.3. Non-additivity of the entropy
The entropy in the thermodynamic limit is an extensive property; however, this may not be so at a microscopic scale since the range of interactions is comparable to the size of the system. To verify if the entropy is additive or not in our nanocrystal model, we carried out the following numeric experiment -see Fig. 4a-: we took a nanocrystal B of side
Similarly, the energy quanta M of the original nanocrystal were distributed evenly among the 8 new nanocrystals. The question is whether the initial and final entropies are equal; and the answer is no. We find that S(B) > 8S(A). Figure 4b shows the difference in entropy δS = S(B) - 8S(A) as a function of the size of the system. It is clear that the difference increases as
Fitting numerically the straight line shown in Fig. 5, results in a slope of approximately 12 which can be understood easily. The difference in the number of oscillators δNosc between the large cube and the 8 small ones is:
The non-additivity of entropy is generated by the partition of the initial cube into 8 smaller cubes, which creates a larger surface and, therefore, increases the number of surface atoms with fewer oscillators associated than the inner atoms.
2.4. Microcanonical and canonical ensembles
The calculations performed so far have followed the microcanonical ensemble. As stated above, this is the proper method for analyzing systems with a few hundred atoms. If a system of this size is connected to an energy reservoir, their interaction would get the nanocrystal out of equilibrium. Besides, it is usually claimed in any statistical mechanics course that the results obtained using either the canonical or microcanonical ensembles coincide in the thermodynamic limit. The question that arises is, then, whether this also happens when working with very few particles.
In our notation, ZFS
is the partition function of the nanocrystal of side nL
and with
In Eqs. (8) and (9), the superscript “can” indicates that the specific heats are obtained using the canonical ensemble -in contrast to our work so far following the microcanonical ensemble. By comparing (8) and (9), it can be seen that the surface introduces a correction of 1/nL to the specific heat.
Do the specific heat values computed with both ensembles coincide? The answer is in Fig. 5 and it shows a very good agreement between both results, even for the small system (nL = 3). The specific heats graphed have been calculated within the microcanonical (superscript mic) and canonical (superscript can) ensembles, for the FP and FS cases. The energy per atom (M/NA ), which is the natural variable in the microcanonical ensemble, was used as independent variable. For the case of the specific heats (8) and (9), Eq. (6) allows to replace the temperature with the energy per atom. For the case of nL = 10, there is no discernible difference between the results from the canonical and microcanonical ensembles and has, therefore, not been illustrated.
3. Conclusions
As in Ref. 5, the first point to highlight is the strength of the results in statistical mechanics: even with systems formed by 27 elements, the thermostatistical properties are not essentially different from the results found in the thermodynamic limit. In this sense, the usual approximations of statistical physics are reassuring; one may trust them.
Related to this point, the surface is clearly important for the properties of the studied systems. For the system with 27 atoms, when there is no consideration for the surface, the specific heat differs by approximately 1% from the thermodynamic value, while there is no difference for the system with 1000 atoms. However, when the surface is considered, the difference is ∼30% in the first case and ∼10% in the second. This means that the presence of the surface has a greater impact than doing statistical mechanics with few particles.
It is clear for what sizes one has to use the exact microcanonical calculation and when the usual thermodynamic results can be employed. For systems formed by dozens of particles, it is justified to proceed in the meticulous way followed in this article. For systems with several thousand atoms, usual thermodynamic results are good enough.
The entropy of the nanocrystal is non-additive, apparently contrary to a basic principle of thermodynamics. However, this is not so; the non-additivity can be explained by the presence of the surface, and it can be clearly shown that it is recovered in the thermodynamic limit.
Finally, our work confirms that the results obtained in the canonical and microcanonical ensembles are equal for systems with NA > 103. This is relevant because calculations are always easier in the canonical ensemble. Thus, people working with nano-systems can use the calculation techniques of the canonical ensemble even when the problem under analysis is isolated and calls for the microcanonical ensemble.
In a nutshell: one can use statistical mechanics even for systems with tens of particles but should be aware of the effects of the surface.
Financial support from Universidad Nacional de Cuyo through Proyecto SECTyP 06/M072 is acknowledged. One of the authors (ENM) thanks Diego Molina for helping in the last stage of this project.