enhanced absorption and photoluminescence from dye .../_publications/...2018/10/15 · of the film...

Enhanced absorption and photoluminescence from dye-containing thin polymer film on plasmonic array

SHUNSUKE MURAI,1,2 SAHO OKA,1 SHAIMAA I. AZZAM,3 ALEXANDER V. KILDISHEV,3 SATOSHI ISHII,4 AND KATSUHISA TANAKA

1 1Department of Material Chemistry, Graduate School of Engineering, Kyoto University, Nishikyo-ku, Kyoto, 615-8510, Japan 2PRESTO, Japan Science and Technology Agency (JST), Kawaguchi, Saitama 332-0012, Japan 3School of Electrical and Computer Engineering, Birck Nanotechnology Center and Purdue Quantum Center, Purdue University, West Lafayette, IN 47907, USA 4WPI International Center for Material Nanoarchitectonics, National Institute for Materials Science, Namiki, Tsukuba, Ibaraki 305-0044, Japan *[email protected]

Abstract: Thin films containing light emitters act as light-to-light converters that absorb the incident light and emit luminescence. This well-known phenomenon is photoluminescence (PL). When a photoluminescent film is notably thinner than the absorption length of emitters, it exhibits weak absorption of incident light. The absorption can be increased by depositing the thin film on a plasmonic array of metallic nanocylinders arranged with a specific periodicity. The array couples the incident light into the thin film, facilitating the plasmon-enhanced absorption by the emitters in the film. In this study, we demonstrate both experimentally and numerically the plasmon-enhanced absorption of a rhodamine 6G-containing film that is thinner than its absorption length using a periodic array of Al nanocylinders. The experimental results demonstrate that the spectrally integrated PL intensity is increased up to 3.78 times. In addition to enhanced absorption, the array is also found to diffract the PL into a direction determined by the periodicity, thereby facilitating the multiplied enhancement of PL. The combination of the two factors yields a PL intensity enhanced up to 10 times at a specific angle and wavelength. Numerical simulations combining the carrier kinetics with full-wave electromagnetics in the time-domain support the experimental observations.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Interaction of quantum emitters with metallic nanostructures has been a key topic in plasmonics [1–7]. The effect of a metallic nanostructure on photoluminescence (PL) is best illustrated by the three optical parameters: absorption, quantum yield (ratio between the radiative and non-radiative decay rates), and outcoupling [6,8–11]. All the preceding works follow this scheme, although it is not often described explicitly. The metallic nanostructures reported in the previous works include optical bow-ties [12–14], Yagi-Uda antennas [15,16], plasmonic metasurfaces [17–20], and plasmonic arrays of nanoholes [21–26] /nanocylinders [27–30]. Among them, plasmonic nanocylinder arrays refer to the metallic nanocylinders arranged periodically with lattice constants on the wavelength scale. In such arrays, the plasmonic-photonic hybrid modes, i.e., lattice plasmons, represented by the coexistence of both localized surface plasmon polariton (LSPP) and diffraction modes [31–35] facilitate coherent and strong interaction between the incident light and the array. The array was used in conjunction with an emitter film for enhanced PL intensity [8]. The main contribution of the plasmonic array to PL enhancement originates from outcoupling whereby a fraction of PL that would otherwise be confined within the emitter film in the absence of the array is scattered out of the film into a specific direction predefined by the periodicity. As the array can enhance the

Vol. 27, No. 4 | 18 Feb 2019 | OPTICS EXPRESS 5083

#348160 Journal © 2019

https://doi.org/10.1364/OE.27.005083 Received 15 Oct 2018; revised 27 Dec 2018; accepted 1 Jan 2019; published 11 Feb 2019

https://doi.org/10.1364/OA_License_v1

https://crossmark.crossref.org/dialog/?doi=10.1364/OE.27.005083&domain=pdf&date_stamp=2019-02-11

outcoupling of the generated PL towards free space (light extraction), plasmonic arrays are especially useful to enhance the PL of high quantum yields emitters [36–39].

In this study, we examine the absorption enhancement contributed by the plasmonic array. The effect of absorption enhancement is especially critical for the films thinner than the absorption length (i.e., the distance over which light propagates before its intensity becomes 1/e of its original value). In such thinner films, the plasmonic array can act as an absorption enhancer: It diffracts the excitation light into the plane of the array thereby enlarging the effective thickness of the film and increasing the absorption. The absorption enhancement is very useful for efficient light-to-light conversion in solar luminescent converters [40,41] and light-to-carrier conversion in solar cells [42–44]. In our experiments, we deposited a poly(methyl methacrylate) (PMMA) film containing rhodamine 6G (R6G) dye, referred to as PMMA + R6G film, which was thinner than its absorption length on a plasmonic array. R6G was particularly chosen because of its high quantum yield and broad excitation band in the visible range that allowed us to examine the effect of absorption enhancement on PL using a range of excitation wavelength. Furthermore, we validated the experimental observations with multiphysics numerical simulations coupling the carrier kinetics of quantum emitters with full-wave electromagnetics in the time-domain. Specifically, the rate equations of carrier kinetics (derived from a density matrix formalism) are combined with a classical Maxwell solver in a joint finite-difference time-domain (FDTD) multiphysics framework using the auxiliary differential equations (ADE) approach [45–47]. The correlation between absorption enhancement and PL intensity enhancement was clarified both by the experiment and the simulation.

2. Experimental

2.1 Sample preparation

