Abstract

Polymer films containing plasmonic nanostructures are of increasing interest for development of responsive energy, sensing, and therapeutic systems. The present work evaluates heat dissipated from power absorbed by resonant gold (Au) nanoparticles (NP) with negligible Rayleigh scattering cross sections randomly dispersed in polydimethylsiloxane (PDMS) films. Finite element analysis (FEA) of heat transport was coordinated with characterization of resonant absorption by Mie theory and coupled dipole approximation (CDA). At AuNP particle separation greater than resonant wavelength, correspondence was observed between measured and CDA-predicted optical absorption and FEA-derived power dissipation. At AuNP particle separation less than resonant wavelength, measured extinction increased relative to predicted values, while FEA-derived power dissipation remained comparable to CDA-predicted power absorption before lagging observed extinguished power at higher AuNP content and resulting particle separation. Effects of isolated particles, for example, scattering, and particle-particle interactions, for example, multiple scattering, aggregation on observed optothermal activity were evaluated. These complementary approaches to distinguish contributions to resonant heat dissipation from isolated particle absorption and interparticle interactions support design and adaptive control of thermoplasmonic materials for a variety of implementations.

1. Introduction

Thermal damping of resonant nanoparticles (NPs) dispersed in optically transparent polymer could affect potential implementations in biomedical therapeutics [1, 2], solar cells [36], optical interconnects [79], sensing [1012], and chemical separation [13]. Optical damping by subwavelength Mie scatterers [14] in the Rayleigh regime [15] whose scattering cross sections are negligible [16] is due to absorption [17]. Resonant power absorption of isolated nanoparticles [16] and colloid suspensions [18] in homogeneous dielectric environments is describable using Beer-Lambert linearization [19, 20] of Mie theory. Interactions between isolated subwavelength scatterers and their effects on absorption require more extensive characterization, such as dipole approximation [21, 22], finite difference time domain, or T-matrix. Additionally, interference to forward scattering, for example, reflection, refraction, or diffraction, due to an obstacle, [23] dielectric interface [24], applied field [25, 26], or other heterogeneity, may also impact optical extinction [2729], thermal dynamics [30], and temperature dispersion [3133]. Evaluation of heat dissipated from power absorbed in nanocomposites to distinguish contributions from isolated scatterers, interacting scatterers, and heterogeneity is important to advance understanding and guide design of thermally responsive nanocomposite materials.

Mie’s solution to Maxwell’s equations elegantly determines scattering and absorption of an electromagnetic plane wave by an isolated, homogenous sphere [34]. The solution, which depends on particle geometry, composition, surrounding dielectric medium, and angle of incidence, allows calculation of efficiencies and cross sections of scattering and absorption as well as intensity distributions. The coupled dipole approximation (CDA) to Maxwell’s equations incorporates Mie particle polarizability and also accounts for interparticle interactions [21, 22]. Multiplying the extinction cross section from homogeneous nanocomposite media containing spherical particles determined via Mie or CDA by the number of particles yields the overall optical response [35]. This is valid when absorbance is linearly proportional to concentration, as in Beer-Lambert law [17]. Beer-Lambert’s law, , relates spectral absorbance, , to AuNP concentration, , and optical pathlength, , using a material- and wavelength-specific absorbance coefficient, [17, 36]. Deviations from Mie and Beer-Lambert Law arise due to particle interactions, for example, multiple scattering, or strongly absorbing media [37, 38]. For nanoparticles in the Rayleigh regime with appreciable absorption cross section, damping of resonant absorption yields a strong photothermal response influenced by incident intensity, NP morphology [39, 40], and host media [4143].

Heating by absorptive nanoparticles has been examined widely, but accounting for incident and dissipated energy by all mechanisms, that is, calibrating heat dissipated per unit absorption, is rare. Pulsed or continuous resonant irradiation of isolated NP or their assemblies is sufficient to reshape NP [44] and melt (evaporate) surrounding solids (liquids) [42, 45, 46]. The first comprehensive microscale calorimetric description of resonant absorption by NP in colloidal suspension analyzed dissipation via convection, conduction, and radiation [36, 47]. It has since been applied to quantitate resonant absorptive heating of NP deposited randomly on ceramic [41] and polymer [31] substrates and in multiphase [46], phase-change, and open [42] systems. Compact analytic expressions are useful to intuit dissipation and characterize thermal dynamics from geometric and thermodynamic nanocomposite features for effectively one-dimensional dissipation [30]. However, computational analysis remains necessary for systems with multidimensional absorption and dissipation. Correlating resonant absorption due to isolated particles by Mie theory and to particle interaction by CDA with heat dissipation by finite element analysis (FEA) for NPs dispersed randomly in multidimensional systems would increase understanding of these systems to support their integration into opto- and bioelectronic devices [4851].

