Study of Structural and Melting Properties of Gold Nanorods

Rida Essajai

#### Abstract

MD simulations combined with the embedded-atom method have been applied to study the structural and melting properties of gold nanorods (AuNRs) of different sizes. The simulation results for the actual structure of AuNRs obtained after energy minimization processes revealed that the AuNRs with largest cohesive energies tend to be structurally more stable than those with smallest ones. Then, it was found that each actual structure of AuNR is classified as an irregular structure composed of a crystalline gold core covered by an amorphous gold shell. In addition, the results showed that the melting of the AuNR surface is an inhomogeneous, gradually occurring process. Besides, it was established that the premelting ratio is inversely correlated with the AuNR size, indicating that the premelting phenomenon is more pronounced in large NP sizes than in small ones.

Keywords: gold nanorods, size effect, thermodynamic properties, melting temperature, premelting temperature, MD simulations

#### 1. Introduction

Gold-based nanorods have attracted and continue to attract the attention of a vast amount of scientists from all over the world thanks to their fundamental and pragmatic significance. Due to their exciting properties that are found to be absent in corresponding bulk counterpart, the gold nanorods (AuNRs) are useful nanoobjects in many applications [1–14]. Especially, considerable attention has been paid over the past decades to AuNRs, because of their great importance in fabricating the new generation of nanomolecular and molecular electronics as mentioned in Ref. [15]. Besides, AuNRs can be added to base fluids in order to better improve the effective thermal properties of the nanofluids [16]. It is worth noting here that the technological properties of molecular electronic devices and the thermal property enhancement of nanofluids strongly depend on both the state of AuNR surface and temperature at which their structures change. Nonetheless, the temperature-dependent surface structure of single metal nanoparticles was ignored by several authors, where they did not take into account the surface premelting stage during the analysis of the melting process in the nanoparticles [17–19]. Therefore, it is necessary to understand the AuNR melting properties before fabricating new nanodevices based on them or suspending them in an energetic system (fluid). However, to our knowledge, the investigation on the melting properties of these nanoobjects is still lacking in the theoretical and experimental reports and, thus, needs intensive studies. To this end, the MD simulations combined with the

embedded-atom method (EAM) will be applied to study the melting behaviors of AuNRs. For this purpose, a series of results will be presented in the present work that is organized as follows: In Section 2, there is a short description of the computational method and the calculation procedure. Section 3 deals with the simulation results of AuNRs. Finally, the main conclusions arising from this work are summarized in Section 4.

where the factor <sup>1</sup>

Study of Structural and Melting Properties of Gold Nanorods

DOI: http://dx.doi.org/10.5772/intechopen.84682

temperature from 0 to 2300 K using NVT.

couples, rij ¼ rj ! � ri � � !�

density is given by

constant.

average.

59

3. Results and discussions

<sup>2</sup> is included in order to avoid double computation of the atom

ρ<sup>i</sup> rij � � (2)

� is the scalar distance, and ρ<sup>i</sup> is the host electron density

induced by all encirclement atoms j at the position of atom i. The host electron

<sup>ρ</sup><sup>i</sup> <sup>¼</sup> <sup>X</sup> i6¼ j

In order to investigate the melting properties of AuNRs of different sizes, a series of MD simulations were carried out. First, all models considered in this work are equilibrated over 500,000 time steps in the NVT ensemble at zero temperature using the velocity Verlet method with a fixed time step of 1 fs and neglecting the periodic boundary conditions (PBCs) in three dimensions. Then after energy minimization process, the actual model systems at 0 K were obtained and characterized using two structural techniques, namely, the coordination number (CN) analysis and the common neighbor analysis (CNA) method. The final point of the simulations is that all the simulated equilibrium samples are heated by increasing the

In order to predict the melting temperatures of AuNRs, the variation of thermodynamic properties (including the potential energy per atom (U) and the heat capacity (C)) with respect to temperature was analyzed. The heat capacity is

BT<sup>2</sup> <sup>U</sup><sup>2</sup> � � � h i <sup>U</sup> <sup>2</sup> � � <sup>þ</sup>

where T is the temperature, N is the number of atoms, and KB is the Boltzmann

Besides the thermodynamic properties, the melting temperature of NPs is also explored by analyzing the variation in other physical quantities in terms of the temperature [29]. In this context, the Lindemann index is regarded as one of the most important of these quantities. Additionally, this index is a tool to examine the NP in the layer form in order to determine from where the NP melting begins. The

