<figure><img src="image_1.png" /></figure>

HAL
open science

# Microscale Temperature Shaping Using Spatial Light Modulation on Gold Nanoparticles

Ljiljana Durdevic, Hadrien M L Robert, Benoit Wattellier, Serge Monneret, Guillaume Baffou

► To cite this version:

Ljiljana Durdevic, Hadrien M L Robert, Benoit Wattellier, Serge Monneret, Guillaume Baffou. Microscale Temperature Shaping Using Spatial Light Modulation on Gold Nanoparticles. Scientific Reports, 2019, 9 (1), 10.1038/s41598-019-40382-3. hal-02073966

HAL Id: hal-02073966

https://hal.science/hal-02073966

Submitted on 24 Nov 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_3.png" /></figure>

OPEN

# Microscale Temperature Shaping Using Spatial Light Modulation on Gold Nanoparticles

Ljiljana Durdevic<fnref n="1" />, Hadrien M. L. Robert<fnref n="1" />, Benoit Wattellier², Serge MonneretD&G<fnref n="1" /> & Guillaume Baffou<fnref n="1" />

Received: 19 October 2018

Accepted: 12 February 2019

Published online: 15 March 2019

Heating on the microscale using focused lasers gave rise to recent applications, e.g., in biomedicine, biology and microfluidics, especially using gold nanoparticles as efficient nanoabsorbers of light. However, such an approach naturally leads to nonuniform, Gaussian-like temperature distributions due to the diffusive nature of heat. Here, we report on an experimental means to generate arbitrary distributions of temperature profiles on the micrometric scale (e.g. uniform, linear, parabolic, etc) consisting in illuminating a uniform gold nanoparticle distribution on a planar substrate using spatially contrasted laser beams, shaped using a spatial light modulator (SLM). We explain how to compute the light pattern and the SLM interferogram to achieve the desired temperature distribution, and demonstrate the approach by carrying out temperature measurements using quantitative wavefront sensing.

Heating over a microscale area is becoming an important concept with the development of nano- and microtechnologies. In particular, heating gold nanoparticles by light absorption is at the basis of a more and more active field of research named thermoplasmonics<fnref n="1" />,², addressing problems in biology³⁻⁵, biomedicine⁶⁻⁸, microscale fluid dynamics⁹,¹⁰, phase transitions (bubbles)¹¹⁻¹⁵, thermophoresis¹⁶,¹⁷ or chemistry¹⁸⁻²⁰.

The task of measuring the temperature on the microscale is now well-mastered by numerous optical microscopy techniques<fnref n="1" />,²¹, usually based on fluorescence measurements, more rarely label-free²²⁻²⁵, sometimes even in three dimensions²⁶, and with a diffraction-limited spatial resolution. However, controlling the temperature spatial distribution is, for the time being, overlooked. Light-heating of gold nanoparticle distributions usually results in non-uniform, Gaussian-like profiles, due to the diffusive nature of heat²⁷. Yet, for some applications, the temperature spatial profile can play an important role. For instance, any temperature gradient can generate important thermophoretic forces on colloids dispersed in liquids¹⁶,¹⁷,²⁸, or non-uniform temperatures may be detrimental when working with biological systems, sensitive to temperature variations of typically 0.5 K⁴. For such applications, it would be important to accurately monitor the gradients or control the uniformity of the microscale temperature profile.

In a previous work, we achieved microscale temperature shaping at will by illuminating sophisticated and non-uniform nanoparticle distributions, made by e-beam lithography, using uniform laser beams²⁹. Albeit effective, this approach suffers from a lack of flexibility: each lithographed area is associated with a predefined temperature profile that could neither be dynamically changed, nor moved to any other area of interest. This may be a stringent constraint for instance in applications involving living cells in culture, whose position can never be predicted in advance.

In this article, we introduce an experimental procedure to dynamically generate distributions of temperature profiles on the micrometric scale with arbitrary complexity by illuminating a uniform gold nanoparticle distribution using spatially contrasted laser beams, shaped using a spatial light modulator (SLM). We explain hereinafter how to compute the light pattern and the SLM interferogram to achieve the desired temperature distributions, explain the benefits of using gold nanoparticles and demonstrate the approach by carrying out temperature measurements using quantitative phase imaging. A final part is dedicated to an illustration of the interest of this approach for the field of thermophoresis of colloids.