We prepared a square array of Al nanocylinders with a period of a = 400 nm that sustained hybrid modes in the visible range of wavelength. This array was designed to show hybrid modes at the spectral ranges of both excitation and emission of R6G, as will be experimentally shown in section 4.1. Al thin films were prepared on SiO2 glass substrates using a sputtering method. Al nanocylinders were fabricated using nanoimprint lithography in conjunction with reactive ion etching (RIE). The procedure is described in detail as follows. First, a resist was deposited on the Al thin film. Then, the surface of the resist was nanostructured using nanoimprint techniques (EntreTM3, Obducat) replicating the surface morphology of an Si mold. The Al thin film with a patterned resist on the top was then structured by RIE (RIE-101iPH, Samco). The array structures were examined using a scanning electron microscopy (SU8000, Hitachi).

A PMMA + R6G film containing 0.8 wt% R6G was spin-coated on the array. The thickness of the film was estimated with a surface stylus profiler (ASIQ KLA Tencor). The dielectric function of the thin film was examined using spectroscopic ellipsometry (FE-5000, Otsuka Electronics Co.) over a wavelength range of 300~800 nm (see Fig. 1).

2.2 Optical characterization

Zeroth-order optical transmission spectra of the samples were measured as functions of the angle of incidence θin. For the sake of measurement, we used an s-polarized collimated beam from a halogen lamp with a beam diameter of ~0.5 mm. The sample was mounted on a computer-controlled rotation stage. The absolute zeroth-order transmission (T) was obtained by normalizing the transmission of incident light through the sample to that through a glass substrate. The specular reflection (R) was measured at θin = 5°, and the absorptance (A) at θin = 5° was deduced from the relation: A = 1 – T – R.

The PL spectra were measured by illuminating the sample with a laser diode (LD). Three different LDs with λex = 445, 473, and 532 nm were used. The LD output (s-polarized) was


incident fromcollected fromthe sample surAl nanocylindλ = 570 nm(Quantaurus-T(temporal reso

Fig. 1fit o

( )ε ω

0ω =

3. Numerica

Fig. 2basedx-y plboundperiodthicknnm, reof R6

m the substrate m the opposite srface. A film oder array was um was measuTau, Hamamatolution of < 1 n

. Refractive indexof a Lorentzi

(2

inf 0) f= +ε ω ω

( )2.33 ,96 eV=

al simulation

2. (a) Sketch of theon FDTD Maxwe

lane (top view) of dary conditions ardicity a = 400 nm.nesses of the PMMespectively. (b) SkG.

side at θin = 5side by a fiber-

of PMMA + R6sed as a referenured using a tsu Photonics) ns) with λex = 4

and extinction coian model to

)2 2

0 i− +ω ω γω

and (0.0722 =γ

ns

e unit cell of simulell solver with em

f the structure. There applied in x- a The bottom pane

MA + R6G film, t,ketch of the four-lev

5° unless other-coupled spectr6G prepared onnce to evaluate

time-correlatequipped with

470 nm.

efficient of the PMthe spectrosc

, with inf

2=ε

)V .e

lated structure usembedded quantum e

e black circle indiand y-directions tl shows the x-z pla and the SiO2 glasvel energy diagram

rwise noted, anrometer at an an a flat SiO2 glae the PL enhancted single-phoh a pulsed ligh

MMA + R6G film copic ellipsome

2.3185, 6.8f =

ed for a 3D multipemitters. The top icates the Al nanoto simulate the sqane (side view) of ss substrate, tsub, wm simulating the el

nd the PL speangle θem from nass substrate wcement. The PLoton countinght-emitting diod

extracted from theetry data, i.e.

( )48755 ,10 eV

−×

physics frameworkpanel displays the

ocylinder. Periodicquare lattice with

f the structure. Thewere 700 and 2000lectronic transition

ctra were normal to ithout the

L decay at g module de (LED)

e ,

k e c h e 0 n


We used a semi-classical model in which the dye molecules were treated quantum-mechanically, and the electromagnetic wave was treated classically in the time domain [42–44]. The unit cell of the model is shown in Fig. 2(a). The top panel shows the top view of the structure. The nanocylinder has a diameter of 150 nm and a height of 150 nm. The array has a 400 nm period in both x and y directions. Periodic boundary conditions are applied in the lateral directions to simulate the periodicity of the array. The bottom panel shows the cross-section of the structure consisting of an SiO2 glass substrate and an Al nanocylinder covered with a 700-nm-thick PMMA + R6G film. The thickness obtained from the surface stylus profiler measurement varied from place to place on the film, and we used a representative value of 700 nm. We used the value of dispersive permittivity of Al taken from ellipsometry measurements. A constant refractive index of 1.45 was used for SiO2 glass and 1.00 for air. The incident light was linearly-polarized with the electric field parallel to the y-direction, while θin was defined in the z-x plane. The absorptance A was simulated by using the calculated values of the transmittance T and reflectance R in the relation A = 1– T– R. For the calculation of A, the PMMA + R6G film was treated as an effective medium with a dispersive permittivity whose value was taken from the fit to the ellipsometry measurement data (see Fig. 1). We simulated the absorption spectra at θin = 5° for the PMMA + R6G films on the Al nanocylinder array and that on the flat glass substrate.

Since spontaneous emission is a stochastic process, a direct implementation of it by a time-domain modeling approach is associated with a significant computational load. To simplify the complexity of the problem and to overcome the computational cost, we phenomenologically incorporated the spontaneous emission into our FDTD solvers by introducing randomly-distributed and randomly-oriented classical quantum emitters (electric dipoles) with a Lorentzian dispersion centered at an emission frequency [45].