The present work used FEA to account for heat dissipation from ca. 1 mm thick PDMS films in which 16 nm AuNP with negligible Rayleigh scattering cross sections were randomly dispersed. Absorbed resonant power was characterized by Mie and CDA solutions. Power dissipation corresponding to measured and FEA-computed temperature profiles at thermal equilibrium was compared with estimates of extinguished power from Mie and CDA. Heat dissipated per unit absorption in AuNP-PDMS was consistent with reported optoplasmonic efficiencies for 15 to 18 nm AuNPs [52]. The compact description developed herein could be useful to guide intuition, design, and development of responsive plasmonic energy materials, sensors, MEMS, and therapeutics for heat-sensitive applications.

2. Materials and Methods

2.1. AuNP-PDMS Thin Film Fabrication

Briefly [29], polydimethylsiloxane (PDMS) films containing randomly dispersed AuNP were fabricated by mixing 1 mg/mL isopropanol suspension of 16 nm ± 2.4 nm diameter AuNP into uncured PDMS at increasing volumes to achieve AuNP concentrations ranging from 0.234 to 3.52 × 1012 NP/cm3. Resulting mixtures were degassed, poured in polystyrene sample boxes, covered, and cured in an oven at 60°C for 24 h.

2.2. Optical Characterization

A light microscope (Eclipse LV100, Nikon Instruments, Melville, NY, USA) integrated with a spectrometer (Shamrock 303, Andor Technology, Belfast, UK) was used to measured AuNP-PDMS spectral response. Amplitudes at 532 nm, the excitation wavelength for subsequent thermal experiments, were determined as the difference between 532 nm extinction and extinction measured off-resonance (800 nm) for each sample. This accounted for contributions to extinction from the polymer matrix.

2.3. Mie Solution

Results from Mie solutions for isolated particles were obtained using a publicly available resource, http://nanocomposix.com/pages/tools. Simulations used a 16 nm AuNP in a refractive index corresponding to PDMS (1.42). For these conditions, resonant absorption cross sections for 16 nm AuNP are ~99% of total extinction. Mie cross sections are reported in nm2 as a function of wavelength.

2.4. Coupled Dipole Approximation

The CDA [21, 53] treats the nanoparticle in an ensemble [54, 55] as a single dipole whose polarization is directly proportional to the local electric field [56] where is a frequency-dependent polarizability and is the sum of incident irradiation () and fields scattered from other particles in the lattice (). Particle polarizability can be obtained analytically [55], computationally [56, 57], or by series approximations for higher order modes [58]. Effects of changing the angle of incident radiation can be evaluated [59]. The present work solved the CDA with a user-defined array by matrix inversion [51, 60]. Each dipole pair was calculated for a finite number of dipoles and superposed to determine the polarization vector, , at each dipole. Simulations in the present work were performed for 16 nm diameter spheres in a media refractive index for PDMS (1.42). A square 150 × 150 grid was used (90,601 dipoles) at a grating constant of double the Wigner-Seitz radius; , where is the particle radius, is the density of gold (19.3 g/cm3), is gold mass per cubic centimeter of host media, is the media volume, and is the number of particles.

2.5. Thermal Characterization

Resonant laser irradiation of AuNP-PDMS was performed with a 532 nm laser (MXL-FN-532, CNI, Changchun, CN) with an intensity of ca. 25 mW. An infrared camera (ICI 7320, P-Series, Beaumont, TX, USA) captured thermal images of the AuNP-PDMS at 5 Hz during a 10-second ambient capture, 2-minute-and-50-second heating (laser irradiation), and 3-minute cooling period. Tweezers were used to position AuNP-PDMS films such that the incident laser spot (~3 mm) was towards the top of the film, distal from the tweezer tip. Experimental apparatus was enclosed during data collection to minimize effects of forced convection. Each trial was performed in duplicate. Reported values were from a single representative trial.