<footnote n="1">¹Institut Fresnel, CNRS, Aix Marseille Univ, Centrale Marseille, Marseille, France. ²PHASICS S.A., Parc technologique de Saint Aubin, Route de l'Orme des Merisiers, 91190, Saint Aubin, France. Correspondence and requests for materials should be addressed to G.B. (email: guillaume.baffou@fresnel.fr)</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** (a) Scheme of the microscope setup. A spatial light modulator (SLM) is used to create arbitrary laser beam profiles at the sample (S) plane, via a diaphragm (D), to remove the zero-order diffraction spot) and a dichroic beam splitter (DBS). A Köhler illumination (KI) is implemented on top of the sample. A QSI wavefront sensor, consisting of sCMOS camera assembled to a reimaging system and coupled with a two-dimensional grating (commonly named a Modified Hartmann Mask, MHM) are used to map the wavefront distortion created by the temperature gradient resulting from the laser heating of the gold nanoparticles at the sample plane. (b) Schematic of the sample, which consists of a water layer sandwiched between two glass coverslips. In our study, the nanoparticles are gold nanorods. They are deposited on top of the lower coverslip, and are meant to be heated by a laser illumination through the objective lens (40×). (c) SEM image of the gold nanorods. (d) Typical extinction spectrum of the gold nanorod samples used in this study.</figcaption></figure>

## Results and Discussion

A single home-made microscope was used to both heat substrates of gold nanoparticles, and measure the resulting temperature distribution, on the microscale. The microscope configuration is sketched in Fig. 1a. A spatial light modulator (SLM, HSP256-1064 Boulder Nonlinear Systems) was used to shape the profile (see next subsection) of a laser beam (Ti:Sapph laser, $\lambda = 800$ nm) that was sent to the sample using a dichroic beam splitter (DBS). The association of a Köhler illumination (Thorlabs LED $\lambda = 625 \pm 12$ nm), a modified Hartmann mask (MHM, Phasics SA), a reimaging system (Sid4-Element, Phasics) and a sCMOS camera (Zyla, Andor) were used to acquire the temperature distribution of the sample resulting from the heating of the gold nanoparticles. This label-free technique, called TIQSI for temperature imaging using quadriwave lateral shearing interferometry, is described in a previous publication³⁰. This technique does not only enable the mapping of the temperature but also of the heat source density (power per unit area).

The sample is depicted in Fig. 1b. It was composed of a water layer sandwiched between two glass coverslips (1 in × 1 in). Depending on the experiments presented in this article, the water layer thickness could vary from a few microns to 1 mm. On the bottom coverslip, a uniform layer of gold nanorods ($46 \pm 4$ nm in length and $20 \pm 2$ nm in diameter) featuring a plasmonic resonance around the laser wavelength (800 nm) was deposited over the whole area of the coverslip (see ref.⁴ for details regarding the fabrication of this kind of sample).

When such a gold nanoparticle substrate is illuminated using a uniform laser beam profile, the heat source density is perfectly uniform due to the uniform distribution of nanoparticles. Importantly, despite the discrete and nanometric nature of the sources of heat (the nanoparticles), the temperature profile does not consist of an assembly of nanoscale hot spots. On the contrary, the temperature distribution is smooth and continuous over the heated area, i.e., over 10 s of microns, due to collective thermal effects²⁷. Thus, the heated area can be considered as a microscale hot plate delivering a uniform heat power density (HPD, power per unit area). In such a condition, even if the heat power per unit area is uniform (i.e., even if the laser intensity is uniform over a given area), the temperature is not uniform. The center of the heated area is hotter than the boundary due to thermal diffusion²⁷,²⁹, as illustrated by Fig. 2. In 2014, we introduced a technique aimed at deviating from this Gaussian-like shape to obtain almost any temperature distribution on the microscale²⁹. The technique consisted in shaping the HPD by illuminating non-uniform gold nanoparticle distributions. As explained above, although this technique is effective, it does not enable a *dynamic* shaping of the temperature profile. Once lithographed, the specific nanoparticle distribution results in a fixed temperature profile at a fixed location of the sample. Here, we propose another approach to spatially shape the HPD, which consists in using a *uniform* nanoparticle distribution and a *non-uniform* laser beam. Because the HPD is now produced by a laser beam, the spatial resolution of the HPD is limited by the diffraction (typically 500 nm), while it was limited by the capabilities of e-beam lithography in our previous study (typically 100 nm). For this reason, this new approach yields a slightly poorer spatial resolution compared to the lithography approach. However, this limited spatial resolution may not be a problem for many applications and this light-shaping approach enables one to dynamically position and shape microscale
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Heating of a uniform carpet of nanoparticles using a circular and uniform laser beam characterized by TIQSI. (a) Heat source density (scale bar 47 µm, the diameter of the laser beam). (b) Associated temperature distribution, featuring a non-uniform temperature distribution.</figcaption></figure>

