Melting point depression

Melting point depression is a term referring to the phenomenon of reduction of the melting point of a material with reduction of its size. This phenomenon is very prominent in nanoscale materials which melt at temperatures hundreds of degrees lower than bulk materials.

Introduction

The melting temperature of a bulk material is not dependent on its size. However as the dimensions of a material decrease towards the atomic scale, the melting temperature scales with the material dimensions. The decrease in melting temperature can be on the order of tens to hundreds of degrees for metals with nanometer dimensions. [1-4]

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.

This article focuses on nanoparticles because researchers have compiled a large amount of size-dependent melting data for near spherical nanoparticles. [1-4] Nanoparticles are easiest to study due their ease of fabrication and simplified conditions for theoretical modeling. The melting temperature of a nanoparticle decreases sharply as the particle reaches critical diameter, usually < 50 nm for common engineering metals. [1,2,5] Figure 1 shows the shape of a typical melting curve for a metal nanoparticle as a function of its diameter.

Melting point depression is a very important issue for applications involving nanoparticles, as it decreases the functional range of the solid phase. Nanoparticles are currently used or proposed for prominent roles in catalyst, sensor, medicinal, optical, magnetic, thermal, electronic, and alternative energy applications. [6] . Nanoparticles must be in the solid state to function at elevated temperatures in several of these applications.

History

Takagi first observed melting point depression of several types of metal nanoparticles in 1954. [4] A variable intensity electron beam from a transmission electron microscope melted metal nanoparticles in early experiments. Diffraction patterns changed from characteristic crystalline patterns to liquid patterns as the small particles melted, allowing Takagi to estimate the melting temperature from the electron beam energy.

Measurement Techniques

Two techniques allow measurement of the melting point of nanoparticle. As described above, the electron beam of the TEM can be used to melt nanoparticles. [4,7] The melting temperature is estimated from the beam intensity, while changes in the diffraction conditions to indicate phase transition from solid to liquid. This method allows direct viewing of nanoparticles as they melt, making it possible to test and characterize samples with a wider distribution of particle sizes. The TEM limits the pressure range at which melting point depression can be tested.

More recently, researchers developed nanocalorimeters that directly measure the enthalpy and melting temperature of nanoparticles. [5] Nanocalorimeters provide the same data as bulk calorimeters, however additional calculations must account for the presence of the substrate supporting the particles. A narrow size distribution of nanoparticles is required since the procedure does not allow users to view the sample during the melting process. There is no way to characterize the exact size of melted particles during experiment

Physics of Melting Point Depression

Nanoparticles have a much greater surface to volume ratio than bulk materials. The increased surface to volume ratio means surface atoms have a much greater effect on chemical and physical properties of a nanoparticle. 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 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 according to Equation 1. [8]

$E=E\_B(1-frac\{d\}\{D\})$

Where: D=nanoparticle size
::d=atomic size
::E_b=cohesive energy of bulk

As Equation 1 shows, the effective cohesive energy of a nanoparticle approaches that of the bulk material as the material extends beyond atomic size range (D>>d).

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 potential. The Lennard-Jones pair-potential shown in Figure 2 models the cohesive energy between atoms as a function of separation distance.

The cohesive energy of an atom is directly related to the thermal energy required to free the atom from the solid. According to Lindemann’s criterion, the melting temperature of a material is proportional to it’s cohesive energy, a_v (T_M=Ca_v). [9] Since atoms near the surface have fewer bonds and reduced cohesive energy, they require less energy to free from the solid phase. Melting point depression of high surface to volume ratio materials results from this effect. For the same reason, surfaces of bulk materials can melt at lower temperatures than the bulk material. [10]

The theoretical size-dependent melting point of a material can be calculated through classical thermodynamic analysis. The result is the Gibbs-Thomson relation shown in Equation 2. [2]