In the model, shown in Fig. 2(b), we assumed a four-level scheme of R6G, with the pumping transition terminating at a quasi-stable short-lived state above the metastable state [48]. The rate equations are as follows:

3 3 30

32 30

1N N

t tτ ω∂ ∂

= − + ⋅∂ ∂

PE

(1i)

32 2 21

32 21 21

1NN N

t tτ τ ω∂ ∂

= − + ⋅∂ ∂

PE

(1ii)

1 2 1 21

21 10 21

1N N N

t tτ τ ω∂ ∂

= − − ⋅∂ ∂

PE

(1iii)

0 301

10 30

1N N

t tτ ω∂ ∂

= − ⋅∂ ∂

PE

(1iv)

where i is the population level density in the ith level, γij = 1/τij = 1/τr, ij + 1/τnr, ij where τij is the total lifetime between the two levels, Pij is the induced transition polarization, ωij is the transition angular frequency given by ωij = 2πc/λij, and E is the electric field. The oscillator equation that describes the temporal evolution of the polarization densities induced into the materials due to the ijth transition is:

( )2

ij ij 2ij ij ij ij j i2

k N Ntt

ω ω∂ ∂

+ Δ + = −∂∂

P PP E (2)

Where ijωΔ is the transition line width, dominated by the de-phasing time T2, 1

ij 2,ij~ 2T −Δω .

De-phasing times are usually in the fs range making the contributions from other factors


negligible. Coupling coefficients ijk are proportional to the radiative decay rates and are

defined as 3 2ij 0 r,ij h ij6 /k cπε γ ε ω= , where hε is the host material dielectric constant. The

induced polarizations are then coupled to the classical treatment of electromagnetic fields via the Maxwell curl equation:

0 h ijijt

ε ε ∂∇× = +∂ E

H P (3)

Positions of dipoles, orientation, as well as their polarization were randomly chosen during simulation. In the excitation process, the active medium was pumped from the substrate side with an angle of 5° in the z-x plane, and the emission was collected from the other side of the sample. Three pump wavelengths were used for the excitation source: λex = 445, 473, and 532 nm. The random dipoles were used to stimulate the radiation emission with a wavelength centered at 589 nm. Numerous simulations were carried out each time with different dipole properties, and all the results were averaged. That randomization and averaging were used to eliminate the coherence in the numerical model of spontaneous emission. The main parameters of numerical simulations are outlined in Table 1.

Table 1. Parameter list for numerical simulations.

Radiative transfer lifetimes

τ10 = 50 fs, τ21 = 3 ns, τ30 = 1 ns, τ32 = 50 fs

Dephasing times Τ2,30 = 19.2356 fs Τ2,21 = 22.786 fs

Transition angular frequencies

ω30 = 2πc/λ30, ω21 = 2πc/λ21, λ30 = 532 nm, λ21 = 589 nm

Number of independent runs

~100

Dipole source

( ) ( )2dip

0 0

dip 2 dip a dip

0

dip

0si 50 fsn e ;xp ,

t tE N t E t= =

−−

σωσ

dip0t is generated randomly with normal distribution;

dip tot

0 totdip tot

0

, 0.rand 7 ps.5

tt t t= =

σ

Dipole amplitudes are proportional to the population of the excited quasi-stable state, i.e., second level, λ21 = 589 nm.

Pump source 2pump

0 0

pump pump pum

pump pum

p

b 0

p

0 0

b

84 fs, 0.25 ps,

={445, 47

2

3, 532}nm

s

.

in exp ,t t

tE E t−

= = =−

σσ

λ

π

λ

4. Results and discussion

4.1 Optical absorption

Figure 3(a) shows the SEM image of the array used. The Al nanocylinders (150 nm in diameter) are arranged in a square lattice with a period of a = 400 nm. The SiO2 glass substrate is flat, and the horizontal bright and dark lines are due to the charge build-up of the insulating glass substrate.

Figures 3(b) and 3(c) show the optical transmittance (T) as a function of θin for the array before and after the deposition of PMMA + R6G film. Also shown are the conditions for the in-plane diffraction. The in-plane diffraction condition is referred to as Rayleigh anomaly, which satisfies the following relation for a square lattice when the light is incident parallel to the x-direction [49],


2 2

2 2out|| in 1 2

2 2,k k m m

a a

π π = + + (4)

where out|| 2k nπ λ= and ( )in in 2 sink nπ λ θ= are the wave vectors of the scattered light and

the incident light, respectively, n is the refractive index of the surrounding medium, a is the periodicity, and m1 and m2 are the diffraction orders, respectively. We shall refer to the magnitude of k as k. Before the film deposition, a broad dip appears at around λ = 580 nm, partly modulated by the diffraction. The modulation is much stronger for the array after the film deposition (see Fig. 3(c)) because of the mitigation of refractive index mismatch between the up- and substrates of the array by the film deposition. Figure 3(c) also shows several features that are unassignable to the in-plane diffractions. They are most probably due to excitation of quasi-waveguide modes: the PMMA + R6G thin film acts as a slab resonator to trap the incident light. Figure 3(d) compares T at θin = 0°. For the array before the film deposition (top panel), asymmetric shape of the dip with a steep edge at λ = 580 nm clearly shows that the LSPP of Al nanocylinders is modulated by the in-plane diffraction. After the deposition (bottom), the main dip redshifts to λ = 630 nm because of the increase in the refractive index of the medium around Al nanocylinders. At shorter wavelengths, several features newly appear. Absorption of R6G should appear at around λ = 532 nm, but it is weaker compared to the other resonances associated with the array (see Fig. 4 for detail) and is indistinguishable to the other resonances.

