Abstract
The recently developed ωFQFμ model (ACS Photonics, 9, 3,025–3,034) is extended to bimetallic nanoparticles, such as nanoalloys and core-shell systems. The method finds its grounds in basic physical concepts, such as Drude conduction theory, electrostatics, interband transitions, and quantum tunneling. The approach, which is parametrized on ab initio simulations of Ag-Au nanoalloys, is challenged against complex Ag-Au nanostructures (spheres, nanorods, and core-shell nanoparticles). Remarkable agreement with available experimental data is found, thus demonstrating the reliability of the newly developed approach.
1 Introduction
Metal nanoparticles (NPs) exhibit unique optical properties, which are mainly due to the formation of surface plasmons, i.e., collective excitations of conductive electrons. (Moskovits, 1985; Nie and Emory, 1997; Maier, 2007; ; ; Santhosh et al., 2016). At the plasmon resonance frequency (PRF) such surface plasmons exhibit a resonant behavior, which gives rise to a huge enhancement of the electric field in the proximity of the NP surface. This phenomenon is exploited to increase the detection limit of common analytical techniques. (; Maier et al., 2003; Muehlschlegel et al., 2005; Lim et al., 2010; ; Neuman et al., 2018). PRF can be tuned by varying the NP shape, dimension, and chemical composition (Ag, Au, Al, …). As an alternative, bimetallic nanoalloys or core-shell NPs can be used. Such systems are constituted by two different metal elements, e.g., Ag-Au, Au-Cu, Cu-Ag (Xiang et al., 2008; Wang et al., 2009; ; ; ; Ma et al., 2020; ). In this case, the PRF can be tuned not only by modifying the aforementioned variables (size, shape, …) but also the relative concentration of the two metals.
Many theoretical approaches have been developed to describe the plasmonic properties of single metal NPs(; ; ; Myroshnychenko et al., 2008; ; ; ; ; ; ; ; ; ; ), among which is worth mentioning quantum hydrodynamic models, which are able to describe both nonlocal and electron spill-out effects (; Raza et al., 2015; ). However, the most widely exploited methods are based on classical physics, (; ; ; ; Pérez-González et al., 2010; ; ; ), such as the Mie Theory (Mie, 1908), the finite difference time domain (FDTD) (Shuford et al., 2006), the Discrete Dipole Approximation (DDA) () or exploit a continuum representation of the NP by means of the Boundary Element Method (BEM) (; ; ; ). All the classical approaches mentioned above rely on defining a suitable permittivity function ɛ. Clearly, in the case of bimetallic systems, ɛ-dependent methods are not flexible enough to treat a priori any chemical composition and need to resort to experimentally measured ɛ values for specific alloy concentrations or to approximations, which may not be physically justified (such as a linear combination of the ɛ values of the two metals) (Peña-Rodríguez and Pal, 2011a; Verbruggen et al., 2013; Ma et al., 2015; Putra et al., 2017; ; Newmai et al., 2022).
Atomistic approaches appear as the most natural choice to overcome this problem. Full ab initio methods can be exploited, however, their prohibitive computational cost limits their application to relatively small systems (; López Lozano et al., 2013; ; ; Olobardi et al., 2019; ; ). For this reason, classical atomistic approaches have been developed, (; ; ; Zakomirnyi et al., 2019; Zakomirnyi et al., 2020), however they have only marginally been applied to the description of bimetallic NPs (Sørensen et al., 2021).
In this work, we extend ωFQFμ () to simulate the optical properties of bimetallic NPs, with special emphasis on Ag-Au nanoalloys and core-shell systems. ωFQFμ is a fully atomistic, classical, approach, which assigns each atom of the metal nanoalloy with an electric charge and an electric dipole moment, which vary as a response to an externally applied electric field. The theoretical foundations of ωFQFμ lay in Drude’s theory of conduction, classical electrodynamics, interband transitions, and quantum tunneling. In fact, the Drude mechanism of charge exchange is modeled through the equation of motion of the electric charges (; ; ; ; ), while interband contributions are taken into account by means of effective interband polarizability () introduced to mimic d-shell polarizability (). Also, quantum tunneling effects, which play a crucial role in nanojunctions and nanoaggregates (; ), are described by tuning charge exchange through a phenomenological function. Therefore, ωFQFμ conceptually differs from other approaches which are based on classical physics only.
The manuscript is organized as follows. In the next section, the formulation of ωFQFμ for bimetallic systems is presented. After a brief section presenting the computational details, the method is tested against the reproduction of the plasmonic properties of alloyed spherical NPs, alloyed nanorods, and spherical core-shell systems. Conclusions and a sketch of the future perspectives of the approach end the manuscript.
2 Materials and methods
2.1 Theoretical model
In this section, we extend ωFQFμ () to describe the optical properties of bimetallic systems. To this end, we follow the derivation reported in () and (). In ωFQFμ, the fluctuations of the electronic density of a plasmonic substrate composed of N atoms emerging as a response to the external electric field are represented by a set of N discrete fluctuating complex-valued charges qi and N fluctuating complex-valued dipoles μi, which are located at atomic positions. The equation of motion of charges q originates from the continuity equation () (atomic units are used throughout the paper):where ni is the electron density on atom i, is the momentum of the electron associated with atom i averaged over all the trajectories towards atom j. is the unit vector connecting i and j. Aij is the effective area dividing ith and jth atoms, and its value modulates the charge exchange between each atom pair. By following (; ; ), Aij is approximated as an atomic parameter (Ai).The electron momentum pi can be estimated by means of the Drude model of conductance () as follows:where Ei(t) is the electric field acting on the ith atom and τi is the relaxation time associated with scattering events. Assuming the external uniform monochromatic electric field Eext(t) = Eext(ω) exp (−iωt) oscillates at frequency ω, we can reformulate Eqs 1, 2 in the frequency domain, i.e.:Following the original derivation reported in (), we can write:where ϕi is the chemical potential of atom i and lij = lji is the distance between ith and jth atoms. By this, Eq. 3 becomes:where quantum tunneling effects are expressed in terms of a Fermi-like function fij (lij), which reads:where is the equilibrium distance between atoms i and j, and the parameters dij and sij determine the shape of the damping function. Notice that the inclusion of a phenomenological description of quantum tunneling is needed to correctly describe the optical response of plasmonic subnanometer junctions and hot-spots (; ; Scholl et al., 2013; ; ; ). Indeed, in the case of plasmonic nanoaggregates, charge transfer plasmons (CT) may occur. In ωFQFμ, the charge exchange among different NPs is governed by Drude and tunneling mechanisms, and the total system charge is conserved (see Eq. 5), allowing for a physically consistent description of such plasmonic modes. To make the notation more compact, the following terms are introduced ():Dividing Eq. 5 by wi(ω), we obtain:where the chemical potential ϕi can be written as:The first and the second terms in Eq. 11 are the electric potential generated by the charges and the dipoles of the system, mediated by the interaction tensors Tqq and Tqμ (see () for their definitions). is the electric potential associated with the external electric field. By plugging Eq. 11 into Eq. 10, the linear equations ruling charge evolution are obtained, i.e.:where the following matrices are introduced:The electric dipoles μi, which are introduced to properly model interband contributions, are obtained by solving the following set of linear equations:where αi(ω) is the frequency-dependent polarizability of the ith atom, related to interband terms, which multiplies the total electric field acting on the ith dipole, composed by the external field , the field generated by the fluctuating charges, , and the field generated by the fluctuating dipoles, . If the system is homogeneous, αi is calculated as reported in ():where is the interband contribution to the frequency-dependent permittivity function and nA is the number density of atoms of the chemical element A. ɛIB is calculated by subtracting the Drude term from the frequency-dependent permittivity.
If the system is heterogeneous, the link between the macroscopic, bulk dynamical permittivity ɛIB(ω) and the microscopic, atomistic frequency-dependent interband polarizability αi(ω) is not straightforward. Let us focus on bimetallic systems, composed of two chemical elements A and B. We assume the interband polarizability of the ith atom belonging to the A moiety to be a function of the local composition of the system, i.e., , where NA,i and NB,i are the number of nearest neighbors of A and B type, respectively. The simplest approach is to model as the weighted average of and . Therefore, the interband frequency-dependent polarizability of the ith atom of species A in the alloy composed by materials A and B is defined as follows:Notice that, possible alternative approaches, e.g., to resort to arithmetic averages have been investigated. However, the weighted harmonic mean indeed gives the best numerical results (see Supplementary Figure S1 and Sec. S1 in the Supplementary Materials – SM).
Starting from Eq. 17, the final expression for the electric dipoles becomes ():where the field generated by the multipoles is mediated by the interaction tensors Tμq and Tμμ, defined according to ().
The coupled charge-dipole equations can finally be written in a compact notation as follows:where ZIB is a diagonal matrix, of which the elements are defined according to Eq. 19. Notice that in the case of homogenous systems, i.e., A = B, the standard ωFQFμ equations presented by are recovered. From the values of charges and dipoles, the complex polarizability can be computed, from which the absorption cross section σabs can be calculated:In the previous equations d is the total complex dipole moment, i runs over NP atoms, k represents x, y, z positions of the ith atom, and l runs over x, y, z directions of the external electric field. c is the speed of light and is the imaginary part of the complex polarizability .
2.2 Computational details
The developed ωFQFμ approach is challenged to reproduce the optical response of several Ag-Au nanostructures. In particular, we first consider Ag-Au spherical nanoalloys with diameter D = 5.2 nm (4,347 atoms) and Ag-Au nanorods with D = 2.5 nm and length L = 10 nm (2,560 atoms), which are generated by using Atomic Simulation Environment (ASE) Python module v. 3.17 (). A lattice constant of 4.08 Å () and a Face-Centered Cubic (FCC) packing are exploited. For both spheres and nanorods, ten alloy compositions are considered by starting from pure Ag structures and increasing the percentage of Au atoms, which randomly replace Ag atoms (from 0% to 100%, with a constant step of 10%). To gain statistical significance, for each Au percentage, ten nanoalloy structures are generated by randomly replacing the proper fraction of Ag with Au, similarly to the strategy followed in previous studies (Sørensen et al., 2021). As a further example, core-shell spherical systems, which are constituted by an inner Au sphere (D = 5.0 nm, 3,851 atoms) surrounded by an outer Ag shell (D = 6.25 nm, 3,698 atoms), are investigated.
ωFQFμ equations are solved by using a stand-alone Fortran 95 package. ωFQFμ parameters (see Eqs 5, 6, 18) for Ag and Au atoms are taken from (). In particular, Drude and interband parameters are recovered from Ag and Au experimental permittivity functions (). The parameters dij = 12.0 and sij = 0.95 entering Eq. 6 are obtained by fitting ωFQFμ results for i = Ag and j = Au (and vice versa) on reference time-dependent density functional theory + tight binding (TD-DFT + TB) () absorption spectra of Ag-Au clusters (see Sec. S2 in the SM).
3 Results and discussion
3.1 Spherical nanoparticles and nanorods
The modeling of the dependence of the absorption cross-section σabs of a spherical NP on the Au concentration has received much attention in the literature (Papavassiliou, 1976; Link et al., 1999; Rioux et al., 2014; Ristig et al., 2015). In particular, a simple linear combination of Ag and Au permittivity functions combined with the Mie theory cannot reproduce the linear behavior of the plasmon red-shift upon increasing Au concentration which is experimentally observed (Link et al., 1999). This issue can be solved by following the strategy proposed by Rioux et al. (2014), who have developed an analytical and rigorous model to predict a composition-dependent complex dielectric function of Ag-Au alloys, based on critical point analysis of the band-structure of Ag and Au. By coupling such modeling to the Mie theory, the correct trend is followed for spherical nanoalloys. It is worth noting that such a strategy might be exploited within the more refined, yet continuum, BEM method, thus overcoming the limitations of this approach in describing the plasmonic properties of nanoalloys. However, while the methodology proposed by Rioux et al. (2014) is rigorous from the theoretical point of view, its extension to different alloy compositions, and to alloys of different chemical nature rather than Au and Ag is far from trivial, due to the high (about 30) number of parameters that would need to be fitted to obtain the composition dependent dielectric function.
Computed ωFQFμ values for σabs of a spherical NP with different Au concentrations are reported in (see Figure 1). In all cases, spectra are characterized by a main peak (associated with the dipolar mode, see Supplementary Figure S2 in the SM), which redshifts and lowers in intensity as the Au concentration increases. By plotting the PRF (in nm) as a function of the percentage of Au (see Figure 1C), a linear trend is observed (R2 ∼ 1.00). This behavior perfectly follows Vegard’s law, which reads:
FIGURE 1
where x is the percentage of Au atoms. Such a linear dependence of the absorption wavelength on Au percentage is reported by many experimental works (Papavassiliou, 1976; Link et al., 1999; Rioux et al., 2014; Ristig et al., 2015). The slope (2.06 nm/% Au) of the line fitted on ωFQFμ results is in good agreement with experimental data (Link et al., 1999) (∼ 1.35 nm/% Au). Furthermore, by moving from pure Ag to pure Au NPs, the intensity of the plasmon band exponentially decreases (R2 = 0.99, see Figure 1D), thus perfectly reproducing experimental findings (Link et al., 1999).
We note that, remarkably, other classical atomistic approaches may fail at reproducing the experimentally reported linear trend in wavelength and the exponentially decreasing intensity (Sørensen et al., 2021). Differently, ωFQFμ can correctly reproduce all the plasmonic properties of these structures, which is a key result of this paper.
We now move to discuss the plasmonic properties of Ag-Au alloy nanorods. Also in this case, we study the dependence of the absorption frequency on the Au concentration, which varies from 0% (pure Ag nanorod) to 100% (pure Au nanorod), with a constant step of 10%. ωFQFμ absorption cross sections, which are reported in Figure 2B, are clearly dominated by an intense peak which redshifts and decreases in intensity as the gold concentration increases. This peak is associated with the longitudinal dipolar mode (see also Supplementary Figure S3 in the SM). From a comparison with spherical NPs (see Figure 1B), we note that the dipolar peak moves to lower energies, independently of the chemical composition of the alloy nanorod. Also, in this case, additional bands can be appreciated in the region between 350 and 500 nm. There, octupolar plasmons arise (), together with transversal dipolar modes. The associated bands associated are graphically depicted as an inset in Figure 2B, and redshift and lower in intensity as the Au concentration increases.
FIGURE 2
The PRFs (in nm) of both the longitudinal and transversal dipolar modes as a function of the Au percentage are graphically depicted in Figure 2C. The transversal mode follows the linear trend (circles, R2 = 0.99) predicted by Vegard’s law (see Eq. 24), while for the longitudinal plasmon, an evident deviation from the linear regime is observed, especially for Au concentrations larger than 40%. Our findings are in agreement with the experimental measurements reported by , where for the transversal peak a linear trend is measured. Also, the experimental slope (1.47 nm/%Au) correlates significantly well with our calculations (1.92 nm/%Au). Differently, the longitudinal absorption wavelength is reported to be rather independent of the chemical composition for Au percentages 40%, while a slight blueshift is expected by increasing the Ag concentration . Note that a DDA description is able to reproduce the composition dependence of the longitudinal mode but fails at describing the linear trend of the transversal peak (). This is probably related to the exploited linear combination of Ag and Au permittivity functions. Remarkably, ωFQFμ correctly models both behaviors.
These findings can be explained by considering that the longitudinal dipolar mode falls at energies that are reasonably far from interband transition regions of both Ag and Au. In this regime, Drude free electrons play a dominant role (Link et al., 1999; ) and the good agreement with experiments reported by both DDA () and ωFQFμ show that the permittivity function in this region can be reasonably approximated as a linear combination of Ag and Au Drude contributions. Differently, the transversal plasmon mode absorbs at much higher energies, where interband effects cannot be neglected (). This proves the reliability of ωFQFμ, which, differently from other models (), describes the plasmonic resonance by means of both electric charges and dipoles, modeling the intraband and interband mechanisms, respectively.
3.2 Au@Ag core-shell NPs
ωFQFμ is not limited to the description of nanoalloys, but can generally describe bimetallic nanostructures. To showcase its applicability, we consider pure core-shell (Au@Ag-0%) spherical NPs with a diameter D of 6.25 nm (3,698 atoms, see Figure 3A), constituted by a Au core (DAu = 5 nm), and an external Ag shell (external DAg = 6.25 nm, DAu/DAg = 0.8). In addition, to show the flexibility of ωFQFμ, the effect of alloying such structures, which has been experimentally realized (Ristig et al., 2015; ), is also taken into account, by randomly replacing 18% of the core Au atoms with Ag, and viceversa for the Ag shell (Au@Ag-18%, see Figure 3A).
FIGURE 3
Computed ωFQFμ absorption cross sections in the visible range for both core-shell systems are reported in Figure 3B, together with reference spectra for Ag and Au spherical NPs with D = 6.25 nm. Similarly to Figure 1B, Ag and Au reference absorption spectra are characterized by a main peak, located at about 360 nm and 570 nm for Ag and Au, respectively. The spectrum of the pure Au@Ag-0% core-shell is instead dominated by two peaks A and B, in agreement with (Peña-Rodríguez and Pal, 2011a; ; ), which fall at about 356 nm (A) and 517 nm (B). To deeply investigate the nature of the plasmons associated with these two bands, we graphically report in Figure 3C, left panel, the electron densities at the excitation energies for pure Au@Ag-0% core-shell (on a slice in the xy plane). As it can be appreciated, both peaks are associated with a global dipolar plasmon. However, they differ for the electron distribution in both the Au core and the interface between the two regions (see yellow dashed line in Figure 3C). Indeed, for peak A a huge charge accumulation at the interface is reported, while for peak B, the boundary between the Ag shell and the Au core is only partially marked. This is due to the fact that our model allows for charge exchange between the two layers, which is for instance not accounted for by other classical descriptions (Szántó et al., 2021). However, such a charge transfer is limited by the different tunneling barriers between the two metals, as introduced in Eq. 6, and by the different chemical nature (chemical hardness, polarizability, Drude parameters, …). As a consequence, charge accumulation at the Ag-Au interface is expected.
Similarly, the Au@Ag-18% absorption cross section (see Figure 3B) is characterized by two main peaks A′ and B′, which are blue- and red-shifted as compared to their counterparts in Au@Ag-0% system. Such a behavior can be justified by considering that by increasing the alloying degree up to 50%, a fully alloyed Ag-Au spherical NP is obtained. Therefore, in this situation, a single band is expected at about 450 nm (see also Figure 1B). The plasmonic nature of the electron densities computed at the A′ and B′ frequencies is graphically reported in Figure 3C, right panel. A global dipolar plasmon is observed. However, by also comparing with Au@Ag-0%, charge accumulation at the core-shell boundary is still present, but less pronounced, in particular for B′. This confirms the above speculation: indeed, by increasing the alloying degree, a decrease in the potential barrier at the interface is obtained.
To conclude the discussion on Au@Ag core-shell systems, we compare our results with experimental absorption cross sections reported by Peña-Rodríguez and Pal (2011b) (see inset in Figure 3B), who studied the plasmonic response of Au@Ag core-shell spherical NPs with D = 44.2 nm and characterized by DAu/DAg = 0.81 (comparable with our simulations). As it can be noticed, the experimental spectrum is characterized by two main peaks at about 400 nm and 497 nm. This is in very good agreement with our calculations, thus further demonstrating the reliability of our model. In fact, only a slight discrepancy in relative intensities, probably due to different experimental conditions (size of the NP, solvent, …), is observed.
4 Conclusion
ωFQFμ has been extended to the description of bimetallic metal nanoalloys and core-shell NPs, and has been applied to spherical alloyed NPs and nanorods, for which alternative methods proposed in the literature fail, due to incorrect modeling of interband contributions. Remarkably, ωFQFμ correctly reproduces experimental trends, also in the case of core-shell systems; this is a direct consequence of the atomistic nature of the approach, which permits a physically consistent picture of the local environment of each specific atom.
ωFQFμ is general enough to describe any kind of bimetallic system; for instance, it can be applied to complex geometrical arrangements, such as subnanometer junctions, which form when two NPs approach each other. Such complex nanostructures require appropriate treatment of interband transitions and reliable modeling of quantum tunneling effects. The effect of alloying in subnanometer junctions has been only marginally investigated, due to the lack of theoretical approaches to correctly describe the plasmonic properties of alloys and generic bimetallic systems. ωFQFμ can indeed be applied to these systems and this investigation will be the topic of future studies. Also, ωFQFμ can be extended to generic multi-metallic nanostructures, through the generalization of Eq. 19 followed by a specification of the quantities entering Eq. 6.
As a final technological perspective, relevant for instance in sensing applications, the approach can be potentially applied to compute spectral properties of molecules adsorbed on bimetallic systems, by generalizing previous studies on homogenous substrates (Payton et al., 2013; ). In this way, Surface-Enhanced spectroscopic signals can be simulated by retaining an atomistic description of the plasmonic substrate, which plays a crucial role in bimetallic systems, as it is demonstrated in this work.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
LN, PL, and TG developed the theoretical methodology. LN, PL, and PG ran ωFQFμ calculations. TG, LN, and PL implemented the method into the stand-alone Fortran95 code for ωFQFμ calculations. LN, PL, PG, LB, and TG analyzed data and wrote the original draft manuscript. TG and CC discussed and supervised the whole project. CC acquired the funding and finalized the manuscript (final version). All authors contributed to the article and approved the submitted version.
Funding
This work has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement N. 818064). CC acknowledges funding from PNRR MUR project PE0000023-NQSTI.
Acknowledgments
We gratefully acknowledge the Center for High-Performance Computing (CHPC) at SNS for providing the computational infrastructure.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphot.2023.1199598/full#supplementary-material
References
1
AizpuruaJ.HanarpP.SutherlandD.KällM.BryantG. W.García de AbajoF. J. (2003). Optical properties of gold nanorings. Phys. Rev. Lett.90, 057401. 10.1103/physrevlett.90.057401
2
AnkerJ. N.HallW. P.LyandresO.ShahN. C.ZhaoJ.Van DuyneR. P. (2008). Biosensing with plasmonic nanosensors. Nat. Mater.7, 442–453. 10.1038/nmat2162
3
Asadi-AghbolaghiN.RugerR.JamshidiZ.VisscherL. (2020). Td-dft+ tb: An efficient and fast approach for quantum plasmonic excitations. J. Phys. Chem. C124, 7946–7955. 10.1021/acs.jpcc.0c00979
4
AtwaterH. A.PolmanA. (2010). Plasmonics for improved photovoltaic devices. Nat. Mater.9, 205–213. 10.1038/nmat2629
5
AwadaC.DabC.GrimaldiM. G.AlshoaibiA.RuffinoF. (2021). High optical enhancement in au/ag alloys and porous au using surface-enhanced Raman spectroscopy technique. Sci. Rep.11, 4714. 10.1038/s41598-021-84093-0
6
BadeW. (1957). Drude-model calculation of dispersion forces. i. general theory. J. Chem. Phys.27, 1280–1284. 10.1063/1.1743991
7
BaghramyanH. M.Della SalaF.CiracìC. (2021). Laplacian-level quantum hydrodynamic theory for plasmonics. Phys. Rev. X11, 011049. 10.1103/physrevx.11.011049
8
BarcaroG.BroyerM.DuranteN.FortunelliA.StenerM. (2011). Alloying effects on the optical properties of ag–au nanoclusters from tddft calculations. J. Phys. Chem. C115, 24085–24091. 10.1021/jp2087219
9
BarcaroG.SementaL.FortunelliA.StenerM. (2014). Optical properties of pt and ag–pt nanoclusters from tddft calculations: Plasmon suppression by pt poisoning. J. Phys. Chem. C118, 28101–28108. 10.1021/jp508824w
10
BarcaroG.SementaL.FortunelliA.StenerM. (2015). Optical properties of nanoalloys. Phys. Chem. Chem. Phys.17, 27952–27967. 10.1039/c5cp00498e
11
BlommaertsN.VanrompayH.NutiS.LenaertsS.BalsS.VerbruggenS. W. (2019). Unraveling structural information of turkevich synthesized plasmonic gold–silver bimetallic nanoparticles. Small15, 1902791. 10.1002/smll.201902791
12
BokH.-M.ShufordK. L.KimS.KimS. K.ParkS. (2009). Multiple surface plasmon modes for gold/silver alloy nanorods. Langmuir25, 5266–5270. 10.1021/la803900w
13
BonattiL.GilG.GiovanniniT.CorniS.CappelliC. (2020). Plasmonic resonances of metal nanoparticles: Atomistic vs. continuum approaches. Front. Chem.8, 340. 10.3389/fchem.2020.00340
14
BonattiL.NicoliL.GiovanniniT.CappelliC. (2022). In silico design of graphene plasmonic hot-spots. Nanoscale Adv.4, 2294–2302. 10.1039/d2na00088a
15
CaoM.LiuQ.ChenM.ChenL.YangD.HuH.et al (2018). Fully alloying auag nanorods in a photothermal nano-oven: Superior plasmonic property and enhanced chemical stability. ACS Omega3, 18623–18629. 10.1021/acsomega.8b03020
16
ChenY.WuH.LiZ.WangP.YangL.FangY. (2012). The study of surface plasmon in au/ag core/shell compound nanoparticles. Plasmonics7, 509–513. 10.1007/s11468-012-9336-6
17
ChenX.MooreJ. E.ZekariasM.JensenL. (2015). Atomistic electrodynamics simulations of bare and ligand-coated nanoparticles in the quantum size regime. Nat. Commun.6, 8921. 10.1038/ncomms9921
18
CiraciC.Della SalaF. (2016). Quantum hydrodynamic theory for plasmonics: Impact of the electron density tail. Phys. Rev. B93, 205405. 10.1103/physrevb.93.205405
19
CiracìC.HillR.MockJ.UrzhumovY.Fernández-DomínguezA.MaierS.et al (2012). Probing the ultimate limits of plasmonic enhancement. Science337, 1072–1074. 10.1126/science.1224823
20
CiraciC.PendryJ. B.SmithD. R. (2013). Hydrodynamic model for plasmonics: A macroscopic approach to a microscopic problem. ChemPhysChem14, 1109–1116. 10.1002/cphc.201200992
21
CocciaE.FregoniJ.GuidoC.MarsiliM.PipoloS.CorniS. (2020). Hybrid theoretical models for molecular nanoplasmonics. J. Chem. Phys.153, 200901. 10.1063/5.0027935
22
CorniS.TomasiJ. (2001). Enhanced response properties of a chromophore physisorbed on a metal particle. J. Chem. Phys.114, 3739–3751. 10.1063/1.1342241
23
CorniS.TomasiJ. (2002a). Excitation energies of a molecule close to a metal surface. J. Chem. Phys.117, 7266–7278. 10.1063/1.1507579
24
CorniS.TomasiJ. (2002b). Surface enhanced Raman scattering from a single molecule adsorbed on a metal particle aggregate: A theoretical study. J. Chem. Phys.116, 1156–1164. 10.1063/1.1428349
25
DanielisN.VegaL.FronzoniG.StenerM.BruixA.NeymanK. M. (2021). Agpd, aupd, and aupt nanoalloys with ag-or au-rich compositions: Modeling chemical ordering and optical properties. J. Phys. Chem. C125, 17372–17384. 10.1021/acs.jpcc.1c04222
26
Della SalaF.PachterR.SukharevM. (2022). Advances in modeling plasmonic systems. J. Chem. Phys.157, 190401. 10.1063/5.0130790
27
Della SalaF. (2022). Orbital-free methods for plasmonics: Linear response. J. Chem. Phys.157, 104101. 10.1063/5.0100797
28
DraineB. T.FlatauP. J. (1994). Discrete-dipole approximation for scattering calculations. J. Opt. Soc. Am. A11, 1491–1499. 10.1364/josaa.11.001491
29
EstebanR.BorisovA. G.NordlanderP.AizpuruaJ. (2012). Bridging quantum and classical plasmonics with a quantum-corrected model. Nat. Commun.3, 825. 10.1038/ncomms1806
30
EstebanR.ZugarramurdiA.ZhangP.NordlanderP.García-VidalF. J.BorisovA. G.et al (2015). A classical treatment of optical tunneling in plasmonic gaps: Extending the quantum corrected model to practical situations. Faraday Discuss.178, 151–183. 10.1039/c4fd00196f
31
EtchegoinP. G.Le RuE.MeyerM. (2006). An analytic model for the optical properties of gold. J. Chem. Phys.125, 164705. 10.1063/1.2360270
32
GaoC.HuY.WangM.ChiM.YinY. (2014). Fully alloyed ag/au nanospheres: Combining the plasmonic property of ag with the stability of au. J. Am. Chem. Soc.136, 7474–7479. 10.1021/ja502890c
33
García de AbajoF. J.HowieA. (2002). Retarded field calculation of electron energy loss in inhomogeneous dielectrics. Phys. Rev. B65, 115418. 10.1103/physrevb.65.115418
34
GianniniV.Fernández-DomínguezA. I.HeckS. C.MaierS. A. (2011). Plasmonic nanoantennas: Fundamentals and their use in controlling the radiative properties of nanoemitters. Chem. Rev.111, 3888–3912. 10.1021/cr1002672
35
GiannoneG.SmigaS.D’AgostinoS.FabianoE.Della SalaF. (2021). Plasmon couplings from subsystem time-dependent density functional theory. J. Phys. Chem. A125, 7246–7259. 10.1021/acs.jpca.1c05384
36
GiovanniniT.PuglisiA.AmbrosettiM.CappelliC. (2019a). Polarizable qm/mm approach with fluctuating charges and fluctuating dipoles: The qm/fqfμ model. J. Chem. Theory Comput.15, 2233–2245. 10.1021/acs.jctc.8b01149
37
GiovanniniT.RosaM.CorniS.CappelliC. (2019b). A classical picture of subnanometer junctions: An atomistic drude approach to nanoplasmonics. Nanoscale11, 6004–6015. 10.1039/c8nr09134j
38
GiovanniniT.BonattiL.PoliniM.CappelliC. (2020). Graphene plasmonics: Fully atomistic approach for realistic structures. J. Phys. Chem. Lett.11, 7595–7602. 10.1021/acs.jpclett.0c02051
39
GiovanniniT.BonattiL.LafioscaP.NicoliL.CastagnolaM.IllobreP. G.et al (2022). Do we really need quantum mechanics to describe plasmonic properties of metal nanostructures?ACS Photonics9, 3025–3034. 10.1021/acsphotonics.2c00761
40
GongC.LeiteM. S. (2016). Noble metal alloys for plasmonics. ACS Photonics3, 507–513. 10.1021/acsphotonics.5b00586
41
HalasN. J.LalS.ChangW.-S.LinkS.NordlanderP. (2011). Plasmons in strongly coupled metallic nanostructures. Chem. Rev.111, 3913–3961. 10.1021/cr200061k
42
HaoF.NehlC. L.HafnerJ. H.NordlanderP. (2007). Plasmon resonances of a gold nanostar. Nano Lett.7, 729–732. 10.1021/nl062969c
43
HaynesW. (2014). CRC handbook of chemistry and physics. Boca Raton: CRC Press.
44
HohenesterU.TrüglerA. (2012). MNPBEM – a Matlab toolbox for the simulation of plasmonic nanoparticles. Comput. Phys. Commun.183, 370–381. 10.1016/j.cpc.2011.09.009
45
HohenesterU. (2015). Quantum corrected model for plasmonic nanoparticles: A boundary element method implementation. Phys. Rev. B91, 205436. 10.1103/physrevb.91.205436
46
HuangJ.ZhuY.LiuC.ZhaoY.LiuZ.HedhiliM. N.et al (2015). Fabricating a homogeneously alloyed auag shell on au nanorods to achieve strong, stable, and tunable surface plasmon resonances. Small11, 5214–5221. 10.1002/smll.201501220
47
JensenL. L.JensenL. (2008). Electrostatic interaction model for the calculation of the polarizability of large noble metal nanoclusters. J. Phys. Chem. C112, 15697–15703. 10.1021/jp804116z
48
JensenL. L.JensenL. (2009). Atomistic electrodynamics model for optical properties of silver nanoclusters. J. Phys. Chem. C113, 15182–15190. 10.1021/jp904956f
49
KneippK.WangY.KneippH.PerelmanL. T.ItzkanI.DasariR. R.et al (1997). Single molecule detection using surface-enhanced Raman scattering (sers). Phys. Rev. Lett.78, 1667–1670. 10.1103/physrevlett.78.1667
50
KreibigU.VollmerM.KreibigU.VollmerM. (1995). “Theoretical considerations,” in Optical properties of metal clusters, 13–201.
51
KuddahM.PutraM.DjuhanaD. (2020). “The incident electrical field angle effect in localized surface plasmon resonance (lspr) of bimetallic ag-au nanorod using mnpbem simulation,” in IOP Conference Series: Materials Science and Engineering (IOP Publishing), 012058.
52
LafioscaP.GiovanniniT.BenziM.CappelliC. (2021). Going beyond the limits of classical atomistic modeling of plasmonic nanostructures. J. Phys. Chem. C125, 23848–23863. 10.1021/acs.jpcc.1c04716
53
LafioscaP.NicoliL.BonattiL.GiovanniniT.CorniS.CappelliC. (2023). Qm/classical modeling of surface enhanced Raman scattering based on atomistic electromagnetic models. J. Chem. Theory Comput.10.1021/acs.jctc.3c00177
54
LarsenA. H.MortensenJ. J.BlomqvistJ.CastelliI. E.ChristensenR.DułakM.et al (2017). The atomic simulation environment—A python library for working with atoms. J. Phys. Condens. Matter29, 273002. 10.1088/1361-648x/aa680e
55
LiebschA. (1993). Surface-plasmon dispersion and size dependence of mie resonance: Silver versus simple metals. Phys. Rev. B48, 11317–11328. 10.1103/physrevb.48.11317
56
LimD.-K.JeonK.-S.KimH. M.NamJ.-M.SuhY. D. (2010). Nanogap-engineerable Raman-active nanodumbbells for single-molecule detection. Nat. Mater.9, 60–67. 10.1038/nmat2596
57
LinkS.WangZ. L.El-SayedM. (1999). Alloy formation of gold-silver nanoparticles and the dependence of the plasmon absorption on their composition. J. Phys. Chem. B103, 3529–3533. 10.1021/jp990387w
58
López LozanoX.MottetC.WeisskerH.-C. (2013). Effect of alloying on the optical properties of ag–au nanoparticles. J. Phys. Chem. C117, 3062–3068. 10.1021/jp309957y
59
MaY.-W.ZhangL.-H.WuZ.-W.YiM.-F.ZhangJ.JianG.-S. (2015). The study of tunable local surface plasmon resonances on au-ag and ag-au core-shell alloy nanostructure particles with dda method. Plasmonics10, 1791–1800. 10.1007/s11468-015-9997-z
60
MaJ.LiuX.WangR.ZhangJ.JiangP.WangY.et al (2020). Bimetallic core–shell nanostars with tunable surface plasmon resonance for surface-enhanced Raman scattering. ACS Appl. Nano Mater.3, 10885–10894. 10.1021/acsanm.0c02144
61
MaierS. A.KikP. G.AtwaterH. A.MeltzerS.HarelE.KoelB. E.et al (2003). Local detection of electromagnetic energy transport below the diffraction limit in metal nanoparticle plasmon waveguides. Nat. Mater.2, 229–232. 10.1038/nmat852
62
MaierS. A. (2007). Plasmonics: Fundamentals and applications. Springer Science & Business Media.
63
MieG. (1908). Beiträge zur optik trüber medien, speziell kolloidaler metallösungen. Ann. Phys.-Berlin330, 377–445. 10.1002/andp.19083300302
64
MoskovitsM. (1985). Surface-enhanced spectroscopy. Rev. Mod. Phys.57, 783–826. 10.1103/revmodphys.57.783
65
MuehlschlegelP.EislerH.-J.MartinO. J.HechtB.PohlD. (2005). Resonant optical antennas. Science308, 1607–1609. 10.1126/science.1111886
66
MyroshnychenkoV.Rodríguez-FernándezJ.Pastoriza-SantosI.FunstonA. M.NovoC.MulvaneyP.et al (2008). Modelling the optical response of gold nanoparticles. Chem. Soc. Rev.37, 1792–1805. 10.1039/b711486a
67
NeumanT.EstebanR.CasanovaD.García-VidalF. J.AizpuruaJ. (2018). Coupling of molecular emitters and plasmonic cavities beyond the point-dipole approximation. Nano Lett.18, 2358–2364. 10.1021/acs.nanolett.7b05297
68
NewmaiM. B.VermaM.DahiyaA.KumarP. S. (2022). Monomer driven growth of catalytically active agau plasmonic nanoalloys. J. Phys. Chem. Solids161, 110371. 10.1016/j.jpcs.2021.110371
69
NieS.EmoryS. R. (1997). Probing single molecules and single nanoparticles by surface-enhanced Raman scattering. Science275, 1102–1106. 10.1126/science.275.5303.1102
70
OlobardiS.VegaL.FortunelliA.StenerM.VinesF.NeymanK. M. (2019). Optical properties and chemical ordering of ag–pt nanoalloys: A computational study. J. Phys. Chem. C123, 25482–25491. 10.1021/acs.jpcc.9b07382
71
PapavassiliouG. C. (1976). Surface plasmons in small au-ag alloy particles. J. Phys. F. Met. Phys.6, L103–L105. 10.1088/0305-4608/6/4/004
72
PaytonJ. L.MortonS. M.MooreJ. E.JensenL. (2013). A hybrid atomistic electrodynamics–quantum mechanical approach for simulating surface-enhanced Raman scattering. Acc. Chem. Res.47, 88–99. 10.1021/ar400075r
73
Peña-RodríguezO.PalU. (2011a). Au@ ag core–shell nanoparticles: Efficient all-plasmonic fano-resonance generators. Nanoscale3, 3609–3612. 10.1039/c1nr10625b
74
Peña-RodríguezO.PalU. (2011b). Enhanced plasmonic behavior of bimetallic (ag-au) multilayered spheres. Nanoscale Res. Lett.6, 279. 10.1186/1556-276X-6-279
75
Pérez-GonzálezO.ZabalaN.BorisovA.HalasN.NordlanderP.AizpuruaJ. (2010). Optical spectroscopy of conductive junctions in plasmonic cavities. Nano Lett.10, 3090–3095. 10.1021/nl1017173
76
PutraM.DjuhanaD.FauziaV.HarmokoA.ImawanC. (2017). “A numerical study of the sensitivity of surface plasmon resonance bimetallic silver-gold alloys using boundary element method,” in IOP Conference Series: Materials Science and Engineering (IOP Publishing), 012016.
77
RazaS.BozhevolnyiS. I.WubsM.MortensenN. A. (2015). Nonlocal optical response in metallic nanostructures. J. Phys. Condens. Matter27, 183204. 10.1088/0953-8984/27/18/183204
78
RiouxD.VallièresS.BesnerS.MuñozP.MazurE.MeunierM. (2014). An analytic model for the dielectric function of au, ag, and their alloys. Adv. Opt. Mater.2, 176–182. 10.1002/adom.201300457
79
RistigS.PrymakO.LozaK.GocylaM.Meyer-ZaikaW.HeggenM.et al (2015). Nanostructure of wet-chemically prepared, polymer-stabilized silver–gold nanoalloys (6 nm) over the entire composition range. J. Mater. Chem. B3, 4654–4662. 10.1039/c5tb00644a
80
SanthoshK.BittonO.ChuntonovL.HaranG. (2016). Vacuum rabi splitting in a plasmonic cavity at the single quantum emitter limit. Nat. Commun.7, 11823. 10.1038/ncomms11823
81
SchollJ. A.García-EtxarriA.KohA. L.DionneJ. A. (2013). Observation of quantum tunneling between two plasmonic nanoparticles. Nano Lett.13, 564–569. 10.1021/nl304078v
82
ShufordK. L.RatnerM. A.GrayS. K.SchatzG. C. (2006). Finite-difference time-domain studies of light transmission through nanohole structures. Appl. Phys. B84, 11–18. 10.1007/s00340-006-2218-x
83
SørensenL. K.UtyushevA. D.ZakomirnyiV. I.ÅgrenH. (2021). Atomistic description of plasmonic generation in alloys and core shell nanoparticles. Phys. Chem. Chem. Phys.23, 173–185. 10.1039/d0cp04854b
84
SzántóG.CsarnovicsI.BonyárA. (2021). Numerical investigation of the refractive index sensitivity of au/ag core-shell nanostructures for sensing applications. Sens. Bio-Sens.32, 100414. 10.1016/j.sbsr.2021.100414
85
VerbruggenS. W.KeulemansM.MartensJ. A.LenaertsS. (2013). Predicting the surface plasmon resonance wavelength of gold–silver alloy nanoparticles. J. Phys. Chem. C117, 19142–19145. 10.1021/jp4070856
86
WangC.PengS.ChanR.SunS. (2009). Synthesis of auag alloy nanoparticles from core/shell-structured ag/au. Small5, 567–570. 10.1002/smll.200801169
87
XiangY.WuX.LiuD.LiZ.ChuW.FengL.et al (2008). Gold nanorod-seeded growth of silver nanostructures: From homogeneous coating to anisotropic coating. Langmuir24, 3465–3470. 10.1021/la702999c
88
ZakomirnyiV. I.RinkeviciusZ.BaryshnikovG. V.SørensenL. K.ÅgrenH. (2019). Extended discrete interaction model: Plasmonic excitations of silver nanoparticles. J. Phys. Chem. C123, 28867–28880. 10.1021/acs.jpcc.9b07410
89
ZakomirnyiV. I.RasskazovI. L.SørensenL. K.CarneyP. S.RinkeviciusZ.ÅgrenH. (2020). Plasmonic nano-shells: Atomistic discrete interaction versus classic electrodynamics models. Phys. Chem. Chem. Phys.22, 13467–13473. 10.1039/d0cp02248a
Summary
Keywords
atomistic, alloys, core-shell, plasmonics, gold, silver, bimetallic
Citation
Nicoli L, Lafiosca P, Grobas Illobre P, Bonatti L, Giovannini T and Cappelli C (2023) Fully atomistic modeling of plasmonic bimetallic nanoparticles: nanoalloys and core-shell systems. Front. Photonics 4:1199598. doi: 10.3389/fphot.2023.1199598
Received
03 April 2023
Accepted
07 June 2023
Published
16 June 2023
Volume
4 - 2023
Edited by
Stefania D'Agostino, National Research Council (CNR), Italy
Reviewed by
Cristian Ciracì, Italian Institute of Technology (IIT), Italy
Maxim A. Yurkin, Institute of Chemical Kinetics and Combustion (RAS), Russia
Updates
Copyright
© 2023 Nicoli, Lafiosca, Grobas Illobre, Bonatti, Giovannini and Cappelli.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tommaso Giovannini, tommaso.giovannini@sns.it; Chiara Cappelli, chiara.cappelli@sns.it
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.