Au-Cu alloys are functional materials with nonlinear optical applications. However, the optical properties of such alloys are difficult to predict due to the random mixing of materials. In this paper, we present a quasi-chemical model to simulate the optical properties of Au-Cu alloy systems based on the mixing of Gibbs free energy. This model is also able to predict the position of the surface plasmon resonance peaks for Au-Cu alloy nanoparticles. The model can be applied to predict the optical properties of alloy systems in the fields of plasmonics and nanophotonics.
Keywords:Dielectric properties; Gibbs free energy; Nanostructures; Metals; Optical properties
Metal nanoparticles (NPs) have attracted much research interest due to their unusual chemical and physical properties, such as catalytic activity, novel electronics, optics, and magnetic properties, and they have potential applications in solar cells and biosensors [1-7].
Alloy nanoparticle systems have been found to exhibit optical limiting properties due to surface plasmon resonance and have been used in biodiagnostic applications [8,9]. Alloy nanoparticles are materials used to tune the position of surface plasmon resonance, and thus help to produce materials for use in nonlinear optical applications [10-14]. Au-Cu alloy system is a completely dissoluble alloy. The position of surface plasmon resonance for Au NPs is about 520 nm. The position of surface plasmon resonance for Cu NPs is 570 ~ 580 nm . At low temperatures, Au, Au3Cu, AuCu, AuCu3, and Cu exist and order easily in Au-Cu alloys system. The prediction of the optical properties of such alloy systems is desirable if they are to be used in the design of optical devices. However, the optical properties of alloy systems are difficult to predict because of the random mixing of materials.
The quasi-chemical method is a statistical approach for predicting the short-range-order of Au-Cu alloys system according to Gibbs free energy. While the optical properties of Au-Cu alloys can be computed by the quasi-chemical model based on the energy potential between the electric field and induced dipole, few works have attempted to do this.
In this study, we thus simulate the optical properties of Au and Cu using a quasi-chemical model, based on the energy potential between the electric field and induced dipole. We then used this quasi-chemical method to modify the statistics for the short-range-order of Au-Cu alloy system. Then the optical properties are simulated by combining the Gibbs free energy and electric potential energy. The light extinction of nanoparticles is calculated by using Mie theory. The results show that the model is suitable for predicting the position of surface plasmon resonance peaks.
Au-Cu alloy system refers to a solid solution. Properties of a regular solution are best examined based on the concept of excess function . The excess value of an extensive thermodynamic solution property is simply the difference between its actual value and the value that it would have if the solution were ideal based on the Gibbs free energy of the solution,
In which G is the molar Gibbs free energy of the solution, Gid is the molar Gibbs free energy that the solution would have if it were ideal, and GXS is the excess molar Gibbs free energy of the solution. Because the two components have equal molar volumes and do not exhibit a change in molar volume when mixed, their regular solution behavior can be understood by the application of a statistical mixing model, i.e., a quasi-chemical model.
The energy of the solution is the sum of its interatomic bond energies. Consider 1 mol of a mixed crystal containing NA atoms of A and NB atoms of B such that
where NO is Avogadro’s number. The mixed crystal, or solid solution, contains three types of atomic bond: A-A bonds, B-B bonds, and A-B bonds. A-A bonds the energy of each of which is UAA, B-B bonds the energy of each of which is UBB, A-B bonds the energy of each of which is UAB. If in the solution, there are PAA A-A bonds, PBB B-B bonds, and PAB A-B bonds, the energy of the solution U is obtained as the linear combination
and the problem of calculating U becomes one of calculating the values of PAA, PBB, and PAB. Thus,
The change in volume is negligible. Since ΔVM = 0,
Ideal mixing requires the condition UAB = UAA = UBB. If ΔHM = 0, the mixing of the NA atoms with the NB atoms of B is random.
The quasi-chemical model is a statistical mixing model in Gibbs free energy. According to Equations 5 and 6, the mixing Gibbs free energy will be presented. In the ‘Results and discussion’ section, the dipole energy in Gibbs free energy was utilized to consider the optical properties with different frequencies of incident light.
Results and discussion
The probability that a neighboring pair of sites contains an A-B pair is 2XAXB, an A-A pair is XA2, and B-B pair is XB2, and The quasi-chemical model is a statistical mixing model that describes the mixing cluster. The difference in Gibbs energy is presented as follows:
Combining Equations 5 and 6 with Equation 7 produces the following:
Because PAB = 2XAXB,
The Gibbs free energy of the solution is as follows:
where is the induced dipole moment of metamaterial, is the induced dipole moment of material A, is the induced dipole moment of material B, and is the induced dipole moment due to the interaction of materials A and B.
The Gibbs energy was subtracted when applying an electric field from that without applying one, as follows:
The dielectric function of the mixed material includes the interaction term and independent terms ϵAXA2 and ϵBXB2. When is assumed to be an experience constant, Λ, the dielectric function of mixing material is reduced to the following form:
The Newton formula  is used to apply these concepts to the clustered material. The dielectric functions refer to these clusters and the embedding matrix,
The quasi-chemical model describes the mixing of clusters in which the interaction term is approximated by Newton formula mixing. Combining the probability of neighboring pairs with the Newton formula, the optical model of the regular solution is as follows:
The effective dielectric complex of the alloy is presented in Figure 1.
Figure 1. Effective dielectric complex of the alloy. (a) Real part, ϵr. (b) Imaginary part, ϵi, of the dielectric complex of Au-Cu alloy.
According to Mie theory [18,19], the resonances denoted as surface plasmon were relative with the onset of the quantum size and shape effects of Au NPs. There is one SPR band for metal NPs, and this is shown as follows [20,21]:
where ϵh is the dielectric constant of the host medium embedding Au NPs, ϵm is the dielectric constant of Au NPs, f is the volume fraction of Au NPs, ϵi is the total dielectric constant, and Γi is a set of three parameters defined along the principal axes of the particle characterizing its shape. Γ1 + Γ2 + Γ3 = 1 and the other parameters range from 0 to 1. The frequencies of the surface plasmon of nonspherical metal NPs have two or three bands, depending on their shape. The extinction coefficients of alloy metal NPs with different sizes and environments are presented in Figures 2, 3, 4.
Figure 2. Extinction of Au-Cu alloy nanoparticles. Extinction of Au-Cu alloy nanoparticles (10 nm) when (a)n = 1, (b)n = 1.4, and (c)n = 1.8 (Qabs is the extinction coefficient).
Figure 3. Extinction of different sized NPs. (a) Au, (b) Au3Cu, (c) AuCu, (d) AuCu3, and (e) Cu alloy nanoparticles (n = 1; Qabs is the extinction coefficient).
Figure 4. Extinction of different refractive index. (a) Au, (b) Au3Cu, (c) AuCu, (d) AuCu3, and (e) Cu alloy nanoparticles.
The quasi-chemical model is used to calculate the optical properties of Au-Cu alloys. The real part of the dielectric complex is negative for Au-Cu alloy system. The imaginary part of dielectric constant for Au-Cu alloy system shows the peaks that appear in range from 430 to 520 nm due to the electronic transition between the d band and sp band. The real and imaginary parts of the dielectric complex for Au-Cu alloys system are as shown in Figure 1a,b, respectively.
We use Mie theory to predict the spectrum and position of surface plasmon resonance. Figure 2b shows the extinction of a 10-nm diameter Au-Cu nanoparticle in different refractive index surroundings. For n = 1.4, the surface plasmon resonance peaks are 532, 538, 561, 567, and 578 nm for Au, Au3Cu, AuCu, AuCu3, and Cu, respectively, and these results which are in agreement with those of other experimental results .
The extinction spectra of Au-Cu bimetallic nanoparticle with size effect are presented in Figure 3. As the size of nanoparticles increase, the peak of surface plasmon resonance red-shifts. When the size is less than 50 nm, the size effect becomes more significant. The higher the ratio of Cu to Au of is, the more the surface plasmon resonance red-shifts. When the size is greater than 50 nm, the size effect is less significant due to the small increase in the cross section. The refractive index effect is shown in Figure 4. As the refractive index increases, the surface resonance peak will red-shift and become increasingly sharp. Based on this, it is possible to predict the surface plasmon resonance peaks of regular solution alloys, such as Au-Cu, Cu-Ag, Ag-Cu, and Au-Cu-Ag systems.
In this work we used the quasi-chemical model to compute the optical properties of Au-Cu alloy system. The results show that it is possible to use this approach to predict the positions of surface plasmon resonance peaks. This model is thus a useful tool in the development of for future applications of alloy nanoparticles for plasmonics and nanophotonics.
The authors declare that they have no competing interests.
YHS and WLW contribute in writing and model setting in all these works. Both authors read and approved the final manuscript.
YHS is an assistant professor and WLW is a student in the Department of Materials Science and Engineering in National Cheng Kung University, Taiwan.
This work was financially supported by the National Science Council of Taiwan (nos. 100-2218-E-259-003-MY3 and 102-2221-E-006-293-MY3) which is gratefully acknowledged. This research was, in part, supported by the Ministry of Education, Taiwan, Republic of China and the Aim for the Top University Project of the National Cheng Kung University (NCKU).
Wustholz KL, Henry AI, McMahon JM, Freeman RG, Valley N, Piotti ME, Natan MJ, Schatz GC, Van Duyne RP: Structure-activity relationships in gold nanoparticle dimers and trimers for surface-enhanced Raman spectroscopy.
Appl Phys Lett 2012, 101:243901. Publisher Full Text
J Phys Chem C 2013, 117:13676. Publisher Full Text
Light Sci Appl 2012, 1:e14. Publisher Full Text
Mater Today 2013, 16:133. Publisher Full Text
Appl Phys Lett 2013, 103:063703. Publisher Full Text
J Amer Chem Soc 2001, 123:7961. Publisher Full Text
Mater Sci Eng B 2005, 117:173. Publisher Full Text
J Am Chem Soc 2000, 122:11005. Publisher Full Text
J Phys Chem B 2000, 104:6133. Publisher Full Text
Phys Rev B 2000, 62:13160. Publisher Full Text
J Phys Chem B 2000, 104:564. Publisher Full Text
Phys Rev B 1975, 11:1732. Publisher Full Text
J Mater Chem 2003, 13:1789. Publisher Full Text