j <k

where L ið Þ is the ith layer, rjk is the distance between atoms k and j, NL ið Þ is number of atoms of the ith layer, and the bracket hi represents the ensemble

Although the initial structures of gold nanorods can be generated by using LAMMPS, these original structures are not the nanorods with equilibrium states (i.e., metastables). Therefore, it is necessarily to get the actual model systems (i.e.,

Using the MD simulation, the temporal evolution of each nanostructure passes through several states until it reaches a more stable one, corresponding to minimum

the real calculation models) before studying the structural analysis.

3

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi rjk<sup>2</sup> � � � rjk � �<sup>2</sup> <sup>q</sup>

<sup>2</sup> (3)

rjk � � (4)

obtained from the ensemble averages of U and its square, as [28]

<sup>¼</sup> <sup>1</sup> NK<sup>2</sup>

C KB

Lindemann index of each layer, δL ið Þ, is represented as [30]

<sup>δ</sup>L ið Þ <sup>¼</sup> <sup>2</sup>

NL ið Þ NL ið Þ � <sup>1</sup> � �<sup>X</sup>

#### 2. Computational methodology

All MD simulations were performed with the LAMMPS package [20], and the atomistic visualizations were carried out by Ovito [21]. First of all, eight different model systems of rod-shaped nanoparticles were constructed using LAMMPS, where the ratio between their diameters to their lengths is more than 1:10. (Example of structure can be seen in Figure 1.) The detailed parameters of all rod-shaped nanoparticles are listed in Table 1. Furthermore, embedded-atom method (EAM) potential was adopted to describe the interatomic interaction that has been developed from the density functional theory (DFT) by Daw and Baskes [22]. Moreover, this potential has been successfully used by our group in a variety of computational disciplines [23–26]. The EAM potential used here was parameterized by Grochola et al. [27] to describe the gold–gold interactions by adjusting its parameters with experiments and ab initio calculations.

The total energy given by the EAM is written as a sum of an embedding function Fi <sup>ρ</sup><sup>i</sup> ð Þ and pair potential <sup>ϕ</sup><sup>i</sup>, <sup>j</sup> rij � � as

$$E = \sum\_{i} F\_{i}(\overline{\rho\_{i}}) + \frac{1}{2} \sum\_{i,j} \phi\_{i,j}(r\_{\vec{\eta}}) \tag{1}$$

Figure 1.

Structure of rod-shaped Au nanoparticle generated by lamps (D = 6.3 nm, Lz = 61.86 nm, D:lz = 0.1, and 111147 atoms). The atomic structures are visualized by Ovito.


Table 1.

Diameter D, length lz, and total number atoms (N) for the AuNRs.

Study of Structural and Melting Properties of Gold Nanorods DOI: http://dx.doi.org/10.5772/intechopen.84682

embedded-atom method (EAM) will be applied to study the melting behaviors of AuNRs. For this purpose, a series of results will be presented in the present work that is organized as follows: In Section 2, there is a short description of the computational method and the calculation procedure. Section 3 deals with the simulation results of AuNRs. Finally, the main conclusions arising from this work are summa-

All MD simulations were performed with the LAMMPS package [20], and the atomistic visualizations were carried out by Ovito [21]. First of all, eight different model systems of rod-shaped nanoparticles were constructed using LAMMPS, where the ratio between their diameters to their lengths is more than 1:10. (Example of structure can be seen in Figure 1.) The detailed parameters of all rod-shaped nanoparticles are listed in Table 1. Furthermore, embedded-atom method (EAM) potential was adopted to describe the interatomic interaction that has been developed from the density functional theory (DFT) by Daw and Baskes [22]. Moreover, this potential has been successfully used by our group in a variety of computational disciplines [23–26]. The EAM potential used here was parameterized by Grochola et al. [27] to describe the gold–gold interactions by adjusting its parameters with

The total energy given by the EAM is written as a sum of an embedding function

ϕ<sup>i</sup>, <sup>j</sup> rij

� � (1)

Fi <sup>ρ</sup><sup>i</sup> ð Þþ <sup>1</sup> 2 X i, j

