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 𝑇=300 K, and 𝑛𝑖 remains small enough so as to have a negligible effect on the operation of most of the devices until π‘‡β‰ˆ1000 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 astotalsurfaceenergy=freesurfaceenergy+interfacialenergy+substrateenergy.(1)

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Ξ“=2πœ‹π‘…π»π›Ύπ‘ ξ€·π›Ύ+πœ‹π»(2π‘…βˆ’π»)π‘ π‘βˆ’π›Ύπ‘ξ€Έ+Γ𝑏,(2) 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πœ•Ξ“+πœ•π‘…πœ•Ξ“πœ•π»=0.(3) To find the relationship between 𝑅 and 𝐻, substituting (2) in (3), we get𝐻=Δ𝛾𝑠𝑏𝛾𝑠𝑅,(4) where Δ𝛾𝑠𝑏 is a parameter called wetting or spreading parameter which is given by Δ𝛾𝑠𝑏=π›Ύπ‘βˆ’π›Ύπ‘ βˆ’π›Ύπ‘ π‘.(5)

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, Δ𝛾𝑠𝑏>0, then there is a possibility of wetting the substrate, and if it is negative, that is, Δ𝛾𝑠𝑏<0, 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, 𝐻=βˆ’Ξ”π›Ύπ‘ π‘π›Ύπ‘ π‘….(6) Substituting the value of 𝐻 (height of the spherical nanoparticle on the substrate) in (3), we getΞ“βˆ—=2πœ‹π›Ύπ‘ ξ‚΅π‘Žπ‘…π‘ ξ‚Άπ‘Ž2+Γ𝑏.(7)

π‘…βˆ—π‘  is the corresponding radius of curvature of supported spherical nanoparticle, and it is given byπ‘…βˆ—π‘ =12ξ‚„1/3ξ‚Έ2𝛾𝑠Δ𝛾𝑠𝑏2/3π‘Žξ€Ί3+Δ𝛾𝑠𝑏/𝛾𝑠1/3,(8) 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π‘…βˆ—π‘ =πΏπΊβŽ›βŽœβŽœβŽœβŽπ›Ώξ€·1βˆ’4𝛾2π‘ πœ‹2ξ€Έ3ξ‚ƒπ‘Ÿξ€·4𝛾2π‘ πœ‹2ξ€Έ+ξ€·π›Ύπ‘βˆ’π›Ύπ‘ βˆ’π›Ύπ‘ π‘ξ€Έ2ξ‚„3⎞⎟⎟⎟⎠,(9) 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π‘…βˆ—π‘™=πΏπΊβŽ›βŽœβŽœβŽœβŽπ›Ώξ€·1βˆ’4𝛾2π‘™πœ‹2ξ€Έ3ξ‚ƒπ‘Ÿξ€·4𝛾2π‘™πœ‹2ξ€Έ+ξ€·π›Ύπ‘βˆ’π›Ύπ‘™βˆ’π›Ύπ‘™π‘ξ€Έ2ξ‚„3⎞⎟⎟⎟⎠,(10) 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]:𝐹𝑠=43πœ‹π‘Ž3πœŒπ‘“π‘ +4πœ‹π‘Ž3πœŒπ›Ύπ‘ π‘…βˆ—π‘ ,𝐹𝑙=43πœ‹π‘Ž3πœŒπ‘“π‘™+4πœ‹π‘Ž3πœŒπ›Ύπ‘™π‘…βˆ—π‘™,(11) 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Δ𝐹=πΉπ‘ βˆ’πΉπ‘™,4Δ𝐹=3πœ‹π‘Ž3πœŒξ‚Έξ€·π‘“π‘ βˆ’π‘“π‘™ξ€Έξ‚΅π›Ύ+3π‘ π‘…βˆ—π‘ βˆ’π›Ύπ‘™π‘…βˆ—π‘™.ξ‚Άξ‚Ή(12) It is already known thatπ‘“π‘™βˆ’π‘“π‘ ξ‚΅π‘‡=𝜌𝐿1βˆ’π‘‡π‘ξ‚Ά,(13) 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, Δ𝐹=0. Substituting the above two equations, we getπ‘‡π‘š=𝑇𝑐21βˆ’ξ‚΅π›ΎπœŒπΏπ‘ π‘…βˆ—π‘ βˆ’π›Ύπ‘™π‘…βˆ—π‘™.ξ‚Άξ‚Ά(14) If π‘…βˆ—π‘ =π‘…βˆ—π‘™=π‘…βˆ—, that is, if no surface melting takes place then, π‘‡π‘š=𝑇freeπ‘š(π‘…βˆ—). 𝑇freeπ‘š 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:π‘‡π‘š=𝑇freeπ‘šξ‚΅21βˆ’ξ‚΅π›ΎπœŒπΏπ‘ π‘…βˆ—π‘ βˆ’π›Ύπ‘™π‘…βˆ—π‘™.ξ‚Άξ‚Ά(15)

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𝑇GN2=𝑇1βˆ’πœŒπ‘ πΏπΊξƒ©π›Ύπ‘ π‘…π‘ βˆ’π›Ύπ‘™π‘…π‘™ξ‚΅πœŒπ‘ πœŒπ‘™ξ‚Ά2/3ξƒͺξƒ­,(16) where 𝑇GN 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𝛿=π»βˆ’β„Ž=π‘…βˆ’π‘Ÿ.(17) The total energy of the surface-melted particle is then a function of 𝛿 and is given as𝐹(𝛿)=𝑉𝑠𝑓(𝛿)π‘ βˆ’π‘“π‘™ξ€Έ+𝑉𝑙𝑓𝑙+Ξ“(𝛿),(18) 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):𝑐Helixlength=𝑐×,helixlength𝐢=helixlengthΓ—helixlength𝑐=2πœ‹π‘…βˆ—,π‘…βˆ—=(2πœ‹π‘…βˆ—)2+𝐻2Γ—βˆš(2πœ‹π‘…)2+𝐻2,2πœ‹Γ—2πœ‹π‘…2πœ‹π‘…βˆ—=(2πœ‹π‘…βˆ—)2+𝐻2Γ—ξƒŽ(2πœ‹π‘…)2+𝐻2,𝑅2πœ‹π‘…βˆ—=(πœ‹/2)2𝑅2+(𝐻/4)2(πœ‹/2)2𝑅,(A.1) 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π‘…βˆ—π‘ =𝑅4𝛾2π‘ πœ‹2ξ€Έ+ξ€·π›Ύπ‘βˆ’π›Ύπ‘ βˆ’π›Ύπ‘ π‘ξ€Έ2ξ‚„4𝛾2π‘ πœ‹2.(A.2)

Similarly, the corresponding radius of curvature of the nanoparticle with a liquid film surrounding the solid is calculated to beπ‘…βˆ—π‘™=𝑅4𝛾2π‘™πœ‹2ξ€Έ+ξ€·π›Ύπ‘βˆ’π›Ύπ‘ βˆ’π›Ύπ‘™π‘ξ€Έ2ξ‚„4𝛾2π‘™πœ‹2.(A.3) where  𝛾𝑙  is the free surface energy of the liquid-vapor interface and is the interfacial energy between the liquid and substrate.