Fig. 3. (a) Top-view SEM image of the array. Scale bar = 500 nm. The x-, y-, and z-axes used in this study are also defined. The inset is the experimental configuration: the incident light is polarized along the y-direction (s-polarized), and θin was varied in the z-x plane, with the azimuth angle being fixed. (b), (c) Wavelength and incident angle dependence of zeroth-order transmission for the array without (b) and with (c) the PMMA + R6G thin film. The dashed curves are the diffraction conditions (Eq. (4)) with different orders. (d) Cuts of transmittance at θin = 0 ° for the array without (top) and with (bottom) the PMMA + R6G thin film.

Figure 4(a) shows absorptance A at the angle of incidence θin = 5°, obtained from the relation A = 1 – T – R, where T and R represent transmittance and reflectance, respectively. The spectrum for the PMMA + R6G film on a flat glass (AR6G, represented by a grey area) shows the absorption peak at λ = 530 nm, which is typical of R6G molecules. Using the AR6G and the


obtained film thickness (t = 700 nm), the absorption length (a) is calculated using the relation,

1 – A = exp(–t/a). The values of a are larger than the film thickness at the three excitation

wavelengths used in this study: As listed in Table 2, we obtained a = 2690, 2690, and 2090 nm

for λ = 445, 473, and 532 nm, respectively. Thus, the absorption is found to be moderate indicating the scope for further enhancement. The absorption of the film deposited on the array of Al nanocylinders (Aarray + R6G) is obviously larger than AR6G. The spectral shape appears irregular owing to the excitations of LSPP and light diffraction of different orders in the array. Here it is noted that the array sample shows light scattering and A = 1 – T – R does not hold in a rigorous sense for the specular components measured in this experiment. Here we used the term A for the sake of simplicity.

Figure 4(b) shows the enhancement of absorptance by the array, defined as Aarray + R6G/AR6G. Aarray + R6G/AR6G is found to be 2.32, 1.81, and 1.04 for λ = 445, 473, and 532 nm, respectively. It is noted that at θin = 5° and λ = 532 nm, there is nearly no change in absorption as inferred from the value of Aarray + R6G/AR6G ~1. The increase in Aarray + R6G compared to AR6G suggests that a part of the incident light, which would have otherwise transmitted through the film, is scattered by the array within the PMMA + R6G film and eventually absorbed by the R6G dye molecules and/or the Al nanocylinders.

Fig. 4. (a) Absorptance, A, calculated by the relation A = 1 – Transmittance – Reflectance for s-polarized light at θin = 5°. The solid line represents A for the PMMA + R6G film on the array (Aarray + R6G), while the grey area is for the same film on the flat substrate (AR6G). The arrows show the in-plane diffraction conditions. (b) Absorptance enhancement by the array, Aarray + R6G / AR6G.

Table 2. Absorption lengths at the three wavelengths in the reference film calculated by the relation 1 – A = exp(–t/a) where t is the film thickness for s-polarized light at θin = 5°.

Wavelength (nm) AR6G Absorption length (nm)

445 0.201 2690

473 0.205 2690

532 0.260 2090


4.2 PL enhancement

Fig. 5. (a) Emission angle θem dependence of PL enhancement for λex = 445 (top), 473 (middle), and 532 nm (bottom), respectively. θem was varied from 0 to 30 °, with θin being fixed to 5°. PL enhancement was defined as the spectra from the sample (Iarray + R6G) divided by that of the reference (IR6G), i.e., Iarray + R6G/IR6G. The inset shows the configuration where θem is varied in z-x-plane. (b) PL spectra at θin = 5° and collected at θem = 0°(see the inset). The dots and grey area represent the spectra from the sample and the reference, respectively. λex = 445 (top panel), 473 (middle), and 532 nm (bottom). The sharp peak in the bottom panel denoted by an arrow is the excitation laser line. (c) PL enhancement, defined as Iarray + R6G/IR6G for λex = 445, 473, and 532 nm, respectively (left axis) and zeroth-order T at θin = 0° (right). (d) Comparison between Aarray +

R6G/AR6G at λex and Iarray + R6G/IR6G averaged spectrally over λ = 545 to 650 nm.

Figure 5 summarizes the PL properties of the PMMA + R6G film on the array under the excitation with three LDs of different wavelengths: λex = 445, 473, and 532 nm. Figure 5(a) shows the PL enhancement, defined as the PL from the PMMA + R6G film on the array (Iarray +

R6G) normalized to that from the reference film on a flat substrate (IR6G), i.e., Iarray + R6G/IR6G, as a


function of emission angle, θem (see the inset for the definition). Also shown are the in-plane diffraction conditions as in Figs. 3(b) and 3(c). At θem = 0°, notable enhancement is observed, which blueshifts with the increase in θem following the (0, ± 1) diffraction order. Another line of enhancement along the (−1,0) order is also notable. This similarity in angular profiles between the PL enhancement and the diffraction indicates the outcoupling effect.

Figure 5(b) compares the PL spectra at θem = 0°. The dotted lines and the grey area represent Iarray + R6G and IR6G, respectively. The reference shows a broad PL spectrum centered at λ = 563 nm, which is typical of R6G. In addition to the peak of this broad spectra, there are extra PL peaks observed at approximately λ = 576 and 585 nm for Iarray + R6G. It is noted that these peaks are associated to diffraction; the (−1,0), (0, ± 1), and (1,0) diffraction orders are degenerate at λ = 580 nm and θem = 0°.