Structure of rod-shaped Au nanoparticle generated by lamps (D = 6.3 nm, Lz = 61.86 nm, D:lz = 0.1, and

Nanoparticles NP1 NP2 NP3 NP4 NP5 NP6 NP7 NP8 D (nm) 2.03 3.84 4.64 5.07 6.3 7.32 8.42 9.44 Lz (nm) 19.43 30.43 41.25 50.55 61.86 72.43 83.34 93.04 N (atoms) 8907 13,907 42,375 74,164 111,147 149,133 194,139 243,880

� � as

<sup>E</sup> <sup>¼</sup> <sup>X</sup> i

rized in Section 4.

Nanorods and Nanocomposites

2. Computational methodology

experiments and ab initio calculations.

111147 atoms). The atomic structures are visualized by Ovito.

Diameter D, length lz, and total number atoms (N) for the AuNRs.

Fi ρ<sup>i</sup> ð Þ and pair potential ϕ<sup>i</sup>, <sup>j</sup> rij

Figure 1.

Table 1.

58

where the factor <sup>1</sup> <sup>2</sup> is included in order to avoid double computation of the atom couples, rij ¼ rj ! � ri � � !� � is the scalar distance, and ρ<sup>i</sup> is the host electron density induced by all encirclement atoms j at the position of atom i. The host electron density is given by

$$\overline{\rho}\_{i} = \sum\_{i \neq j} \rho\_{i}(r\_{\vec{\eta}}) \tag{2}$$

In order to investigate the melting properties of AuNRs of different sizes, a series of MD simulations were carried out. First, all models considered in this work are equilibrated over 500,000 time steps in the NVT ensemble at zero temperature using the velocity Verlet method with a fixed time step of 1 fs and neglecting the periodic boundary conditions (PBCs) in three dimensions. Then after energy minimization process, the actual model systems at 0 K were obtained and characterized using two structural techniques, namely, the coordination number (CN) analysis and the common neighbor analysis (CNA) method. The final point of the simulations is that all the simulated equilibrium samples are heated by increasing the temperature from 0 to 2300 K using NVT.

In order to predict the melting temperatures of AuNRs, the variation of thermodynamic properties (including the potential energy per atom (U) and the heat capacity (C)) with respect to temperature was analyzed. The heat capacity is obtained from the ensemble averages of U and its square, as [28]

$$\frac{\mathbf{C}}{K\_B} = \frac{1}{N\mathbf{K}\_B^2 T^2} \left( \left< \mathbf{U}^2 \right> - \left< \mathbf{U} \right>^2 \right) + \frac{\mathbf{3}}{2} \tag{3}$$

where T is the temperature, N is the number of atoms, and KB is the Boltzmann constant.

Besides the thermodynamic properties, the melting temperature of NPs is also explored by analyzing the variation in other physical quantities in terms of the temperature [29]. In this context, the Lindemann index is regarded as one of the most important of these quantities. Additionally, this index is a tool to examine the NP in the layer form in order to determine from where the NP melting begins. The Lindemann index of each layer, δL ið Þ, is represented as [30]

$$\delta\_{L(i)} = \frac{2}{N\_{L(i)} \left( N\_{L(i)} - 1 \right)} \sum\_{j - \left< r\_{jk} \right>^2}}{\left< r\_{jk} \right>} \tag{4}$$

where L ið Þ is the ith layer, rjk is the distance between atoms k and j, NL ið Þ is number of atoms of the ith layer, and the bracket hi represents the ensemble average.

#### 3. Results and discussions

Although the initial structures of gold nanorods can be generated by using LAMMPS, these original structures are not the nanorods with equilibrium states (i.e., metastables). Therefore, it is necessarily to get the actual model systems (i.e., the real calculation models) before studying the structural analysis.

Using the MD simulation, the temporal evolution of each nanostructure passes through several states until it reaches a more stable one, corresponding to minimum total energy, but remains metastable, resulting from its surface atoms always wanting to aggregate to reduce surface energy as was indicated in Ref. [5]. After this pretreatment stage, the structure stability of AuNRs was investigated by the sizedependent cohesive energy (the absolute value of total energy) of the AuNR, and the results were plotted in Figure 2.

As expected, it is found that the calculated cohesive energies of the particles are lower than that of the bulk Au which was obtained previously [23, 31]. It was also established that the cohesive energy of AuNR is very sensitive to the size. The cohesive energy value increases with increasing AuNR size. This behavior seems to indicate that the AuNRs of higher energies tend to be structurally more stable and it can be explained in terms of reducing the influence of surface free energy (i.e., energy of all of the dangling bonds in the surface) resulting from the decrease in the dangling bonds to total bonds ratio when increasing the size.

