1. Introduction
The field of plasmonic photocatalysis has experienced a great development in the recent years [
1,
2,
3] as a clean technique for solar energy conversion. It has arisen as an alternative to conventional semiconductor catalysis in order to overcome some of its drawbacks [
2,
3]. Semiconductor photocatalysis uses semiconductors to absorb light and create electron-hole pairs (
e-h) that can initiate reduction and/or oxidation chemical reactions. However, during this process, two main difficulties arise: low photocatalytic efficiency due to the
e-h recombination, and the lack of high-performance photocatalytic materials with band gaps that can absorb visible light. For instance, high performing semiconductor photocatalysts, like TiO
or CeO
, absorb light in the near-ultraviolet hampering solar-light driven photocatalytic reactions. The inclusion of metallic nanoparticles into the semiconductor enable new paths to circumvent these problems. At the metal-semiconductor interface, a Schottky junction appears, creating an internal electric field that forces both
e and
h into different directions inhibiting their recombination. In addition, metallic nanoparticles (NPs) support localized surface plasmon resonances (LSPRs) that also bring some benefits. Although common metals used in plasmonic photocatalysis such as gold (Au) or silver (Ag) present LSPR visibly, they are able to drive a photocatalytic response of large-band gap semiconductors (e.g., TiO
or CeO
) through a charge-transfer mechanism. It consists of the transfer of photo-excited electrons from the metal into the conduction band of the semiconductor [
2,
3]. LSPRs also create very strong electromagnetic fields confined in subwavelength dimensions at the vicinity of the nanoparticles, which can increase the
e-h pair generation [
2]. Moreover, the excitation of LSPRs in metals inherently produces heat due to Joule losses. This temperature increment increases redox reaction rates, mass transfer and polarizes nonpolar molecules for better absorption.
A lot of effort has been invested in the synthesis and design of metal-semiconductor photocatalysts. This includes choosing both the semiconductor and metal, and the configuration of the systems (size, shape, interface design, ...) [
2,
4]. Nanostructuring the metal-semiconductor system into nanoparticles has been widely studied as a way to boost the photocatalytic performance by increasing the surface-to-volume ratios. The synthesis of large variety of metal-semiconductor NPs has been reported in the literature e.g., metal-semiconductor core-shell NPs or semiconductor NPs with their surfaces decorated with metal inclusions [
4,
5]. Specifically, core-shell nanostruxctures present special advantages for photocatalytic reactions such as protecting the metal and enlarging the active sites [
4]. These NP systems are usually supported in a host matrix or solvent. To maximize solar absorption, it is critical to accurately model the reflectance/absorptance/transmittance of the host-matrix with metal-semiconductor NPs as inclusions. In principle, these systems, assuming that the inclusions behave as electric dipoles, could be modeled using, for example, classical Maxwell Garnett effective medium theory (EMT) [
6]. Maxwell Garnett EMT allows the calculation of the effective dielectric function of the host matrix with metal-semiconductor NPs inclusions taking as inputs the dielectric function of both the host and the inclusions. At this point, the problem arises regarding:
which is the dielectric function of the metal-semiconductor inclusion? In this contribution, we propose a new procedure to calculate the effective dielectric function of metal-semiconductor core-shell NPs. In this way, we enable the use of Maxwell Garnett EMT or equivalent EMTs for the modeling of state-of-the-art plasmonic photocatalytic systems in which metal-semiconductor core-shell photocatalysts are supported in a matrix of arbitrary dielectric function [
7,
8]. Note that similar configurations are also widely used in solar energy harvesting applications [
9,
10,
11]. This procedure is illustrated in
Figure 1: the proposed approach is used to calculate the effective dielectric function of the core-shell inclusions
. Then, the problem is reduced to the one of a matrix with spherical inclusions randomly distributed. This system is commonly solved by applying Maxwell Garnett EMT or equivalent EMTs; taking as input the effective dielectric function of the core-shell
and the host matrix
, the effective dielectric function of the whole system
is calculated. Therefore, the initial system (i.e., host material with core-shell inclusions) is reduced to a homogeneous medium with dielectric function
.
The paper is organized as follows: first, we perform a detailed study on the effective medium theories used by other authors in the literature to model core-shell NPs. Secondly, we present the basis of the new EMT for core-shell NPs developed in this work. In the following section, and in order to validate our method, we compare the exact extinction efficiencies of different metal-semiconductor core-shell NPs calculated through Mie theory with that obtained by using the effective dielectric functions in Mie theory for a homogeneous sphere. Afterwards, we give an example about the applicability of the proposed method developed here: the calculation of the reflectance of a composite with core-shell type inclusions. Finally, the main conclusions of this research are presented.
2. Effective Medium Theories
Although the electromagnetic response of spherical core-shell NPs has an analytical solution given by Mie theory [
12], it is not trivial to obtain an effective dielectric function of this system. Mie theory for a coated sphere gives exact analytical expressions for the absorption/scattering/extinction efficiencies and scattered electric and magnetic fields. However, it does not provide an expression for the effective dielectric function of a coated sphere. Therefore, effective medium theories (EMTs) give an alternative to deal with this problem. EMTs allow us to obtain the dielectric function of an equivalent sphere with the same size and electromagnetic behavior as a given spherical core-shell NP. This process is illustrated in
Figure 2.
Although in the literature several EMTs have been used to model the electromagnetic response of core-shell NPs [
11,
13,
14,
15,
16,
17], here we present a new approach that overcomes some of the limitations of the existing EMTs. Specifically, we compare the proposed approach with the weighted average, Bruggeman, Maxwell Garnett, Internal Homogenization and Mie theory based Maxwell Garnett.
In the following subsections, we will review the aforemention EMTs, their limitations and their applicability to the proposed problem: calculating the effective dielectric function of a metal-semiconductor core-shell spherical nanoparticles.
2.1. Weighted Average (WA)
One of the simplest approaches to get the effective dielectric function (
) of a core-shell nanoparticle is given by the weighted average of the dielectric functions of both the core (
) and the shell (
),
where
is the fraction of volume that occupies the core. When
(core-material particle),
. Conversely, when
(shell-material particle),
. This approach has been used by Kuzma et al. [
13] to model the effective dielectric function of Ag/Ag
O core-shell nano-spheres.
2.2. Bruggeman
Bruggeman’s effective medium theory [
12,
18] is proposed for non-homegenous media in which two materials are present, the main one—the matrix—hosting the secondary one—the inclusion. Both materials, however, are given the same importance, i.e., it is invariant to an interchange between the host matrix and the inclusions. In this approximation, the effective dielectric function
is given by
where
and
are the dielectric functions of the host (in our case the shell) and the inclusions (the core), respectively.
f is the filling factor of the inclusions. If in Equation (
2), we now consider the filling fraction of the host (
) instead of the one of the inclusion (
), we obtain an equivalent expression. This means that which component is labelled as inclusion or host does not matter, both are interchangeable. Bruggeman effective medium theory has been used by Knight et al. [
14] to approximate the effective dielectric function of oxidized Al disks.
2.3. Maxwell Garnett (MG)
The Maxwell Garnett theory [
6,
12] considers the case of inclusions randomly spread in a continuous matrix. This theory is based on the Clausius–Mossotti relation, which relates the polarizability (
) with the dielectric function (
). The effective dielectric function of a suspension of small spheres in a host material is given by
where
and
are the dielectric functions of host and inclusions materials, and
f is the volume fraction of the embedded particles. Because a core-shell nanoparticle can be considered as an inclusion (core with dielectric function
) in a host medium (shell with dielectric function
) where
and
, its effective dielectric function can be expressed as
Again,
f is the filling factor that in the case of a spherical nano-shell is given by
. Notice that when
(core-material particle)
. On the contrary, when
(shell-material particle)
. Many authors have used this EMT to model the response of core-shell NPs [
11,
16,
17].
2.4. Internal Homogenization (IH)
A more complex approach was introduced by Chettiar and Enghetta [
15] through the concept of internal homogenization. Through this process the polarizability of the equivalent sphere is equated to that of a core-shell in the electrostatic approximation [
12]. According to this, the effective dielectric function for a core-shell NP in vacuum is given by
where
f is the previously defined filling factor.
By rearranging Equation (
5), the effective dielectric function can be expressed as the one given by an MG EMT, Equation (
4), assuming a core embedded in a shell material medium with a filling factor of
f. This equation that is independent of the particle size is only valid in the regime in which the size of the particle is much smaller than the incident wavelength (electrostatic approximation).
2.5. Mie Theory Based Maxwell Garnett (MMG)
A size dependent extension of the MG model, Equation (
3), was proposed by Doyle [
19,
20] for a suspension of metallic spheres. In this case, the polarizability of the core embedded in shell material is given in terms of the Mie coefficient
[
12] rather than the Clausius–Mossoti relation,
This approach considers the size effects on the polarizability of the inclusion—core—in the host material—shell. In addition, by considering the polarizability of the core embedded in shell material, the validity of this model is restricted to the limit .
3. New Approach
All the previous EMTs, already used in the literature, are either limited by the size of the particles (electrostatic approximation) or by the condition . In this section, we present a new model aimed to overcome these limitations by considering the exact dipolar polarizability of the core-shell sphere in terms of the dipolar electric coefficient from the Mie expanssion , and introducing radiating and depolarization effects in the polarizability of the effective sphere.
As stated above, we will first consider the exact polarizability of the core-shell nanoparticle (
) in terms of its Mie dipolar electric coefficient
[
21],
By doing this, size effects are being considered and the
condition is removed.
On the other hand, size effects in the effective dielectric function of the equivalent sphere can be taken into account by introducing dynamic depolarization and radiation damping factors [
22]:
where
x is the size paramenter given by
being
R the size of the particle
. The depolarization factor considered here,
, is the one derived by Meier and Wokaun [
22] by using a self-consistent derivation of the polarization of a sphere. It is worth mentioning that other expressions can be found for this factor in the literature calculated from a power-series expansion of the Mie coefficient
. For instance, Meier and Wokaun using a power expansion of the Mie coefficient obtained a depolarization factor of
. Other authors like Schebarchov et al. [
23] and Tzarouchis et al. [
24] through a Taylor and Padé expansions, respectively, have expressed this term scaling
. Regarding the radiation damping,
, all authors have reported the same expression.
Thus, by considering the exact polarizability of the core-shell, Equation (
7), and the depolarization and radiation damping factors in the polarizability of the equivalent sphere, Equation (
8), we can express the effective dielectric function of an spherical nano-shell as
This approach is more of an internal homogenization process [
15] than a Maxwell Garnett effective medium theory [
19]. This new approach is a step forward in the methodology proposed by Chettiar and Engheta [
15] given that our approach considers size effects in the polarizability of both effective sphere and coated sphere. For these reasons, from now on, we will refer to this proposed theory as Extended Internal Homogenization (Ext. IH).
4. Results
To evaluate the validity of the EMTs, we will compare the extinction efficiency spectra (
) of metal-semiconductor core-shell spherical nanoparticles calculated with Mie theory for a coated sphere [
12] (exact solution) with that obtained by using the effective dielectric function calculated through each of the aforementioned EMTs in Mie theory for a homogeneous sphere [
12]. For this study, we have chosen as plasmonic metal Au and three common semiconductors widely used in semiconductor photocatalysis: Cu
O, CeO
and TiO
. A broad variety of Au-core Cu
O/CeO
/TiO
-shell nanoparticles have been reported in the literature for different photocatalytic applications [
5,
25,
26]. A clear example of the utility of these core-shell nanostructures is the use Au-TiO
core-shell NPs in the water-splitting reaction. The efficient performance of this chemical reaction, in which water molecules dissociated into oxygen and hydrogen (2H
O → 2H
+ O
), would have a great impact in the hydrogen economy [
27]. Whereas TiO
NPs can drive the water-splitting reaction by themselves when illuminated with UV light, when using Au-TiO
, the reaction can be driven by visible light, for instance that coming from the sun [
28,
29].
Figure 3 shows the dielectric function (
) of the chosen materials. The optical constants have been taken from different sources of the literature: Au [
30], Cu
O [
31], CeO
[
32] and TiO
[
33].
Figure 4,
Figure 5 and
Figure 6 show the efficiency spectra (
) of Au-CeO
, Au-Cu
O and Au-TiO
core-shell NPs with external radius
10, 20 and 30 nm and varying core radii
. These sizes have chosen in order to work in the dipolar regime. The exact
spectra calculated with Mie theory for a coated sphere are compared with those obtained by using the effective dielectric functions in Mie theory for a homogeneous sphere. The considered effective dielectric functions have been calculated with the following EMTs: weighted average (WA), Bruggeman (Brugg), Maxwell Garnett (MG) (which leads to the same expression as the IH proposed by Chettiar–Engheta), Mie theory based Maxwell Garnett (MMG) and the proposed approach (Ext. IH).
Figure 4,
Figure 5 and
Figure 6 show that Maxwell Garnett (MG) and the proposed method (Ext. IH) give very similar results and both are the ones that better reproduce the exact solution given by Mie theory for a coated sphere. Conversely, Bruggeman EMT (Brug) is the one that worst describes the plasmonic response of the core-shell NPs. The reason arises from the basis of this theory i.e., both matrix and inclusions materials are interchangeable. In the system under study, interchanging the core and the shell materials lead to systems which, in principle, have a completely different plasmonic response. Those effective dielectric functions calculated by a weighted average EMT (WA) do not reproduce the exact solution. Although a plasmonic peak appears for the largest
, it is considerably blue-shifted and higher than those of the exact solution. The origin of these discrepancies comes from the construction of (
1). When
(almost 100% shell material in the NP) and
(almost 100% core material in the NP), it is expected to reach the convergence with the exact solution since
and
. However, for the intermediate values of
, and considering the high contrast of
(metallic material) and
(dielectric material), the obtained average is far from that of the considered material because none of the underlying physics of the system are included in Equation (
1). In the case of the Mie theory based Maxwell Garnett (MMG), although it presents a good agreement with small values of
, it fails to reproduce the values of
. All these problems emerge from the basis of this model. It considers the core embedded in an infinite matrix of shell material. When the metallic core is small, under illumination, the evanescent field exited by the LSPRs remains attached to it and dissipates within the dielectric shell. Therefore, the core can be considered to be completely embedded in shell material. However, as the shell thickness decreases, part of the evanescent field reaches beyond the outer shell surface. It is in this regime where MMG fails because the core cannot be considered to be hosted longer in an infinite shell [
34]. In general, for all the EMTs under study, the predictions are worse as the size of the core-shell NP (
) increases because quadrupolar effects start to emerge. For the same reason, the predictions are better for the smallest core sizes.
Since MG and Ext. IH EMTs fit the exact solution best, we have performed a more systematic comparison of the performance of both EMTs. For that reason,
Figure 7 shows the spectral shift (
, red right axis) of the
peak predicted through MG and Ext. IH EMTs with respect to the exact solution as a function of the core radius
. Indeed, the accurate prediction of the spectral position of the localized surface plasmon resonances is crucial in the design of core-shell nanoparticles for their application in photocatalysis. The precise selection of the resonant wavelength will determine the efficiency in which the catalyst absorbs light and the coupling of the enhanced electromagnetic field with the band gap of the semiconductor shell. Therefore, low values of
are desirable. In addition, in
Figure 7, the predicted values of the
peaks (blue left axis) along with the exact value are plotted. Again, Au-CeO
(top row), Au-Cu
O (middle row) and Au-TiO
core-shell NPs (bottom row) are considered. In this case, the total radii of the NPs are
nm (left column),
nm (middle column) and
nm (right column). The calculations have been performed for core sizes ranging in
nm.
Concerning the values of the peaks, the difference between the exact value and the predicted one increases as both the core and total radii become larger. Specifically, the predicted values are lower than the exact ones, those being given by Ext. IH closer to the exact solution. The only exception is produced for the largest of the and 35 nm particles, where the value of MG is closer to that of Ext. IH. This is due to the fact that, for , the effective dielectric function approaches the exact solution. This issue also appears when comparing the value of of both EMTs. In the case of MG, it can be seen how, for and , the values are very close to that of the exact solution while the discrepancy is larger for intermediate values of . On the contrary, for the Ext. IH, monotonically increases with the value of . However, the values of are smaller than for MG in all the considered cases. For example, for and 25 nm is always smaller than ≈ 1 nm while MG shows values of up to ≈ 3 and 7.5 nm, respectively. For the biggest NP size and the largest , both MG and Ext. IH converge. However, for smaller values of , nm for Ext. IH while for MG reaches up to 15 nm.
Therefore, taking into account all these results, we can conclude that the proposed approach shows a better performance in predicting the electromagnetic behavior of metal-semiconductor core-shells than conventional effective medium theories such as weighted average [
13], Bruggeman [
18], Maxwell Garnett [
6] and Mie theory based Maxwell Garnett [
20].
5. Application: Reflectance of a Core-Shell Inclusion Type Composite
As an example of the applicability of the results in this research, here we will compare the reflectance at normal incidence of a system consisting of a homogeneous slab with metal-semiconductor core-shell inclusions calculated using different methods. The slab (i.e., host medium) has been considered to have a thickness of
nm and a dielectric function
(water). The core-shell inclusions have an Au core of radius
= 10 nm, and the shell has a thickness of 10 nm and dielectric function
(TiO
). The filling fraction of inclusions in the slab is
. A scheme of this system is sketched in
Figure 8a. Note that there is not an analytical solution to calculate the reflectance of these types of systems.
The reflectances have been calculated using three different methodologies:
Using FDTD Lumerical software (Lumerical Inc., Vancouver, Canada) (numerical solution of Maxwell’s equations), we have built a unit cell periodically replicated in the plane perpendicular to the wave vector of the illuminating beam (
k in
Figure 8a) that reproduces the characteristics of the system under study.
Using the Transfer Matrix Method [
35] taking as inputs the effective dielectric constant of the slab with the core-shell inclusions and its thickness (
nm). The procedure followed is illustrated in
Figure 9. (
i) The initial system is reduced to one in which the core-shell inclusions are transformed into homogeneous spheres with dielectric function
. The dielectric function of these homogeneous spheres (
) has been calculated following the proposed approach (Equation (
9)). (
) Using Maxwell Garnett, taking as inputs the dielectric function of the host medium (water,
), the effective dielectric function of the inclusions calculated using Equation (
9) (
), and the filling fraction (
), we have calculated the effective medium of the whole system
. Finally, using the Transfer Matrix Method (TMM) [
35], we have calculated the reflectance of a homogeneous slab with dielectric function
and thickness
nm at normal incidence.
Following the same recipe of the previous point, but, instead of using the proposed approach to calculate the effective dielectric function of the core-shell with Equation (
9), we have used Maxwell Garnett (Equation (
4)), a widely used approach to model core-shell NPs [
11,
16,
17]. This EMT has shown the second best performance in our study. It is worth mentioning that this equation is equivalent to that proposed by Chettiar–Enghetta for the internal homogenization of coated spheres (Equation (
5)) [
15].
The result of this comparison is shown in
Figure 8b. Yellow dots represent the numerical solution given by the FDTD simulations. The reflectances obtained by TMM using the proposed approach and Maxwell Garnett to calculate the effective dielectric function of the inclusions are represented in blue and red, respectively. The numerical solution presents a characteristic dip at 663 nm. While using Maxwell Garnett to calculate the effective dielectric function of the core-shell completely fails to reproduce this dip, our approach clearly reproduces it. More precisely, in the latter case, the dip is produced at 657 nm (showing a 5 nm deviation from the numerical solution). The discrepancy between the exact solution and that derived by using our approach to calculate the effective dielectric function of the inclusion as input for the Maxwell Garnett of the composite comes from the fact that Maxwell-Garnett assumes inclusions in the dipolar regime. As the size of the inclusions increases, and high order multipoles start to contribute, the inaccuracy of the effective dielectric function also increases. This size effects are enhanced when the particles are surrounded by a medium denser than vacuum since the wavelength inside the host matrix is shorter, making the particle effectively bigger. Therefore, by using our proposed approach, we can get a very good estimation of the reflectance on this type of systems without the need for using advanced electromagnetic simulations software.