temperature profiles. This is a requisite for several applications, e.g., in microscale thermophoresis of colloids or cell biology, where the location and motion of colloids or cells on the sample cannot be predicted.

For a heat source located at the interface between two semi-infinite media separated by a planar interface, the temperature everywhere in the universe reads

$$T(\mathbf{r}) = [G \otimes p](\mathbf{r}) \qquad (1)$$

where $G(\mathbf{r}) = 1/(4\pi\bar{\kappa}|\mathbf{r}|)$ is the thermal Green's function and $p(\mathbf{r})$ the HPD (power per unit area at the planar interface) with $\bar{\kappa} = (\kappa_1 + \kappa_2)/2$ is the average of the two thermal conductivities of the two semi-infinite media. This Green's function can be used when the liquid layer is much thicker than the size of the laser beam impinging on the nanoparticles, so that the liquid layer thickness can be considered as infinite. If the water layer is reduced, then the Green's function has to be modified accordingly. It cannot be considered anymore as scaling as $1/r$ in any direction. One has to consider the three-layer Green's function for accurate computation$^{26,31}$. This convolution between $G$ and $p$ can be numerically inverted to obtain the function $p$ corresponding to a desired $T_{\text{target}}$ distribution throughout the interface$^{29}$:

$$p(\mathbf{r}) = \mathcal{F}^{-1}[\mathcal{F}(T_{\text{target}}(\mathbf{r}))/\mathcal{F}(G(\mathbf{r}))] \qquad (2)$$

where $\mathcal{F}$ represents the Fourier transform. As explained above, the approach we introduce consists in reproducing this HPD using a spatially contrasted laser beam on a uniform absorbing layer, composed here of gold nanoparticles. Hence, the laser beam intensity profile $I_{\text{laser}}$ (power per unit area of the laser beam, at the interface) has to reproduce the HPD $p$ to within a scalar multiplication factor. This factor equals the absorbance $A$ of the sample:

$$I_{\text{laser}}(\mathbf{r}) = p(\mathbf{r})/A \qquad (3)$$

The profile of laser beam can be adjusted using different approaches. We chose to use a spatial light modulator (SLM) to achieve this task. The SLM we used consisted of a 2D-matrix of liquid crystals applying a phase shift to a laser beam in reflection. When placed within the optical path of a light beam entering the objective lens of a microscope, it can be used to modify the laser intensity by only modifying its phase in the Fourier plane, i.e., at the entrance pupil of the objective lens. This is why the SLM is optically conjugated with the objective lens location in Fig. 1. To compute the phase profile $\Phi(\mathbf{r})$ to be applied in the Fourier space to generate a given light intensity profile $I_{\text{laser}}(\mathbf{r})$ in the sample plane, we used the Gerchberg-Saxton phase-retrieval algorithm (GSA, see further on for a more detailed discussion)$^{32}$. This iterative algorithm enables the determination of a phase profile $\Phi(\mathbf{r})$ such that

$$I_{\text{laser}}(\mathbf{r}) = \mathcal{F}[I_0(\mathbf{r}) \exp(i\Phi(\mathbf{r}))] \qquad (4)$$

where $I_0(\mathbf{r})$ is the light intensity profile impinging on the SLM.

Figure 3 summarizes the different steps (algorithms and experiments) to achieve temperature shaping and control using spatial light modulation (TS-SLM). The first step consists in choosing a 2-dimensional temperature profile over an area of interest. Figure 3a considers the particular example of a uniform temperature increase over a circular area. Note that some singular temperature profiles cannot be created. For instance, specifying negative temperature variations won't obviously be possible. This would generate a calculated HPD with negative values. Also, too steep gradients may also generate negative HPD in some pixels. However, the specification of uniform
COMPUTATIONS

<figure><img src="image_6.png" /></figure>

