In nucleation theories, the work of formation of a nucleus is often denoted by W = ΔG. This convention misleads that the nucleation should be considered in the isothermal-isobaric system. However, the pressure in the system with a nucleus is no longer uniform due to Laplace’s equation. Instead, the chemical potential is uniform throughout the system for the critical nucleus. Therefore, one can consider the nucleation in the grand ensemble properly. Accordingly, W is found to be the grand potential difference and the interfacial tension is also turned to be an interfacial excess grand potential. This treatment is not entirely new; however, to explicitly treat in the grand potential formalism is for the first time. We have successfully given an overwhelmingly clear description.
The work of formation of a nucleus is often written as. It leads one to understand the work of formation of the critical nucleus as a difference of the Gibbs energy. The meaning of the form of the work of formation of a critical nucleus (Equation (4) in the text) becomes, however, clear straightforwardly if we deal the system including a critical nucleus as an isothermal-isochoric open system. The treatment as an isothermal-isobaric closed system brings confusions. The concept of the Gibbs dividing surface is more clearly understood in the isothermal-isochoric open system. As will be stated in the text, the treatments of an isothermal-isochoric open system appeared in literatures already. In this paper, we will give a clearer and direct statement in the grand potential formalism for nucleation, aiming at helping researchers who are not specialists in thermodynamics. In other words, by describing with definite terminologies we will put forward understandings—some terminologies will be for the first time used definitely in this paper.
Gibbs established the interfacial thermodynamic formula for the work of formation of a critical nucleus in 1870s [
Throughout this paper we restrict ourselves to the case of spherical interfaces for simplicity and for the sake of avoiding complexity for better understanding. For example, two principal curvatures appear in general; this may bring confusion. Also, for the same sake we limit ourselves to unary cases. Also, for the same sake we omit the structure of both two phases; if at least one of the coexisting phases is crystalline, the interfacial tension becomes, strictly speaking, crystallographic orientation dependent.
One of purposes of the thermodynamics of nucleation is to calculate the reversible work of formation of a critical nucleus of a stable phase in an undercooled parent phase [
where is the difference between the chemical potentials of the nucleating phase (phase) and the parent phase (phase). The direct interpretation of Equation (1) is as follows. Limiting ourselves to the case that the molecular volumes (volumes per molecule) of the and phases are equal1, let us denote the molecular volume. Hereafter, the subscript indicates the molecular quantities. Then, the quantity is defined as the number of molecule consisting the nucleus, which is equal to with being the volume of the nucleus. The first term in Equation (1) is the volume term, which is the reversible work associated with the transformation from the phase to the phase of molecules. The second term in Equation (1) is the surface term, which is the reversible work to form a surface of area. Here, is the radius of the nucleus; the rigorous definition of will be given later. Remembering that the chemical potential is equal to the molecular Gibbs energy, the expression of seems at apparent appropriate. The question arises whether the expression of Equation (1) is only valid for the case that no volume change is associated with the phase transition or not. Exact expression for the reversible work was already given and the approximation which reduces the exact expression to Equation (1) was derived [
To review the Gibbs’ formalism for evaluating W is not only heuristic but also ingredient for understanding the thermodynamic “ensemble” appropriate for the system of nucleation. In other words, due to this one can find why the grand potential formalism is appropriate; that is, constant condition is imposed. It is sufficient to limit ourselves to the unary case; formulation for the multi-component system is seen, for example, in a previous paper [
Consider a spherical nucleus of the phase in an undercooled phase of the chemical potential at the temperature The chemical potential and temperature are regarded as those of the reservoir. Along with the isothermal condition, for the critical nucleus one can regard a cluster of the phase is in equilibrium with the phase with respect to the material transport. One can select as independent variables specifying the total system with being the volume of the total system. The following is the procedure of the calculation of the work of formation of a critical nucleus.
1) The pressure of the phase is determined by the equilibrium equation with respect to the materials transport, i.e.
2) Presuming the interfacial tension as known, the radius is determined by Laplace’s equation,
where is the pressure of the phase corresponding to.
3) The work of formation of the critical nucleus of radius R is calculated by
We note that eliminating using Laplace’s equation (Equation (3)), Equation (4) is rewritten into
We should note that the work of formation of a critical nucleus consists of two terms; as has been mentioned the first term is the volume term and, in tern, the second term is the surface term. The first term in Equation (4) is the work to replace the phase of volume with the phase. The second term, , is understood as the work associated with the formation of area of the surface free energy per unit area. In other words, in writing the work of formation of the critical nucleus we divide the process of nucleus formation into two. One is to form a hypothetical nucleus of the phase possessing the bulk properties throughout the entire volume in the parent phase. The other is regarded to that to form a actual structure of the interface.
For the first one of the two works of formation of a critical nucleus, the mathematical surface of radius is a key concept. This surface is called the Gibbs dividing surface. Owing to introducing the dividing surface one can divide the work of formation of a nucleus into two. The volume term is the work of formation of a hypothetical cluster as illustrated in
curved interface possessing the mechanical tension. Therefore, the interfacial free energy is called the interfacial tension.
Sometimes very unnatural variables are specified [2,11]. That is, the internal energy, the entropy, and the amount of substances are selected as independent variable. The mass as well as the number of molecule can be employed as the amount of substances. Nevertheless, Nishioka [11,13] derived a correct conclusion that is equal to the superficial grand potential through an entangled argument.
As pointed in Section 1.2 the chemical potential throughout the system is uniform. Along with the fact that the system is considered as isothermal, it is appropriate to select the temperature and the chemical potential as independent variables. In this case, because at least one extensive variable is necessary for complete description, the total system volume must be, in general, selected as one of the independent variables. We note that the uniformity of the chemical potential was already pointed out [
As mentioned above the temperature and the chemical potential are uniform throughout the system. One can regard that the system is exposed to the isobaric reserver because if the chemical potential and the temperature are kept constant, the corresponding pressure, which is a function of and, is also constant. In
account the change of the total volume in calculation of the work of formation of a nucleus for the former case. This is somewhat complicated. Hence, it is convenient to treat the system as isothermal-isochoric open one. Of course, two ways of description are both correct. The reversible work calculated as the Gibbs energy difference should coincide to that calculated as the grand potential difference. Indeed, a consideration with confusions led to the correct answer [
At least in Japan, a thermodynamics class does not teach the grand potential systematically. One can, however, obtain isochoric open system by Legendre transformation of the isothermal-isochoric closed system, i.e., the independent variable is transformed from the amount of substances to the chemical potential to obtain this system [
where is the Gibbs energy. To reach to the last expression we have used the definition. One may be familiar with this form in the grand canonical ensemble (ensemble) through the bridging relation in this ensemble [
By virtue of the last expression of Equation (6), we obtain the volume term of the work of formation of a critical nucleus, as the grand potential difference between the system including the hypothetical nucleus and the homogeneous phase, as
where and are the pressures of respective bulk phases; even though there is no bulk part of the phase in reality such as for a small nucleus, the pressure is well defined (through Equation (2)). Due to the positive interfacial tension between the and phases, the pressure of the phase inside the dividing surface is greater than (thermodynamic derivation of this relation will be given in Section 2.2). In this way, we have the first term in Equation (4), which is negative and corresponding to the volume bulk term in Equation (1).
As argued up to now, we know that the work of formation of a critical nucleus is composed of the volume term, which is corresponding to the first term in Equation (1) and given by Equation (7), and the surface term, which is corresponding to the second term in Equation (1). If the equilibrium with respect to the materials transport holds between the parent phase and the nucleus, the pressure inside the nucleus, , is obtained by solving
which corresponds to Equation (2) and consistent to the isothermal open system (ensemble). Because the phase is metastable and the phase is the stable phase; that is,
holds, one can derive. Recalling the GibbsDuhem relation, we draw schematically the chemical potentials as functions of the pressure in
In this way, the negativity of the volume term is understood. The criterion for the dividing surface has not been given yet. The surface term, in general, take a form [7,9,11]
Here, denotes that this coefficient depends on the criterion for the dividing surface. The surface of tension is defined by. Only for this choice of, the coefficient coincide with the interfacial tension. In other words, the surface term consist of, in general, the interfacial area dependent term and the curvature dependent term. The surface of tension is defined for which the curvature dependent term vanishes. We note that takes the minimum for the surface of tension [
In this way, we have obtained Equation (4) for the work of formation of a critical nucleus. We give a note here. The work for the formation of the critical nucleus takes, however, the same value if the physical condition is unchanged; that is, it is not dependent on the criterion of the dividing surface. Therefrom, one can derive the relation between the general and the interfacial tension. This was done by Kondo [
Noting and, let us solve the equation that the derivative with respect to of Equation (4) vanishes. By a simple calculation we have Laplace’s equation (Equation (3)). This is a mechanical balance equation. Namely, in a case that two phases are coexisting via an interface of a curvature radius with an interrfacial tension, the force acting from the inside of the interface due to the pressure balances with the composed force of the force due to the outside pressure and that due to the interfacial tension (corresponding to). The quantity defined as the interfacial free energy per unit area of the interface is, if one chooses the surface of tension as the dividing surface, coincides with the mechanical interfacial tension. Readers can readily confirm the coincidence between the unit of the energy per area and the tension.
Now, let us derive the form of the first term in Equation (1), following Nishioka and Kusaka [
which is nothing other than the Gibbs-Duhem relation for the isothermal case. We consider a case that an incompressible phase nucleus is nucleated in the phase. Let us integrate Equation (11) for the phase for from to.
Eliminating in Equation (4) using the equation derived by dividing Equation (12), we have an equation corresponding to Equation (1):
where
To reach to the last expression, Equation (8) has been used. One can integrate Equation (11) for the phase to obtain the form of Equation (1) in a case that the phase is incompressible. This is, however, not the present concern. It should be noted that for a case that no volume change is associated with the phase transition, a form far form Equation (1) is obtained [
In this section, we derive the Gibbs adsorption isotherm
where represents the chemical potential of the materials reservoir, which is equal to, and is the superficial number density per unit area of the interface, sometimes referred to as the excess number density or the interfacial adsorption quantity. A rigorous definition of will be given later.
We define the system as a spherical cone as illustrated in
For the hypothetical system, because of the bulk properties, the following fundamental equations (Gibbs relations) hold for two parts of the system:
Here, according to a convention is used to represent the internal energy. This equation is rewritten in terms of the grand potentials as
Those equations hold for both systems with the solid angle and the entire sphere. In those expressions
and we should note that and are independent variables.
Let us denote quantities for the entire spherical system by symbols with a superscript and those for the system with the solid angle by symbols without a superscript. For a while, let us consider again a general dividing surface. Denoting the contribution due to the nucleus by…, the fundamental equation
holds. Here, and because the is an independent variable,. Let us rewrite Equation (20) using and
. Because (from), we have
Here, we express the contribution of the nucleus by introducing the coefficient defined by
as previously done [2,3,5,9,11,13]. In those previous papers, the expression in the square brackets was given.
Differentiating and using Equation (20), we have
By comparing Equations (21) and (23), we obtain
In previous papers [2,3,5,9,11,13], the last expression was given, despite that the mid expression is conceptually meaningful. This equation is the equation obtained from the relation on the basis of the fact that when the solid angle is multiplied by, the grand potential
is transformed as
(Euler’s theorem). We note that Nishioka [
In Equation (21), existence of is due to the nucleus. Therefore, one can write
(pay attention on the independent variables). The first two terms are of the hypothetical system defined in Section 3.1. The last two terms are for forming interfacial structure after the formation of the hypothetical system. As mentioned above, we note that a term depending on the derivative of the curvature radius, , appears. This term, also as mentioned above, vanishes if the surface of tension is taken as the dividing surface.
Let us go forward the argument by taking the surface of tension as the dividing surface. Using the equation obtained by putting in Equation (25), we rewrite Equation (21) into
The fundamental equation for the hypothetical system is just the addition of both of Equation (17):
Subtracting Equation (27) from Equation (26), we have the fundamental equation for the superficial grand potential:
where and are, respectively, the superficial entropy and the superficial number of molecules. In this equation has been eliminated because the state of the interface is independent of the selection of; in other words, has been fixed at the position.
Euler’s relation obtained from the fact that is transformed as when is multiplied by as is
To derive this equation, one can use the same method to derive Equation (24). From Equation (29), the interfacial tension is revealed to be the superficial grand potential per unit area of the interface. Introducing the superficial quantities per unit area of the interface, , and, we have
The last expressions in Equations (29) and (30) have already be given in previous papers [3,4,5,9,11-13,15,20, 21,31]. In those papers, except for [12,20,21,31]—Rusanov et al. [
A general way to obtain the Gibbs-Duhem relation is to take differential of Euler’s relation and subtract the fundamental equation. For the interface, the same procedure is possible; we can have the Gibbs-Duhem relation for the interface
by taking differential of Equation (29) and subtract the first equation of Equation (28) and dividing by. We can, also, obtain Equation (31) by direct differentiation of Equation (30) and using the fundamental equation for. From Equation (31) we have Equation (15) or . This is the Gibbs adsorption isotherm.
We have given a grand potential formalism for the interfacial thermodynamics. It is revealed that the work of formation of a critical nucleus is equal to the grand potential difference. This makes a point of view clearer overwhelmingly than regarding the work of formation of the nucleus as the Gibbs energy difference. Also, the interfacial tension is revealed to be defined as the superficial grand potential per unit area of the interface. Although equivalent form was given previously [3-5,9,11, 13], this paper has explicitly closed up the grand potential property for the first time.
This paper is base on a lecture [in Japanese] at the 35th research meeting on the crystal growth (Toronkai) held by Japanese Association for Crystal Growth on Sept. 7-9, 2011 at Shimotsuma, Japan.