SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series EMN Meeting.

Open Access Nano Express

Surface plasmon resonance of Au-Cu bimetallic nanoparticles predicted by a quasi-chemical model

Yen-Hsun Su12* and Wen-Lin Wang1

Author Affiliations

1 Department of Materials Science and Engineering, National Cheng Kung University, Tainan, Taiwan

2 Advanced Optoelectronic Technology Center, National Cheng Kung University, Tainan, Taiwan

For all author emails, please log on.

Nanoscale Research Letters 2013, 8:408  doi:10.1186/1556-276X-8-408

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/8/1/408

Received:17 July 2013
Accepted:28 August 2013
Published:2 October 2013

© 2013 Su and Wang; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


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.

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 [15]. 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.



Regular solution

Au-Cu alloy system refers to a solid solution. Properties of a regular solution are best examined based on the concept of excess function [16]. 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,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M1">View MathML</a>


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.

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

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M2">View MathML</a>


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

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M3">View MathML</a>


and the problem of calculating U becomes one of calculating the values of PAA, PBB, and PAB. Thus,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M4">View MathML</a>


The change in volume is negligible. Since ΔVM = 0,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M5">View MathML</a>


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.

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M6">View MathML</a>


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 <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M7">View MathML</a> The quasi-chemical model is a statistical mixing model that describes the mixing cluster. The difference in Gibbs energy is presented as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M8">View MathML</a>


Combining Equations 5 and 6 with Equation 7 produces the following:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M9">View MathML</a>


Because PAB = 2XAXB,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M10">View MathML</a>


The Gibbs free energy of the solution is as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M11">View MathML</a>


After applying the electric field <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M12">View MathML</a>,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M13">View MathML</a>


where <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M14">View MathML</a> is the induced dipole moment of metamaterial, <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M15">View MathML</a> is the induced dipole moment of material A, <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M16">View MathML</a> is the induced dipole moment of material B, and <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M17">View MathML</a> 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:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M18">View MathML</a>


Because <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M19">View MathML</a>, the above equation can be rewritten as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M20">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M21">View MathML</a>


The dielectric function of the mixed material includes the interaction term <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M22">View MathML</a> and independent terms ϵAXA2 and ϵBXB2. When <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M23">View MathML</a> is assumed to be an experience constant, Λ, the dielectric function of mixing material is reduced to the following form:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M24">View MathML</a>


The Newton formula [17] is used to apply these concepts to the clustered material. The dielectric functions refer to these clusters and the embedding matrix,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M25">View MathML</a>


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:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M26">View MathML</a>


The effective dielectric complex of the alloy is presented in Figure  1.

thumbnailFigure 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]:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/408/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/408/mathml/M27">View MathML</a>


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.

thumbnailFigure 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).

thumbnailFigure 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).

thumbnailFigure 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 [22].

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.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

YHS and WLW contribute in writing and model setting in all these works. Both authors read and approved the final manuscript.

Authors’ information

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).


  1. Banholzer MJ, Osberg KD, Li S, Mangelson BF, Schatz GC, Mirkin CA: Silver-based nanodisk codes.

    ACS Nano 2010, 4:5446. PubMed Abstract | Publisher Full Text OpenURL

  2. 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.

    J Am Chem Soc 2010, 132:10903. PubMed Abstract | Publisher Full Text OpenURL

  3. Zhang XL, Song JF, Li XB, Feng J, Sun HB Sun : Optical Tamm states enhanced broad-band absorption of organic solar cells.

    Appl Phys Lett 2012, 101:243901. Publisher Full Text OpenURL

  4. Sen A, Lin CJ, Kaun CC: Single-molecule conductance through chiral gold nanotubes.

    J Phys Chem C 2013, 117:13676. Publisher Full Text OpenURL

  5. Su YH, Ke YF, Cai SL, Yao QY: Surface plasmon resonance of layer-by-layer gold nanoparticles induced photoelectric current in environmentally-friendly plasmon-sensitized solar cell.

    Light Sci Appl 2012, 1:e14. Publisher Full Text OpenURL

  6. Stratakis E, Kymakis E: Nanoparticle-based plasmonic organic photovoltaic devices.

    Mater Today 2013, 16:133. Publisher Full Text OpenURL

  7. Su YH, Hsu CY, Chang CC, Tu SL, Shen YH: Ultra-thin titanium nanolayers for plasmon-assisted enhancement of bioluminescence of chloroplast in biological light emitting devices.

    Appl Phys Lett 2013, 103:063703. Publisher Full Text OpenURL

  8. Cao YW, Jin R, Mirkin CA: DNA-modified core-shell Ag/Au nanoparticles.

    J Amer Chem Soc 2001, 123:7961. Publisher Full Text OpenURL

  9. Nair AS, Suryannarayanan V, Pradeep T, Thomas J, Anija M, Philip R: AuxAgy@ZrO2 core - shell nanoparticles: synthesis, characterization, reactivity and optical limiting.

    Mater Sci Eng B 2005, 117:173. Publisher Full Text OpenURL

  10. Ispasoiu RG, Balogh L, Varnavski OP, Tomalia DA, Goodson T: Large optical limiting from novel metal-dendrimer nanocomposite materials.

    J Am Chem Soc 2000, 122:11005. Publisher Full Text OpenURL

  11. Sun WF, Dai Q, Worden JG, Huo Q: Optical limiting of a covalently bonded gold nanoparticle/polylysine hybrid material.

    J Phys Chem B 2005, 109:20854. PubMed Abstract | Publisher Full Text OpenURL

  12. Wang G, Sun WF: Optical limiting of gold nanoparticle aggregates induced by electrolytes.

    J Phys Chem B 2006, 110:20901. PubMed Abstract | Publisher Full Text OpenURL

  13. Francois L, Mostafavi M, Belloni J, Delouis JF, Delaire J, Feneyrou P: Optical limitation induced by gold clusters.

    J Phys Chem B 2000, 104:6133. Publisher Full Text OpenURL

  14. Philip R, Kumar GR, Sandhyarani N, Pradeep T: Picosecond optical nonlinearity in monolayer-protected gold, silver, and gold-silver alloy nanoclusters.

    Phys Rev B 2000, 62:13160. Publisher Full Text OpenURL

  15. Yeh YH, Yeh MS, Lee YP, Yeh CS: Formation of Cu nanoparticles from CuO powder by laser ablation in 2-propanol.

    Chem Lett 1998, 11:1183. OpenURL

  16. Gaskell DR: Introduction to the Thermodynamics of Materials. 5th edition. New York: Taylor & Francis; 2008:249. OpenURL

  17. Theiss W:

    Optische Eigenschaften Inhomogener Materialien. Dissertation, RWTH Aachen. 1989. OpenURL

  18. Templeton AC, Pietron JJ, Murray RW, Mulvaney P: Solvent refractive index and core charge influences on the surface plasmon absorbance of alkanethiolate monolayer-protected gold clusters.

    J Phys Chem B 2000, 104:564. Publisher Full Text OpenURL

  19. Fuchs R: Theory of the optical properties of ionic crystal cubes.

    Phys Rev B 1975, 11:1732. Publisher Full Text OpenURL

  20. Prasad PN: Nanophotonics. Hoboken: Wiley; 2004:130. OpenURL

  21. Fox M: Optical Properties of Solids. London: Oxford University Press; 2001:151. OpenURL

  22. Kim MJ, Na HJ, Lee KC, Yoo EA, Lee M: Preparation and characterization of Au-Ag and Au-Cu alloy nanoparticles in chloroform.

    J Mater Chem 2003, 13:1789. Publisher Full Text OpenURL