2.6. Finite Element Analysis (FEA)

Briefly [32], finite element analysis (FEA) used the Heat Transfer in Solids module in Version 5.2a of COMSOL Multiphysics (COMSOL, Stockholm, SE). Thermal equilibrium was simulated with an applied heat source that corresponded to laser irradiation for comparison with measured steady-state temperatures. Radiative, conductive, and convective cooling boundary conditions were used to estimate the rate at which heat was transferred to the surrounding environment (air). Thermal diffusivity controlled the developed temperature profile within the nanocomposite media. Measured two-dimensional AuNP-PDMS temperature profiles at the first instant of cooling were projected into the three-dimensional COMSOL geometry for comparison with steady-state simulated heating results.

The AuNP-PDMS films were assigned physical dimensions corresponding to measured height, length, and thickness. Each film had slightly different dimensions but was approximately 5 × 5 × 1 mm. Values of density and specific heat capacity for each AuNP-PDMS film were estimated as weighted averages of respective values for Au and PDMS based on mass fraction of Au in the film. The thermal conductivity of AuNP-PDMS was assumed to be that of Au-free PDMS, as low mass fraction metallic dispersions in polymer have a negligible effect on bulk thermal conductivity [61]. Density, specific heat capacity, and thermal conductivity values for PDMS used were 970 kg/m3, 1460 kJ/kg K, and 0.16 W/m K, respectively [62]. For Au, values used were 19320 kg/m3, 128 J/kg K, and 317.9 W/mK [63], respectively.

Plasmonic heating from laser irradiation was represented by a Gaussian volumetric heat source centered at the measured laser spot center with a standard deviation of , where was the measured laser spot radius of ca. 1.5 mm. The full width at half maximum (FWHM) of this Gaussian distribution was 1.8 mm. Based on a raw incident laser intensity of 25 mW and this Gaussian function, the laser power incident on the sample, , was estimated by centering the Gaussian function at the apparent laser spot center based on measured thermal results. This resulted in slightly suppressed powers ranging from 20 to 24 mW; that is, outer portions of the incident beam were off the sample. Spectral extinguished powers () were obtained from , where was spectral extinction in absorbance units. Thermally dissipated power value was the fitting parameter in FEA used to match model-derived equilibrium temperature distributions with temperature distributions measured by the infrared camera at steady-state. Its value was adjusted until resulting FEA equilibrium maximum and minimum temperatures within the laser spot were within 0.1°C of measured steady-state values.

3. Results and Discussion

3.1. Optical Response Increased as Particle Separation Decreased below Resonant Wavelength

Measured optical responses of AuNP-PDMS films increased with NP concentration consistent with Mie and CDA results at low concentrations. However, when interparticle separation exceeded resonant wavelength, measured responses increased above calculated extinction. Figure 1(a) compares measured (solid), Mie (dotted), and CDA (dashed) spectra from PDMS films containing 16 nm AuNP at 1.17 × 1012 NP/cm3 (0.005 mass-percent AuNP; blue) and 2.34 × 1012 NP/cm3 (0.01 mass-percent AuNP; red). Spectral units are absorbance coefficient, calculated as optical extinction in absorbance units (AU) at 532 nm divided by sample thickness (mm), relative to off-resonance extinction at 800 nm. At 1.17 × 1012 NP/cm3, the maximum extinction per mm for measured, Mie, and CDA values was 0.121, 0.140, and 0.120, respectively. At 2.34 × 1012 NP/cm3, the maximum measured extinction per thickness was 0.313, thirty percent on average above Mie and CDA estimates of 0.252 and 0.23, respectively. Whether the 2.6-fold increase in optical response from a 2-fold rise in AuNP concentration translated to increased heat dissipation motivated further study.