In order to confirm the outcoupling effect, we plot Iarray + R6G/IR6G at θem = 0° and T at θin = 0° (Fig. 5(c)). The spectral positions of the dips in T (at λ = 572 and 585 nm) correspond to the extra peaks in Iarray + R6G/IR6G thereby clearly indicating the outcoupling effect: The dips in T mean that the light is coupled inside the thin film for the specific condition of incidence. Conversely, when the PL is generated inside the film, a fraction of the PL that possesses the corresponding λ is coupled outside the film at the corresponding direction, i.e., the array acts as a spatial and spectral filter for PL.

Although the spectral shapes are similar, the degree of enhancement depends notably on λex. While the largest value of Iarray + R6G/IR6G, as large as 10 times at λ = 585 nm (see Fig. 5(c)), is found for λex = 445 nm, the smallest value is found for λex = 532 nm. This variation comes from the absorption enhancement at each λex. Figure 5(d) compares the Iarray + R6G/IR6G averaged over λ = 545 to 650 nm and the Aarray + R6G/AR6G at the three λex. It is found that Iarray + R6G/IR6G remarkably scales with Aarray + R6G/AR6G. As it is also inferred from Fig. 4(b) that Aarray + R6G/AR6G = 1.04 (close to unity) for λex = 532 nm, it is evident that the PL enhancement observed for this wavelength is not owing to the absorption enhancement but owing to the outcoupling aspect. For λex = 445 and 473 nm, the integrated values of Iarray + R6G/IR6G are 3.78 and 3.20 times, respectively, implying that the enhancement originated from both outcoupling and absorption. As the outcoupling depends only on the PL process and does not depend on the excitation conditions, the difference in the values of Iarray + R6G/IR6G between the three conditions is attributed to absorption enhancement. Compared to λex = 532 nm where Iarray + R6G/IR6G = 1.39 and absorption enhancement is negligibly small as explained above, the Iarray + R6G/IR6G for λex = 445 nm (473 nm) is 2.72 (2.30) times larger while Aarray + R6G/AR6G being larger by 2.23 (1.74) times. The slight differences between Iarray + R6G/IR6G and Aarray + R6G/AR6G indicate that the fraction of Aarray + R6G absorbed by R6G dyes varies with λex; not only the R6G dye molecules, but the metallic nanocylinders can also absorb.

Fig. 6. PL decay curves of the PMMA + R6G film on the array (red curve) and that on the glass substrate(black). The excitation wavelength was 470 nm.

We investigated the other enhancement factor, quantum yield, using PL decay measurement. Figure 6 shows the PL decay curves excited with a pulsed LED with λex = 470


nm for the PMMA + R6G films deposited on the Al nanocylinder array (red line), and on the glass substrate (black line) as a reference. The PL lifetimes (τ) for the film on the array and reference are τarray + R6G = 2.15 ns and τR6G = 2.34 ns, respectively. We speculated the decrease in τarray + R6G equally occurs for all the three excitation wavelengths used, because the PL lifetime only depends on the density of optical states in the initial and final states of the transition and is not dependent on the excitation wavelength. The decrease in τ is only 8%, which is small with respect to the spectrally-integrated PL intensity enhancement up to 3.78 times. We speculated that the quantum yield modification plays a minor contribution to the PL enhancement in the present system.

4.3 Comparison of experimental results with numerical simulations

Fig. 7. Simulated absorptance spectra of PMMA + R6G films (thickness = 700 nm) at θin = 5° on the Al nanocylinder array (Aarray + R6G) (denoted as a red line) and that on the flat glass substrate (AR6G) (grey area)

Figure 7 compares the simulated absorption spectra for the PMMA + R6G films on the Al array (Aarray + R6G) and on the flat substrate (AR6G), where the dispersive permittivity of PMMA + R6G film was taken from measurement data of spectroscopic ellipsometry (see Fig. 1). Absorption of the PMMA + R6G film on the flat substrate, AR6G, shows an absorption band centered at λ = 530 nm, which is typical for the R6G dye. Unlike the experimental spectrum (Fig. 4), the simulated absorption spectrum for the case of the flat substrate is found to exhibit a fringe-like pattern indicating the effect of Fresnel reflectance at the surfaces. This fringe is not observed in the experiment owing to the surface corrugation. In contrast, Aarray + R6G exhibits a considerable increase in the magnitude and also the additional peaks. This simulation qualitatively reproduces the experimental absorption spectra in Fig. 4(a), including the spectral positions of the in-plane diffractions; the discrepancy may be attributable to the film surface morphology. It is also noted that the simulation contains both specular and diffuse components, while the experiment contains specular component only. This difference also causes the discrepancy. A comparison between the curves of Aarray + R6G and AR6G clearly manifests the increased absorption by the R6G molecules attributable to the functioning of the PMMA + R6G film as a slab waveguide that receives the incident light scattered by the Al nanocylinder array. At each excitation wavelength, we calculated the simulated absorption enhancement Aarray + R6G/AR6G that was found to be 8.1, 18.7, and 2.6 at λ = 445, 473, and 532 nm, respectively.


Fig. 8. Simulated PL spectra at θin = 5° and θem = 0°. The solid line and grey area represent the spectra from the sample (Iarray + R6G) and the reference (IR6G), respectively. λex = 445 (top panel), 473 (middle), and 532 nm (bottom).