MEASUREMENTS

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Experimental procedure of the TS-SLM technique. (a) Choice of a targeted microscale temperature profile. (b) Numerical calculation of the required heat power density to achieve this temperature profile. (c) Associated laser beam profile to achieve this heat power density. (d) Phase of the light beam at the entrance of the objective lens to achieve the laser beam profile (e) at the sample plane. (e,f) Experimental measurements of the temperature and heat power density using TIQSI. (g) Image of the light scattering by the gold nanoparticle layer to render the profile of the light beam at the sample plane. Images (e-g) have to be compared with their analog images (a-c).</figcaption></figure>

temperature profiles, no matter the complexity of the shape, necessarily gives rise to valid HPD. The HPD (Fig. 3b) has to be calculated by deconvolution (see Eq. (2)). The core of the technique consists in generating this HPD by illuminating a uniform absorbing layer (here composed of gold nanoparticles) by a light pattern (Fig. 3c) that reproduces the HPD profile. The light pattern irradiance can be determined from the HPD provided the absorption of the sample is known (see Eq. (3)), which is usually not the case. This quantity is difficult to measure. One usually rather measures the extinction, which arises not only from absorption but also from scattering. But determining A is not a requisite in our approach. Besides, we did not do it in our study. By arbitrarily setting the laser power, the temperature distribution will eventually feature the appropriate profile, up to a constant factor, and the amplitude of the temperature increase could be adjusted a posteriori by playing with the laser power. The implementation of the temperature microscopy technique is then mandatory to monitor the overall temperature increase and adjust the laser power accordingly to achieve the appropriate temperature increase. To obtain the targeted light pattern, one option could be to use an array of micromirrors (digital light processor, DLP) but we preferred to opt for a spatial light modulator in this study. It more elegantly redistributes the light intensity in the Fourier space instead of switching pixels on and off and loosening a lot of laser intensity. A DLP would do the job, provided a laser of a few watts is used. The phase profile to be applied to the SLM to generate the appropriate light pattern was calculated using the Gerchberg–Saxton algorithm (GSA)³². Because no simple inversion procedure exists, GSA is based on an iterative process. The code is simple (no more than 10-line long, see Suppl. Info.), but some refinement is often implemented, in particular to remove the occurrence of speckle, inherent to this approach (see the speckle in Fig. 3g)³³,³⁴. But in our case, thermal diffusion smoothes any small non-uniformities of the heat generation on the microscale, so that any procedure to damp the speckle is not mandatory. It is, however, mandatory to remove the zero-order spot, that would create a heating spot in the middle of the image. For this purpose, a linear phase gradient is applied to the phase image calculated by the GSA in order to laterally shift the light pattern on the sample plane away from the zero-order of diffraction. Then, a diaphragm (D) is positioned at an intermediary image plane (see Fig. 1a) to remove this light spot and only keep the light pattern of interest.

Once the appropriate light pattern is projected onto the sample, the resulting temperature field has to be controlled. For this purpose, many temperature microscopies have been developed, most of them based on optical measurements, usually involving fluorescence (see Chap. 4 of ref. ¹). In this study, we used QSI as it is label-free: it does not involve the addition of any molecular temperature probe since the measurement is based on temperature-induced variations of the refractive index of the liquid layer³⁰. Moreover, this technique enables not
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Thermophoresis-driven motion of 28 beads. (a) Linear temperature gradient produced by TS-SLM. (b) Dark-field image of beads dispersed in the liquid layer. (c) Trajectory of the beads with no heating during 14 seconds (dark crosses) exhibiting a Brownian motion, and upon heating (brown crosses) undergoing a linear motion following the temperature gradient.</figcaption></figure>

only the measurement of the temperature, but also of the HPD, quantitatively. Such measurements are presented in Fig. 3e,f, and they are in agreement with the targeted temperature and HPD (see Fig. 3a,b). QSI also enables the acquisition of intensity images (not only phase), as displayed in Fig. 3g. In that case, the optical filter was changed in front of the camera to remove the light from the Köhler illumination and let the laser pass.

Gold nanoparticles have been our preferred choice for this study. Although any two-dimensional, uniform absorbing layer should be effective as well, our choice was driven by the multiple benefits provided by gold nano-particles: First, when using sufficiently small nanoparticles (typically below 40 nm), it is possible to make absorption the only interaction process existing between light and nanoparticles: (i) Scattering can be made negligible because it scales as the volume squared, while the absorption scales with the volume of the nanoparticle. Then, no reflection of the laser beam occurs like what would happen with a metal layer for instance, acting as a mirror. Thus, the transduction of the light energy can be optimized, with no light energy lost in unwanted side-effects. In our case, gold nanorods 46 × 20 nm in size feature an absorption/scattering ratio of 7 at resonance.