This work defined extinction per NP to compare extinction in PDMS films containing 16 nm AuNP at concentrations ranging from 0.234 to 3.52 × 1012 NP/cm3 [29]. Measured extinction per NP was calculated as , where was the extinction magnitude in absorbance units (AU), at 532 nm, was the concentration in NP/cm3, was the spectrum-averaged refractive index of PDMS (1.42), and was film thickness in cm. Extinction per NP was estimated from Mie theory as using Beer-Lambert absorbance yielding , where was the Mie theory peak extinction cross section in cm2, 2.78 × 10−12 cm2, computed from an online tool (http://nanocomposix.com/pages/tools), and was a conversion factor between log-bases. Extinction per NP calculated from estimated Beer-Lambert absorption is thus independent of AuNP concentration. The output of CDA was extinction efficiency which was converted to an extinction cross section by multiplication by NP geometric cross section; that is, , where was sphere radius of 8 nm. Extinction per NP was then estimated from CDA as using this cross section, 2.44 × 10−12 cm2.

CDA was used to evaluate effects of interparticle interactions on extinction per NP. CDA results were computed using a square particle lattice with an interparticle separation distance double the Wigner-Seitz radius (upper-axis Figure 1(b)). The Wigner-Seitz radius, , for a given particle concentration corresponds to the mean spherical volume of medium per particle. It is defined as , where is particle radius, is gold density (19.3 g/cm3), and is gold mass per cubic centimeter of PDMS. It was introduced in condensed matter physics to describe the density of a system [64] and estimates distance between NP in 3D dispersions [29, 32, 36].

Measured extinction per NP in AuNP-PDMS films was 9% lower on average than Mie and 2% higher on average than CDA estimates at AuNP concentrations for which the Wigner-Seitz radius, , was greater than localized surface plasmon resonance wavelength,  nm; that is, . For , measured extinction per NP increased up to 1.4-fold relative to Mie and CDA estimates at the highest AuNP concentration. Figure 1(b) shows measured resonant extinction per NP for 16 nm AuNP with negligible scattering cross section randomly dispersed in PDMS films (blue circles). Corresponding results from the Mie (red dotted) and CDA (green dash) solutions are shown.

The top horizontal axis in Figure 1(b) shows that Wigner-Seitz radii decreased from 1006 nm to 408 nm as corresponding NP concentration on the bottom horizontal axis increased. Table 1 correlates AuNP concentration to Wigner-Seitz radius. At 1.76 × 1012 NP/cm3, the first sample with Wigner-Seitz radii below the LSPR wavelength, extinction per NP measured in AuNP-PDMS films surpassed the corresponding Mie estimate (8.91 versus 8.45 × 10−13 cm2/NP). Measured extinction per NP surpassed CDA estimates at 1.17 × 1012 NP/cm3 (7.65 versus 7.41 × 10−13 cm2/NP). As AuNP content increased, measured values diverged upward further from constant estimates, rising to a value of 12.0 × 10−13 cm2/NP, 1.4- and 1.6-fold above values predicted using Mie and CDA at 3.52 × 1012 NP/cm3. Depending on the mechanism producing it, an increase in optical extinction could benefit solar cells, as plasmon-induced electric fields and light scattering are reported to improve efficiency [6].

The observed increase in optical response for subwavelength particles randomly dispersed to interparticle distances less than the resonant wavelength could arise due to effects from isolated particles [15], particle-particle interactions [23], or optical interference, for example, dielectric heterogeneity [65]. The Mie results in Figure 1 correspond to resonant absorption from an isolated 16 nm Au nanosphere in the Rayleigh regime [15] whose scattering cross section is negligible [16]. Extinction per isolated AuNP in PDMS is constant since particle composition and dielectric environment remain unchanged. Optical effects of anomalous nonresonant scatterers in the AuNP-polymer matrix were removed by calculating resonant absorption by difference from off-resonant absorption at 800 nm. Control experiments were performed with AuNP-PDMS films containing 76 nm Au nanospheres with appreciable scattering cross section [29]. These films exhibited optical responses up to 1.3-fold greater than Mie estimates at times the , attributable to multiple scattering. Figure 1(b) indicates extinction per NP values consistent with Mie and CDA estimates at . These results indicate scattering from isolated particles appears is unlikely from PDMS films containing 16 nm AuNP. Thus, isolated particle effects due to scattering from resonant or nonresonant scatterers appear unlikely to contribute significantly to increased optical response at or near resonance .

Considering particle-particle interactions, the CDA results in Figure 1 correspond to collective extinction from one such 16 nm Au nanosphere interacting with others in a lattice. Extinction per NP decreases in CDA relative to Mie due to particle-particle interactions. Extinction per interacting AuNP in lattices at values of examined here remained constant since separation of particles with negligible scattering cross section remains sufficiently large, that is, in the far-field, such that negligible effect on the LSPR resonance accrues. Larger particles that resulted in appreciable scattering could accrue multiple scattering. Multiple scattering has been modeled by radiative transport and diffusive theory, typically for opaque samples that lack ballistic light transmission [66, 67]. However, negligible scattering cross sections for 16 nm AuNPs minimize effects of multiple scattering. On the other hand, near-field interactions at particle separation much less than values, that is, on the order of AuNP size (hot spots), have been reported to enhance resonance extinction [68, 69]. In this work, the 16 nm AuNPs studied were coated with ~10 nm polyvinylpyrrolidone (PVP) layer that diminished likelihood of such near-field interactions. This coating was observed to eliminated AuNP aggregation that was previously found after similar film fabrication using uncoated AuNP [31, 32]. In addition, LSPR wavelengths of AuNP-PDMS films blue-shifted as AuNP content increased from 542 to 532 nm rather than red-shifting as might be expected from near-field interaction. Moreover, optical responses of 16 nm AuNPs at up to 103 the concentrations evaluated here were within 12 and 17% of values predicted by CDA and Mie theory, respectively, in films sufficiently thin to preclude optical interference. Scanning electron microscope images of such films did not exhibit close particle interactions. Furthermore, random occurrence of such near-field interactions would not be expected to result in the observed linear divergence which was duplicated reproducibly at three particle concentrations. Therefore, particle-particle interactions appeared unlikely to contribute to increased optical response for .

Considering optical interference, reproducible extinction per NP values consistent with Mie and CDA estimates at in Figure 1(b) indicate a prevailing optical interference arises only when interparticle distance decreases to values less than resonant wavelengths. The random distribution of the dispersed AuNP precludes optical interference from ordered scatterers, for example, lattice diffraction [51, 57, 58, 60, 68, 69]. Use of Maxwell Garnett effective medium theory [70, 71] to calculate change in effective dielectric constant due to increasing AuNP content yields extinction per NP estimates consistent with Mie and CDA. This absolves interference, for example, refraction, from dielectric heterogeneity. Geometric interference alone, for example, reflection at the PDMS-air interface [72], does not produce the increase at interparticle separation greater than >. Remaining possibilities include diffraction in scattering systems, reported for coupled waveguides and photonic crystals [73], which could be measured using dynamic scattering, or time-averaged intensity of scattered light [74]. Control experiments with 76 nm AuNP with appreciable scattering cross section demonstrated that divergence of optical rays, for example, by multiple scattering or diffraction, convolved with internal reflection could increase observed extinction per AuNP relative to Mie-calculated values [29]. Thus, increased optical response for could be attributed to a geometric optical effect at the dielectric interface between polymer and air convolved with interference to incident radiation produced by adjacent resonant AuNP cross sections for .

3.2. Heat Dissipation in AuNP-PDMS Films with Decreasing Interparticle Separation

Dissipation of power absorbed by AuNP dispersed in PDMS films was characterized by finite element analysis (FEA) using relations for heat conduction, natural convection, and radiation. Figure 2(a) compares measured (top) and simulated (bottom) temperature profiles for the 1.17 × 1012 NPs/cm3 AuNP-PDMS film. A resonant absorbed power was introduced in order to calculate equilibrium temperature profiles that matched measured thermal profiles. Values for this power at consecutively higher AuNP contents increased from 1.2 to 11.1 mW as shown in Table 1, which summarizes optically extinguished powers and FEA-fitted heat dissipated powers. Measured, Mie, and CDA-predicted optically extinguished powers () were calculated as , where was 532 nm spectral extinction in absorbance units, based on measured, Mie, and CDA theory, respectively, and was incident power of 532 nm irradiation (19–24 mW) for each AuNP-PDMS film. Spectral extinction was estimated from Mie theory and CDA using Beer-Lambert absorbance [35] as reported in Section 3.1.

CDA-predicted optically extinguished power corresponded to FEA-fitted heat dissipated power derived from measured thermal response across the AuNP-PDMS films. Figure 2(b) plots FEA-derived thermal power dissipated (-axis) versus optically extinguished power (-axis) based on a priori CDA prediction (green hollow triangles) and measured optically extinguished power (filled blue circles). The blue dashed line indicates one-to-one correspondence. Overall agreement between CDA-predicted and FEA-derived powers at , that is, AuNP ≤ 1.17 × 1012 NP/cm3, was consistent with ~98% agreement between CDA-derived and measured extinction per NP as shown in Figure 1(b). Incorporation of particle-particle interactions in CDA improved agreement with measured optical responses relative to Mie results, which reflected only isolated particle effects and are thus excluded from Figure 2(b).

For , FEA-derived heat dissipated power lagged measured power values. The observed increase in extinction at these particle spacings appeared not to correspond to increased power dissipation. The optically extinguished power increased linearly with increasing incident power up to 25 mW, as shown in Figure 2(b) inset for the 3.52 × 1012 NPs/cm3 AuNP-PDMS film. This suggests that neither aggregate thermal effects on particle and media dielectric values nor saturation measurably alters measured optical extinction relative to values obtained spectroscopically across the power and temperature ranges examined. Suppressed heat dissipation relative to optically extinguished power appeared attributable in part to edge effects. Temperature distributions of heated samples exceeded ambient values at film boundaries, particularly for off-center laser incidence. The heat transfer coefficients that were employed underpredicted heat transfer from the film’s narrow vertical surface (~5 mm × ~1 mm) due to unaccounting for inward fluid flow induced near the edge. This effect has been reported to increase heat transfer rates by up to double for narrow, vertical surfaces [75]. As a result, modeled and measured temperature values diverged near film edges. Edge effects on geometric optics near film edges may have contributed to less power absorbed than that accounted for in measured extinction. Samples in which incident irradiation was at edges could have been affected by geometric optics, for example, reflection and refraction. FEA-derived heat dissipation values were subject to user specification of absorbed power distributions in the FEA model which was important for characterizing sensitivity of incident beam distribution on resulting FEA-fitted heat source powers. The Gaussian power distribution used for all data shown enabled heat dissipated powers that reached an average temperature within the laser spot that matched measured results within 0.1°C. In contrast, uniformly distributing power over the 1.5 mm radius laser spot broadened the temperature profile with maximum (minimum) temperatures below (above) measured values. Decreasing power distribution radius to 0.5 mm resulted in maximum (minimum) temperatures above (below) measurements.

4. Conclusions

This work compared absorbed power derived from measured and predicted optical extinction with dissipated power obtained fitting measured temperature distributions to a finite element model in PDMS films containing 16 nm AuNPs with negligible scattering cross sections. Measured optical extinction corresponded to values predicted using CDA which accounted for particle-particle interactions at particle separation greater than resonance wavelength. Measured extinction exceeded CDA predictions at particle separation below resonance wavelength. Dissipated power estimated from measured temperature profiles corresponded to CDA-predicted optical extinction values at all particle separation studied. Spectroscopically measured extinction values corresponded to measured power extinguished up to 25 mW. But increases in spectroscopically measured extinction values relative to CDA-predicted extinction at particle separation less than resonance wavelength appeared to not result in corresponding increases in heat dissipation. Prediction of optical absorbance by Mie theory and coupled dipole approximation accounted for effects of isolated particle and particle-particle interactions, respectively, but not for interference with incidence radiation or dielectric heterogeneity, that is, geometric optical effects. Considered effects of near-field interactions and multiple scattering on measured optical responses did not appear significant for the particles studied. Coordinated examination of isolated particle, particle-particle, and dielectric heterogeneity effect on optothermal responses comparing experiment and theory is important to guide intuition, design, and development of responsive plasmonic energy materials for a variety of applications.

Disclosure

Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported in part by NSF CBET-1134222, NSF EEC-1260301, NSF Graduate Research Fellowship awarded to Jeremy R. Dunklin, the University of Arkansas Foundation, and the Walton Family Charitable Foundation. The authors graciously acknowledge Caitlyn Chambers for thermal data collection, Keith Berry, Jr. for conversion of temperature maps into three dimensions, and Gregory T. Forcherio for helpful discussions.