Abstract
The study of variation of the size and shape effect on the melting property of gallium nitride nanoparticles with their spherical and cylindrical geometrical feature is theoretically explored. A numerical thermodynamical model has been devoted for the study. A comparative investigation is made between the two shapes, at the range of ~3βnm dia. The cylindrical GaN nanoparticles, whose melting point has been reported to decrease with decreasing particle radius, become larger than spherical-shaped nanoparticles. The melting temperature obtained in the present study is in line with the function of radius of curvature.
1. Introduction
It has been well established both experimentally and theoretically that the melting temperature () of nanoparticles depends on the particle size [1β11]. Pawlow in 1909 developed a thermodynamical model that predicts a point depression of nanoparticles with the particle size. An attempt to confirm this experimentally has been made first by Pawlow [12] in 1910. Subsequently, other researchers have investigated the variation of melting temperature with particle size [13β17].
A numerical thermodynamical model is implemented for our exploration of the temperature distribution with respect to other size of the spherical and cylindrical nanoparticles of gallium nitride (GaN). In the last decade, GaN has been investigated intensively, both experimentally and theoretically [18β22]. It is a wide bandgap semiconductor that exists in both wurtzite and zinc blende crystal structure. All devices, currently constructed for applications with the use of wurtzite structure. Because of the wide bandgap, the intrinsic carrier concentration of GaN is essentially Zero at βK, and remains small enough so as to have a negligible effect on the operation of most of the devices until βK. This property makes this wide bandgap material suitable for use in high-temperature environments. GaN is a very hard, mechanically stable wide bandgap semiconductor material with high heat capacity and thermal conductivity [18]. GaN can be doped with silicon (Si) or with oxygen to n-type and with magnesium (Mg) to p-type [19]; however, the Si and Mg atoms change the way GaN crystals grow, introducing tensile stresses and making them brittle [20]. GaN compounds also tend to have a high spatial defect frequency, on the order of a hundred million to ten billion defects per square centimeter. The crystalline quality GaN led to the discovery of p-type GaN [19], p-n junction blue/UV-LEDs [19], and room-temperature-stimulated emission. The very high breakdown voltage [21] high electron mobility, and saturation velocity of GaN has also made it an ideal candidate for high-power and high-temperature microwave applications, as evidenced by its high Johnsonβs figure of merit. Moreover, GaN-based MOSFET and MESFET transistors also offer many advantages in high-power electronics, especially in automotive and electric car applications [22]. Nanotubes of GaN are proposed for applications in nanoscale electronics, optoelectronics, and biochemical-sensing applications. They are also useful in military electronics such as active electronically scanned array radars. A GaN nanoparticle in a BN matrix by nitridation with urea has a lot of applications. Bionanotechnology is the use of biomolecules for applications in nanotechnology, including use of viruses.
The interest has been heightened recently at the nanoscale, because nanostructures are pervasive in nature and in modern industry. We depict that the variation of the properties of GaN material with its geometrical feature size has a long history because of its importance in many fields. In physics and chemistry, the effect of particle size on melting has been discussed since 1900s, and this effect is not restricted to any particular material; rather, it is observed in a variety of materials from metals and alloys, and the typical size range over which the melting temperature undergoes a large change in the range 5β100βnm [23β31]. Many phenomena in solid-state physics and materials science also exhibit size dependence. For example, the elastic constants of Ag and Pb nanowires of diameter 30βnm are nearly twice those of the bulk metals [32]. Such increase in stiffness cannot be explained by structural modifications of the materials as the nanoscale reduction in the size of the solids also results in a change of their failure mode.
Thus, when the size of brittle calcium carbonate particles is reduced to a critical value of 850βnm [33, 34], the particles behave as if they were ductile. Size dependence of the melting temperature at nanoscale has enormous implications in the production of nanocrystals and in the thermal stability of quantum dots. A large data has accumulated on this size dependence, and a number of theoretical models have been proposed to explain it. Thermodynamically, the melting temperature of nanoparticles has been described by three models: (1) the homogeneous melting and growth, (2) the liquid shell nucleation, and (3) the liquid nucleation and growth. All the three models predict a size-dependent melting temperature.
The melting temperature of the nanoparticles will be different in different shapes when considering the radius of curvature of nanoparticles; especially the particle size is mentioned with respect to its radius of curvature. Therefore, an indispensable course of action to take the particle size into consideration is when one develops the models for the melting temperature of nanoparticles. In the present work, the radius of curvature is considered to account for the particle shape difference and develop the model for the cylindrical nanoparticles. According to the relation between the melting temperature and radius of curvature of the nanoparticles, an expression (15) for the size- and shape-dependent melting temperature of nanoparticles is developed. The theoretical prediction of this expression for the melting temperature of GaN nanoparticles is compared between spherical and cylindrical shapes.
2. Model of Calculation
2.1. Surface Energy Calculation
It is known that the total surface energy involved in any heterogeneous nucleation is given as
It has been assumed that the geometry of the heterogeneous nucleation is that of the cylindrical-shaped particle. So the equation can be written as where are the radius and the height of the cylindrical nanoparticle on the substrate, is the surface energy of the solid vapor interface, is the surface energy of the bare substrate, is the interfacial energy between the solid and the substrate and is the total energy of the substrate.
Now, (2) can be minimized by using the condition To find the relationship between and , substituting (2) in (3), we get where is a parameter called wetting or spreading parameter which is given by
The magnitude of the wetting or spreading parameter () is used to determine whether the surface melting takes place or not in the given substrate. If the magnitude of the spreading parameter is positive, that is, , then there is a possibility of wetting the substrate, and if it is negative, that is, , then the surface melting takes place.
In the case of nanofilms, the free surface energy of the solid is always higher than the energy density of the bare substrate, thus,ββββis always less than zero and it is negative. From (5), it is seen clearly that if is negative, then is also negative, and this is not possible. Thus, to keep positive, a negative sign is included in (4), that is, Substituting the value of (height of the spherical nanoparticle on the substrate) in (3), we get
is the corresponding radius of curvature of supported spherical nanoparticle, and it is given by where is the equilibrium surface energy, and is the corresponding radius of curvature of the supported gallium nitride (GaN) solid spherical nanoparticles [35] and it is given by where is latent heat of melting material of GaN, is fitting parameter, and it is a measure of GaN solid shell.
Similarly, the corresponding radius of curvature of the GaN nanoparticle with liquid film surrounding the solid [35] is given as whereββββis the free surface energy of the liquid-vapors interface, and is the interfacial energy between the liquid and substrate. The basic formula for the radius of curvature of cylindrical nanoparticle which is used in the above expression has been attached in the Appendix [35].
2.2. Melting Point Calculation
The expression for the melting of nanoparticle can be derived using free energies density of solid and liquid. The free energy density expression for the solid () and the corresponding liquid ()ββare given as follows [36]: where is the radius of the melting of nanoparticle, is the surface energy of solid, and is the surface energy of liquid. Now the difference in free energies () is calculated to be It is already known that where is the density of the materials, is the latent heat of the material, is the bulk melting temperature, and is the melting point of the nanomaterial as a function of size. At melting point the difference in free energies is zero, that is, . Substituting the above two equations, we get If , that is, if no surface melting takes place then, (). is the melting point of nanofilms with free interfacial energy between liquid and substrate. Thus, the expression for melting point for nanofilms () as a function of radius of the particle is given as follows:
Using the expression above, the variation of the melting point with respect to the GaN particle can be plotted.
Similarly, the expression for the melting point GaN nanofilms as a function of the radius of the particle [35] is given as where is the temperature of nanoparticle of GaN, ββis the Bulk melting temperature of GaN, Density of GaN solid phase, is the density of GaN liquid phase, is the latent heat of melting GaN material, is the surface energy of solid GaN, is the surface energy of liquid GaN, is the radius of curvature of supported GaN nanoparticle, and is the radius of curvature of nanoparticle with liquid film.
2.3. Surface Melting Calculations
The schematic diagram for the surface melting is as shown in Figure 1. A quasiliquid film of radius is formed on the core solid particle of radius . The difference between and gives the thickness of the liquid film.
At the onset of surface melting, we assume that the geometry is close to that of the solid particle in a quasiliquid film of radius of curvature formed on the core solid particle of radius of curvature . The difference between and gives the thickness of the liquid film. The solid particle is in its equilibrium geometry and it is initially wet by a molten layer with uniform thickness. The radius of curvature of solid particle isββββand height isββ [35].
Liquid film thickness is given as The total energy of the surface-melted particle is then a function of and is given as where is the volume of the solid, is the volume of the liquid, is the thickness-dependent surface energy.
3. Results and Discussion
The crystal shape and the supporting substrate influence the size dependence of melting point, where the cylindrical shape is considered. The mentioned gallium nitride (GaN) nanoparticles are extremely small in 1β100βnm range. There is a considerable difference in the calculated values of the melting point as a function of size of the particle. It is known that the radius of curvature varies according to the shapes of the nanoparticle which decide the size of the particle.
The main difference between (8) and other expression (9) for the size-dependent melting temperature is that the radius of curvature of spherical and cylindrical shape is considered in (8) and (9), respectively, and derived using a separate helical method which is given in the Appendix. Based on expression (8), corresponding melting point for GaN is obtained in expression (15).
Using the expressions (8) and (15), the size of the supported spherical nanoparticle is calculated and also the variation of the melting point, which is listed in Table 1.
Similarly, using the expressions (9) and (16), the size of the supported cylindrical nanoparticle is calculated and also the variation of the melting point, which is listed in Table 2. The thickness of the quasiliquid film is plotted as a function of the size of the supported nanoparticles.
3.1. Melting Point versus Radius of the Particle
Using the data for the gallium nitride (GaN) nanofilms in the expression that is derived for melting point, the variation of the melting point with respect to the size of the supported nanoparticles can be tabulated. Using Tables 1 and 2, the melting point can be plotted against the radius of the supported spherical and cylindrical nanoparticles, respectively.
From Figures 2 and 3, it can be clearly seen that with decrease in size of the supported nanoparticles, the melting point reduces. It is also seen that there is a variation in the melting point for free nanoparticles and the nanofilms. Thus, the effect of substrate is determined.
An important part of our present work is a comparative study of GaN film between spherical and cylindrical nanoparticles has been done, and the theoretical calculated values are tabulated in Table 3, comparative Figure 4 plotted by the tabulated values shows the variation of melting point of GaN nanofilms against the size of the spherical and cylindrical nanoparticles.
In this, it is clearly shown that the melting temperature decreases when the particle size reduces little more in the case of cylindrical than the spherical nanoparticle. At nanoscales, particles exhibit many thermophysical features distinct from those found at microscales. As the size decreases beyond a critical value, due to the increase in the surface-to-volume ratio, the melting temperature derivates from the bulk values and becomes a size-dependent property. This phenomenon has been studied experimentally by means of transmission electron diffraction by Wronski for nanosized tin [37]. The melting point agrees reasonably well with the predictions based on classical theories, which show a nonlinear relationship with the reciprocal of the particle size. The melting temperature of gallium nitride (GaN) becomes size dependent for particles smaller than 10βnm, reducing from 1618βK range for a 5βnm particle.
Pawlow [1] improved the Gibbs-Thompson model by considering the equilibrium of a liquid spherical drop with both a solid spherical particle of the same material and its vapor. This model leads to the following Well-Lenon Equation (15). The suggestion that the melting point of a particle should depend on its size is implicit in the work Thomson [38].
It is well known that the melting temperature of Au (1064βK) decreases when particle dimensions are reduced to the nanoscale. Therefore, at ~3βnm diameter, Au particle can melt at temperature ~500βK [23, 39, 40]. Similarly, the melting temperature of B4C (2450βK) lowered to ~764βK range with spherical-shaped and ~495βK ranges with cylindrical shaped nanoparticles at ~3βnm dia [35]. In the present work, the melting temperature of GaN (2770βK) decreases significantly when the particle dimensions are reduced to the nanoscale and we got that an ~3βnm diameter of GaN spherical shaped particle melted at temperature ~1747βK range, but the same particle with an ~3βnm diameter melted at the temperature ~1553βK range, since the particle shape is considered in cylindrical shape.
3.2. Melting Point versus Inverse of the Radius of the Particle
In Figures 5 and 6 there are two regions namely, logarithmic increases at low values of radius and an exponential increase at the higher values of radius. A linear figure (Figures 5 and 6) between melting point and the inverse of the radius of the particle is drawn.
3.3. Thickness of the Liquid Film versus Radius of the Particle
From Figure 2 it can be seen clearly that with increased in the radius of the supported nanoparticles, it increases up to certain limit after which it is saturated. The melting temperature is a function of the particle radius evaluated from the calculations and the analytical formulation reported in Buffat and Borel [23]. In that, the calculation of latent heat is a function of the particle radius, which is valid for small particle with spherical shape, and it is the radius of curvature varies according to the shape of the particle. The obtained expression (9) because is for the corresponding radius of curvature of the GaN solid cylindrical nanoparticle. Similarly, the expression (16) obtained for the melting point of GaN nanofilms is a function of radius of the particle.
Melting-point depression is most evident in nanowires, nanotubes, and nanoparticles, which all melt at lower temperatures than bulk amounts of the same material. Changes in melting point occur because nanoscale materials have a much larger surface-to-volume ratio than bulk materials, drastically altering their thermodynamic and thermal properties. The decrease in melting temperature can be on the order of tens to hundreds of degrees for metals with nanometer dimensions [41β44]. Surface atoms bind in the solid phase with less cohesive energy because they have fewer neighboring atoms in close proximity compared to atoms in the bulk of the solid. Each chemical bond an atom shares with a neighboring atom provides cohesive energy, so atoms with fewer bonds and neighboring atoms have lower cohesive energy. The average cohesive energy per atom of a nanoparticle has been theoretically calculated as a function of particle size [45].
Atoms located at or near the surface of the nanoparticle have reduced cohesive energy due to a reduced number of cohesive bonds. An atom experiences an attractive force with all nearby atoms according to the Lennard-Jones potentials [46]. The cohesive energy of an atom is directly related to the thermal energy required to free the atom from the solid. Since atoms near the surface have fewer bonds and reduced cohesive energy, they require less energy to be free from the solid phase. Melting point depression of high surface-to-volume ratio materials results from this effect [47]. The liquid drop model (LDM) represents the binding energy of a nanoparticle as a function of the free energies of the volume and surface [46]. The liquid shell nucleation model (LSN) predicts that a surface layer of atoms melts prior to the bulk of the particle [48]. The melting temperature of a nanoparticle is a function of its radius of curvature according to the LSN. The bond-order-length-strength (BOLS) model calculates the melting temperature for individual atoms from the sum of their cohesive bonds. As a result, the BOLS predicts the surface layers of a nanoparticle melt at lower temperatures than the bulk of the nanoparticle [49]. The lowered coordination number changes the equilibrium bond length between atoms near the surface of the nanoparticle. The integrated cohesive energy for surface atoms is much lower than bulk atoms due to the reduced coordination number and overall decrease in cohesive energy.
Nanoparticle shape impacts the melting point of a nanoparticle. Facets, edges, and deviations from a perfect sphere all change the magnitude of melting point depression [46]. These shape changes affect the surface-to-volume ratio, which affects the cohesive energy and thermal properties of a nanostructure. The cohesive energy of nanocrystals, which can be determined by experiment [50] and computed by different theoretical models such as the SE model [51], the BE model [52], the Lennard-Jones potential model [53], Jiangβs model [54], the liquid-drop model [55β57], and the bond-OLS model [58, 59], is regarded as directly related to the nature of the thermal stability of nanocrystals [60].
In observing the evolution of numerical thermodynamical model yields the calculated data, which is compared to the melting point of spherical nanoparticle, the cylindrical nanoparticle has low melting point, since the surface area of cylindrical nanoparticle (169.56βnm2 at ~3βnm radius range) is larger than the spherical (113.04βnm2 at ~3βnm radius range). The number of atoms present on the surface will have only fewer neighboring atoms. Hence, on the larger surface, the more number of atoms will have fewer neighboring atoms in close proximity compared to atoms in the bulk of the nanomaterial. Each chemical bond of an atom shared with a neighboring atom provides cohesive energy, so atom with fewer bonds and neighboring atoms have lower cohesive energy; they require less energy to be free from the solid phase. Hence, melting point depression of high surface (cylindrical shape) is reduced more than the lower surface (spherical shape) nanoparticles. It should be mentioned that the radius of curvature with respect to the shape of the nanoparticle and as the present work describes the melting temperature effect on the shape difference between the spherical nanoparticles and cylindrical nanoparticles. This same property of GaN nanowires with triangular cross-section has been studied using molecular dynamics simulation, in that it was reported that the melting temperature of the GaN nanowires increases with increasing cross-sectional area to a saturation value [61]. The βapproximatelyβ is stressed here due to the fact that some different shape of the nanoparticle may have different radius of curvature which decides the particle size. However, present calculation shows that the present radius of curvature expression is enough for predicting the shape-dependent size effect on the temperature of nanoparticle. This method could be used to study the nanoparticles, and it can be experimentally determined by measuring the particle shape with radius of curvature, which is subjected to future experiments.
It is reported that the atomic radius of metallic nanoparticles contracts with decreasing their particle size, which means that the atomic radius will change a little if the particle size is small. The thickness of the liquid film against the radius of the supported nanoparticle is also done for GaN which has been given in the tabulation (Tables 4 and 5) and in Figures 7 and 8.
The regulation of nanotechnology has made a growing debate related to the human health and safety risks associated with nanotechnology [62]. An important property of the material, the melting temperature of GaN, was predicted by top-down approach through classical thermodynamics [63]. Hence, the phenomenon of size and shape effect on melting temperature of nanoparticles is utilized in nanotechnology field. The expression for the size-dependent melting temperature of nanoparticles in present work is derived from their size-dependent radius of curvature according to the shape of the nanoparticles. It should be mentioned that, according to the shape of the nanoparticles, the radius of curvature varies. And it should be noted that the particle size-dependent on its radius of curvature. Similar melting phenomena have been reported for Au nanoparticles [19], and the melting temperature of prism-shaped nanoparticles [26] and especially for B4C nanoparticles [45] had been analyzed. It has been argued that substrate-particle interaction plays a significant role in the melting behavior of nanoparticles. It has also been shown that the extrapolated bulk melting temperature is lower than the experimental value. It has also been noted by different researchers that the bulk melting temperature cannot be extrapolated from the nanoscale results [64]. In our present work, we showed that the melting temperature of GaN decreased when the size of the nanoparticle reduced. Moreover this reduction is little more in cylindrical nanoparticle when compared with spherical nanoparticle.
4. Conclusion
The present exploration showed, the size- and shape-dependent temperature of GaN nanoparticles, using the numerical thermodynamical model, where spherical and cylindrical shapes of the nanoparticle are considered with respect to its radius of curvature. It is shown that the present results of the melting temperature of GaN nanoparticles consistently decrease with the size of the nanoparticle. The cylindrical GaN nanoparticles, whose melting point has been decreased with decreasing particle radius, become larger than spherical-shaped nanoparticle. Further, it is found that the particle shape can affect the melting temperature of nanoparticles, and this effect on the melting temperature becomes larger with decreasing particle size. Because melting temperature is a very important parameter, the model developed in the present investigation may have potential applications in the temperature-related phenomena of nanoparticles.
Appendix
The following formula for radius of curvature of a cylinder (by Helix method) (see Figure 9): where : Height required for helix to complete one revolution about the cylinder, : Circumference of cylinder, : Radius of cylinder, : Circumference of baluster circle, : Radius of curvature of cylinder.
By using the above expression, the radius of curvature of solid cylindrical nanoparticle is obtained as
Similarly, the corresponding radius of curvature of the nanoparticle with a liquid film surrounding the solid is calculated to be whereββββis the free surface energy of the liquid-vapor interface and is the interfacial energy between the liquid and substrate.