Figure 8 shows the simulated PL spectra excited at the three different wavelengths. The solid line and the grey area represent the PL from the PMMA + R6G film on the array (Iarray +

R6G) and the reference film on the flat substrate (IR6G), respectively. The reference shows a broad PL spectrum centered at λ = 583 ~596 nm, which is typical of R6G. In addition to this broad peak, extra PL peaks at approximately λ = 600 ~612 nm are observed for the spectra of PMMA + R6G film on the arrays. The PL spectra reproduce the experimental results qualitatively although the simulated PL shows a single additional peak at λ = 600 ~612 nm and the experimental PL shows a main peak at λ = 585 nm with a subpeak at λ = 576 nm. The discrepancy in spectral positions, as well as the number of the peaks, may be attributable to the fitting parameters adopted, including λ21, the number density of the dipole and the geometrical parameters (dimensions) of the Al nanocylinders. A better fit could be achieved by finely tuning these parameters. We were not keen on further optimization because the simulated results already reproduced the optical phenomenon in line with the theoretical expectations. We evaluated the PL intensity enhancement Iarray + R6G/IR6G averaged over λ = 545 to 650 nm from the spectra in Fig. 8, and it was found to be 2.7, 2.8, and 2.2 times for λex = 445, 473 and 532 nm, respectively.

Fig. 9. Comparison between the simulated Aarray + R6G/AR6G (left axis) and Iarray + R6G/IR6G averaged spectrally over λ = 545 to 650 nm (right).


Figure 9 compares the absorption enhancement and the PL enhancement. Except for the difference in the order of the simulated and the experimental results (see Fig. 5(d)) of PL enhancement which is attributed to the difference in the surface flatness between the sample and the model as explained above, the correlation between the absorption and the PL enhancements is found to be qualitatively reproduced. The discrepancy in the magnitude of Aarray + R6G/AR6G between the experiment and the simulation also comes from the surface corrugation. The flat surface of the simulated structure results in the Fresnel-like spectral fringe for the reference film that underestimate AR6G, which results in the larger Aarray + R6G/AR6G for the simulation. The simulation confirms that the variation of PL enhancement with the excitation wavelength originates from the dependency of the absorption enhancement at the wavelength.

5. Summary

In summary, we evaluated the absorption enhancement of the light emitter films thinner than the absorption length deposited on a plasmonic array of Al nanocylinders, both experimentally and numerically. For this purpose, a PMMA + R6G film was deposited on the square array of Al nanocylinders with a 400 nm periodicity. The array was designed to support hybrid modes in a spectral range that covered the excitation and emission wavelength of R6G, thus; it can influence both processes via absorption enhancement and outcoupling, respectively. The PL enhancement at the three excitation wavelengths (λex = 445, 473, and 532 nm) was found to scale with the absorption enhancement. The results helped to establish that the array scattered the excitation light that was eventually reabsorbed by the R6G molecules in the film, thereby enhancing the PL intensity. The comparison to other materials, including plasmonic Ag and nonplasmonic TiO2 and Si, would give further insight into the effect of LSPPs on the absorption enhancement. Our study also brings out the importance of quantum emitters with a small Stokes shift, such as the R6G dye used in this work, as a key utility in controlling both the plasmon-enhanced absorption and PL outcoupling by a single array.

Funding

Grant-in-Aid for Scientic Research (B, 16H04217 and Exploratory, 17K19176) by MEXT;

Precursory Research for Embryonic Science and Technology (PRESTO, JPMJPR131B) from JST; the Asahi Glass Foundation; Nippon Sheet Glass Foundation for Materials Science and Engineering; Nanotech CUPAL (SM and SI); Defense Advanced Research Projects Agency (DARPA/DSO) Extreme Optics and Imaging (EXTREME) Program, (HR00111720032, A.V.K and S.I.A).

Acknowledgment

A part of this work was supported by Kyoto University Nano Technology Hub and NIMS Nanofabrication Platform in the “Nanotechnology Platform Project”, Japan.

References

1. H. Morawitz, “Self-coupling of a two-level system by a mirror,” Phys. Rev. 187(5), 1792–1796 (1969). 2. K. H. Drexhage, “Influence of a dielectric interface on fluorescence decay time,” J. Lumin. 1–2, 693–701 (1970). 3. K. H. Tews, “Zur Variation von Lumineszens-Lebensdauern,” Ann. Phys. (Leipz.) 29(2), 97–120 (1973). 4. T. Okamoto, F. H’Dhili, and S. Kawata, “Towards plasmonic band gap laser,” Appl. Phys. Lett. 85(18), 3968–

3970 (2004). 5. R. R. Chance, A. Prock, and R. Silbey, “Lifetime of an excited molecule near a metal mirror: Energy transfer in the

Eu3+/silver system,” J. Chem. Phys. 60(5), 2184–2185 (1974). 6. K. Okamoto, I. Niki, A. Shvartser, Y. Narukawa, T. Mukai, and A. Scherer, “Surface-plasmon-enhanced light

emitters based on InGaN quantum wells,” Nat. Mater. 3(9), 601–605 (2004). 7. J. R. Lakowicz, K. Ray, M. Chowdhury, H. Szmacinski, Y. Fu, J. Zhang, and K. Nowaczyk, “Plasmon-controlled

fluorescence: a new paradigm in fluorescence spectroscopy,” Analyst (Lond.) 133(10), 1308 (2008). 8. G. Lozano, D. J. Louwers, S. R. K. Rodríguez, S. Murai, O. T. A. Jansen, M. A. Verschuuren, and J. Gómez Rivas,