Within the framework of the local structural characterization of AuNRs, the common neighbor analysis (CNA) technique was adopted. For more details on this method, see Refs. [23, 32, 33]. Once AuNRs in their equilibrium states were obtained (see Figure 3a which displays the section of [100]-oriented equilibrium AuNR at different sizes), the two categories of Au atoms in each NR were calculated by the CNA technique; then their percentage were illustrated in Figure 3b. It was found from Figure 3b that for the smallest AuNR (D = 1.622 and Lz = 15.26 nm), there is a significant percentage of unidentified structures (47.5%), while proportion of fcc structures is identified by 52.5%. It was moreover observed that for all AuNR sizes, the fraction of Au atoms having fcc structures is higher than those possessing unidentified structures. Also, it was found that when the AuNR size increases, the fraction of both the fcc structure and unidentified one are, respectively, an increasing and decreasing function until they reach percentages of 82.5 and 17.5% at size correspondent to D = 8.704 nm and Lz = 86.12 nm (NP8).

layers: the first contains atoms located on the subsurface layer (CN = 11 and 10), the second consists of atoms located on the free-surface layer, corresponding to both low-index facets [i.e., {110} and {100} facets with the CN = 7 and CN = 8 and their perimeter (edge and corner with the CN < 7)]. Apparently, the shell part of the NP is inherently different from the core one, where its atoms possess structure intermediate between the bulk liquid structure and the bulk solid one. It is worth mentioning that the classification of AuNRs in core-shell constructions obtained here supports the previous experimental and simulated studies carried out on single metal nanoparticles [23–26, 34, 35]. Furthermore, it was seen from Figure 4 that for all simulated AuNRs, the fraction of liquid-like atoms located in the shell part is smaller than that of solid-like atoms located in the core part. It was also found that for all actual structures of AuNRs, the liquid-like atoms located in the area of lowindex facets (7 and 8 coordinated atoms) exhibit the highest fraction in the shell part. In addition, when the AuNR size increases, the fraction of atoms' core part increases from 51 to 83%, while the total fraction of shell atoms decreases from 49 to 17%. These findings are consistent with those previously reported by the CNA analysis. The results of the current study indicates that the ratio of low-coordinated

(a) Sections of [100]-oriented AuNRs after energy minimization. (b) The atomic structures are visualized by Ovito. The percentage of each category of atoms in AuNRs is counted by using the CNA method, yellow for fcc

Figure 3.

61

structures and purple for unidentified structures.

Study of Structural and Melting Properties of Gold Nanorods

DOI: http://dx.doi.org/10.5772/intechopen.84682

atoms to full-coordinated atoms is inversely correlated with the AuNR size.