Second, unlike dyes, gold nanoparticles do not photobleach, do not oxidize, and get damaged/reshaped only at temperatures higher than 150 °C. This makes gold nanoparticles samples particularly stable and robust. Third, although gold is an expensive metal, several inexpensive bottom-up fabrication approaches exist to yield uniform distributions of nanoparticles over macroscopic areas. Fourth, gold nanoparticles can feature resonant absorption in the infrared, if they deviate from a spherical geometry, like in this study. This offers the possibility to build a sample that is absorbant in the infrared, but which remains quite transparent in the visible range. This is highly desirable to ease any fluorescence imaging for example.

The only limitation of using gold nanoparticles is that they can reshape when typically exceeding 150 °C to 200 °C (under cw illumination) depending on the heating duration³⁵. Reshaping induces a variation of their absorption efficiency, which would compromise the required uniformity of the sample. To avoid this problem in high-temperature applications, the gold nanoparticle layer may be embedded in a silica matrix³⁶.

In the previous part, we chose to illustrate the principle of the TS-SLM approach on a uniform temperature distribution because it is particularly suited for experiments in chemistry²⁰ and cell biology⁴, where the setting of a precise temperature over the observation area may be mandatory. But another identified relevant temperature profile for other applications consists of *linear temperature gradients*, e.g., for the study of microscale thermophoresis in liquids¹⁶,²⁸,³⁷. This is what we wish to evidence in Fig. 4. In this experiment, silica beads, 1.24 µm in diameter, were dispersed in milli-Q water, with a fluid thickness of 34 µm. Reducing the thickness of the fluid enables the damping of fluid convection, that would cancel the thermophoretic motion of the beads. Figure 4a displays a linear temperature gradient of 0.36 K/µm applied over a rectangular area by TS-SLM. In this experiment, a super-heated state of water (water above 100 °C without boiling) was achieved¹⁵. Figure 4b is a dark-field image of the beads dispersed in water. When no heating is performed, these beads are simply undergoing Brownian motion at ambient temperature (21 °C). When the linear temperature gradient is applied, a strong thermophoretic motion is observed, moving the beads from hot to cold. Video V1, provided in Suppl. Info., shows such a process. During the first 14 s, no heating was applied. Then, heating was applied during 14 s, and the beads exhibited a drift. At this scale, the temperature distribution observed in Fig. 4a establishes on the millisecond scale, and can be considered as instantaneous, compared to the time scale of the physical process under consideration. Figure 4c
represents the motion of the beads of Fig. 4b and in the video. Dark crosses correspond to the no-heating phase, lasting 14 s. No drift in specific direction is observed, only Brownian motion. The brown crosses correspond to the first 14 s of heating, where clear directed motions are observed, along the temperature gradient. The four beads in the centre of the image, where the temperature gradient is linear, enable the quantification of the thermophoretic properties of the beads. In this temperature range (80–120 °C), the measured thermophoretic mobility $D_T$ is 3.5 ± 0.6 µm²/(K·s). Refined measurements could be performed with statistics on a larger number of beads, but this will be the purpose of a future work dedicated to microscale thermophoresis in liquids. The idea here was to illustrate the applicability and the interest of TS-SLM for an original application.

## Conclusions

In summary, we introduced an optical microscopy technique suited to generate microscale temperature patterns with arbitrary shapes, in particular uniform distributions and linear temperature gradients, the most useful cases for envisioned applications. This technique involves a spatial light modulator (to generate the temperature profile) and a wavefront sensing camera (to measure the temperature profile) or any other temperature mapping technique. This approach enables a dynamic control of the temperature profile, even down to the millisecond scale. This experimental scheme can be simply implemented on any regular microscope. Applications are envisioned when an accurate temperature setting over the observation area of interest in mandatory, for instance in microscale thermal-assisted chemistry, fluid dynamics, cell biology and more generally in most applications of thermoplasmonics.

## Methods