$T\_M(d)=T\_\{MB\}(1-frac\{4sigma,\_\{sl\{H\_f\; ho,\_sd\})$

Where: T_MB=Bulk Melting temperature
::σ_sl=solid liquid interface energy
::H_f=Bulk heat of fusion
::ρ_s=density of solid
::d=particle diameter

A normalized Gibbs-Thomson Equation for gold nanoparticles is plotted in Figure 1, and the shape of the curve is in general agreement with those obtained through experiment. [11]

emiconductor/Covalent Nanoparticles

Equation 2 gives the general relation between the melting point of a metal nanoparticle and its diameter. However, recent work indicates the melting point of semiconductor and covalently bonded nanoparticles may have a different dependence on particle size. [12] The covalent character of the bonds changes the melting physics of these materials. Researchers have demonstrated that Equation 3 more accurately models melting point depression in covalently bonded materials. [12]

$T\_M(d)=T\_\{MB\}(1-(frac\{c\}\{d\})^2)$
Where: T_MB=Bulk Melting temperature
::c=materials constant
::d=particle diameter

Equation 3 indicates that melting point depression is less pronounced in covalent nanoparticles due to the quadratic nature of particle size dependence in the melting Equation.

Proposed Mechanisms for Nanoparticle Melting

The specific melting process for nanoparticles is currently unknown. The scientific community currently accepts several mechanisms as possible models of nanoparticle melting. [12] Each of the corresponding models effectively match experimental data for melting of nanoparticles. Three of the four models detailed below derive the melting temperature in a similar form using different approaches based on classical thermodynamics.

Liquid Drop Model

The liquid drop model (LDM) assumes that an entire nanoparticle transitions from solid to liquid at a single temperature. [9] This feature distinguishes the model, as the other models predict melting of the nanoparticle surface prior to the bulk atoms. If the LDM is true, a solid nanoparticle should function over a greater temperature range than other models predict. The LDM assumes that the surface atoms of a nanoparticle dominate the properties of all atoms in the particle. The cohesive energy of the particle is identical for all atoms in the nanoparticle.

The liquid drop model represents the binding energy of a nanoparticles as function of the free energies of the volume and surface. [9] Equation 4 gives the normalized, size dependent melting temperature of a material according the liquid-drop model.

$T\_M(d)=frac\{4T\_\{MB\{H\_fd\}(sigma,\_\{sv\}-sigma,\_\{lv\}(frac\{\; ho,\_s\}\{\; ho,\_l\})^\{2/3\})$

Where: σ_sv=solid-vapor interface energy
::σ_lv=liquid-vapor interface energy
::H_f=Bulk heat of fusion
::ρ_s=density of solid
::ρ_l=density of liquid
::d=diameter of nanoparticle

Liquid Shell Nucleation Model

The liquid shell nucleation model (LSN) predicts that a surface layer of atoms melts prior to the bulk of the particle. [13] The melting temperature of a nanoparticle is a function of its radius of curvature according to the LSN. Large nanoparticles melt at greater temperatures as a result of their larger radius of curvature.

The model calculates melting conditions as a function of two competing order parameters using Landau potentials. One order parameter represents a solid nanoparticle, while the other represents the liquid phase. Each of the order parameters is a function of particle radius.

The parabolic Landau potentials for the liquid and solid phases are calculated at a given temperature, with the lesser Landau potential assumed to be the equilibrium state at any point in the particle. In the temperature range of surface melting, the results show that the Landau curve of the ordered state is favored near the center of the particle while the Landau curve of the disordered state is smaller near the surface of the particle.

The Landau curves intersect at a specific radius from the center of the particle. The distinct intersection of the potentials means the LSN predicts a sharp, unmoving interface between the solid and liquid phases at a given temperature. The exact thickness of the liquid layer at a given temperature is the equilibrium point between the competing Landau potentials.

Equation 5 gives the condition at which an entire nanoparticle melts according to the LSN model. [14]

$T\_M(d)=frac\{4T\_\{MB\{H\_fd\}(frac\{sigma,\_\{sv\{1-frac\{d\_0\}\{d-sigma,\_\{lv\}(1-frac\{\; ho,\_s\}\{\; ho,\_l\}))$
Where: d₀=atomic diameter

Liquid Nucleation and Growth Model

The liquid nucleation and growth model (LNG) treats nanoparticle melting as a surface initiated process. [15] The surface melts initially, and the liquid-solid interface quickly advances through the entire nanoparticle. The LNG defines melting conditions through the Gibbs-Duhem relations, yielding a melting temperature function dependent on the interfacial energies between the solid and liquid phases, volumes and surface areas of each phase, and size of the nanoparticle. The model calculations show that the liquid phase forms at lower temperatures for smaller nanoparticles. Once the liquid phase forms, the free energy conditions quickly change and favor melting. Equation 6 gives the melting conditions for a spherical nanoparticle according to the LNG model. [14]

$T\_M(d)=frac\{2T\_\{MB\{H\_fd\}(sigma,\_\{sl\}-sigma,\_\{lv\}3(sigma,\_\{sv\}-sigma,\_\{lv\}frac\{\; ho,\_s\}\{\; ho,\_l\}))$

Bond Order Length Strength (BOLS) Model

The bond-order-length-strength (BOLS) model employs an atomistic approach to explain melting point depression. [14] The model focuses on the cohesive energy of individual atoms rather than a classical thermodynamic approach. The 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.

The BOLS mechanism states that if one bond breaks the remaining neighbouring ones become shorter and stronger. The cohesive energy, or the sum of bond energy, of the less coordinated atoms determines the thermal stability including meting, evaporating and other phase transition. The lowered CN changes the equilibrium bond length between atoms near the surface of the nanoparticle. The bonds relax towards equilibrium lengths, increasing the cohesive energy per bond between atoms, independet of the exact form of the specific interatomi cpotential. However, 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.

Using a core-shell configuration, the melting point depression of nanoparticles is dominated by the outermost two atomic layers yet atoms in the core interior remain their bulk nature.

The BOLS model and the core-shell structure have been applied to other size dependency of nanostrtcures such as the mechanical strength, chemical and thermal stability, lattice dynamics (optical and accoustic phonons), Photon emmision and absorption, electronic colevel shift and work fnction moduluation, magnetism at various temperatures, and dielectrics due to electron polarization etc. Reproduction of experimental observations in the above mentioned size dependency has been realized. Quantitative information such as the energy level of an isolated atom and the vibration frequency of individual dimer has been obtained by matching the BOLS predictions to the measured size dependency [15] .

Particle Shape

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. [9] These shape changes affect the surface to volume ratio, which affects the cohesive energy and thermal properties of a nanostructure. Equation 7 gives a general shape corrected formula for the theoretical melting point of a nanoparticle based on its size and shape. [9]

$T\_M(d)=T\_\{MB\}(1-frac\{c\}\{zd\})$
Where: c=materials constant
::z=shape parameter of particle

The shape parameter is 1 for sphere and 3/2 for a very long wire, indicating that melting point depression is suppressed in nanowires compared to nanoparticles. Past experimental data show that nanoscale tin platelets melt within a narrow range of 10 ^oC of the bulk melting temperature. [7] The melting point depression of these platelets was suppressed compared to spherical tin nanoparticles. [5]

ubstrate

Several nanoparticle melting simulations theorize that the supporting substrate affects the extent of melting point depression of a nanoparticle. [1,16] These models account for energetic interactions between the substrate materials. A free nanoparticle, as many theoretical models assume, has a different melting temperature (usually lower) than a supported particle due to the absence cohesive energy between the nanoparticle and substrate. However, measurement of the properties of a freestanding nanoparticle remains impossible, so the extent of the interactions cannot be verified through experiment. Ultimately, substrates currently support nanoparticles for all nanoparticle applications, so substrate/nanoparticle interactions are always present and must impact melting point depression.

References

[1] A. Jiang, N. Awasthi, A. N. Kolmogorov, W. Setyawan, A. Borjesson, K. Bolton, A. R. Harutyunyan and S. Curtarolo, Phys. Rev. B, 75, 205426 (2007).

[2] J. Sun and S. L, Simon, Thermochimica Acta, 463, 32 (2007).

[3] A. F. Lopeandia and J. Rodriguez-Viejo, Thermochimica Acta, 461, 82 (2007).

[4] M. Takagi, J Phys. Soc. Jpn., 9, 359 (1954).

[5] S. L. Lai, J. Y. Guo, V. Petrova, G. Rammath and L. H. Allen, Phys. Rev. Lett., 77, 99 (1996).

[6] G. G. Wildgoose, C. E. Banks and R. G. Compton, Small, 2,182 (2006).

[7] G. L. Allen, R. A. Bayles, W. W. Giles and W. A. Jesser, Thin Solid Films, 144, 297 (1986).

[8] W. H. Qi and M. P. Wang, J. Mat. Sci. Lett., 21, 1743 (2002).

[9] K. K. Nanda, S. N. Sahu and S. N. Behera, Phys. Rev. A, 66, 013208 (2002).

[10] J. W. M. Frenken and J. F. van der Veen, Phys. Rev. Lett. 54, 134 (1985).

[11] P. H. Buffat and J. P. Borrel, Phys. Rev. A, 13, 2287 (1976).

[12] H. H. Farrell and C. D. Van Sicien, J. Vac. Sci. and Technol. B, 25, 1441 (2007).

[13] H. Sakai, Surf. Sci., 351, 285 (1996).

[14] C. Q. Sun, Y. Wang, B. K. Tay, S. Li, H. Huang and Y. B. Zhang, J. Phys. Chem. B, 106, 10701 (2002).

[15] C. Q. Sun, Size dependence of nanostructures: impact or bond order deficiency, Prog Solid State Chem 35, 1-159 (2007).

[16] P. R. Couchman and W. A. Jesser, Nature, 269, 481 (1977).

[17] S. P. Huang, D. S. Mainardi and P. B. Balbuena, Surf. Sci., 545, 163 (2003).