“Plasmonics for solid-state lighting: enhanced excitation and directional emission of highly efficient light sources,” Light Sci. Appl. 2(5), e66 (2013).


9. G. Lozano, S. R. K. Rodriguez, M. A. Verschuuren, and J. Gómez Rivas, “Metallic nanostructures for efficient LED lighting,” Light Sci. Appl. 5(6), e16080 (2016).

10. S. Murai, M. A. Verschuuren, G. Lozano, G. Pirruccio, S. R. K. Rodriguez, and J. G. Rivas, “Hybrid plasmonic-photonic modes in diffractive arrays of nanoparticles coupled to light-emitting optical waveguides,” Opt. Express 21(4), 4250–4262 (2013).

11. J. Feng, T. Okamoto, and S. Kawata, “Highly directional emission via coupled surface-plasmon tunneling from electroluminescence in organic light-emitting devices,” Appl. Phys. Lett. 87(24), 241109 (2005).

12. A. Kinkhabwala, Z. F. Yu, S. H. Fan, Y. Avlasevich, K. Mullen, and W. E. Moerner, “Large single-molecule fluorescence enhancements produced by a bowtie nanoantenna,” Nat. Photonics 3(11), 654–657 (2009).

13. O. L. Muskens, V. Giannini, J. A. Sánchez-Gil, and J. Gómez Rivas, “Strong Enhancement of the Radiative Decay Rate of Emitters by Single Plasmonic Nanoantennas,” Nano Lett. 7(9), 2871–2875 (2007).

14. R. M. Bakker, H.-K. Yuan, Z. Liu, V. P. Drachev, A. V. Kildishev, V. M. Shalaev, R. H. Pedersen, S. Gresillon, and A. Boltasseva, “Enhanced localized fluorescence in plasmonic nanoantennae,” Appl. Phys. Lett. 92(4), 043101 (2008).

15. T. Kosako, Y. Kadoya, and H. F. Hofmann, “Directional control of light by a nano-optical Yagi-Uda antenna,” Nat. Photonics 4(5), 312–315 (2010).

16. T. Coenen, E. J. R. Vesseur, A. Polman, and A. F. Koenderink, “Directional emission from plasmonic Yagi-Uda antennas probed by angle-resolved cathodoluminescence spectroscopy,” Nano Lett. 11(9), 3779–3784 (2011).

17. M. Y. Shalaginov, S. Ishii, J. Liu, J. Liu, J. Irudayaraj, A. Lagutchev, A. V. Kildishev, and V. M. Shalaev, “Broadband enhancement of spontaneous emission from nitrogen-vacancy centers in nanodiamonds by hyperbolic metamaterials,” Appl. Phys. Lett. 102(17), 173114 (2013).

18. D. Lu, J. J. Kan, E. E. Fullerton, and Z. Liu, “Enhancing spontaneous emission rates of molecules using nanopatterned multilayer hyperbolic metamaterials,” Nat. Nanotechnol. 9(1), 48–53 (2014).

19. J. S. T. Smalley, F. Vallini, S. A. Montoya, L. Ferrari, S. Shahin, C. T. Riley, B. Kanté, E. E. Fullerton, Z. Liu, and Y. Fainman, “Luminescent hyperbolic metasurfaces,” Nat. Commun. 8, 13793 (2017).

20. K. Tanaka, E. Plum, J. Y. Ou, T. Uchino, and N. I. Zheludev, “Multifold Enhancement of Quantum Dot Luminescence in Plasmonic Metamaterials,” Phys. Rev. Lett. 105(22), 227403 (2010).

21. M. Saboktakin, X. Ye, U. K. Chettiar, N. Engheta, C. B. Murray, and C. R. Kagan, “Plasmonic Enhancement of Nanophosphor Upconversion Luminescence in Au Nanohole Arrays,” ACS Nano 7(8), 7186–7192 (2013).

22. A. G. Brolo, S. C. Kwok, M. G. Moffitt, R. Gordon, J. Riordon, and K. L. Kavanagh, “Enhanced Fluorescence from Arrays of Nanoholes in a Gold Film,” J. Am. Chem. Soc. 127(42), 14936–14941 (2005).

23. C.-H. Lu, C.-C. Lan, Y.-L. Lai, Y.-L. Li, and C.-P. Liu, “Enhancement of Green Emission from InGaN/GaN multiple quantum wells via coupling to surface plasmons in a two-dimensional silver array,” Adv. Funct. Mater. 21(24), 4719–4723 (2011).

24. X. Meng, J. Liu, A. V. Kildishev, and V. M. Shalaev, “Highly directional spaser array for the red wavelength region,” Laser Photonics Rev. 8(6), 896–903 (2014).

25. A. V. Dorofeenko, A. A. Zyablovsky, A. P. Vinogradov, E. S. Andrianov, A. A. Pukhov, and A. A. Lisyansky, “Steady state superradiance of a 2D-spaser array,” Opt. Express 21(12), 14539–14547 (2013).

26. F. van Beijnum, P. J. van Veldhoven, E. J. Geluk, M. J. A. de Dood, G. W. ’t Hooft, and M. P. van Exter, “Surface plasmon lasing observed in metal hole arrays,” Phys. Rev. Lett. 110(20), 206802 (2013).

27. W. Zhou, M. Dridi, J. Y. Suh, C. H. Kim, D. T. Co, M. R. Wasielewski, G. C. Schatz, and T. W. Odom, “Lasing action in strongly coupled plasmonic nanocavity arrays,” Nat. Nanotechnol. 8(7), 506–511 (2013).