**Gold nanoparticle sample fabrication.** The gold nanorods synthesis is adapted from the protocols developed by Nikoobakht et al.38 and Liu et al.39, and detailed in ref.4. The as-prepared gold nanorods were functionalized with PVP. Before the deposition of the PVP-functionalized gold nanorods onto the surface of the coverslip, the latter has been recovered by 3 layers of polyelectrolytes as described in the work of Vial et al.40. The functionalized coverslips were immersed in a solution of PVP-gold nanorods for 3 hours. The deposition was promoted via electrostatic interactions between the positively charged coverslip and the negatively charged gold nanorods.

## References

1. Baffou, G. *Thermoplasmonics*. Heating Metal Nanoparticles with Light (Cambridge University Press, 2017).

2. Baffou, G. & Quidant, R. Thermo-plasmonics: Using metallic nanostructures as nano-sources of heat. *Laser & Photon. Rev.* 7, 171–187 (2013).

3. Li, M., Lohmüller, T. & Feldmann, J. Optical injection of gold nanoparticles into living cells. *Nano Lett.* 15, 770–775 (2015).

4. Robert, H. M. L. et al. Photothermal control of heat-shock protein expression at the single cell level. *Small* 14, 1801910 (2018).

5. Bahadori, A., Oddershede, L. B. & Bendix, P. M. Hot-nanoparticle-mediated fusion of selected cells. *Nano Research* 10, 2034–2045 (2017).

6. Lal, S., Clare, S. E. & Halas, N. J. Nanoshell-enabled photothermal cancer therapy: Impending clinical impact. *Acc. Chem. Res.* 41, 1842 (2008).

7. Paithankar, D. Y. et al. Acne treatment based on selective photothermolysis of sebaceous follicles with topically delivered light-absorbing gold microparticles. *J. Invest. Dermatol.* 135, 1727–1734 (2015).

8. Kharlamov, A. N. et al. Silica-gold nanoparticles for atheroprotective management of plaques: Results of the nanom-fim trial. *Nanoscale* 7, 8003–8015 (2015).

9. Namura, K., Nakajima, K. & Suzuki, M. Quasi-stokeslet induced by thermoplasmonic marangoni effect around a water vapor microbubble. *Sci. Rep.* 7, 45776 (2017).

10. Donner, J. S., Baffou, G., McCloskey, D. & Quidant, R. Plasmon-assisted optofluidics. *ACS Nano* 5, 5457–5462 (2011).

11. Siems, A., Weber, S. A. L., Boneberg, J. & Plech, A. Thermodynamics of nanosecond nanobubble formation at laser-excited metal nanoparticles. *New J. Phys.* 13, 043018 (2011).

12. Hou, L., Yorulmaz, M., Verhart, N. R. & Orrit, M. Explosive formation and dynamics of vapor nanobubbles around a continuously heated gold nanosphere. *New J. Phys.* 17, 013050 (2015).

13. Lombard, J., Biben, T. & Merabia, S. Nanobubbles around plasmonic nanoparticles: Thermodynamic analysis. *Phys. Rev. E* 91, 043007 (2015).

14. Lachaine, R., Boulais, E. & Meunier, M. From thermo- to plasma-mediated ultrafast laser-induced plasmonic nanobubbles. *ACS Photonics* 1, 331–336 (2014).

15. Baffou, G., Polleux, J., Rigneault, H. & Monneret, S. Super-heating and micro-bubble generation around plasmonic nanoparticles under cw illumination. *J. Phys. Chem. C* 118, 4890 (2014).

16. Braun, M., Würger, A. & Cichos, F. Trapping of single nano-objects in dynamic temperature fields. *Phys. Chem. Chem. Phys.* 16, 15207 (2014).

17. Braun, M., Bregulla, A. P., Günther, K., Mertig, M. & Cichos, F. Single molecules trapped by dynamic inhomogeneous temperature. *Nano Lett.* 15, 5499–5505 (2015).

18. Baffou, G. & Quidant, R. Nanoplasmonics for chemistry. *Chem. Soc. Rev.* 43, 3898–3907 (2014).

19. Ni, W., Ba, H., Lutich, A. A., Jäckel, F. & Feldmann, J. Enhancing single-nanoparticle surface-chemistry by plasmonic overheating in an optical trap. *Nano Lett.* 12, 4647–4650 (2012).

20. Robert, H. M. L. et al. Light-assisted solvothermal chemistry using plasmonic nanoparticles. *ACS Omega* 1, 2–8 (2016).

21. Brites, C. D. S. et al. Thermometry at the nanoscale. *Nanoscale* 4, 4799 (2012).