Figure 5 shows the typical trend of potential energy (U) and heat capacity (C) of different AuNR sizes as a function of temperature T during the heating process from 0 to 2300 k. For each AuNR, it was established that the U curve presents a jump at the same temperature that the C curve has a sharp peak corresponding to the complete melting point (Tpm) of the AuNR (observed at 1270.17 K for NP1, 1306.43 for NP2, 1330.34 for NP3, 1349.18 for NP4, 1361.85 for NP5, 1367.07 for NP6, 1372.52 for NP7, and 1377.65 K for NP8). As predicted, the results observed here reveal that the complete melting points of AuNRs are lower than that of bulk Au obtained previously [23, 36]. It was also found that, as the AuNR size increases, the simulated value of Tpm increases. This behavior is in agreement with previous results [37, 38]. Notably, the increase of Tpm could be attributed to the decrease in the ratio of low-coordinated atoms to full-coordinated ones (i.e., the total number

In order to describe in more details the local atomic-level structures in different size of AuNRs, the coordination number (CN) about each individual atom was calculated; then the results are presented in Figure 4. It was shown that each actual structure of the AuNR obtained after energy minimization process reveals the coreshell structure composed by a crystalline gold core (full-coordinated atoms with the CN = 12) covered by an amorphous gold shell. This last part is divided into two

Figure 2. The cohesive energy per atom as a function of the AuNR size (N) after energy minimization.

Study of Structural and Melting Properties of Gold Nanorods DOI: http://dx.doi.org/10.5772/intechopen.84682

#### Figure 3.

total energy, but remains metastable, resulting from its surface atoms always wanting to aggregate to reduce surface energy as was indicated in Ref. [5]. After this pretreatment stage, the structure stability of AuNRs was investigated by the sizedependent cohesive energy (the absolute value of total energy) of the AuNR, and

As expected, it is found that the calculated cohesive energies of the particles are lower than that of the bulk Au which was obtained previously [23, 31]. It was also established that the cohesive energy of AuNR is very sensitive to the size. The cohesive energy value increases with increasing AuNR size. This behavior seems to indicate that the AuNRs of higher energies tend to be structurally more stable and it can be explained in terms of reducing the influence of surface free energy (i.e., energy of all of the dangling bonds in the surface) resulting from the decrease in the

Within the framework of the local structural characterization of AuNRs, the common neighbor analysis (CNA) technique was adopted. For more details on this method, see Refs. [23, 32, 33]. Once AuNRs in their equilibrium states were obtained (see Figure 3a which displays the section of [100]-oriented equilibrium AuNR at different sizes), the two categories of Au atoms in each NR were calculated by the CNA technique; then their percentage were illustrated in Figure 3b. It was found from Figure 3b that for the smallest AuNR (D = 1.622 and Lz = 15.26 nm), there is a significant percentage of unidentified structures (47.5%), while proportion of fcc structures is identified by 52.5%. It was moreover observed that for all AuNR sizes, the fraction of Au atoms having fcc structures is higher than those possessing unidentified structures. Also, it was found that when the AuNR size increases, the fraction of both the fcc structure and unidentified one are, respectively, an increasing and decreasing function until they reach percentages of 82.5 and 17.5% at size correspondent to D = 8.704 nm and Lz = 86.12 nm (NP8).

In order to describe in more details the local atomic-level structures in different

size of AuNRs, the coordination number (CN) about each individual atom was calculated; then the results are presented in Figure 4. It was shown that each actual structure of the AuNR obtained after energy minimization process reveals the coreshell structure composed by a crystalline gold core (full-coordinated atoms with the CN = 12) covered by an amorphous gold shell. This last part is divided into two

The cohesive energy per atom as a function of the AuNR size (N) after energy minimization.

the results were plotted in Figure 2.

Nanorods and Nanocomposites

Figure 2.

60

dangling bonds to total bonds ratio when increasing the size.

(a) Sections of [100]-oriented AuNRs after energy minimization. (b) The atomic structures are visualized by Ovito. The percentage of each category of atoms in AuNRs is counted by using the CNA method, yellow for fcc structures and purple for unidentified structures.

layers: the first contains atoms located on the subsurface layer (CN = 11 and 10), the second consists of atoms located on the free-surface layer, corresponding to both low-index facets [i.e., {110} and {100} facets with the CN = 7 and CN = 8 and their perimeter (edge and corner with the CN < 7)]. Apparently, the shell part of the NP is inherently different from the core one, where its atoms possess structure intermediate between the bulk liquid structure and the bulk solid one. It is worth mentioning that the classification of AuNRs in core-shell constructions obtained here supports the previous experimental and simulated studies carried out on single metal nanoparticles [23–26, 34, 35]. Furthermore, it was seen from Figure 4 that for all simulated AuNRs, the fraction of liquid-like atoms located in the shell part is smaller than that of solid-like atoms located in the core part. It was also found that for all actual structures of AuNRs, the liquid-like atoms located in the area of lowindex facets (7 and 8 coordinated atoms) exhibit the highest fraction in the shell part. In addition, when the AuNR size increases, the fraction of atoms' core part increases from 51 to 83%, while the total fraction of shell atoms decreases from 49 to 17%. These findings are consistent with those previously reported by the CNA analysis. The results of the current study indicates that the ratio of low-coordinated atoms to full-coordinated atoms is inversely correlated with the AuNR size.

Figure 5 shows the typical trend of potential energy (U) and heat capacity (C) of different AuNR sizes as a function of temperature T during the heating process from 0 to 2300 k. For each AuNR, it was established that the U curve presents a jump at the same temperature that the C curve has a sharp peak corresponding to the complete melting point (Tpm) of the AuNR (observed at 1270.17 K for NP1, 1306.43 for NP2, 1330.34 for NP3, 1349.18 for NP4, 1361.85 for NP5, 1367.07 for NP6, 1372.52 for NP7, and 1377.65 K for NP8). As predicted, the results observed here reveal that the complete melting points of AuNRs are lower than that of bulk Au obtained previously [23, 36]. It was also found that, as the AuNR size increases, the simulated value of Tpm increases. This behavior is in agreement with previous results [37, 38]. Notably, the increase of Tpm could be attributed to the decrease in the ratio of low-coordinated atoms to full-coordinated ones (i.e., the total number

#### Figure 4.

The distribution of the coordination numbers of AuNR for different sizes, orange for fcc structures, royal for second layer (CN = 11, 10), olive for Perimeter atoms (CN < 7), and dark cyan for facet atoms (CN = 7, 8).

of bonds in the NP increases); hence, AuNRs with large sizes require more temperature to completely melt than those with small sizes.

lattice positions. Then, δ<sup>L</sup> abruptly jumped, and its factor is many times greater than before, from which the premelting phenomenon occurs at the temperature before that of complete melting; hence some liquid-like atoms start to appear. It is evident that after the sharp rise, δ<sup>L</sup> increased with fluctuations at several temperatures. This latter can be ascribed to the nonhomogeneous nature of the melting of free-surface atoms originating from the coexistence between two types of atoms: those who are moving and those who are remaining strongly constrained in their solid ordering

Variation of the heat capacity and potential energy per atom of AuNR for different sizes, during the heating

Study of Structural and Melting Properties of Gold Nanorods

DOI: http://dx.doi.org/10.5772/intechopen.84682

Figure 5.

process.

63

Now we will turn to Figure 6, which shows the temperature-dependent Lindemann index value of free-surface atoms (δL) for different AuNR sizes. It is easy to see that in the beginning of the heating process, the δ<sup>L</sup> increases slowly and linearly, thus suggesting that all surface atoms merely vibrate around their original

Figure 5.

Variation of the heat capacity and potential energy per atom of AuNR for different sizes, during the heating process.

lattice positions. Then, δ<sup>L</sup> abruptly jumped, and its factor is many times greater than before, from which the premelting phenomenon occurs at the temperature before that of complete melting; hence some liquid-like atoms start to appear. It is evident that after the sharp rise, δ<sup>L</sup> increased with fluctuations at several temperatures. This latter can be ascribed to the nonhomogeneous nature of the melting of free-surface atoms originating from the coexistence between two types of atoms: those who are moving and those who are remaining strongly constrained in their solid ordering

of bonds in the NP increases); hence, AuNRs with large sizes require more temper-