28. P. Törmä and W. L. Barnes, “Strong coupling between surface plasmon polaritons and emitters: a Review,” Rep. Prog. Phys. 78(1), 013901 (2015).

29. R. Guo, S. Derom, A. I. Väkeväinen, R. J. van Dijk-Moes, P. Liljeroth, D. Vanmaekelbergh, and P. Törmä, “Controlling quantum dot emission by plasmonic nanoarrays,” Opt. Express 23(22), 28206–28215 (2015).

30. F. Laux, N. Bonod, and D. Gerard, “Single emitter fluorescence enhancement with surface lattice resonances,” J. Phys. Chem. C 121(24), 13280–13289 (2017).

31. V. G. Kravets, A. V. Kabashin, W. L. Barnes, and A. N. Grigorenko, “Plasmonic surface lattice resonances: A review of properties and applications,” Chem. Rev. 118(12), 5912–5951 (2018).

32. V. G. Kravets, F. Schedin, and A. N. Grigorenko, “Extremely narrow plasmon resonances based on diffraction coupling of localized plasmons in arrays of metallic nanoparticles,” Phys. Rev. Lett. 101(8), 087403 (2008).

33. B. Auguié and W. L. Barnes, “Collective resonances in gold nanoparticle arrays,” Phys. Rev. Lett. 101(14), 143902 (2008).

34. Y. Z. Chu, E. Schonbrun, T. Yang, and K. B. Crozier, “Experimental observation of narrow surface plasmon resonances in gold nanoparticle arrays,” Appl. Phys. Lett. 93(18), 181108 (2008).

35. S. Zou, N. Janel, and G. C. Schatz, “Silver nanoparticle array structures that produce remarkably narrow plasmon lineshapes,” J. Chem. Phys. 120(23), 10871–10875 (2004).

36. P. Anger, P. Bharadwaj, and L. Novotny, “Enhancement and quenching of single-molecule fluorescence,” Phys. Rev. Lett. 96(11), 113002 (2006).

37. Y. Kawachiya, S. Murai, M. Saito, H. Sakamoto, K. Fujita, and K. Tanaka, “Collective plasmonic modes excited in Al nanocylinder arrays in the UV spectral region,” Opt. Express 26(5), 5970–5982 (2018).

38. S. Murai, M. Saito, H. Sakamoto, M. Yamamoto, R. Kamakura, T. Nakanishi, K. Fujita, M. Verschuuren, Y. Hasegawa, and K. Tanaka, “Directional outcoupling of photoluminescence from Eu (III)-complex thin films by plasmonic array,” APL Photonics 2(2), 026104 (2017).


39. R. Kamakura, S. Murai, Y. Yokobayashi, K. Takashima, M. Kuramoto, K. Fujita, and K. Tanaka, “Enhanced photoluminescence and directional white-light generation by plasmonic array,” J. Appl. Phys. 124(21), 213105 (2018).

40. Z. Hosseini, N. Taghavinia, and E. Wei-Guang Diau, “Luminescent spectral conversion to improve the performance of dye-sensitized solar cells,” ChemPhysChem 18(23), 3292–3308 (2017).

41. H. J. Hovel, R. T. Hodgson, and J. M. Woodall, “The effect of fluorescent wavelength shifting on solar cell spectral response,” Sol. Energy Mater. 2(1), 19–29 (1979).

42. H. A. Atwater and A. Polman, “Plasmonics for improved photovoltaic devices,” Nat. Mater. 9(3), 205–213 (2010). 43. V. E. Ferry, L. A. Sweatlock, D. Pacifici, and H. A. Atwater, “Plasmonic nanostructure design for efficient light

coupling into solar cells,” Nano Lett. 8(12), 4391–4397 (2008). 44. Y. Zhao, N. Hoivik, M. N. Akram, and K. Wang, “Study of plasmonics induced optical absorption enhancement of

Au embedded in titanium dioxide nanohole arrays,” Opt. Mater. Express 7(8), 2871–2879 (2017). 45. X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett. 85(1), 70–73 (2000). 46. S. Azzam, Z. Wang, S. Murai, S. Ishii, A. Boltasseva, and A. Kildishev, “Time domain modeling of lasing

dynamics in hyperbolic metamaterials,” Conference on Lasers and Electro-Optics, OSA Technical Digest Online, Optical Society of America (2017) JTu5A.40 .

47. S. I. Azzam, M. Saito, S. Murai, S. Ishii, K. Tanaka, and A. Kildishev, “Enhanced spontaneous emission of quantum emitters in the vicinity of TiN thin films,” in Conference on Lasers and Electro-Optics, OSA Technical Digest (online) (Optical Society of America, 2018), JTh2A.87.

48. M. A. Noginov, G. Zhu, M. Bahoura, C. E. Small, C. Davison, J. Adegoke, V. P. Drachev, P. Nyga, and V. M. Shalaev, “Enhancement of spontaneous and stimulated emission of a rhodamine 6G dye by an Ag aggregate,” Phys. Rev. B Condens. Matter Mater. Phys. 74(18), 184203 (2006).

49. G. Vecchi, V. Giannini, and J. Gómez Rivas, “Surface modes in plasmonic crystals induced by diffractive coupling of nanoantennas,” Phys. Rev. B Condens. Matter Mater. Phys. 80(20), 201401 (2009).


enhanced absorption and photoluminescence from dye .../_publications/...2018/10/15 · of the film...

Documents