22. Setoura, K., Werner, D. & Hashimoto, S. Optical scattering spectral thermometry and refractometry of a single gold nanoparticle under cw laser excitation. *J. Phys. Chem. C* 116, 15458–15466 (2012).

23. Chen, Z. et al. Imaging local heating and thermal diffusion of nanomaterials with plasmonic thermal microscopy. *ACS Nano* 9, 11574–11581 (2015).

24. Xie, X. & Cahill, D. G. Thermometry of plasmonic nanostructures by anti-stokes electronic raman scattering. *Appl. Phys. Lett.* 109, 183104 (2016).

25. Carattino, A., Caldarola, M. & Orrit, M. Gold nanoparticles as absolute nanothermometers. *Nano Lett.* 18, 874–880 (2018).

26. Bon, P. et al. Three-dimensional temperature imaging around a gold microwire. *Appl. Phys. Lett.* 102, 244103 (2013).

27. Baffou, G. et al. Photoinduced heating of nanoparticle arrays. *ACS Nano* 7, 6478–6488 (2013).

28. Würger, A. Thermal non-equilibrium transport in colloids. *Rep. Prog. Phys.* 73, 126601 (2010).

29. Baffou, G. et al. Deterministic temperature shaping using plasmonic nanoparticle assemblies. *Nanoscale* 6, 8984–8989 (2014).

30. Baffou, G. et al. Thermal imaging of nanostructures by quantitative optical phase analysis. *ACS Nano* 6, 2452–2458 (2012).
31. Ramanan, V. S., Muthukumar, M., Gnanasekaran, S., Venkataramana Reddy, M. J. & Emmanuel, B. Green's functions for the laplace equation in a 3-layer medium, boundary element integrals and their application to cathodic protection. *Eng. Anal. Bound. Elem.* **23**, 777–786 (1999).

32. Gerchberg, R. W. & Saxton, W. O. A practical algorithm for the determination of the phase from image and diffraction plane pictures. *Optik* **35**, 237 (1972).

33. Chen, C. Y. et al. Speckle reduction by combination of digital filter and optical suppression in a modified gerchberg-saxton algorithm computer-generated hologram. *Appl. Opt.* **53**, G163–G168 (2014).

34. Chang, C., Xia, J., Yang, L., Yang, Z. & Chen, J. Speckle-suppressed phase-only holographic three-dimensional display based on double-constraint gerchberg-saxton algorithm. *Appl. Opt.* **54**, 6994–7001 (2015).

35. Petrova, H. et al. On the temperature stability of gold nanorods: Comparison between thermal and ultrafast laser-induced heating. *Phys. Chem. Chem. Phys.* **8**, 814–821 (2006).

36. Chen, Y. S. et al. Enhanced thermal stability of silica-coated gold nanorods for photoacoustic imaging and image-guided therapy. *Opt. Express* **18**, 8867–8877 (2010).

37. Fayolle, S., Bickel, T. & Würger, A. Thermophoresis of charged colloidal particles. *Phys. Rev. E* **77**, 041404 (2008).

38. Nikoobakht, B. & El-Sayed, M. A. Preparation and growth mechanism of gold nanorods (nrs) using seed-mediated growth method. *Chem. Mater.* **15**, 1957–1962 (2003).

39. Liu, M. & Guyot-Sionnest, P. Synthesis and optical characterization of au/ag core/shell nanorods. *J. Phys. Chem. B* **108**, 5882–5888 (2004).

40. Vial, S., Pastoriza-Santos, I., Pérez-Juste, J. & Liz-Marzán, L. M. Plasmon coupling in layer-by-layer assembled gold nanorod films. *Langmuir* **23**, 4606–4611 (2007).

## Acknowledgements

This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement no. 772725, project HiPhore). The authors thank Sophie Brasselet and Carolina Rendon Barrasa for their valuable help in this work. BW is a member of the Phasics company, the manufacturer of the wavefront sensing camera used in the reported experiments.

## Author Contributions

L.D. mounted the setup. L.D. and H.M.L.R. performed the experiments and postprocessed the data. H.M.L.R. fabricated the gold nanorod samples. G.B. designed the project, supervised the experiments and wrote the first version of the article. S.M. provided his expertise on quantitative phase imaging and spatial light modulation. All the authors discussed the results and contributed to write the manuscript.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-40382-3.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019