The distribution of the coordination numbers of AuNR for different sizes, orange for fcc structures, royal for second layer (CN = 11, 10), olive for Perimeter atoms (CN < 7), and dark cyan for facet atoms (CN = 7, 8).

Now we will turn to Figure 6, which shows the temperature-dependent Lindemann index value of free-surface atoms (δL) for different AuNR sizes. It is easy to see that in the beginning of the heating process, the δ<sup>L</sup> increases slowly and linearly, thus suggesting that all surface atoms merely vibrate around their original

ature to completely melt than those with small sizes.

Figure 4.

Nanorods and Nanocomposites

62

Figure 6. Temperature-dependent Lindemann index of first layer for AuNRs.

against premelting phenomenon. Besides, it was observed that the fluctuations presenting in δ<sup>L</sup> curve are attenuated progressively with the increase in temperature until completely damped. This may be due to the increased amount of liquid-like atoms until they form a continuous liquid layer. It can be clearly shown from Figure 6 that the higher temperature region where δ<sup>L</sup> increases monotonically is in the liquid-like range of the first layer for AuNR different sizes. This liquid region of δ<sup>L</sup> curve may correspond with the fact that all surface atoms moved cooperatively.

that AuNRs of smaller cohesive energies tend to be structurally less stable than those larger ones. Then, the local structure characterization of the actual structure of AuNRs carried out by using the CN and CNA techniques showed that each AuNR is classified as an irregular structure which is composed of a crystalline gold core that is covered by an amorphous gold shell. Further, it was shown that melting temperatures of the actual structure of AuNRs are below that in bulk state. In addition, it was found that the surface melting in the actual structure of AuNRs is an

inhomogeneous, gradually occurring process. Besides, the result of the sizedependent premelting ratio proved that the premelting phenomenon is less

Group of Semiconductors and Environmental Sensor Technologies- Energy Research Center, Faculty of Science, Mohammed V University, Rabat, Morocco

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

\*Address all correspondence to: rida.essajai@gmail.com

provided the original work is properly cited.

pronounced in small AuNRs than in large ones.

Premelting ratio as a function of the particle size (N).

Study of Structural and Melting Properties of Gold Nanorods

DOI: http://dx.doi.org/10.5772/intechopen.84682

Author details

Rida Essajai

65

Figure 7.

Once again according to Figure 6, the onset of premelting temperature (Tsm) of each actual structure of the AuNR is determined from the temperature where δ<sup>L</sup> curve shows a sharp increase. The results showed that Tsm decreased from 1198.82 K for NP8 to 1142.87 K for NP1, and considering the Tpm already determined from U and C curves (Figure 5), its value is above the Tsm one. It was also found that the Tsmvalue decreased by 177.51 for NP8, 174.17 for NP7, 170.57 for NP6, 165.40 for NP5, 159.25 for NP4, 151.98 for NP3, 139.18 for NP2, and 127.3 for NP1, compared with the Tpm one. On the basis of the above, it was concluded that all AuNRs in the present study are passing through the surface premelting stage during heating process.

Our attention will now be drawn to identifying how the ratio between Tsm and Tpm of AuNR also called the premelting ratio varies in terms of size. One of the main objectives of this investigation is to provide insights into the AuNR size-dependent premelting phenomenon. Hence, this study aims to fill this void in the existing literature. Figure 7 is given to show the behavior of the premelting ratio of AuNR as a function of the size (N). It was found that the premelting ratio value decreases when the NP size increases, thus indicating that the premelting phenomenon is less pronounced in small NP sizes than in large ones.

#### 4. Conclusions

In summary, the structural and melting properties of AuNRs have been investigated by means of MD simulations combined with the embedded-atom method (EAM) potential. The simulation results for the actual structure of AuNRs revealed Study of Structural and Melting Properties of Gold Nanorods DOI: http://dx.doi.org/10.5772/intechopen.84682

Figure 7. Premelting ratio as a function of the particle size (N).

against premelting phenomenon. Besides, it was observed that the fluctuations presenting in δ<sup>L</sup> curve are attenuated progressively with the increase in temperature until completely damped. This may be due to the increased amount of liquid-like atoms until they form a continuous liquid layer. It can be clearly shown from Figure 6 that the higher temperature region where δ<sup>L</sup> increases monotonically is in the liquid-like range of the first layer for AuNR different sizes. This liquid region of δ<sup>L</sup> curve may correspond with the fact that all surface atoms moved cooperatively. Once again according to Figure 6, the onset of premelting temperature (Tsm) of each actual structure of the AuNR is determined from the temperature where δ<sup>L</sup> curve shows a sharp increase. The results showed that Tsm decreased from 1198.82 K for NP8 to 1142.87 K for NP1, and considering the Tpm already determined from U and C curves (Figure 5), its value is above the Tsm one. It was also found that the Tsmvalue decreased by 177.51 for NP8, 174.17 for NP7, 170.57 for NP6, 165.40 for NP5, 159.25 for NP4, 151.98 for NP3, 139.18 for NP2, and 127.3 for NP1, compared with the Tpm one. On the basis of the above, it was concluded that all AuNRs in the present study are passing through the surface premelting stage during heating

Temperature-dependent Lindemann index of first layer for AuNRs.

Our attention will now be drawn to identifying how the ratio between Tsm and Tpm of AuNR also called the premelting ratio varies in terms of size. One of the main objectives of this investigation is to provide insights into the AuNR size-dependent premelting phenomenon. Hence, this study aims to fill this void in the existing literature. Figure 7 is given to show the behavior of the premelting ratio of AuNR as a function of the size (N). It was found that the premelting ratio value decreases when the NP size increases, thus indicating that the premelting phenomenon is less

In summary, the structural and melting properties of AuNRs have been investigated by means of MD simulations combined with the embedded-atom method (EAM) potential. The simulation results for the actual structure of AuNRs revealed

pronounced in small NP sizes than in large ones.

process.

Figure 6.

Nanorods and Nanocomposites

4. Conclusions

64

that AuNRs of smaller cohesive energies tend to be structurally less stable than those larger ones. Then, the local structure characterization of the actual structure of AuNRs carried out by using the CN and CNA techniques showed that each AuNR is classified as an irregular structure which is composed of a crystalline gold core that is covered by an amorphous gold shell. Further, it was shown that melting temperatures of the actual structure of AuNRs are below that in bulk state. In addition, it was found that the surface melting in the actual structure of AuNRs is an inhomogeneous, gradually occurring process. Besides, the result of the sizedependent premelting ratio proved that the premelting phenomenon is less pronounced in small AuNRs than in large ones.
